A model for intracellular actin waves explored by nonlinear local perturbation analysis

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.

Versatile Patterns in the Actin Cortex of Motile Cells: Self-Organized Pulses Can Coexist with Macropinocytic Ring-Shaped Waves

Self-organized patterns in the actin cytoskeleton are essential for eukaryotic cellular life. They are the building blocks of many functional structures that often operate simultaneously to facilitate, for example, nutrient uptake and movement of cells. However, identifying how qualitatively distinct actin patterns can coexist remains a challenge. Using bifurcation theory of a mass conserved activator-inhibitor system, we uncover a generic mechanism of how different actin waves-traveling waves and excitable pulses-organize and simultaneously emerge. Live-cell imaging experiments indeed reveal that narrow, planar, and fast-moving excitable pulses may coexist with ring-shaped macropinocytic actin waves in the cortex of motile amoeboid cells.

Biophysical Models of PAR Cluster Transport by Cortical Flow in C. elegans Early Embryogenesis

The clustering of membrane-bound proteins facilitates their transport by cortical actin flow in early Caenorhabditis elegans embryo cell polarity. PAR-3 clustering is critical for this process, yet the biophysical processes that couple protein clusters to cortical flow remain unknown. We develop a discrete, stochastic agent-based model of protein clustering and test four hypothetical models for how clusters may interact with the flow. Results show that the canonical way to assess transport characteristics from single-particle tracking data used thus far in this area, the Péclet number, is insufficient to distinguish these hypotheses and that all models can account for transport characteristics quantified by this measure. However, using this model, we demonstrate that these different cluster-cortex interactions may be distinguished using a different metric, namely the scalar projection of cluster displacement on to the flow displacement vector. Our results thus provide a testable way to use existing single-particle tracking data to test how endogenous protein clusters may interact with the cortical flow to localize during polarity establishment. To facilitate this investigation, we also develop both improved simulation and semi-analytic methodologies to quantify motion summary statistics (e.g., Péclet number and scalar projection) for these stochastic models as a function of biophysical parameters.

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.

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.

Why a Large-Scale Mode Can Be Essential for Understanding Intracellular Actin Waves

During the last decade, intracellular actin waves have attracted much attention due to their essential role in various cellular functions, ranging from motility to cytokinesis. Experimental methods have advanced significantly and can capture the dynamics of actin waves over a large range of spatio-temporal scales. However, the corresponding coarse-grained theory mostly avoids the full complexity of this multi-scale phenomenon. In this perspective, we focus on a minimal continuum model of activator-inhibitor type and highlight the qualitative role of mass conservation, which is typically overlooked. Specifically, our interest is to connect between the mathematical mechanisms of pattern formation in the presence of a large-scale mode, due to mass conservation, and distinct behaviors of actin waves.

Tuning Cell Motility via Cell Tension with a Mechanochemical Cell Migration Model

Cell migration is orchestrated by a complicated mechanochemical system. However, few cell migration models take into account the coupling between the biochemical network and mechanical factors. Here, we construct a mechanochemical cell migration model to study the cell tension effect on cell migration. Our model incorporates the interactions between Rac-GTP, Rac-GDP, F-actin, myosin, and cell tension, and it is very convenient in capturing the change of cell shape by taking the phase field approach. This model captures the characteristic features of cell polarization, cell shape change, and cell migration modes. It shows that cell tension inhibits migration ability monotonically when cells are applied with persistent external stimuli. On the other hand, if random internal noise is significant, the regulation of cell tension exerts a nonmonotonic effect on cell migration. Because the increase of cell tension hinders the formation of multiple protrusions, migration ability could be maximized at intermediate cell tension under random internal noise. These model predictions are consistent with our single-cell experiments and other experimental results.

Simple Rho GTPase Dynamics Generate a Complex Regulatory Landscape Associated with Cell Shape

Migratory cells exhibit a variety of morphologically distinct responses to their environments that manifest in their cell shape. Some protrude uniformly to increase substrate contacts, others are broadly contractile, some polarize to facilitate migration, and yet others exhibit mixtures of these responses. Prior studies have identified a discrete collection of shapes that the majority of cells display and demonstrated that activity levels of the cytoskeletal regulators Rac1 and RhoA GTPase regulate those shapes. Here, we use computational modeling to assess whether known GTPase dynamics can give rise to a sufficient diversity of spatial signaling states to explain the observed shapes. Results show that the combination of autoactivation and mutually antagonistic cross talk between GTPases, along with the conservative membrane binding, generates a wide array of distinct homogeneous and polarized regulatory phenotypes that arise for fixed model parameters. From a theoretical perspective, these results demonstrate that simple GTPase dynamics can generate complex multistability in which six distinct stable steady states (three homogeneous and three polarized) coexist for a fixed set of parameters, each of which naturally maps to an observed morphology. From a biological perspective, although we do not explicitly model the cytoskeleton or resulting cell morphologies, these results, along with prior literature linking GTPase activity to cell morphology, support the hypothesis that GTPase signaling dynamics can generate the broad morphological characteristics observed in many migratory cell populations. Further, the observed diversity may be the result of cells populating a complex morphological landscape generated by GTPase regulation rather than being the result of intrinsic cell-cell variation. These results demonstrate that Rho GTPases may have a central role in regulating the broad characteristics of cell shape (e.g., expansive, contractile, polarized, etc.) and that shape heterogeneity may be (at least partly) a reflection of the rich signaling dynamics regulating the cytoskeleton rather than intrinsic cell heterogeneity.

