How cells integrate complex stimuli: the effect of feedback from phosphoinositides and cell shape on cell polarization and motility

Microfluidic platform for omics analysis on single cells with diverse morphology and size: A review

BACKGROUND

Microfluidic techniques have emerged as powerful tools in single-cell research, facilitating the exploration of omics information from individual cells. Cell morphology is crucial for gene expression and physiological processes. However, there is currently a lack of integrated analysis of morphology and single-cell omics information. A critical challenge remains: what platform technologies are the best option to decode omics data of cells that are complex in morphology and size?

RESULTS

This review highlights achievements in microfluidic-based single-cell omics and isolation of cells based on morphology, along with other cell sorting methods based on physical characteristics. Various microfluidic platforms for single-cell isolation are systematically presented, showcasing their diversity and adaptability. The discussion focuses on microfluidic devices tailored to the distinct single-cell isolation requirements in plants and animals, emphasizing the significance of considering cell morphology and cell size in optimizing single-cell omics strategies. Simultaneously, it explores the application of microfluidic single-cell sorting technologies to single-cell sequencing, aiming to effectively integrate information about cell shape and size.

SIGNIFICANCE AND NOVELTY

The novelty lies in presenting a comprehensive overview of recent accomplishments in microfluidic-based single-cell omics, emphasizing the integration of different microfluidic platforms and their implications for cell morphology-based isolation. By underscoring the pivotal role of the specialized morphology of different cells in single-cell research, this review provides robust support for delving deeper into the exploration of single-cell omics data.

From actin waves to mechanism and back: How theory aids biological understanding

Actin dynamics in cell motility, division, and phagocytosis is regulated by complex factors with multiple feedback loops, often leading to emergent dynamic patterns in the form of propagating waves of actin polymerization activity that are poorly understood. Many in the actin wave community have attempted to discern the underlying mechanisms using experiments and/or mathematical models and theory. Here, we survey methods and hypotheses for actin waves based on signaling networks, mechano-chemical effects, and transport characteristics, with examples drawn from Dictyostelium discoideum, human neutrophils, Caenorhabditis elegans, and Xenopus laevis oocytes. While experimentalists focus on the details of molecular components, theorists pose a central question of universality: Are there generic, model-independent, underlying principles, or just boundless cell-specific details? We argue that mathematical methods are equally important for understanding the emergence, evolution, and persistence of actin waves and conclude with a few challenges for future studies.

Morphomigrational description as a new approach connecting cell's migration with its morphology

The examination of morphology and migration of cells plays substantial role in understanding the cellular behaviour, being described by plethora of quantitative parameters and models. These descriptions, however, treat cell migration and morphology as independent properties of temporal cell state, while not taking into account their strong interdependence in adherent cells. Here we present the new and simple mathematical parameter called signed morphomigrational angle (sMM angle) that links cell geometry with translocation of cell centroid, considering them as one morphomigrational behaviour. The sMM angle combined with pre-existing quantitative parameters enabled us to build a new tool called morphomigrational description, used to assign the numerical values to several cellular behaviours. Thus, the cellular activities that until now were characterized using verbal description or by complex mathematical models, are described here by a set of numbers. Our tool can be further used in automatic analysis of cell populations as well as in studies focused on cellular response to environmental directional signals.

Geometry-induced patterns through mechanochemical coupling

Intracellular protein patterns regulate a variety of vital cellular processes such as cell division and motility, which often involve dynamic cell-shape changes. These changes in cell shape may in turn affect the dynamics of pattern-forming proteins, hence leading to an intricate feedback loop between cell shape and chemical dynamics. While several computational studies have examined the rich resulting dynamics, the underlying mechanisms are not yet fully understood. To elucidate some of these mechanisms, we explore a conceptual model for cell polarity on a dynamic one-dimensional manifold. Using concepts from differential geometry, we derive the equations governing mass-conserving reaction-diffusion systems on time-evolving manifolds. Analyzing these equations mathematically, we show that dynamic shape changes of the membrane can induce pattern-forming instabilities in parts of the membrane, which we refer to as regional instabilities. Deformations of the local membrane geometry can also (regionally) suppress pattern formation and spatially shift already existing patterns. We explain our findings by applying and generalizing the local equilibria theory of mass-conserving reaction-diffusion systems. This allows us to determine a simple onset criterion for geometry-induced pattern-forming instabilities, which is linked to the phase-space structure of the reaction-diffusion system. The feedback loop between membrane shape deformations and reaction-diffusion dynamics then leads to a surprisingly rich phenomenology of patterns, including oscillations, traveling waves, and standing waves, even if these patterns do not occur in systems with a fixed membrane shape. Our paper reveals that the local conformation of the membrane geometry acts as an important dynamical control parameter for pattern formation in mass-conserving reaction-diffusion systems.

Simulating 3D Cell Shape with the Cellular Potts Model

Computer simulations have become a widely used method for the field of mechanobiology. An important question is whether one can predict the shape and forces of cells as a function of the extracellular environment. Different types of models have been described before to simulate cell and tissue shapes in structured environments. In this chapter, we give a brief overview of commonly used models and then describe the Cellular Potts Model, a lattice-based modelling framework, in more detail. We provide a hands-on guide on how to build a model that simulates the shape of a single cell on a micropattern in three dimensions in different open source software packages using the Cellular Potts framework. A simulation is set up with an initial configuration of generalized cells that change shape and position due to an energy function that incorporates cellular volume and surface area constraints as well as interaction energies between the generalized cells.

Sensing the shape of a cell with reaction diffusion and energy minimization

Some dividing cells sense their shape by becoming polarized along their long axis. Cell polarity is controlled in part by polarity proteins, like Rho GTPases, cycling between active membrane-bound forms and inactive cytosolic forms, modeled as a "wave-pinning" reaction-diffusion process. Does shape sensing emerge from wave pinning? We show that wave pinning senses the cell's long axis. Simulating wave pinning on a curved surface, we find that high-activity domains migrate to peaks and troughs of the surface. For smooth surfaces, a simple rule of minimizing the domain perimeter while keeping its area fixed predicts the final position of the domain and its shape. However, when we introduce roughness to our surfaces, shape sensing can be disrupted, and high-activity domains can become localized to locations other than the global peaks and valleys of the surface. On rough surfaces, the domains of the wave-pinning model are more robust in finding the peaks and troughs than the minimization rule, although both can become trapped in steady states away from the peaks and valleys. We can control the robustness of shape sensing by altering the Rho GTPase diffusivity and the domain size. We also find that the shape-sensing properties of cell polarity models can explain how domains localize to curved regions of deformed cells. Our results help to understand the factors that allow cells to sense their shape-and the limits that membrane roughness can place on this process.

