Epidemiological modeling in StochSS Live!

SUMMARY: We present StochSS Live!, a web-based service for modeling, simulation and analysis of a wide range of mathematical, biological and biochemical systems. Using an epidemiological model of COVID-19, we demonstrate the power of StochSS Live! to enable researchers to quickly develop a deterministic or a discrete stochastic model, infer its parameters and analyze the results. AVAILABILITY AND IMPLEMENTATION: StochSS Live! is freely available at https://live.stochss.org/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Coordinating cell polarization and morphogenesis through mechanical feedback.

Many cellular processes require cell polarization to be maintained as the cell changes shape, grows or moves. Without feedback mechanisms relaying information about cell shape to the polarity molecular machinery, the coordination between cell polarization and morphogenesis, movement or growth would not be possible. Here we theoretically and computationally study the role of a genetically-encoded mechanical feedback (in the Cell Wall Integrity pathway) as a potential coordination mechanism between cell morphogenesis and polarity during budding yeast mating projection growth. We developed a coarse-grained continuum description of the coupled dynamics of cell polarization and morphogenesis as well as 3D stochastic simulations of the molecular polarization machinery in the evolving cell shape. Both theoretical approaches show that in the absence of mechanical feedback (or in the presence of weak feedback), cell polarity cannot be maintained at the projection tip during growth, with the polarization cap wandering off the projection tip, arresting morphogenesis. In contrast, for mechanical feedback strengths above a threshold, cells can robustly maintain cell polarization at the tip and simultaneously sustain mating projection growth. These results indicate that the mechanical feedback encoded in the Cell Wall Integrity pathway can provide important positional information to the molecular machinery in the cell, thereby enabling the coordination of cell polarization and morphogenesis.

An arbitrary Lagrangian Eulerian smoothed particle hydrodynamics (ALE-SPH) method with a boundary volume fraction formulation for fluid-structure interaction.

We present a new weakly-compressible smoothed particle hydrodynamics (SPH) method capable of modeling non-slip fixed and moving wall boundary conditions. The formulation combines a boundary volume fraction (BVF) wall approach with the transport-velocity SPH method. The resulting method, named SPH-BVF, offers detection of arbitrarily shaped solid walls on-the-fly, with small computational overhead due to its local formulation. This simple framework is capable of solving problems that are difficult or infeasible for standard SPH, namely flows subject to large shear stresses or at moderate Reynolds numbers, and mass transfer in deformable boundaries. In addition, the method extends the transport-velocity formulation to reaction-diffusion transport of mass in Newtonian fluids and linear elastic solids, which is common in biological structures. Taken together, the SPH-BVF method provides a good balance of simplicity and versatility, while avoiding some of the standard obstacles associated with SPH: particle penetration at the boundaries, tension instabilities and anisotropic particle alignments, that hamper SPH from being applied to complex problems such as fluid-structure interaction in a biological system.

Parameter uncertainty quantification using surrogate models applied to a spatial model of yeast mating polarization.

A common challenge in systems biology is quantifying the effects of unknown parameters and estimating parameter values from data. For many systems, this task is computationally intractable due to expensive model evaluations and large numbers of parameters. In this work, we investigate a new method for performing sensitivity analysis and parameter estimation of complex biological models using techniques from uncertainty quantification. The primary advance is a significant improvement in computational efficiency from the replacement of model simulation by evaluation of a polynomial surrogate model. We demonstrate the method on two models of mating in budding yeast: a smaller ODE model of the heterotrimeric G-protein cycle, and a larger spatial model of pheromone-induced cell polarization. A small number of model simulations are used to fit the polynomial surrogates, which are then used to calculate global parameter sensitivities. The surrogate models also allow rapid Bayesian inference of the parameters via Markov chain Monte Carlo (MCMC) by eliminating model simulations at each step. Application to the ODE model shows results consistent with published single-point estimates for the model and data, with the added benefit of calculating the correlations between pairs of parameters. On the larger PDE model, the surrogate models allowed convergence for the distribution of 15 parameters, which otherwise would have been computationally prohibitive using simulations at each MCMC step. We inferred parameter distributions that in certain cases peaked at values different from published values, and showed that a wide range of parameters would permit polarization in the model. Strikingly our results suggested different diffusion constants for active versus inactive Cdc42 to achieve good polarization, which is consistent with experimental observations in another yeast species S. pombe.

Mechanical feedback coordinates cell wall expansion and assembly in yeast mating morphogenesis.

