A novel stochastic simulation approach enables exploration of mechanisms for regulating polarity site movement

Particle-based simulations reveal two positive feedback loops allow relocation and stabilization of the polarity site during yeast mating

Many cells adjust the direction of polarized growth or migration in response to external directional cues. The yeast Saccharomyces cerevisiae orient their cell fronts (also called polarity sites) up pheromone gradients in the course of mating. However, the initial polarity site is often not oriented towards the eventual mating partner, and cells relocate the polarity site in an indecisive manner before developing a stable orientation. During this reorientation phase, the polarity site displays erratic assembly-disassembly behavior and moves around the cell cortex. The mechanisms underlying this dynamic behavior remain poorly understood. Particle-based simulations of the core polarity circuit revealed that molecular-level fluctuations are unlikely to overcome the strong positive feedback required for polarization and generate relocating polarity sites. Surprisingly, inclusion of a second pathway that promotes polarity site orientation generated relocating polarity sites with properties similar to those observed experimentally. This pathway forms a second positive feedback loop involving the recruitment of receptors to the cell membrane and couples polarity establishment to gradient sensing. This second positive feedback loop also allows cells to stabilize their polarity site once the site is aligned with the pheromone gradient.

Close agreement between deterministic versus stochastic modeling of first-passage time to vesicle fusion

Ca2+-dependent cell processes, such as neurotransmitter or endocrine vesicle fusion, are inherently stochastic due to large fluctuations in Ca2+ channel gating, Ca2+ diffusion, and Ca2+ binding to buffers and target sensors. However, previous studies revealed closer-than-expected agreement between deterministic and stochastic simulations of Ca2+ diffusion, buffering, and sensing if Ca2+ channel gating is not Ca2+ dependent. To understand this result more fully, we present a comparative study complementing previous work, focusing on Ca2+ dynamics downstream of Ca2+ channel gating. Specifically, we compare deterministic (mean-field/mass-action) and stochastic simulations of vesicle exocytosis latency, quantified by the probability density of the first-passage time (FPT) to the Ca2+-bound state of a vesicle fusion sensor, following a brief Ca2+ current pulse. We show that under physiological constraints, the discrepancy between FPT densities obtained using the two approaches remains small even if as few as ∼50 Ca2+ ions enter per single channel-vesicle release unit. Using a reduced two-compartment model for ease of analysis, we illustrate how this close agreement arises from the smallness of correlations between fluctuations of the reactant molecule numbers, despite the large magnitude of fluctuation amplitudes. This holds if all relevant reactions are heteroreaction between molecules of different species, as is the case for bimolecular Ca2+ binding to buffers and downstream sensor targets. In this case, diffusion and buffering effectively decorrelate the state of the Ca2+ sensor from local Ca2+ fluctuations. Thus, fluctuations in the Ca2+ sensor's state underlying the FPT distribution are only weakly affected by the fluctuations in the local Ca2+ concentration around its average, deterministically computable value.

Chemotropism and Cell-Cell Fusion in Fungi

Fungi exhibit an enormous variety of morphologies, including yeast colonies, hyphal mycelia, and elaborate fruiting bodies. This diversity arises through a combination of polar growth, cell division, and cell fusion. Because fungal cells are nonmotile and surrounded by a protective cell wall that is essential for cell integrity, potential fusion partners must grow toward each other until they touch and then degrade the intervening cell walls without impacting cell integrity. Here, we review recent progress on understanding how fungi overcome these challenges. Extracellular chemoattractants, including small peptide pheromones, mediate communication between potential fusion partners, promoting the local activation of core cell polarity regulators to orient polar growth and cell wall degradation. However, in crowded environments, pheromone gradients can be complex and potentially confusing, raising the question of how cells can effectively find their partners. Recent findings suggest that the cell polarity circuit exhibits searching behavior that can respond to pheromone cues through a remarkably flexible and effective strategy called exploratory polarization.

Orientation of Cell Polarity by Chemical Gradients

Accurate decoding of spatial chemical landscapes is critical for many cell functions. Eukaryotic cells decode local chemical gradients to orient growth or movement in productive directions. Recent work on yeast model systems, whose gradient sensing pathways display much less complexity than those in animal cells, has suggested new paradigms for how these very small cells successfully exploit information in noisy and dynamic pheromone gradients to identify their mates. Pheromone receptors regulate a polarity circuit centered on the conserved Rho-family GTPase, Cdc42. The polarity circuit contains both positive and negative feedback pathways, allowing spontaneous symmetry breaking and also polarity site disassembly and relocation. Cdc42 orients the actin cytoskeleton, leading to focused vesicle traffic that promotes movement of the polarity site and also reshapes the cortical distribution of receptors at the cell surface. In this article, we review the advances from work on yeasts and compare them with the excitable signaling pathways that have been revealed in chemotactic animal cells. Expected final online publication date for the Annual Review of Biophysics, Volume 51 is May 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Three-dimensional stochastic simulation of chemoattractant-mediated excitability in cells

During the last decade, a consensus has emerged that the stochastic triggering of an excitable system drives pseudopod formation and subsequent migration of amoeboid cells. The presence of chemoattractant stimuli alters the threshold for triggering this activity and can bias the direction of migration. Though noise plays an important role in these behaviors, mathematical models have typically ignored its origin and merely introduced it as an external signal into a series of reaction-diffusion equations. Here we consider a more realistic description based on a reaction-diffusion master equation formalism to implement these networks. In this scheme, noise arises naturally from a stochastic description of the various reaction and diffusion terms. Working on a three-dimensional geometry in which separate compartments are divided into a tetrahedral mesh, we implement a modular description of the system, consisting of G-protein coupled receptor signaling (GPCR), a local excitation-global inhibition mechanism (LEGI), and signal transduction excitable network (STEN). Our models implement detailed biochemical descriptions whenever this information is available, such as in the GPCR and G-protein interactions. In contrast, where the biochemical entities are less certain, such as the LEGI mechanism, we consider various possible schemes and highlight the differences between them. Our simulations show that even when the LEGI mechanism displays perfect adaptation in terms of the mean level of proteins, the variance shows a dose-dependence. This differs between the various models considered, suggesting a possible means for determining experimentally among the various potential networks. Overall, our simulations recreate temporal and spatial patterns observed experimentally in both wild-type and perturbed cells, providing further evidence for the excitable system paradigm. Moreover, because of the overall importance and ubiquity of the modules we consider, including GPCR signaling and adaptation, our results will be of interest beyond the field of directed migration.

Though the term noise usually carries negative connotations, it can also contribute positively to the characteristic dynamics of a system. In biological systems, where noise arises from the stochastic interactions between molecules, its study is usually confined to genetic regulatory systems in which copy numbers are small and fluctuations large. However, noise can have important roles when the number of signaling molecules is large. The extension of pseudopods and the subsequent motion of amoeboid cells arises from the noise-induced trigger of an excitable system. Chemoattractant signals bias this triggering thereby directing cell motion. To date, this paradigm has not been tested by mathematical models that account accurately for the noise that arises in the corresponding reactions. In this study, we employ a reaction-diffusion master equation approach to investigate the effects of noise. Using a modular approach and a three-dimensional cell model with specific subdomains attributed to the cell membrane and cortex, we explore the spatiotemporal dynamics of the system. Our simulations recreate many experimentally-observed cell behaviors thereby supporting the biased-excitable network hypothesis.