Correlated random walks inside a cell: actin branching and microtubule dynamics

Correlated random walks (CRW) have been explored in many settings, most notably in the motion of individuals in a swarm or flock. But some subcellular systems such as growth or disassembly of bio-polymers can also be described with similar models and understood using related mathematical methods. Here we consider two examples of growing cytoskeletal elements, actin and microtubules. We use CRW or generalized CRW-like PDEs to model their spatial distributions. In each case, the linear models can be reduced to a Telegrapher's equation. A combination of explicit solutions (in one case) and numerical solutions (in the other) demonstrates that the approach to steady state can be accompanied by (decaying) waves.

Computational modeling of single-cell mechanics and cytoskeletal mechanobiology

Cellular cytoskeletal mechanics plays a major role in many aspects of human health from organ development to wound healing, tissue homeostasis and cancer metastasis. We summarize the state-of-the-art techniques for mathematically modeling cellular stiffness and mechanics and the cytoskeletal components and factors that regulate them. We highlight key experiments that have assisted model parameterization and compare the advantages of different models that have been used to recapitulate these experiments. An overview of feed-forward mechanisms from signaling to cytoskeleton remodeling is provided, followed by a discussion of the rapidly growing niche of encapsulating feedback mechanisms from cytoskeletal and cell mechanics to signaling. We discuss broad areas of advancement that could accelerate research and understanding of cellular mechanobiology. A precise understanding of the molecular mechanisms that affect cell and tissue mechanics and function will underpin innovations in medical device technologies of the future. WIREs Syst Biol Med 2018, 10:e1407. doi: 10.1002/wsbm.1407 This article is categorized under: Models of Systems Properties and Processes > Mechanistic Models Physiology > Mammalian Physiology in Health and Disease Models of Systems Properties and Processes > Cellular Models.

Fronts and waves of actin polymerization in a bistability-based mechanism of circular dorsal ruffles

During macropinocytosis, cells remodel their morphologies for the uptake of extracellular matter. This endocytotic mechanism relies on the collapse and closure of precursory structures, which are propagating actin-based, ring-shaped vertical undulations at the dorsal (top) cell membrane, a.k.a. circular dorsal ruffles (CDRs). As such, CDRs are essential to a range of vital and pathogenic processes alike. Here we show, based on both experimental data and theoretical analysis, that CDRs are propagating fronts of actin polymerization in a bistable system. The theory relies on a novel mass-conserving reaction–diffusion model, which associates the expansion and contraction of waves to distinct counter-propagating front solutions. Moreover, the model predicts that under a change in parameters (for example, biochemical conditions) CDRs may be pinned and fluctuate near the cell boundary or exhibit complex spiral wave dynamics due to a wave instability. We observe both phenomena also in our experiments indicating the conditions for which macropinocytosis is suppressed.

Circular dorsal ruffles (CDRs) are important for the vesicular uptake of extracellular matter, but the basis of their wave dynamics is not understood. Here, the authors propose and experimentally test a bistable reaction-diffusion system, which they show accounts for the typical CDR expansion and shrinkage and for aberrant formation of pinned waves and spirals.

A mathematical model coupling polarity signaling to cell adhesion explains diverse cell migration patterns

Protrusion and retraction of lamellipodia are common features of eukaryotic cell motility. As a cell migrates through its extracellular matrix (ECM), lamellipod growth increases cell-ECM contact area and enhances engagement of integrin receptors, locally amplifying ECM input to internal signaling cascades. In contrast, contraction of lamellipodia results in reduced integrin engagement that dampens the level of ECM-induced signaling. These changes in cell shape are both influenced by, and feed back onto ECM signaling. Motivated by experimental observations on melanoma cells lines (1205Lu and SBcl2) migrating on fibronectin (FN) coated topographic substrates (anisotropic post-density arrays), we probe this interplay between intracellular and ECM signaling. Experimentally, cells exhibited one of three lamellipodial dynamics: persistently polarized, random, or oscillatory, with competing lamellipodia oscillating out of phase (Park et al., 2017). Pharmacological treatments, changes in FN density, and substrate topography all affected the fraction of cells exhibiting these behaviours. We use these observations as constraints to test a sequence of hypotheses for how intracellular (GTPase) and ECM signaling jointly regulate lamellipodial dynamics. The models encoding these hypotheses are predicated on mutually antagonistic Rac-Rho signaling, Rac-mediated protrusion (via activation of Arp2/3 actin nucleation) and Rho-mediated contraction (via ROCK phosphorylation of myosin light chain), which are coupled to ECM signaling that is modulated by protrusion/contraction. By testing each model against experimental observations, we identify how the signaling layers interact to generate the diverse range of cell behaviors, and how various molecular perturbations and changes in ECM signaling modulate the fraction of cells exhibiting each. We identify several factors that play distinct but critical roles in generating the observed dynamic: (1) competition between lamellipodia for shared pools of Rac and Rho, (2) activation of RhoA by ECM signaling, and (3) feedback from lamellipodial growth or contraction to cell-ECM contact area and therefore to the ECM signaling level.