Viscoelasticity, Like Forces, Plays a Role in Mechanotransduction

Viscoelasticity and its alteration in time and space has turned out to act as a key element in fundamental biological processes in living systems, such as morphogenesis and motility. Based on experimental and theoretical findings it can be proposed that viscoelasticity of cells, spheroids and tissues seems to be a collective characteristic that demands macromolecular, intracellular component and intercellular interactions. A major challenge is to couple the alterations in the macroscopic structural or material characteristics of cells, spheroids and tissues, such as cell and tissue phase transitions, to the microscopic interferences of their elements. Therefore, the biophysical technologies need to be improved, advanced and connected to classical biological assays. In this review, the viscoelastic nature of cytoskeletal, extracellular and cellular networks is presented and discussed. Viscoelasticity is conceptualized as a major contributor to cell migration and invasion and it is discussed whether it can serve as a biomarker for the cells’ migratory capacity in several biological contexts. It can be hypothesized that the statistical mechanics of intra- and extracellular networks may be applied in the future as a powerful tool to explore quantitatively the biomechanical foundation of viscoelasticity over a broad range of time and length scales. Finally, the importance of the cellular viscoelasticity is illustrated in identifying and characterizing multiple disorders, such as cancer, tissue injuries, acute or chronic inflammations or fibrotic diseases.

Two Complementary Signaling Pathways Depict Eukaryotic Chemotaxis: A Mechanochemical Coupling Model

Many eukaryotic cells, including neutrophils and Dictyostelium cells, are able to undergo correlated random migration in the absence of directional cues while reacting to shallow gradients of chemoattractants with exquisite precision. Although progress has been made with regard to molecular identities, it remains elusive how molecular mechanics are integrated with cell mechanics to initiate and manipulate cell motility. Here, we propose a two dimensional (2D) cell migration model wherein a multilayered dynamic seesaw mechanism is accompanied by a mechanical strain-based inhibition mechanism. In biology, these two mechanisms can be mapped onto the biochemical feedback between phosphoinositides (PIs) and Rho GTPase and the mechanical interplay between filamin A (FLNa) and FilGAP. Cell migration and the accompanying morphological changes are demonstrated in numerical simulations using a particle-spring model, and the diffusion in the cell membrane are simulations using a one dimensional (1D) finite differences method (FDM). The fine balance established between endogenous signaling and a mechanically governed inactivation scheme ensures the endogenous cycle of self-organizing pseudopods, accounting for the correlated random migration. Furthermore, this model cell manifests directional and adaptable responses to shallow graded signaling, depending on the overwhelming effect of the graded stimuli guidance on strain-based inhibition. Finally, the model cell becomes trapped within an obstacle-ridden spatial region, manifesting a shuttle run for local explorations and can chemotactically “escape”, illustrating again the balance required in the complementary signaling pathways.

Cellular Tango: how extracellular matrix adhesion choreographs Rac-Rho signaling and cell movement

The small GTPases Rac and Rho are known to regulate eukaryotic cell shape, promoting front protrusion (Rac) or rear retraction (Rho) of the cell edge. Such cell deformation changes the contact and adhesion of cell to the extracellular matrix (ECM), while ECM signaling through integrin receptors also affects GTPase activity. We develop and investigate a model for this three-way feedback loop in 1D and 2D spatial domains, as well as in a fully deforming 2D cell shapes with detailed adhesion-bond biophysics. The model consists of reaction-diffusion equations solved numerically with open-source software, Morpheus, and with custom-built cellular Potts model simulations. We find a variety of patterns and cell behaviors, including persistent polarity, flipped front-back cell polarity oscillations, spiral waves, and random protrusion-retraction. We show that the observed spatial patterns depend on the cell shape, and vice versa.

Exploring the influence of cytosolic and membrane FAK activation on YAP/TAZ nuclear translocation

Membrane binding and unbinding dynamics play a crucial role in the biological activity of several nonintegral membrane proteins, which have to be recruited to the membrane to perform their functions. By localizing to the membrane, these proteins are able to induce downstream signal amplification in their respective signaling pathways. Here, we present a 3D computational approach using reaction-diffusion equations to investigate the relation between membrane localization of focal adhesion kinase (FAK), Ras homolog family member A (RhoA), and signal amplification of the YAP/TAZ signaling pathway. Our results show that the theoretical scenarios in which FAK is membrane bound yield robust and amplified YAP/TAZ nuclear translocation signals. Moreover, we predict that the amount of YAP/TAZ nuclear translocation increases with cell spreading, confirming the experimental findings in the literature. In summary, our in silico predictions show that when the cell membrane interaction area with the underlying substrate increases, for example, through cell spreading, this leads to more encounters between membrane-bound signaling partners and downstream signal amplification. Because membrane activation is a motif common to many signaling pathways, this study has important implications for understanding the design principles of signaling networks.

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.

An optogenetic model reveals cell shape regulation through FAK and fascin

Cell shape regulation is important, but the mechanisms that govern shape are not fully understood, in part due to limited experimental models in which cell shape changes and underlying molecular processes can be rapidly and non-invasively monitored in real time. Here, we used an optogenetic tool to activate RhoA in the middle of mononucleated macrophages to induce contraction, resulting in a side with the nucleus that retained its shape and a non-nucleated side that was unable to maintain its shape and collapsed. In cells overexpressing focal adhesion kinase (FAK; also known as PTK2), the non-nucleated side exhibited a wide flat morphology and was similar in adhesion area to the nucleated side. In cells overexpressing fascin, an actin-bundling protein, the non-nucleated side assumed a spherical shape and was similar in height to the nucleated side. This effect of fascin was also observed in fibroblasts even without inducing furrow formation. Based on these results, we conclude that FAK and fascin work together to maintain cell shape by regulating adhesion area and height, respectively, in different cell types. This article has an associated First Person interview with the first author of the paper.