The shaping of individual cells requires a tight coordination of cell mechanics and growth. However, it is unclear how information about the mechanical state of the wall is relayed to the molecular processes building it, thereby enabling the coordination of cell wall expansion and assembly during morphogenesis. Combining theoretical and experimental approaches, we show that a mechanical feedback coordinating cell wall assembly and expansion is essential to sustain mating projection growth in budding yeast (Saccharomyces cerevisiae). Our theoretical results indicate that the mechanical feedback provided by the Cell Wall Integrity pathway, with cell wall stress sensors Wsc1 and Mid2 increasingly activating membrane-localized cell wall synthases Fks1/2 upon faster cell wall expansion, stabilizes mating projection growth without affecting cell shape. Experimental perturbation of the osmotic pressure and cell wall mechanics, as well as compromising the mechanical feedback through genetic deletion of the stress sensors, leads to cellular phenotypes that support the theoretical predictions. Our results indicate that while the existence of mechanical feedback is essential to stabilize mating projection growth, the shape and size of the cell are insensitive to the feedback.

The effect of cell geometry on polarization in budding yeast.

The localization (or polarization) of proteins on the membrane during the mating of budding yeast (Saccharomyces cerevisiae) is an important model system for understanding simple pattern formation within cells. While there are many existing mathematical models of polarization, for both budding and mating, there are still many aspects of this process that are not well understood. In this paper we set out to elucidate the effect that the geometry of the cell can have on the dynamics of certain models of polarization. Specifically, we look at several spatial stochastic models of Cdc42 polarization that have been adapted from published models, on a variety of tip-shaped geometries, to replicate the shape change that occurs during the growth of the mating projection. We show here that there is a complex interplay between the dynamics of polarization and the shape of the cell. Our results show that while models of polarization can generate a stable polarization cap, its localization at the tip of mating projections is unstable, with the polarization cap drifting away from the tip of the projection in a geometry dependent manner. We also compare predictions from our computational results to experiments that observe cells with projections of varying lengths, and track the stability of the polarization cap. Lastly, we examine one model of actin polarization and show that it is unlikely, at least for the models studied here, that actin dynamics and vesicle traffic are able to overcome this effect of geometry.

Modelling of Yeast Mating Reveals Robustness Strategies for Cell-Cell Interactions.

Mating of budding yeast cells is a model system for studying cell-cell interactions. Haploid yeast cells secrete mating pheromones that are sensed by the partner which responds by growing a mating projection toward the source. The two projections meet and fuse to form the diploid. Successful mating relies on precise coordination of dynamic extracellular signals, signaling pathways, and cell shape changes in a noisy background. It remains elusive how cells mate accurately and efficiently in a natural multi-cell environment. Here we present the first stochastic model of multiple mating cells whose morphologies are driven by pheromone gradients and intracellular signals. Our novel computational framework encompassed a moving boundary method for modeling both a-cells and α-cells and their cell shape changes, the extracellular diffusion of mating pheromones dynamically coupled with cell polarization, and both external and internal noise. Quantification of mating efficiency was developed and tested for different model parameters. Computer simulations revealed important robustness strategies for mating in the presence of noise. These strategies included the polarized secretion of pheromone, the presence of the α-factor protease Bar1, and the regulation of sensing sensitivity; all were consistent with data in the literature. In addition, we investigated mating discrimination, the ability of an a-cell to distinguish between α-cells either making or not making α-factor, and mating competition, in which multiple a-cells compete to mate with one α-cell. Our simulations were consistent with previous experimental results. Moreover, we performed a combination of simulations and experiments to estimate the diffusion rate of the pheromone a-factor. In summary, we constructed a framework for simulating yeast mating with multiple cells in a noisy environment, and used this framework to reproduce mating behaviors and to identify strategies for robust cell-cell interactions.

A framework for discrete stochastic simulation on 3D moving boundary domains.

We have developed a method for modeling spatial stochastic biochemical reactions in complex, three-dimensional, and time-dependent domains using the reaction-diffusion master equation formalism. In particular, we look to address the fully coupled problems that arise in systems biology where the shape and mechanical properties of a cell are determined by the state of the biochemistry and vice versa. To validate our method and characterize the error involved, we compare our results for a carefully constructed test problem to those of a microscale implementation. We demonstrate the effectiveness of our method by simulating a model of polarization and shmoo formation during the mating of yeast. The method is generally applicable to problems in systems biology where biochemistry and mechanics are coupled, and spatial stochastic effects are critical.

Spatial stochastic dynamics enable robust cell polarization.

