Protocol for FRAP-based estimation of nuclear import and export rates in single yeast cells.

Here, we present a protocol for estimating nuclear transport parameters in single cells. We describe steps for performing four consecutive fluorescence recovery after photobleaching experiments, fitting the obtained data to an ordinary differential equations model, and statistical analysis of the fittings using a specialized R package. This protocol permits the estimation of import and export rates, nuclear or cytosolic fixed fractions, and total number of molecules. For complete details on the use and execution of this protocol, please refer to Durrieu et al.1.

rcell2: Microscopy-Based Cytometry in R.

This article describes a method for quantifying various cellular features (e.g., volume, curvature, total and sub-cellular fluorescence localization) of individual cells from sets of microscope images, and for tracking them over time-course microscopy experiments. One purposely defocused transmission image (sometimes referred to as bright-field or BF) is used to segment the image and locate each cell. Fluorescence images (one for each of the color channels or z-stacks to be analyzed) may be acquired by conventional wide-field epifluorescence or confocal microscopy. This method uses a set of R packages called rcell2. Relative to the original release of Rcell (Bush et al., 2012), the updated version bundles, into a single software suite, the image-processing capabilities of Cell-ID, offers new data analysis tools for cytometry, and relies on the widely used data analysis and visualization tools of the statistical programming framework R. © 2023 Wiley Periodicals LLC. Basic Protocol: Extracting quantitative information from single cells Support Protocol 1: Obtaining and installing Cell-ID and R Support Protocol 2: Preparing cells for imaging.

Modulation of transcription factor dynamics allows versatile information transmission.

Cells detect changes in their environment and generate responses, often involving changes in gene expression. In this paper we use information theory and a simple transcription model to analyze whether the resulting gene expression serves to identify extracellular stimuli and assess their intensity when they are encoded in the amplitude, duration or frequency of pulses of a transcription factor's nuclear concentration (or activation state). We find, for all cases, that about three ranges of input strengths can be distinguished and that maximum information transmission occurs for fast and high activation threshold promoters. The three input modulation modes differ in the sensitivity to changes in the promoters parameters. Frequency modulation is the most sensitive and duration modulation, the least. This is key for signal identification: there are promoter parameters that yield a relatively high information transmission for duration or amplitude modulation and a much smaller value for frequency modulation. The reverse situation cannot be found with a single promoter transcription model. Thus, pulses of transcription factors can selectively activate the "frequency-tuned" promoter while prolonged nuclear accumulation would activate promoters of all three modes simultaneously. Frequency modulation is therefore highly selective and better suited than the other encoding modes for signal identification without requiring other mediators of the transduction process.

Characterization of cell-to-cell variation in nuclear transport rates and identification of its sources.

Nuclear transport is an essential part of eukaryotic cell function. Here, we present scFRAP, a model-assisted fluorescent recovery after photobleaching (FRAP)- based method to determine nuclear import and export rates independently in individual live cells. To overcome the inherent noise of single-cell measurements, we performed sequential FRAPs on the same cell. We found large cell-to-cell variation in transport rates within isogenic yeast populations. For passive transport, the variability in NPC number might explain most of the variability. Using this approach, we studied mother-daughter cell asymmetry in the active nuclear shuttling of the transcription factor Ace2, which is specifically concentrated in daughter cell nuclei in early G1. Rather than reduced export in the daughter cell, as previously hypothesized, we found that this asymmetry is mainly due to an increased import in daughters. These results shed light on cell-to-cell variation in cellular dynamics and its sources.

Akt Is S-Palmitoylated: A New Layer of Regulation for Akt.

The protein kinase Akt/PKB participates in a great variety of processes, including translation, cell proliferation and survival, as well as malignant transformation and viral infection. In the last few years, novel Akt posttranslational modifications have been found. However, how these modification patterns affect Akt subcellular localization, target specificity and, in general, function is not thoroughly understood. Here, we postulate and experimentally demonstrate by acyl-biotin exchange (ABE) assay and 3H-palmitate metabolic labeling that Akt is S-palmitoylated, a modification related to protein sorting throughout subcellular membranes. Mutating cysteine 344 into serine blocked Akt S-palmitoylation and diminished its phosphorylation at two key sites, T308 and T450. Particularly, we show that palmitoylation-deficient Akt increases its recruitment to cytoplasmic structures that colocalize with lysosomes, a process stimulated during autophagy. Finally, we found that cysteine 344 in Akt1 is important for proper its function, since Akt1-C344S was unable to support adipocyte cell differentiation in vitro. These results add an unexpected new layer to the already complex Akt molecular code, improving our understanding of cell decision-making mechanisms such as cell survival, differentiation and death.

A new tool to sense pH changes at the neuromuscular junction synaptic cleft.