More than just a barrier: using physical models to couple membrane shape to cell function

The correct execution of many cellular processes, such as division and motility, requires the cell to adopt a specific shape. Physically, these shapes are determined by the interplay of the plasma membrane and internal cellular driving factors. While the plasma membrane defines the boundary of the cell, processes inside the cell can result in the generation of forces that deform the membrane. These processes include protein binding, the assembly of protein superstructures, and the growth and contraction of cytoskeletal networks. Due to the complexity of the cell, relating observed membrane deformations back to internal processes is a challenging problem. Here, we review cell shape changes in endocytosis, cell adhesion, cell migration and cell division and discuss how by modeling membrane deformations we can investigate the inner working principles of the cell.

Droplet-Based Membranous Soft Materials

Inspired by the structure and functionality of natural cellular tissues, droplet interface bilayer (DIB)-based materials strategically combine model membrane assembly techniques and droplet microfluidics. These structures have shown promising results in applications ranging from biological computing to chemical microrobots. This Feature Article briefly explores recent advances in the areas of construction, manipulation, and functionalization of DIB networks; discusses their unique mechanics; and focuses on the contributions of our lab in the advancement of this platform. We also reflect on some of the limitations facing DIB-based materials and how they might be addressed, highlighting promising applications made possible through the refinement of the material concept.

On the influence of cell shape on dynamic reaction-diffusion polarization patterns

The distribution of signaling molecules following mechanical or chemical stimulation of a cell defines cell polarization, with regions of high active Cdc42 at the front and low active Cdc42 at the rear. As reaction-diffusion phenomena between signaling molecules, such as Rho GTPases, define the gradient dynamics, we hypothesize that the cell shape influences the maintenance of the “front-to-back” cell polarization patterns. We investigated the influence of cell shape on the Cdc42 patterns using an established computational polarization model. Our simulation results showed that not only cell shape but also Cdc42 and Rho-related (in)activation parameter values affected the distribution of active Cdc42. Despite an initial Cdc42 gradient, the in silico results showed that the maximal Cdc42 concentration shifts in the opposite direction, a phenomenon we propose to call “reverse polarization”. Additional in silico analyses indicated that “reverse polarization” only occurred in a particular parameter value space that resulted in a balance between inactivation and activation of Rho GTPases. Future work should focus on a mathematical description of the underpinnings of reverse polarization, in combination with experimental validation using, for example, dedicated FRET-probes to spatiotemporally track Rho GTPase patterns in migrating cells. In summary, the findings of this study enhance our understanding of the role of cell shape in intracellular signaling.

Spots, stripes, and spiral waves in models for static and motile cells : GTPase patterns in cells

The polarization and motility of eukaryotic cells depends on assembly and contraction of the actin cytoskeleton and its regulation by proteins called GTPases. The activity of GTPases causes assembly of filamentous actin (by GTPases Cdc42, Rac), resulting in protrusion of the cell edge. Mathematical models for GTPase dynamics address the spontaneous formation of patterns and nonuniform spatial distributions of such proteins in the cell. Here we revisit the wave-pinning model for GTPase-induced cell polarization, together with a number of extensions proposed in the literature. These include introduction of sources and sinks of active and inactive GTPase (by the group of A. Champneys), and negative feedback from F-actin to GTPase activity. We discuss these extensions singly and in combination, in 1D, and 2D static domains. We then show how the patterns that form (spots, waves, and spirals) interact with cell boundaries to create a variety of interesting and dynamic cell shapes and motion.

Using reactive dissipative particle dynamics to understand local shape manipulation of polymer vesicles

Biological cells have long been of interest to researchers due to their capacity to actively control their shape. Accordingly, there is significant interest in generating simplified synthetic protocells that can alter their shape based on an externally or internally generated stimulus. To date, most progress has been made towards controlling the global shape of a protocell, whereas less is known about generating a local shape change. Here, we seek to better understand the possible mechanisms for producing local morphological changes in a popular protocell system, the block copolymer vesicle. Accordingly, we have combined Dissipative Particle Dynamics (DPD) and the Split Reactive Brownian Dynamics algorithm (SRBD) to produce a simulation tool that is capable of modeling the dynamics of self-assembled polymer structures as they undergo chemical reactions. Using this Reactive DPD or RDPD method, we investigate local morphological change driven by either the microinjection of a stimulus or an enzymatically-produced stimulus. We find that sub-vesicle-scale morphological change can be induced by either a solvent stimulus that swells the vesicle membrane, or by a reactant stimulus that alters the chemistry of the block polymer in the membrane corona. Notably, the latter method results in a more persistent local deformation than the former, which we attribute to the slower diffusion of polymer chains relative to the solvent. We quantify this deformation and show that it can be modulated by altering the interaction parameter of the parts of the polymer chain that are affected by the stimulus.

Bridging from single to collective cell migration: A review of models and links to experiments

Mathematical and computational models can assist in gaining an understanding of cell behavior at many levels of organization. Here, we review models in the literature that focus on eukaryotic cell motility at 3 size scales: intracellular signaling that regulates cell shape and movement, single cell motility, and collective cell behavior from a few cells to tissues. We survey recent literature to summarize distinct computational methods (phase-field, polygonal, Cellular Potts, and spherical cells). We discuss models that bridge between levels of organization, and describe levels of detail, both biochemical and geometric, included in the models. We also highlight links between models and experiments. We find that models that span the 3 levels are still in the minority.

Membrane Tension Can Enhance Adaptation to Maintain Polarity of Migrating Cells