Although cell polarity is an essential feature of living cells, it is far from being well-understood. Using a combination of computational modeling and biological experiments we closely examine an important prototype of cell polarity: the pheromone-induced formation of the yeast polarisome. Focusing on the role of noise and spatial heterogeneity, we develop and investigate two mechanistic spatial models of polarisome formation, one deterministic and the other stochastic, and compare the contrasting predictions of these two models against experimental phenotypes of wild-type and mutant cells. We find that the stochastic model can more robustly reproduce two fundamental characteristics observed in wild-type cells: a highly polarized phenotype via a mechanism that we refer to as spatial stochastic amplification, and the ability of the polarisome to track a moving pheromone input. Moreover, we find that only the stochastic model can simultaneously reproduce these characteristics of the wild-type phenotype and the multi-polarisome phenotype of a deletion mutant of the scaffolding protein Spa2. Significantly, our analysis also demonstrates that higher levels of stochastic noise results in increased robustness of polarization to parameter variation. Furthermore, our work suggests a novel role for a polarisome protein in the stabilization of actin cables. These findings elucidate the intricate role of spatial stochastic effects in cell polarity, giving support to a cellular model where noise and spatial heterogeneity combine to achieve robust biological function.

A positive feedback loop involving the Spa2 SHD domain contributes to focal polarization.

Focal polarization is necessary for finely arranged cell-cell interactions. The yeast mating projection, with its punctate polarisome, is a good model system for this process. We explored the critical role of the polarisome scaffold protein Spa2 during yeast mating with a hypothesis motivated by mathematical modeling and tested by in vivo experiments. Our simulations predicted that two positive feedback loops generate focal polarization, including a novel feedback pathway involving the N-terminal domain of Spa2. We characterized the latter using loss-of-function and gain-of-function mutants. The N-terminal region contains a Spa2 Homology Domain (SHD) which is conserved from yeast to humans, and when mutated largely reproduced the spa2Δ phenotype. Our work together with published data show that the SHD domain recruits Msb3/4 that stimulates Sec4-mediated transport of Bud6 to the polarisome. There, Bud6 activates Bni1-catalyzed actin cable formation, recruiting more Spa2 and completing the positive feedback loop. We demonstrate that disrupting this loop at any point results in morphological defects. Gain-of-function perturbations partially restored focal polarization in a spa2 loss-of-function mutant without restoring localization of upstream components, thus supporting the pathway order. Thus, we have collected data consistent with a novel positive feedback loop that contributes to focal polarization during pheromone-induced polarization in yeast.

The effects of replacing Sst2 with the heterologous RGS4 on polarization and mating in yeast.

RGS proteins stimulate the deactivation of heterotrimeric G-proteins. The yeast RGS protein Sst2 is regulated at both the transcriptional and posttranscriptional levels. We replaced the SST2 gene with the distantly related human RGS4 gene, which consists of the catalytic domain and an N-terminal membrane attachment peptide, and replaced the native promoter (P(SST2)) with the heterologous tetracycline-repressible promoter (P(TET)). We then measured the effect of the substitutions on pheromone sensitivity, mating, and polarization. Although the pheromone sensitivity was essentially normal, there were differences in mating and polarization. In particular, the RGS4-substituted strains did not form multiple mating projections at high levels of alpha-factor, but instead formed a single malformed projection, which frequently gave rise to a bud. We provide evidence that this phenotype arose because unlike Sst2, RGS4 did not localize to the projection. We use mathematical modeling to argue that localization of Sst2 to the projection prevents excess G-protein activation during the pheromone response. In addition, modeling and experiments demonstrate that the dose of Sst2 influences the frequency of mating projection formation.

Validation data for a hybrid smoothed dissipative particle dynamics (SDPD) spatial stochastic simulation algorithm (sSSA) method.

We present the validation of the hybrid sSSA-SDPD method for advection-diffusion-reaction problems coupled to discrete biochemical systems, as presented in the publication "A hybrid smoothed dissipative particle dynamics (SDPD) spatial stochastic simulation algorithm (sSSA) for advection-diffusion-reaction problems" (Drawert et al., 2019). We validate 1D diffusion, and 2D diffusion cases against their analytical solutions. We present graphs and tables of data showing the error in the simulation method.

Control analysis of bacterial chemotaxis signaling.