Synaptic transmission triggers transient acidification of the synaptic cleft. Recently, it has been shown that pH affects the opening of postsynaptic channels and therefore the production of tools that allow to study these behaviors should result of paramount value. We fused α-bungarotoxin, a neurotoxin derived from the snake Bungarus multicinctus that binds irreversibly to the acetylcholine receptor extracellular domain, to the pH sensitive GFP Super Ecliptic pHluorin, and efficiently expressed it in Pichia pastoris. This sensor allows synaptic changes in pH to be measured without the need of incorporating transgenes into animal cells. Here, we show that incubation of the mouse levator auris muscle with a solution containing this recombinant protein is enough to fluorescently label the endplate post synaptic membrane. Furthermore, we could physiologically alter and measure post synaptic pH by evaluating changes in the fluorescent signal of pHluorin molecules bound to acetylcholine receptors. In fact, using this tool we were able to detect a drop in 0.01 to 0.05 pH units in the vicinity of the acetylcholine receptors following vesicle exocytosis triggered by nerve electrical stimulation. Further experiments will allow to learn the precise changes in pH during and after synaptic activation.

GPCR receptor phosphorylation and endocytosis are not necessary to switch polarized growth between internal cues during pheromone response in S. cerevisiae.

Chemotactic/chemotropic cells follow accurately the direction of gradients of regulatory molecules. Many G-protein-coupled receptors (GPCR) function as chemoattractant receptors to guide polarized responses. In "a" mating type yeast, the GPCR Ste2 senses the α-cell's pheromone. Previously, phosphorylation and trafficking of this receptor have been implicated in the process of gradient sensing, where cells dynamically correct growth. Correction is often necessary since yeast have intrinsic polarity sites that interfere with a correct initial gradient decoding. We have recently showed that when actively dividing (not in G1) yeast are exposed to a uniform pheromone concentration, they initiate a pheromone-induced polarization next to the mother-daughter cytokinesis site. Then, they reorient their growth to the intrinsic polarity site. Here, to study if Ste2 phosphorylation and internalization are involved in this process, we generated receptor variants combining three types of mutated signals for the first time: phosphorylation, ubiquitylation and the NPFX1,2D Sla1-binding motif. We first characterized their effect on endocytosis and found that these processes regulate internalization in a more complex manner than previously shown. Interestingly, we showed that receptor phosphorylation can drive internalization independently of ubiquitylation and the NPFX1,2D motif. When tested in our assays, cells expressing either phosphorylation or endocytosis-deficient receptors were able to switch away from the cytokinesis site to find the intrinsic polarity site as efficiently as their WT counterparts. Thus, we conclude that these processes are not necessary for the reorientation of polarization.

Mitotic and pheromone-specific intrinsic polarization cues interfere with gradient sensing in Saccharomyces cerevisiae.

Polarity decisions are central to many processes, including mitosis and chemotropism. In Saccharomyces cerevisiae, budding and mating projection (MP) formation use an overlapping system of cortical landmarks that converges on the small G protein Cdc42. However, pheromone-gradient sensing must override the Rsr1-dependent internal polarity cues used for budding. Using this model system, we asked what happens when intrinsic and extrinsic spatial cues are not aligned. Is there competition, or collaboration? By live-cell microscopy and microfluidics techniques, we uncovered three previously overlooked features of this signaling system. First, the cytokinesis-associated polarization patch serves as a polarity landmark independently of all known cues. Second, the Rax1-Rax2 complex functions as a pheromone-promoted polarity cue in the distal pole of the cells. Third, internal cues remain active during pheromone-gradient tracking and can interfere with this process, biasing the location of MPs. Yeast defective in internal-cue utilization align significantly better than wild type with artificially generated pheromone gradients.

Properties of cell signaling pathways and gene expression systems operating far from steady-state.

Ligand-receptor systems, covalent modification cycles, and transcriptional networks are basic units of signaling systems and their steady-state properties are well understood. However, the behavior of such systems before steady-state is poorly characterized. Here, we analyzed the properties of input-output curves for each of these systems as they approach steady-state. In ligand-receptor systems, the EC50 (concentration of the ligand that occupies 50% of the receptors) is higher before the system reaches steady-state. Based on this behavior, we have previously defined PRESS (for pre-equilibrium sensing and signaling), a general "systems level" mechanism cells may use to overcome input saturation. Originally, we showed that, given a step stimulation, PRESS operates when the kinetics of ligand-receptor binding are slower than the downstream signaling steps. Now, we show that, provided the input increases slowly, it is not essential for the ligand binding reaction itself to be slow. In addition, we demonstrate that covalent modification cycles and gene expression systems may also operate in PRESS mode. Thus, nearly all biochemical processes may operate in PRESS mode, suggesting that this mechanism may be ubiquitous in cell signaling systems.