Spatiotemporal Patterns of Noise-Driven Confined Actin Waves in Living Cells

Cells utilize waves of polymerizing actin to reshape their morphologies, which is central to physiological and pathological processes alike. Here, we force dorsal actin waves to propagate on one-dimensional domains with periodic boundary conditions, which results in striking spatiotemporal patterns with a clear signature of noise-driven dynamics. We show that these patterns can be very closely reproduced with a noise-driven active medium at coherence resonance.

Local perturbation analysis: a computational tool for biophysical reaction-diffusion models

Diffusion and interaction of molecular regulators in cells is often modeled using reaction-diffusion partial differential equations. Analysis of such models and exploration of their parameter space is challenging, particularly for systems of high dimensionality. Here, we present a relatively simple and straightforward analysis, the local perturbation analysis, that reveals how parameter variations affect model behavior. This computational tool, which greatly aids exploration of the behavior of a model, exploits a structural feature common to many cellular regulatory systems: regulators are typically either bound to a membrane or freely diffusing in the interior of the cell. Using well-documented, readily available bifurcation software, the local perturbation analysis tracks the approximate early evolution of an arbitrarily large perturbation of a homogeneous steady state. In doing so, it provides a bifurcation diagram that concisely describes various regimes of the model's behavior, reducing the need for exhaustive simulations to explore parameter space. We explain the method and provide detailed step-by-step guides to its use and application.

Dynamics of actin waves on patterned substrates: a quantitative analysis of circular dorsal ruffles

Circular Dorsal Ruffles (CDRs) have been known for decades, but the mechanism that organizes these actin waves remains unclear. In this article we systematically analyze the dynamics of CDRs on fibroblasts with respect to characteristics of current models of actin waves. We studied CDRs on heterogeneously shaped cells and on cells that we forced into disk-like morphology. We show that CDRs exhibit phenomena such as periodic cycles of formation, spiral patterns, and mutual wave annihilations that are in accord with an active medium description of CDRs. On cells of controlled morphologies, CDRs exhibit extremely regular patterns of repeated wave formation and propagation, whereas on random-shaped cells the dynamics seem to be dominated by the limited availability of a reactive species. We show that theoretical models of reaction-diffusion type incorporating conserved species capture partially the behavior we observe in our data.

Actin and PIP3 waves in giant cells reveal the inherent length scale of an excited state

The membrane and actin cortex of a motile cell can autonomously differentiate into two states, one typical of the front, the other of the tail. On the substrate-attached surface of Dictyostelium discoideum cells, dynamic patterns of front-like and tail-like states are generated that are well suited to monitor transitions between these states. To image large-scale pattern dynamics independently of boundary effects, we produced giant cells by electric-pulse-induced cell fusion. In these cells, actin waves are coupled to the front and back of phosphatidylinositol (3,4,5)-trisphosphate (PIP3)-rich bands that have a finite width. These composite waves propagate across the plasma membrane of the giant cells with undiminished velocity. After any disturbance, the bands of PIP3 return to their intrinsic width. Upon collision, the waves locally annihilate each other and change direction; at the cell border they are either extinguished or reflected. Accordingly, expanding areas of progressing PIP3 synthesis become unstable beyond a critical radius, their center switching from a front-like to a tail-like state. Our data suggest that PIP3 patterns in normal-sized cells are segments of the self-organizing patterns that evolve in giant cells.

An efficient, nonlinear stability analysis for detecting pattern formation in reaction diffusion systems

Reaction diffusion systems are often used to study pattern formation in biological systems. However, most methods for understanding their behavior are challenging and can rarely be applied to complex systems common in biological applications. I present a relatively simple and efficient, nonlinear stability technique that greatly aids such analysis when rates of diffusion are substantially different. This technique reduces a system of reaction diffusion equations to a system of ordinary differential equations tracking the evolution of a large amplitude, spatially localized perturbation of a homogeneous steady state. Stability properties of this system, determined using standard bifurcation techniques and software, describe both linear and nonlinear patterning regimes of the reaction diffusion system. I describe the class of systems this method can be applied to and demonstrate its application. Analysis of Schnakenberg and substrate inhibition models is performed to demonstrate the methods capabilities in simplified settings and show that even these simple models have nonlinear patterning regimes not previously detected. The real power of this technique, however, is its simplicity and applicability to larger complex systems where other nonlinear methods become intractable. This is demonstrated through analysis of a chemotaxis regulatory network comprised of interacting proteins and phospholipids. In each case, predictions of this method are verified against results of numerical simulation, linear stability, asymptotic, and/or full PDE bifurcation analyses.