Migratory cells are known to adapt to environments that contain wide-ranging levels of chemoattractant. Although biochemical models of adaptation have been previously proposed, here, we discuss a different mechanism based on mechanosensing, in which the interaction between biochemical signaling and cell tension facilitates adaptation. We describe and analyze a model of mechanochemical-based adaptation coupling a mechanics-based physical model of cell tension coupled with the wave-pinning reaction-diffusion model for Rac GTPase activity. The mathematical analysis of this model, simulations of a simplified one-dimensional cell geometry, and two-dimensional finite element simulations of deforming cells reveal that as a cell protrudes under the influence of high stimulation levels, tension-mediated inhibition of Rac signaling causes the cell to polarize even when initially overstimulated. Specifically, tension-mediated inhibition of Rac activation, which has been experimentally observed in recent years, facilitates this adaptation by countering the high levels of environmental stimulation. These results demonstrate how tension-related mechanosensing may provide an alternative (and potentially complementary) mechanism for cell adaptation.

A Stochastic Modelling Framework for Single Cell Migration: Coupling Contractility and Focal Adhesions

The interaction of the actin cytoskeleton with cell–substrate adhesions is necessary for cell migration. While the trajectories of motile cells have a stochastic character, investigations of cell motility mechanisms rarely elaborate on the origins of the observed randomness. Here, guided by a few fundamental attributes of cell motility, I construct a minimal stochastic cell migration model from ground-up. The resulting model couples a deterministic actomyosin contractility mechanism with stochastic cell–substrate adhesion kinetics, and yields a well-defined piecewise deterministic process. Numerical simulations reproduce several experimentally observed results, including anomalous diffusion, tactic migration and contact guidance. This work provides a basis for the development of cell–cell collision and population migration models.

Mechanical interplay between cell shape and actin cytoskeleton organization

We investigate the mechanical interplay between the spatial organization of the actin cytoskeleton and the shape of animal cells adhering on micropillar arrays. Using a combination of analytical work, computer simulations and in vitro experiments, we demonstrate that the orientation of the stress fibers strongly influences the geometry of the cell edge. In the presence of a uniformly aligned cytoskeleton, the cell edge can be well approximated by elliptical arcs, whose eccentricity reflects the degree of anisotropy of the cell's internal stresses. Upon modeling the actin cytoskeleton as a nematic liquid crystal, we further show that the geometry of the cell edge feeds back on the organization of the stress fibers by altering the length scale at which these are confined. This feedback mechanism is controlled by a dimensionless number, the anchoring number, representing the relative weight of surface-anchoring and bulk-aligning torques. Our model allows to predict both cellular shape and the internal structure of the actin cytoskeleton and is in good quantitative agreement with experiments on fibroblastoid (GDβ1, GDβ3) and epithelioid (GEβ1, GEβ3) cells.

Influence of myosin activity and mechanical impact on keratocyte polarization

In cell migration, polarization is the process by which a stationary cell breaks symmetry and initiates motion. Although a lot is known about the mechanisms involved in cell polarization, the role played by myosin contraction remains unclear. In addition, cell polarization by mechanical impact has received little attention. Here, we study the influence of myosin activity on cell polarization and the initiation of motion induced by mechanical cues using a computational model for keratocytes. The model accounts for cell deformation, the dynamics of myosin and the signaling protein RhoA (a member of the Rho GTPases family), as well as the forces acting on the actomyosin network. Our results show that the attainment of a steady polarized state depends on the strength of myosin down- or up-regulation and that myosin upregulation favors cell polarization. Our results also confirm the existence of a threshold level for cell polarization, which is determined by the level of polarization of the Rho GTPases at the time the external stimuli vanish. In all, this paper shows that capturing the interactions between the signaling proteins (Rho GTPases for keratocytes) and the compounds of the motile machinery in a moving cell is crucial to study cell polarization.

CompuCell3D Simulations Reproduce Mesenchymal Cell Migration on Flat Substrates

Mesenchymal cell crawling is a critical process in normal development, in tissue function, and in many diseases. Quantitatively predictive numerical simulations of cell crawling thus have multiple scientific, medical, and technological applications. However, we still lack a low-computational-cost approach to simulate mesenchymal three-dimensional (3D) cell crawling. Here, we develop a computationally tractable 3D model (implemented as a simulation in the CompuCell3D simulation environment) of mesenchymal cells crawling on a two-dimensional substrate. The Fürth equation, the usual characterization of mean-squared displacement (MSD) curves for migrating cells, describes a motion in which, for increasing time intervals, cell movement transitions from a ballistic to a diffusive regime. Recent experiments have shown that for very short time intervals, cells exhibit an additional fast diffusive regime. Our simulations' MSD curves reproduce the three experimentally observed temporal regimes, with fast diffusion for short time intervals, slow diffusion for long time intervals, and intermediate time -interval-ballistic motion. The resulting parameterization of the trajectories for both experiments and simulations allows the definition of time- and length scales that translate between computational and laboratory units. Rescaling by these scales allows direct quantitative comparisons among MSD curves and between velocity autocorrelation functions from experiments and simulations. Although our simulations replicate experimentally observed spontaneous symmetry breaking, short-timescale diffusive motion, and spontaneous cell-motion reorientation, their computational cost is low, allowing their use in multiscale virtual-tissue simulations. Comparisons between experimental and simulated cell motion support the hypothesis that short-time actomyosin dynamics affects longer-time cell motility. The success of the base cell-migration simulation model suggests its future application in more complex situations, including chemotaxis, migration through complex 3D matrices, and collective cell motion.

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.

Phospholipids of the Plasma Membrane - Regulators or Consequence of Cell Polarity?

Cell polarity is a key feature of many eukaryotic cells, including neurons, epithelia, endothelia and asymmetrically dividing stem cells. Apart from the specific localization of proteins to distinct domains of the plasma membrane, most of these cells exhibit an asymmetric distribution of phospholipids within the plasma membrane too. Notably, research over the last years has revealed that many known conserved regulators of apical-basal polarity in epithelial cells are capable of binding to phospholipids, which in turn regulate the localization and to some extent the function of these proteins. Conversely, phospholipid-modifying enzymes are recruited and controlled by polarity regulators, demonstrating an elaborated balance between asymmetrically localized proteins and phospholipids, which are enriched in certain (micro)domains of the plasma membrane. In this review, we will focus on our current understanding of apical-basal polarity and the implication of phospholipids within the plasma membrane during the cell polarization of epithelia and migrating cells.

Quasi-periodic migration of single cells on short microlanes