Heat-stress triggers MAPK crosstalk to turn on the hyperosmotic response pathway.

Cells make decisions based on a combination of external and internal signals. In yeast, the high osmolarity response (HOG) is a mitogen-activated protein kinase (MAPK) pathway that responds to a variety of stimuli, and it is central to the general stress response. Here we studied the effect of heat-stress (HS) on HOG. Using live-cell reporters and genetics, we show that HS promotes Hog1 phosphorylation and Hog1-dependent gene expression, exclusively via the Sln1 phosphorelay branch, and that the strength of the activation is larger in yeast adapted to high external osmolarity. HS stimulation of HOG is indirect. First, we show that HS causes glycerol loss, necessary for HOG activation. Preventing glycerol efflux by deleting the glyceroporin FPS1 or its regulators RGC1 and ASK10/RGC2, or by increasing external glycerol, greatly reduced HOG activation. Second, we found that HOG stimulation by HS depended on the operation of a second MAPK pathway, the cell-wall integrity (CWI), a well-known mediator of HS, since inactivating Pkc1 or deleting the MAPK SLT2 greatly reduced HOG activation. Our data suggest that the main role of the CWI in this process is to stimulate glycerol loss. We found that in yeast expressing the constitutively open channel mutant (Fps1-Δ11), HOG activity was independent of Slt2. In summary, we suggest that HS causes a reduction in turgor due to the loss of glycerol and the accompanying water, and that this is what actually stimulates HOG. Thus, taken together, our findings highlight a central role for Fps1, and the metabolism of glycerol, in the communication between the yeast MAPK pathways, essential for survival and reproduction in changing environments.

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.

CDK and MAPK Synergistically Regulate Signaling Dynamics via a Shared Multi-site Phosphorylation Region on the Scaffold Protein Ste5.

We report an unanticipated system of joint regulation by cyclin-dependent kinase (CDK) and mitogen-activated protein kinase (MAPK), involving collaborative multi-site phosphorylation of a single substrate. In budding yeast, the protein Ste5 controls signaling through a G1 arrest pathway. Upon cell-cycle entry, CDK inhibits Ste5 via multiple phosphorylation sites, disrupting its membrane association. Using quantitative time-lapse microscopy, we examined Ste5 membrane recruitment dynamics at different cell-cycle stages. Surprisingly, in S phase, where Ste5 recruitment should be blocked, we observed an initial recruitment followed by a steep drop-off. This delayed inhibition revealed a requirement for both CDK activity and negative feedback from the pathway MAPK Fus3. Mutagenesis, mass spectrometry, and electrophoretic analyses suggest that the CDK and MAPK modify shared sites, which are most extensively phosphorylated when both kinases are active and able to bind their docking sites on Ste5. Such collaborative phosphorylation can broaden regulatory inputs and diversify output dynamics of signaling pathways.

Ultrasensitivity in signaling cascades revisited: Linking local and global ultrasensitivity estimations.

Ultrasensitive response motifs, capable of converting graded stimuli into binary responses, are well-conserved in signal transduction networks. Although it has been shown that a cascade arrangement of multiple ultrasensitive modules can enhance the system's ultrasensitivity, how a given combination of layers affects a cascade's ultrasensitivity remains an open question for the general case. Here, we introduce a methodology that allows us to determine the presence of sequestration effects and to quantify the relative contribution of each module to the overall cascade's ultrasensitivity. The proposed analysis framework provides a natural link between global and local ultrasensitivity descriptors and it is particularly well-suited to characterize and understand mathematical models used to study real biological systems. As a case study, we have considered three mathematical models introduced by O'Shaughnessy et al. to study a tunable synthetic MAPK cascade, and we show how our methodology can help modelers better understand alternative models.

Yeast GPCR signaling reflects the fraction of occupied receptors, not the number.

According to receptor theory, the effect of a ligand depends on the amount of agonist-receptor complex. Therefore, changes in receptor abundance should have quantitative effects. However, the response to pheromone in Saccharomyces cerevisiae is robust (unaltered) to increases or reductions in the abundance of the G-protein-coupled receptor (GPCR), Ste2, responding instead to the fraction of occupied receptor. We found experimentally that this robustness originates during G-protein activation. We developed a complete mathematical model of this step, which suggested the ability to compute fractional occupancy depends on the physical interaction between the inhibitory regulator of G-protein signaling (RGS), Sst2, and the receptor. Accordingly, replacing Sst2 by the heterologous hsRGS4, incapable of interacting with the receptor, abolished robustness. Conversely, forcing hsRGS4:Ste2 interaction restored robustness. Taken together with other results of our work, we conclude that this GPCR pathway computes fractional occupancy because ligand-bound GPCR-RGS complexes stimulate signaling while unoccupied complexes actively inhibit it. In eukaryotes, many RGSs bind to specific GPCRs, suggesting these complexes with opposing activities also detect fraction occupancy by a ratiometric measurement. Such complexes operate as push-pull devices, which we have recently described.