Bacteria such as Escherichia coli demonstrate the remarkable ability to migrate up gradients of attractants and down gradients of repellents in a rapid and sensitive fashion. They employ a temporal sensing strategy in which they estimate the concentration of ligand at different time points and continue moving in the same direction if the concentration is increasing in time, and randomly reorient if the concentration is decreasing in time. The key to success is accurate sensing of ligand levels in the presence of extracellular and intracellular disturbances. Research from a control theory perspective has begun to characterize the robustness of the bacterial chemotaxis signal transduction system to these perturbations. Modeling and theory can describe the optimal performance of such a sensor and how it can be achieved, thereby illuminating the design of the network. This chapter describes some basic principles of control theory relevant to the analysis of this sensing system, including sensitivity analysis, Bode plots, integral feedback control, and noise filters (i.e., Kalman filters).

Optimal noise filtering in the chemotactic response of Escherichia coli.

Information-carrying signals in the real world are often obscured by noise. A challenge for any system is to filter the signal from the corrupting noise. This task is particularly acute for the signal transduction network that mediates bacterial chemotaxis, because the signals are subtle, the noise arising from stochastic fluctuations is substantial, and the system is effectively acting as a differentiator which amplifies noise. Here, we investigated the filtering properties of this biological system. Through simulation, we first show that the cutoff frequency has a dramatic effect on the chemotactic efficiency of the cell. Then, using a mathematical model to describe the signal, noise, and system, we formulated and solved an optimal filtering problem to determine the cutoff frequency that bests separates the low-frequency signal from the high-frequency noise. There was good agreement between the theory, simulations, and published experimental data. Finally, we propose that an elegant implementation of the optimal filter in combination with a differentiator can be achieved via an integral control system. This paper furnishes a simple quantitative framework for interpreting many of the key notions about bacterial chemotaxis, and, more generally, it highlights the constraints on biological systems imposed by noise.

Some protein interaction data do not exhibit power law statistics.

It has been claimed that protein-protein interaction (PPI) networks are scale-free, and that identifying high-degree "hub" proteins reveals important features of PPI networks. In this paper, we evaluate the claims that PPI node degree sequences follow a power law, a necessary condition for networks to be scale-free. We provide two PPI network examples which clearly do not have power laws when analyzed correctly, and thus at least these PPI networks are not scale-free. We also show that these PPI networks do appear to have power laws according to methods that have become standard in the existing literature. We explain the source of this error using numerically generated data from analytic formulas, where there are no sampling or noise ambiguities.

Alternative cell polarity behaviours arise from changes in G-protein spatial dynamics.

Yeast cells form a single mating projection when exposed to mating pheromone, a classic example of cell polarity. Prolonged treatment with pheromone or specific mutations results in alternative cell polarity behaviours. The authors performed mathematical modelling to investigate these unusual cell morphologies from the perspective of balancing spatial amplification (i.e. positive feedback that localises components) with spatial tracking (i.e. negative feedback that allows sensing of gradient). First, they used generic models of cell polarity to explore different cell polarity behaviours that arose from changes in the model spatial dynamics. By exploring the positive and negative feedback loops in each stage of a two-stage model, they simulated a variety of cell morphologies including single bending projections, single straight projections, periodic multiple projections and simultaneous double projections. In the second half of the study, they used a two-stage mechanistic model of yeast cell polarity focusing on G-protein signalling to integrate the modelling results more closely with the authors' previously published experimental observations. In summary, the combination of modelling and experiments describes how yeast cells exhibit a diversity of cell morphologies arising from two-stage G-protein signalling dynamics modulated by positive and negative feedbacks.

Yeast G-proteins mediate directional sensing and polarization behaviors in response to changes in pheromone gradient direction.

Yeast cells polarize by projecting up mating pheromone gradients, a classic cell polarity behavior. However, these chemical gradients may shift direction. We examine how yeast cells sense and respond to a 180(o) switch in the direction of microfluidically generated pheromone gradients. We identify two behaviors: at low concentrations of α-factor, the initial projection grows by bending, whereas at high concentrations, cells form a second projection toward the new source. Mutations that increase heterotrimeric G-protein activity expand the bending-growth morphology to high concentrations; mutations that increase Cdc42 activity result in second projections at low concentrations. Gradient-sensing projection bending requires interaction between Gßγ and Cdc24, whereas gradient-nonsensing projection extension is stimulated by Bem1 and hyperactivated Cdc42. Of interest, a mutation in Gα affects both bending and extension. Finally, we find a genetic perturbation that exhibits both behaviors. Overexpression of the formin Bni1, a component of the polarisome, makes both bending-growth projections and second projections at low and high α-factor concentrations, suggesting a role for Bni1 downstream of the heterotrimeric G-protein and Cdc42 during gradient sensing and response. Thus we demonstrate that G-proteins modulate in a ligand-dependent manner two fundamental cell-polarity behaviors in response to gradient directional change.