Cell migration on microlanes represents a suitable and simple platform for the exploration of the molecular mechanisms underlying cell cytoskeleton dynamics. Here, we report on the quasi-periodic movement of cells confined in stripe-shaped microlanes. We observe persistent polarized cell shapes and directed pole-to-pole motion within the microlanes. Cells depolarize at one end of a given microlane, followed by delayed repolarization towards the opposite end. We analyze cell motility via the spatial velocity distribution, the velocity frequency spectrum and the reversal time as a measure for depolarization and spontaneous repolarization of cells at the microlane ends. The frequent encounters of a boundary in the stripe geometry provides a robust framework for quantitative investigations of the cytoskeleton protrusion and repolarization dynamics. In a first advance to rigorously test physical models of cell migration, we find that the statistics of the cell migration is recapitulated by a Cellular Potts model with a minimal description of cytoskeleton dynamics. Using LifeAct-GFP transfected cells and microlanes with differently shaped ends, we show that the local deformation of the leading cell edge in response to the tip geometry can locally either amplify or quench actin polymerization, while leaving the average reversal times unaffected.

Topotaxis of active Brownian particles

Recent experimental studies have demonstrated that cellular motion can be directed by topographical gradients, such as those resulting from spatial variations in the features of a micropatterned substrate. This phenomenon, known as topotaxis, has been extensively studied for topographical gradients at the subcellular scale, but can also occur in the presence of a spatially varying density of cell-sized features. Such a large-scale topotaxis has recently been observed in highly motile cells that persistently crawl within an array of obstacles with smoothly varying lattice spacing. We introduce a toy model of large-scale topotaxis, based on active Brownian particles. Using numerical simulations and analytical arguments, we demonstrate that topographical gradients introduce a spatial modulation of the particles' persistence, leading to directed motion toward regions of higher persistence. Our results demonstrate that persistent motion alone is sufficient to drive large-scale topotaxis and could serve as a starting point for more detailed studies on self-propelled particles and cells.

Cell Size, Mechanical Tension, and GTPase Signaling in the Single Cell

Cell polarization requires redistribution of specific proteins to the nascent front and back of a eukaryotic cell. Among these proteins are Rac and Rho, members of the small GTPase family that regulate the actin cytoskeleton. Rac promotes actin assembly and protrusion of the front edge, whereas Rho activates myosin-driven contraction at the back. Mathematical models of cell polarization at many levels of detail have appeared. One of the simplest based on "wave-pinning" consists of a pair of reaction-diffusion equations for a single GTPase. Mathematical analysis of wave-pinning so far is largely restricted to static domains in one spatial dimension. Here, we extend the analysis to cells that change in size, showing that both shrinking and growing cells can lose polarity. We further consider the feedback between mechanical tension, GTPase activation, and cell deformation in both static, growing, shrinking, and moving cells. Special cases (spatially uniform cell chemistry, the absence or presence of mechanical feedback) are analyzed, and the full model is explored by simulations in 1D. We find a variety of novel behaviors, including "dilution-induced" oscillations of Rac activity and cell size, as well as gain or loss of polarization and motility in the model cell.

From energy to cellular forces in the Cellular Potts Model: An algorithmic approach

Single and collective cell dynamics, cell shape changes, and cell migration can be conveniently represented by the Cellular Potts Model, a computational platform based on minimization of a Hamiltonian. Using the fact that a force field is easily derived from a scalar energy (F = −∇H), we develop a simple algorithm to associate effective forces with cell shapes in the CPM. We predict the traction forces exerted by single cells of various shapes and sizes on a 2D substrate. While CPM forces are specified directly from the Hamiltonian on the cell perimeter, we approximate the force field inside the cell domain using interpolation, and refine the results with smoothing. Predicted forces compare favorably with experimentally measured cellular traction forces. We show that a CPM model with internal signaling (such as Rho-GTPase-related contractility) can be associated with retraction-protrusion forces that accompany cell shape changes and migration. We adapt the computations to multicellular systems, showing, for example, the forces that a pair of swirling cells exert on one another, demonstrating that our algorithm works equally well for interacting cells. Finally, we show forces exerted by cells on one another in classic cell-sorting experiments.

Cells exert forces on their surroundings and on one another. In simulations of cell shape using the Cellular Potts Model (CPM), the dynamics of deforming cell shapes is traditionally represented by an energy-minimization method. We use this CPM energy, the Hamiltonian, to derive and visualize the corresponding forces exerted by the cells. We use the fact that force is the negative gradient of energy to assign forces to the CPM cell edges, and then extend the results to approximate interior forces by interpolation. We show that this method works for single as well as multiple interacting model cells, both static and motile. Finally, we show favorable comparison between predicted forces and real forces measured experimentally.

Bridging the gap between single-cell migration and collective dynamics

Motivated by the wealth of experimental data recently available, we present a cellular-automaton-based modeling framework focussing on high-level cell functions and their concerted effect on cellular migration patterns. Specifically, we formulate a coarse-grained description of cell polarity through self-regulated actin organization and its response to mechanical cues. Furthermore, we address the impact of cell adhesion on collective migration in cell cohorts. The model faithfully reproduces typical cell shapes and movements down to the level of single cells, yet allows for the efficient simulation of confluent tissues. In confined circular geometries, we find that specific properties of individual cells (polarizability; contractility) influence the emerging collective motion of small cell cohorts. Finally, we study the properties of expanding cellular monolayers (front morphology; stress and velocity distributions) at the level of extended tissues.

Experimental and computational analyses reveal that environmental restrictions shape HIV-1 spread in 3D cultures

Pathogens face varying microenvironments in vivo, but suitable experimental systems and analysis tools to dissect how three-dimensional (3D) tissue environments impact pathogen spread are lacking. Here we develop an Integrative method to Study Pathogen spread by Experiment and Computation within Tissue-like 3D cultures (INSPECT-3D), combining quantification of pathogen replication with imaging to study single-cell and cell population dynamics. We apply INSPECT-3D to analyze HIV-1 spread between primary human CD4 T-lymphocytes using collagen as tissue-like 3D-scaffold. Measurements of virus replication, infectivity, diffusion, cellular motility and interactions are combined by mathematical analyses into an integrated spatial infection model to estimate parameters governing HIV-1 spread. This reveals that environmental restrictions limit infection by cell-free virions but promote cell-associated HIV-1 transmission. Experimental validation identifies cell motility and density as essential determinants of efficacy and mode of HIV-1 spread in 3D. INSPECT-3D represents an adaptable method for quantitative time-resolved analyses of 3D pathogen spread.