Push-Pull and Feedback Mechanisms Can Align Signaling System Outputs with Inputs.

Many cell signaling systems, including the yeast pheromone response system, exhibit "dose-response alignment" (DoRA), in which output of one or more downstream steps closely matches the fraction of occupied receptors. DoRA can improve the fidelity of transmitted dose information. Here, we searched systematically for biochemical network topologies that produced DoRA. Most networks, including many containing feedback and feedforward loops, could not produce DoRA. However, networks including "push-pull" mechanisms, in which the active form of a signaling species stimulates downstream activity and the nominally inactive form reduces downstream activity, enabled perfect DoRA. Networks containing feedbacks enabled DoRA, but only if they also compared feedback to input and adjusted output to match. Our results establish push-pull as a non-feedback mechanism to align output with variable input and maximize information transfer in signaling systems. They also suggest genetic approaches to determine whether particular signaling systems use feedback or push-pull control.

Nonlinear mixed-effects modelling for single cell estimation: when, why, and how to use it.

BACKGROUND: Studies of cell-to-cell variation have in recent years grown in interest, due to improved bioanalytical techniques which facilitates determination of small changes with high uncertainty. Like much high-quality data, single-cell data is best analysed using a systems biology approach. The most common systems biology approach to single-cell data is the standard two-stage (STS) approach. In STS, data from each cell is analysed in a separate sub-problem, meaning that only data from the same cell is used to calculate the parameter values within that cell. Because only parts of the data are considered, problems with parameter unidentifiability are exaggerated in STS. In contrast, a related approach to data analysis has been developed for the studies of patient-to-patient variations. This approach, called nonlinear mixed-effects modelling (NLME), makes use of all data, when estimating the patient-specific parameters. NLME would therefore be advantageous compared to STS also for the study of cell-to-cell variation. However, no such systematic evaluation of the two approaches exists. RESULTS: Herein, such a systematic comparison between STS and NLME has been performed. Different examples, both linear and nonlinear, and both simulated and real experimental data, have been examined. With informative data, there is no significant difference in the results for either parameter or noise estimation. However, when data becomes uninformative, NLME is significantly superior to STS. These results hold independently of whether the loss of information is due to a low signal-to-noise ratio, too few data points, or a bad input signal. The improvement is shown to come from both the consideration of a joint likelihood (JLH) function, describing all parameters and data, and from an a priori postulated form of the population parameters. Finally, we provide a small tutorial that shows how to use NLME for single-cell analysis, using the free and user-friendly software Monolix. CONCLUSIONS: When considering uninformative single-cell data, NLME yields more accurate parameter and noise estimates, compared to more traditional approaches, such as STS and JLH.

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.

Impact of upstream and downstream constraints on a signaling module's ultrasensitivity.

Much work has been done on the study of the biochemical mechanisms that result in ultrasensitive behavior of simple biochemical modules. However, in a living cell, such modules are embedded in a bigger network that constrains the range of inputs that the module will receive as well as the range of the module's outputs that network will be able to detect. Here, we studied how the effective ultrasensitivity of a modular system is affected by these restrictions. We use a simple setup to explore to what extent the dynamic range spanned by upstream and downstream components of an ultrasensitive module impact on the effective sensitivity of the system. Interestingly, we found for some ultrasensitive motifs that dynamic range limitations imposed by downstream components can produce effective sensitivities much larger than that of the original module when considered in isolation.

Utilization of extracellular information before ligand-receptor binding reaches equilibrium expands and shifts the input dynamic range.

Cell signaling systems sense and respond to ligands that bind cell surface receptors. These systems often respond to changes in the concentration of extracellular ligand more rapidly than the ligand equilibrates with its receptor. We demonstrate, by modeling and experiment, a general "systems level" mechanism cells use to take advantage of the information present in the early signal, before receptor binding reaches a new steady state. This mechanism, pre-equilibrium sensing and signaling (PRESS), operates in signaling systems in which the kinetics of ligand-receptor binding are slower than the downstream signaling steps, and it typically involves transient activation of a downstream step. In the systems where it operates, PRESS expands and shifts the input dynamic range, allowing cells to make different responses to ligand concentrations so high as to be otherwise indistinguishable. Specifically, we show that PRESS applies to the yeast directional polarization in response to pheromone gradients. Consideration of preexisting kinetic data for ligand-receptor interactions suggests that PRESS operates in many cell signaling systems throughout biology. The same mechanism may also operate at other levels in signaling systems in which a slow activation step couples to a faster downstream step.