Signaling regulated endocytosis and exocytosis lead to mating pheromone concentration dependent morphologies in yeast.

Polarized cell morphogenesis requires actin cytoskeleton rearrangement for polarized transport of proteins, organelles and secretory vesicles, which fundamentally underlies cell differentiation and behavior. During yeast mating, Saccharomyces cerevisiae responds to extracellular pheromone gradients by extending polarized projections, which are likely maintained through vesicle transport to (exocytosis) and from (endocytosis) the membrane. We experimentally demonstrate that the projection morphology is pheromone concentration-dependent, and propose the underlying mechanism through mathematical modeling. The inclusion of membrane flux and dynamically evolving cell boundary into our yeast mating signaling model shows good agreement with experimental measurements, and provides a plausible explanation for pheromone-induced cell morphology.

Mathematical analysis of steady-state solutions in compartment and continuum models of cell polarization.

Cell polarization, in which substances previously uniformly distributed become asymmetric due to external or/and internal stimulation, is a fundamental process underlying cell mobility, cell division, and other polarized functions. The yeast cell S. cerevisiae has been a model system to study cell polarization. During mating, yeast cells sense shallow external spatial gradients and respond by creating steeper internal gradients of protein aligned with the external cue. The complex spatial dynamics during yeast mating polarization consists of positive feedback, degradation, global negative feedback control, and cooperative effects in protein synthesis. Understanding such complex regulations and interactions is critical to studying many important characteristics in cell polarization including signal amplification, tracking dynamic signals, and potential trade-off between achieving both objectives in a robust fashion. In this paper, we study some of these questions by analyzing several models with different spatial complexity: two compartments, three compartments, and continuum in space. The step-wise approach allows detailed characterization of properties of the steady state of the system, providing more insights for biological regulations during cell polarization. For cases without membrane diffusion, our study reveals that increasing the number of spatial compartments results in an increase in the number of steady-state solutions, in particular, the number of stable steady-state solutions, with the continuum models possessing infinitely many steady-state solutions. Through both analysis and simulations, we find that stronger positive feedback, reduced diffusion, and a shallower ligand gradient all result in more steady-state solutions, although most of these are not optimally aligned with the gradient. We explore in the different settings the relationship between the number of steady-state solutions and the extent and accuracy of the polarization. Taken together these results furnish a detailed description of the factors that influence the tradeoff between a single correctly aligned but poorly polarized stable steady-state solution versus multiple more highly polarized stable steady-state solutions that may be incorrectly aligned with the external gradient.

Noise filtering tradeoffs in spatial gradient sensing and cell polarization response.

BACKGROUND: Cells sense chemical spatial gradients and respond by polarizing internal components. This process can be disrupted by gradient noise caused by fluctuations in chemical concentration. RESULTS: We investigated how external gradient noise affects spatial sensing and response focusing on noise-filtering and the resultant tradeoffs. First, using a coarse-grained mathematical model of gradient-sensing and cell polarity, we characterized three negative consequences of noise: Inhibition of the extent of polarization, degradation of directional accuracy, and production of a noisy output polarization. Next, we explored filtering strategies and discovered that a combination of positive feedback, multiple signaling stages, and time-averaging produced good results. There was an important tradeoff, however, because filtering resulted in slower polarization. Simulations demonstrated that a two-stage filter-amplifier resulted in a balanced outcome. Then, we analyzed the effect of noise on a mechanistic model of yeast cell polarization in response to gradients of mating pheromone. This analysis showed that yeast cells likely also combine the above three filtering mechanisms into a filter-amplifier structure to achieve impressive spatial-noise tolerance, but with the consequence of a slow response time. Further investigation of the amplifier architecture revealed two positive feedback loops, a fast inner and a slow outer, both of which contributed to noise-tolerant polarization. This model also made specific predictions about how orientation performance depended upon the ratio between the gradient slope (signal) and the noise variance. To test these predictions, we performed microfluidics experiments measuring the ability of yeast cells to orient to shallow gradients of mating pheromone. The results of these experiments agreed well with the modeling predictions, demonstrating that yeast cells can sense gradients shallower than 0.1% µm-1, approximately a single receptor-ligand molecule difference between front and back, on par with motile eukaryotic cells. CONCLUSIONS: Spatial noise impedes the extent, accuracy, and smoothness of cell polarization. A combined filtering strategy implemented by a filter-amplifier architecture with slow dynamics was effective. Modeling and experimental data suggest that yeast cells employ these elaborate mechanisms to filter gradient noise resulting in a slow but relatively accurate polarization response.