Here, using an integrative experimental and computational approach, Imle et al. show how cell motility and density affect HIV cell-associated transmission in a three-dimensional tissue-like culture system of CD4⁺ T cells and collagen, and how different collagen matrices restrict infection by cell-free virions.

Guided by curvature: shaping cells by coupling curved membrane proteins and cytoskeletal forces

Eukaryote cells have flexible membranes that allow them to have a variety of dynamical shapes. The shapes of the cells serve important biological functions, both for cells within an intact tissue, and during embryogenesis and cellular motility. How cells control their shapes and the structures that they form on their surface has been a subject of intensive biological research, exposing the building blocks that cells use to deform their membranes. These processes have also drawn the interest of theoretical physicists, aiming to develop models based on physics, chemistry and nonlinear dynamics. Such models explore quantitatively different possible mechanisms that the cells can employ to initiate the spontaneous formation of shapes and patterns on their membranes. We review here theoretical work where one such class of mechanisms was investigated: the coupling between curved membrane proteins, and the cytoskeletal forces that they recruit. Theory indicates that this coupling gives rise to a rich variety of membrane shapes and dynamics, while experiments indicate that this mechanism appears to drive many cellular shape changes.This article is part of the theme issue 'Self-organization in cell biology'.

Small GTPase patterning: How to stabilise cluster coexistence

Many biological processes have to occur at specific locations on the cell membrane. These locations are often specified by the localised activity of small GTPase proteins. Some processes require the formation of a single cluster of active GTPase, also called unipolar polarisation (here “polarisation”), whereas others need multiple coexisting clusters. Moreover, sometimes the pattern of GTPase clusters is dynamically regulated after its formation. This raises the question how the same interacting protein components can produce such a rich variety of naturally occurring patterns. Most currently used models for GTPase-based patterning inherently yield polarisation. Such models may at best yield transient coexistence of at most a few clusters, and hence fail to explain several important biological phenomena. These existing models are all based on mass conservation of total GTPase and some form of direct or indirect positive feedback. Here, we show that either of two biologically plausible modifications can yield stable coexistence: including explicit GTPase turnover, i.e., breaking mass conservation, or negative feedback by activation of an inhibitor like a GAP. Since we start from two different polarising models our findings seem independent of the precise self-activation mechanism. By studying the net GTPase flows among clusters, we provide insight into how these mechanisms operate. Our coexistence models also allow for dynamical regulation of the final pattern, which we illustrate with examples of pollen tube growth and the branching of fungal hyphae. Together, these results provide a better understanding of how cells can tune a single system to generate a wide variety of biologically relevant patterns.

Epithelial vertex models with active biochemical regulation of contractility can explain organized collective cell motility

Collective motions of groups of cells are observed in many biological settings such as embryo development, tissue formation, and cancer metastasis. To effectively model collective cell movement, it is important to incorporate cell specific features such as cell size, cell shape, and cell mechanics, as well as active behavior of cells such as protrusion and force generation, contractile forces, and active biochemical signaling mechanisms that regulate cell behavior. In this paper, we develop a comprehensive model of collective cell migration in confluent epithelia based on the vertex modeling approach. We develop a method to compute cell-cell viscous friction based on the vertex model and incorporate RhoGTPase regulation of cortical myosin contraction. Global features of collective cell migration are examined by computing the spatial velocity correlation function. As active cell force parameters are varied, we found rich dynamical behavior. Furthermore, we find that cells exhibit nonlinear phenomena such as contractile waves and vortex formation. Together our work highlights the importance of active behavior of cells in generating collective cell movement. The vertex modeling approach is an efficient and versatile approach to rigorously examine cell motion in the epithelium.

Mechanochemical modeling of neutrophil migration based on four signaling layers, integrin dynamics, and substrate stiffness

Directional neutrophil migration during human immune responses is a highly coordinated process regulated by both biochemical and biomechanical environments. In this paper, we developed an integrative mathematical model of neutrophil migration using a lattice Boltzmann-particle method built in-house to solve the moving boundary problem with spatiotemporal regulation of biochemical components. The mechanical features of the cell cortex are modeled by a series of spring-connected nodes representing discrete cell-substrate adhesive sites. The intracellular signaling cascades responsible for cytoskeletal remodeling [e.g., small GTPases, phosphoinositide-3-kinase (PI3K), and phosphatase and tensin homolog] are built based on our previous four-layered signaling model centered on the bidirectional molecular transport mechanism and implemented as reaction-diffusion equations. Focal adhesion dynamics are determined by force-dependent integrin-ligand binding kinetics and integrin recycling and are thus integrated with cell motion. Using numerical simulations, the model reproduces the major features of cell migration in response to uniform and gradient biochemical stimuli based on the quantitative spatiotemporal regulation of signaling molecules, which agree with experimental observations. The existence of multiple types of integrins with different binding kinetics could act as an adaptation mechanism for substrate stiffness. Moreover, cells can perform reversal, U-turn, or lock-on behaviors depending on the steepness of the reversal biochemical signals received. Finally, this model is also applied to predict the responses of mutants in which PTEN is overexpressed or disrupted.

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.

Dynamic seesaw model for rapid signaling responses in eukaryotic chemotaxis

Directed movement of eukaryotic cells toward spatiotemporally varied chemotactic stimuli enables rapid intracellular signaling responses. While macroscopic cellular manifestation is shaped by balancing external stimuli strength with finite internal delays, the organizing principles of the underlying molecular mechanisms remain to be clarified. Here, we developed a novel modeling framework based on a simple seesaw mechanism to elucidate how cells repeatedly reverse polarity. As a key feature of the modeling, the bottom module of bidirectional molecular transport is successively controlled by three upstream modules of signal reception, initial signal processing, and Rho GTPase regulation. Our simulations indicated that an isotropic cell is polarized in response to a graded input signal. By applying a reversal gradient to a chemoattractant signal, lamellipod-specific molecules (i.e. PIP₃ and PI3K) disappear, first from the cell front, and then they redistribute at the opposite side, whereas functional molecules at the rear of the cell (i.e. PIP₂ and PTEN) act oppositely. In particular, the model cell exhibits a seesaw-like spatiotemporal pattern for the establishment of front and rear and interconversion, consistent with those related experimental observations. Increasing the switching frequency of the chemotactic gradient causes the cell to stay in a trapped state, further supporting the proposed dynamics of eukaryotic chemotaxis with the underlying cytoskeletal remodeling.

A PKC-MARCKS-PI3K regulatory module links Ca2+ and PIP3 signals at the leading edge of polarized macrophages

The leukocyte chemosensory pathway detects attractant gradients and directs cell migration to sites of inflammation, infection, tissue damage, and carcinogenesis. Previous studies have revealed that local Ca²⁺ and PIP₃ signals at the leading edge of polarized leukocytes play central roles in positive feedback loop essential to cell polarization and chemotaxis. These prior studies showed that stimulation of the leading edge Ca²⁺ signal can strongly activate PI3K, thereby triggering a larger PIP₃ signal, but did not elucidate the mechanistic link between Ca²⁺ and PIP₃ signaling. A hypothesis explaining this link emerged, postulating that Ca²⁺-activated PKC displaces the MARCKS protein from plasma membrane PIP₂, thereby releasing sequestered PIP₂ to serve as the target and substrate lipid of PI3K in PIP₃ production. In vitro single molecule studies of the reconstituted pathway on lipid bilayers demonstrated the feasibility of this PKC-MARCKS-PI3K regulatory module linking Ca²⁺ and PIP₃ signals in the reconstituted system. The present study tests the model predictions in live macrophages by quantifying the effects of: (a) two pathway activators—PDGF and ATP that stimulate chemoreceptors and Ca²⁺ influx, respectively; and (b) three pathway inhibitors—wortmannin, EGTA, and Go6976 that inhibit PI3K, Ca²⁺ influx, and PKC, respectively; on (c) four leading edge activity sensors—AKT-PH-mRFP, CKAR, MARCKSp-mRFP, and leading edge area that report on PIP₃ density, PKC activity, MARCKS membrane binding, and leading edge expansion/contraction, respectively. The results provide additional evidence that PKC and PI3K are both essential elements of the leading edge positive feedback loop, and strongly support the existence of a PKC-MARCKS-PI3K regulatory module linking the leading edge Ca²⁺ and PIP₃ signals. As predicted, activators stimulate leading edge PKC activity, displacement of MARCKS from the leading edge membrane and increased leading edge PIP₃ levels, while inhibitors trigger the opposite effects. Comparison of the findings for the ameboid chemotaxis of leukocytes with recently published findings for the mesenchymal chemotaxis of fibroblasts suggests that some features of the emerging leukocyte leading edge core pathway (PLC-DAG-Ca²⁺-PKC-MARCKS-PIP₂-PI3K-PIP₃) may well be shared by all chemotaxing eukaryotic cells, while other elements of the leukocyte pathway may be specialized features of these highly optimized, professional gradient-seeking cells. More broadly, the findings suggest a molecular mechanism for the strong links between phospho-MARCKS and many human cancers.

Computational Modeling of the Dynamics of Spatiotemporal Rho GTPase Signaling: A Systematic Review

The Rho family of GTPases are known to play pivotal roles in the regulation of fundamental cellular processes, ranging from cell migration and polarity to wound healing and regulation of actin cytoskeleton. Over the past decades, accumulating experimental work has increasingly mapped out the mechanistic details and interactions between members of the family and their regulators, establishing detailed interaction circuits within the Rho GTPase signaling network. These circuits have served as a vital foundation based on which a multitude of mathematical models have been developed to explain experimental data, gain deeper insights into the biological phenomenon they describe, as well as make new testable predictions and hypotheses. Due to the diverse nature and purpose of these models, they often vary greatly in size, scope, complexity, and formulation. Here, we provide a systematic, categorical, and comprehensive account of the recent modeling studies of Rho family GTPases, with an aim to offer a broad perspective of the field. The modeling limitations and possible future research directions are also discussed.

Mechanisms of Cell Polarization

Cell polarization is a key step in the migration, development, and organization of eukaryotic cells, both at the single cell and multicellular level. Research on the mechanisms that give rise to polarization of a given cell, and organization of polarity within a tissue has led to new understanding across cellular and developmental biology. In this review, we describe some of the history of theoretical and experimental aspects of the field, as well as some interesting questions and challenges for the future.

Crawling and turning in a minimal reaction-diffusion cell motility model: Coupling cell shape and biochemistry

We study a minimal model of a crawling eukaryotic cell with a chemical polarity controlled by a reaction-diffusion mechanism describing Rho GTPase dynamics. The size, shape, and speed of the cell emerge from the combination of the chemical polarity, which controls the locations where actin polymerization occurs, and the physical properties of the cell, including its membrane tension. We find in our model both highly persistent trajectories, in which the cell crawls in a straight line, and turning trajectories, where the cell transitions from crawling in a line to crawling in a circle. We discuss the controlling variables for this turning instability and argue that turning arises from a coupling between the reaction-diffusion mechanism and the shape of the cell. This emphasizes the surprising features that can arise from simple links between cell mechanics and biochemistry. Our results suggest that similar instabilities may be present in a broad class of biochemical descriptions of cell polarity.

Free boundary problem for cell protrusion formations: theoretical and numerical aspects

In this paper, a free boundary problem for cell protrusion formation is studied theoretically and numerically. The cell membrane is precisely described thanks to a level set function, whose motion is due to specific signalling pathways. The aim is to model the chemical interactions between the cell and its environment, in the process of invadopodia or pseudopodia formation. The model consists of Laplace equation with Dirichlet condition inside the cell coupled to Laplace equation with Neumann condition in the outer domain. The actin polymerization is accounted for as the gradient of the inner signal, which drives the motion of the interface. We prove the well-posedness of our free boundary problem under a sign condition on the datum. This criterion ensures the consistency of the model, and provides conditions to focus on for any enrichment of the model. We then propose a new first order Cartesian finite-difference method to solve the problem. We eventually exhibit the main biological features that can be accounted for by the model: the formation of thin and elongated protrusions as for invadopodia, or larger protrusion as for pseudopodia, depending on the source term in the equation. The model provides the theoretical and numerical grounds for single cell migration modeling, whose formulation is valid in 2D and 3D. In particular, specific chemical reactions that occurred at the cell membrane could be precisely described in forthcoming works.

Regulation of a Coupled MARCKS-PI3K Lipid Kinase Circuit by Calmodulin: Single-Molecule Analysis of a Membrane-Bound Signaling Module

Amoeboid cells that employ chemotaxis to travel up an attractant gradient possess a signaling network assembled on the leading edge of the plasma membrane that senses the gradient and remodels the actin mesh and cell membrane to drive movement in the appropriate direction. In leukocytes such as macrophages and neutrophils, and perhaps in other amoeboid cells as well, the leading edge network includes a positive feedback loop in which the signaling of multiple pathway components is cooperatively coupled. Cytoplasmic Ca²⁺ is a recently recognized component of the feedback loop at the leading edge where it stimulates phosphoinositide-3-kinase (PI3K) and the production of its product signaling lipid phosphatidylinositol 3,4,5-trisphosphate (PIP₃). A previous study implicated Ca²⁺-activated protein kinase C (PKC) and the phosphatidylinositol 4,5-bisphosphate (PIP₂) binding protein MARCKS as two important players in this signaling, because PKC phosphorylation of MARCKS releases free PIP₂ that serves as the membrane binding target and substrate for PI3K. This study asks whether calmodulin (CaM), which is known to directly bind MARCKS, also stimulates PIP₃ production by releasing free PIP₂. Single-molecule fluorescence microscopy is used to quantify the surface density and enzyme activity of key protein components of the hypothesized Ca²⁺-CaM-MARCKS-PIP₂-PI3K-PIP₃ circuit. The findings show that CaM does stimulate PI3K lipid kinase activity by binding MARCKS and displacing it from PIP₂ headgroups, thereby releasing free PIP₂ that recruits active PI3K to the membrane and serves as the substrate for the generation of PIP₃. The resulting CaM-triggered activation of PI3K is complete in seconds and is much faster than PKC-triggered activation, which takes minutes. Overall, the available evidence implicates both PKC and CaM in the coupling of Ca²⁺ and PIP₃ signals and suggests these two different pathways have slow and fast activation kinetics, respectively.

Modeling cell shape and dynamics on micropatterns

Adhesive micropatterns have become a standard tool to study cells under defined conditions. Applications range from controlling the differentiation and fate of single cells to guiding the collective migration of cell sheets. In long-term experiments, single cell normalization is challenged by cell division. For all of these setups, mathematical models predicting cell shape and dynamics can guide pattern design. Here we review recent advances in predicting and explaining cell shape, traction forces and dynamics on micropatterns. Starting with contour models as the simplest approach to explain concave cell shapes, we move on to network and continuum descriptions as examples for static models. To describe dynamic processes, cellular Potts, vertex and phase field models can be used. Different types of model are appropriate to address different biological questions and together, they provide a versatile tool box to predict cell behavior on micropatterns.

Rationalizing Rac1 and RhoA GTPase signaling: A mathematical approach

Precise spatiotemporal dynamics of Rho GTPases are essential for efficient cell migration. Manipulating Rac1 and RhoA signaling is thus a potential intervention strategy to abrogate harmful cell invasion and subsequent metastasis; however GTPase signaling can be extremely complicated due to crosstalk and the multitude of upstream regulators and downstream effectors. Studying Rho GTPase networks in a formal mathematical setting can therefore be of great use. We recently built a predictive model based on Boolean logic which identified a negative feedback loop critical for RhoA and Rac1 activity. Here, we discuss the value and potential pitfalls of different mathematical approaches which have been used to study Rho GTPase dynamics, and highlight the importance of choosing the correct approach given the data available and outputs desired. Overall, a mathematical approach, particularly when combined iteratively with in vitro experiments, can be of great use in deriving new biological insight to further harness the activity of Rho GTPases.

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.

One of the riddles of Nature is how cells interact with one another to create complex cellular networks such as the neural networks in the brain. Forming precise connections between irregularly shaped cells is a challenge for biology. We developed computational methods for simulating these complex cell-cell interactions. We applied these methods to investigate yeast mating in which two yeast cells grow projections that meet and fuse guided by pheromone attractants. The simulations described molecules both inside and outside of the cell, and represented the continually changing shapes of the cells. We found that positioning the secretion and sensing of pheromones at the same location on the cell surface was important. Other key factors for robust mating included secreting a protein that removed excess pheromone from outside of the cell so that the signal would not be too strong. An important advance was being able to simulate as many as five cells in complex mating arrangements. Taken together we used our novel computational methods to describe in greater detail the yeast mating process, and more generally, interactions among cells changing their shapes in response to their neighbors.

Effects of 3D geometries on cellular gradient sensing and polarization

During cell migration, cells become polarized, change their shape, and move in response to various internal and external cues. Cell polarization is defined through the spatio-temporal organization of molecules such as PI3K or small GTPases, and is determined by intracellular signaling networks. It results in directional forces through actin polymerization and myosin contractions. Many existing mathematical models of cell polarization are formulated in terms of reaction-diffusion systems of interacting molecules, and are often defined in one or two spatial dimensions. In this paper, we introduce a 3D reaction-diffusion model of interacting molecules in a single cell, and find that cell geometry has an important role affecting the capability of a cell to polarize, or change polarization when an external signal changes direction. Our results suggest a geometrical argument why more roundish cells can repolarize more effectively than cells which are elongated along the direction of the original stimulus, and thus enable roundish cells to turn faster, as has been observed in experiments. On the other hand, elongated cells preferentially polarize along their main axis even when a gradient stimulus appears from another direction. Furthermore, our 3D model can accurately capture the effect of binding and unbinding of important regulators of cell polarization to and from the cell membrane. This spatial separation of membrane and cytosol, not possible to capture in 1D or 2D models, leads to marked differences of our model from comparable lower-dimensional models.