Analysis of a minimal Rho-GTPase circuit regulating cell shape

Rho of Plants patterning: linking mathematical models and molecular diversity

ROPs (Rho of Plants) are plant specific small GTPases involved in many membrane patterning processes and play important roles in the establishment and communication of cell polarity. These small GTPases can produce a wide variety of patterns, ranging from a single cluster in tip-growing root hairs and pollen tubes to an oriented stripe pattern controlling protoxylem cell wall deposition. For an understanding of what controls these various patterns, models are indispensable. Consequently, many modelling studies on small GTPase patterning exist, often focusing on yeast or animal cells. Multiple patterns occurring in plants, however, require the stable co-existence of multiple active ROP clusters, which does not occur with the most common yeast/animal models. The possibility of such patterns critically depends on the precise model formulation. Additionally, different small GTPases are usually treated interchangeably in models, even though plants possess two types of ROPs with distinct molecular properties, one of which is unique to plants. Furthermore, the shape and even the type of ROP patterns may be affected by the cortical cytoskeleton, and cortex composition and anisotropy differ dramatically between plants and animals. Here, we review insights into ROP patterning from modelling efforts across kingdoms, as well as some outstanding questions arising from these models and recent experimental findings.

A multiscale theory for spreading and migration of adhesion-reinforced mesenchymal cells

We present a chemomechanical whole-cell theory for the spreading and migration dynamics of mesenchymal cells that can actively reinforce their adhesion to an underlying viscoelastic substrate as a function of its stiffness. Our multiscale model couples the adhesion reinforcement effect at the subcellular scale with the nonlinear mechanics of the nucleus-cytoskeletal network complex at the cellular scale to explain the concurrent monotonic area-stiffness and non-monotonic speed-stiffness relationships observed in experiments: we consider that large cell spreading on stiff substrates flattens the nucleus, increasing the viscous drag force on it. The resulting force balance dictates a reduction in the migration speed on stiff substrates. We also reproduce the experimental influence of the substrate viscosity on the cell spreading area and migration speed by elucidating how the viscosity may either maintain adhesion reinforcement or prevent it depending on the substrate stiffness. Additionally, our model captures the experimental directed migration behaviour of the adhesion-reinforced cells along a stiffness gradient, known as durotaxis, as well as up or down a viscosity gradient (viscotaxis or anti-viscotaxis), the cell moving towards an optimal viscosity in either case. Overall, our theory explains the intertwined mechanics of the cell spreading, migration speed and direction in the presence of the molecular adhesion reinforcement mechanism.

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.

A multiscale whole-cell theory for mechanosensitive migration on viscoelastic substrates

Increasing experimental evidence validates that both the elastic stiffness and viscosity of the extracellular matrix regulate mesenchymal cell behavior, such as the rational switch between durotaxis (cell migration to stiffer regions), anti-durotaxis (migration to softer regions), and adurotaxis (stiffness-insensitive migration). To reveal the mechanisms underlying the crossover between these motility regimes, we have developed a multiscale chemomechanical whole-cell theory for mesenchymal migration. Our framework couples the subcellular focal adhesion dynamics at the cell-substrate interface with the cellular cytoskeletal mechanics and the chemical signaling pathways involving Rho GTPase proteins. Upon polarization by the Rho GTPase gradients, our simulated cell migrates by concerted peripheral protrusions and contractions, a hallmark of the mesenchymal mode. The resulting cell dynamics quantitatively reproduces the experimental migration speed as a function of the uniform substrate stiffness and explains the influence of viscosity on the migration efficiency. In the presence of stiffness gradients and absence of chemical polarization, our simulated cell can exhibit durotaxis, anti-durotaxis, and adurotaxis respectively with increasing substrate stiffness or viscosity. The cell moves toward an optimally stiff region from softer regions during durotaxis and from stiffer regions during anti-durotaxis. We show that cell polarization through steep Rho GTPase gradients can reverse the migration direction dictated by the mechanical cues. Overall, our theory demonstrates that opposing durotactic behaviors emerge via the interplay between intracellular signaling and cell-medium mechanical interactions in agreement with experiments, thereby elucidating complex mechanosensing at the single-cell level.

ARL11 correlates with the immunosuppression and poor prognosis in breast cancer: A comprehensive bioinformatics analysis of ARL family members

ADP-ribosylation factor-like protein (ARL) family members (ARLs) may regulate the malignant phenotypes of cancer cells. However, relevant studies on ARLs in breast cancer (BC) are limited. In this research, the expression profiles, genetic variations, and prognostic values of ARLs in BC have been systematically analyzed for the first time using various databases. We find that ARLs are significantly dysregulated in BC according to the TCGA database, which may result from DNA methylation and copy number alteration. Prognostic analysis suggests that ARL11 is the most significant prognostic indicator for BC, and higher ARL11 predicts worse clinical outcomes for BC patients. Further functional enrichment analysis demonstrates that ARL11 enhances the immunosuppression in BC, and dysregulation of ARL11 is significantly associated with immune infiltration in various types of cancer. Our results demonstrate the potential of ARL11 as an immune therapeutic target for BC.

Cell Repolarization: A Bifurcation Study of Spatio-Temporal Perturbations of Polar Cells

The intrinsic polarity of migrating cells is regulated by spatial distributions of protein activity. Those proteins (Rho-family GTPases, such as Rac and Rho) redistribute in response to stimuli, determining the cell front and back. Reaction-diffusion equations with mass conservation and positive feedback have been used to explain initial polarization of a cell. However, the sensitivity of a polar cell to a reversal stimulus has not yet been fully understood. We carry out a PDE bifurcation analysis of two polarity models to investigate routes to repolarization: (1) a single-GTPase ("wave-pinning") model and (2) a mutually antagonistic Rac-Rho model. We find distinct routes to reversal in (1) vs. (2). We show numerical simulations of full PDE solutions for the RD equations, demonstrating agreement with predictions of the bifurcation results. Finally, we show that simulations of the polarity models in deforming 1D model cells are consistent with biological experiments.

Interplay between mechanics and signalling in regulating cell fate

Mechanical signalling affects multiple biological processes during development and in adult organisms, including cell fate transitions, cell migration, morphogenesis and immune responses. Here, we review recent insights into the mechanisms and functions of two main routes of mechanical signalling: outside-in mechanical signalling, such as mechanosensing of substrate properties or shear stresses; and mechanical signalling regulated by the physical properties of the cell surface itself. We discuss examples of how these two classes of mechanical signalling regulate stem cell function, as well as developmental processes in vivo. We also discuss how cell surface mechanics affects intracellular signalling and, in turn, how intracellular signalling controls cell surface mechanics, generating feedback into the regulation of mechanosensing. The cooperation between mechanosensing, intracellular signalling and cell surface mechanics has a profound impact on biological processes. We discuss here our understanding of how these three elements interact to regulate stem cell fate and development.

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.

Modeling cell protrusion predicts how myosin II and actin turnover affect adhesion-based signaling

Orchestration of cell migration is essential for development, tissue regeneration, and the immune response. This dynamic process integrates adhesion, signaling, and cytoskeletal subprocesses across spatial and temporal scales. In mesenchymal cells, adhesion complexes bound to extracellular matrix mediate both biochemical signal transduction and physical interaction with the F-actin cytoskeleton. Here, we present a mathematical model that offers insight into both aspects, considering spatiotemporal dynamics of nascent adhesions, active signaling molecules, mechanical clutching, actin treadmilling, and nonmuscle myosin II contractility. At the core of the model is a positive feedback loop, whereby adhesion-based signaling promotes generation of barbed ends at, and protrusion of, the cell's leading edge, which in turn promotes formation and stabilization of nascent adhesions. The model predicts a switch-like transition and optimality of membrane protrusion, determined by the balance of actin polymerization and retrograde flow, with respect to extracellular matrix density. The model, together with new experimental measurements, explains how protrusion can be modulated by mechanical effects (nonmuscle myosin II contractility and adhesive bond stiffness) and F-actin turnover.

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.

Spatiotemporal development of coexisting wave domains of Rho activity in the cell cortex

The Rho family GTPases are molecular switches that regulate cytoskeletal dynamics and cell movement through a complex spatiotemporal organization of their activity. In Patiria miniata (starfish) oocytes under in vitro experimental conditions (with overexpressed Ect2, induced expression of Δ90 cyclin B, and roscovitine treatment), such activity generates multiple co-existing regions of coherent propagation of actin waves. Here we use computational modeling to investigate the development and properties of such wave domains. The model reveals that the formation of wave domains requires a balance between the activation and inhibition in the Rho signaling motif. Intriguingly, the development of the wave domains is preceded by a stage of low-activity quasi-static patterns, which may not be readily observed in experiments. Spatiotemporal patterns of this stage and the different paths of their destabilization define the behavior of the system in the later high-activity (observable) stage. Accounting for a strong intrinsic noise allowed us to achieve good quantitative agreement between simulated dynamics in different parameter regimes of the model and different wave dynamics in Patiria miniata and wild type Xenopus laevis (frog) data. For quantitative comparison of simulated and experimental results, we developed an automated method of wave domain detection, which revealed a sharp reversal in the process of pattern formation in starfish oocytes. Overall, our findings provide an insight into spatiotemporal regulation of complex and diverse but still computationally reproducible cell-level actin dynamics.

Mass conservative reaction-diffusion systems describing cell polarity

A reaction-diffusion system with mass conservation modeling cell polarity is considered. A range of the parameters is found where the ω-limit set of the solution is spatially homogeneous, containing constant stationary solution as well as possible nontrivial spatially homogeneous orbit.

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.

ARL4C might serve as a prognostic factor and a novel therapeutic target for gastric cancer: bioinformatics analyses and biological experiments

The ADP‐ribosylation factor‐like proteins (ARLs) have been proved to regulate the malignant phenotypes of several cancers. However, the exact role of ARLs in gastric cancer (GC) remains elusive. In this study, we systematically investigate the expression status, interactive relations, potential pathways, genetic variations and clinical values of ARLs in GC. We find that ARLs are significantly dysregulated in GC and involved in various cancer‐related pathways. Subsequently, machine learning models identify ARL4C as one of the two most significant clinical indicators among ARLs for GC. Furthermore, ARL4C silencing remarkably inhibits the growth and metastasis of GC cells both in vitro and in vivo. Moreover, enrichment analysis indicates that ARL4C is highly correlated with TGF‐β1 signalling. Correspondingly, TGF‐β1 treatment dramatically increases ARL4C expression and ARL4C knockdown inhibits the phosphorylation level of Smads, downstream factors of TGF‐β1. Meanwhile, the coexpression of ARL4C and TGF‐β1 worsens the prognosis of GC patients. Our work comprehensively demonstrates the crucial role of ARLs in the carcinogenesis of GC and the specific mechanisms underlying the GC‐promoting effects of TGF‐β1. More importantly, we uncover the great promise of ARL4C‐targeted therapy in improving the efficacy of TGF‐β1 inhibitors for GC patients.

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.

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.

A free boundary mechanobiological model of epithelial tissues

In this study, we couple intracellular signalling and cell-based mechanical properties to develop a novel free boundary mechanobiological model of epithelial tissue dynamics. Mechanobiological coupling is introduced at the cell level in a discrete modelling framework, and new reaction-diffusion equations are derived to describe tissue-level outcomes. The free boundary evolves as a result of the underlying biological mechanisms included in the discrete model. To demonstrate the accuracy of the continuum model, we compare numerical solutions of the discrete and continuum models for two different signalling pathways. First, we study the Rac-Rho pathway where cell- and tissue-level mechanics are directly related to intracellular signalling. Second, we study an activator-inhibitor system which gives rise to spatial and temporal patterning related to Turing patterns. In all cases, the continuum model and free boundary condition accurately reflect the cell-level processes included in the discrete model.

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.

Mechanical Cell Competition in Heterogeneous Epithelial Tissues

Mechanical cell competition is important during tissue development, cancer invasion, and tissue ageing. Heterogeneity plays a key role in practical applications since cancer cells can have different cell stiffness and different proliferation rates than normal cells. To study this phenomenon, we propose a one-dimensional mechanical model of heterogeneous epithelial tissue dynamics that includes cell-length-dependent proliferation and death mechanisms. Proliferation and death are incorporated into the discrete model stochastically and arise as source/sink terms in the corresponding continuum model that we derive. Using the new discrete model and continuum description, we explore several applications including the evolution of homogeneous tissues experiencing proliferation and death, and competition in a heterogeneous setting with a cancerous tissue competing for space with an adjacent normal tissue. This framework allows us to postulate new mechanisms that explain the ability of cancer cells to outcompete healthy cells through mechanical differences rather than an intrinsic proliferative advantage. We advise when the continuum model is beneficial and demonstrate why naively adding source/sink terms to a continuum model without considering the underlying discrete model may lead to incorrect results.

Periodic propagating waves coordinate RhoGTPase network dynamics at the leading and trailing edges during cell migration

Migrating cells need to coordinate distinct leading and trailing edge dynamics but the underlying mechanisms are unclear. Here, we combine experiments and mathematical modeling to elaborate the minimal autonomous biochemical machinery necessary and sufficient for this dynamic coordination and cell movement. RhoA activates Rac1 via DIA and inhibits Rac1 via ROCK, while Rac1 inhibits RhoA through PAK. Our data suggest that in motile, polarized cells, RhoA–ROCK interactions prevail at the rear, whereas RhoA-DIA interactions dominate at the front where Rac1/Rho oscillations drive protrusions and retractions. At the rear, high RhoA and low Rac1 activities are maintained until a wave of oscillatory GTPase activities from the cell front reaches the rear, inducing transient GTPase oscillations and RhoA activity spikes. After the rear retracts, the initial GTPase pattern resumes. Our findings show how periodic, propagating GTPase waves coordinate distinct GTPase patterns at the leading and trailing edge dynamics in moving cells.

Bayesian integrative analysis of epigenomic and transcriptomic data identifies Alzheimer's disease candidate genes and networks

Biomedical research studies have generated large multi-omic datasets to study complex diseases like Alzheimer’s disease (AD). An important aim of these studies is the identification of candidate genes that demonstrate congruent disease-related alterations across the different data types measured by the study. We developed a new method to detect such candidate genes in large multi-omic case-control studies that measure multiple data types in the same set of samples. The method is based on a gene-centric integrative coefficient quantifying to what degree consistent differences are observed in the different data types. For statistical inference, a Bayesian hierarchical model is used to study the distribution of the integrative coefficient. The model employs a conditional autoregressive prior to integrate a functional gene network and to share information between genes known to be functionally related. We applied the method to an AD dataset consisting of histone acetylation, DNA methylation, and RNA transcription data from human cortical tissue samples of 233 subjects, and we detected 816 genes with consistent differences between persons with AD and controls. The findings were validated in protein data and in RNA transcription data from two independent AD studies. Finally, we found three subnetworks of jointly dysregulated genes within the functional gene network which capture three distinct biological processes: myeloid cell differentiation, protein phosphorylation and synaptic signaling. Further investigation of the myeloid network indicated an upregulation of this network in early stages of AD prior to accumulation of hyperphosphorylated tau and suggested that increased CSF1 transcription in astrocytes may contribute to microglial activation in AD. Thus, we developed a method that integrates multiple data types and external knowledge of gene function to detect candidate genes, applied the method to an AD dataset, and identified several disease-related genes and processes demonstrating the usefulness of the integrative approach.

Recent technological advances have led to a new generation of studies that interrogate multiple molecular levels in the same target tissue of a set of subjects, generating complex multi-omic datasets with which to study disease mechanism. These datasets of genetic, epigenomic, transcriptomic, and other data have the potential to reveal novel biological insights; however, integrative analyses remain challenging and require new computational methods. We developed an integrative Bayesian approach to detect genes with consistent differences between case and control samples across multiple data types. The method further integrates prior knowledge about gene function in the form of a gene functional similarity network to improve statistical inference by sharing information between related genes. We applied our method to an Alzheimer’s disease dataset of epigenomic and transcriptomic data and detected and then validated several novel and known candidate genes as well as three major disease-related biological processes. One of these processes reflected microglial activation and included the cytokine CSF1. Single-nucleus data revealed that CSF1 was primarily upregulated in astrocytes, implicating the involvement of this cell type in microglial activation. Hence, we demonstrated that integrative analysis approaches to multi-omic datasets can improve candidate gene detection and thereby generate new insights into complex diseases.

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.

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.

Rac activation is key to cell motility and directionality: An experimental and modelling investigation

Cell migration is a tightly-regulated process that involves protein gradients formed by the Rho family of GTPases, including Rho and Rac. The front (rear) of cells is generally characterized by higher active Rac (Rho) and lower active Rho (Rac) concentrations. Protein clusters, called adhesions, that anchor cells to their external environment have been shown to be dynamic and small (stable and large) at the cell front (rear), forming the force-transmission points necessary for persistent movement. Differences in adhesion sizes and dynamics have been linked to gradients in Rac and Rho activity. Here, we study the effects of Rac activation and gradients in Rac and Rho concentrations and activities on cellular polarity and adhesion size using mathematical and experimental approaches. The former is accomplished by expanding an existing reaction-diffusion model to a 2D domain utilizing stochastic dynamics. The model revealed that a hysteresis between the induced/uninduced states (corresponding to higher/lower Rac concentrations, respectively) along with Rac and Rho activation gradients, generated by chemical cues, were vital for forming polarity. Experimentally, the induced state was generated by increasing the cellular βPIX (a Rac-GEF) level and/or decreasing ROCK (a Rac-GAP effector protein) activity with Y-27632 (a ROCK-inhibitor). In agreement with the simulations, our results showed that cells with elevated RacGTP migrated faster, indicating more robust cellular polarization. However, the directionality of cells was not changed significantly, suggesting that external and/or internal physical or chemical cues were needed. Complementing the faster migration observed, adhesions were smaller, generating the phenotype expected with the induced state.

Spatiotemporal dynamics of a reaction-diffusion model of pollen tube tip growth

A reaction-diffusion model is proposed to describe the mechanisms underlying the spatial distributions of ROP1 and calcium on the pollen tube tip. The model assumes that the plasma membrane ROP1 activates itself through positive feedback loop, while the cytosolic calcium ions inhibit ROP1 via a negative feedback loop. Furthermore it is proposed that lateral movement of molecules on the plasma membrane are depicted by diffusion. It is shown that bistable or oscillatory dynamics could exist even in the non-spatial model, and stationary and oscillatory spatiotemporal patterns are found in the full spatial model which resemble the experimental data of pollen tube tip growth.

Mechanochemical Coupling and Junctional Forces during Collective Cell Migration

Cell migration, a fundamental physiological process in which cells sense and move through their surrounding physical environment, plays a critical role in development and tissue formation, as well as pathological processes, such as cancer metastasis and wound healing. During cell migration, dynamics are governed by the bidirectional interplay between cell-generated mechanical forces and the activity of Rho GTPases, a family of small GTP-binding proteins that regulate actin cytoskeleton assembly and cellular contractility. These interactions are inherently more complex during the collective migration of mechanically coupled cells because of the additional regulation of cell-cell junctional forces. In this study, we adapted a recent minimal modeling framework to simulate the interactions between mechanochemical signaling in individual cells and interactions with cell-cell junctional forces during collective cell migration. We find that migration of individual cells depends on the feedback between mechanical tension and Rho GTPase activity in a biphasic manner. During collective cell migration, waves of Rho GTPase activity mediate mechanical contraction/extension and thus synchronization throughout the tissue. Further, cell-cell junctional forces exhibit distinct spatial patterns during collective cell migration, with larger forces near the leading edge. Larger junctional force magnitudes are associated with faster collective cell migration and larger tissue size. Simulations of heterogeneous tissue migration exhibit a complex dependence on the properties of both leading and trailing cells. Computational predictions demonstrate that collective cell migration depends on both the emergent dynamics and interactions between cellular-level Rho GTPase activity and contractility and multicellular-level junctional forces.

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.

Paxillin phosphorylation at serine 273 and its effects on Rac, Rho and adhesion dynamics

Focal adhesions are protein complexes that anchor cells to the extracellular matrix. During migration, the growth and disassembly of these structures are spatiotemporally regulated, with new adhesions forming at the leading edge of the cell and mature adhesions disassembling at the rear. Signalling proteins and structural cytoskeletal components tightly regulate adhesion dynamics. Paxillin, an adaptor protein within adhesions, is one of these proteins. Its phosphorylation at serine 273 (S273) is crucial for maintaining fast adhesion assembly and disassembly. Paxillin is known to bind to a GIT1-βPIX-PAK1 complex, which increases the local activation of the small GTPase Rac. To understand quantitatively the behaviour of this system and how it relates to adhesion assembly/disassembly, we developed a mathematical model describing the dynamics of the small GTPases Rac and Rho as determined by paxillin S273 phosphorylation. Our model revealed that the system possesses bistability, where switching between uninduced (active Rho) and induced (active Rac) states can occur through a change in rate of paxillin phosphorylation or PAK1 activation. The bistable switch is characterized by the presence of memory, minimal change in the levels of active Rac and Rho within the induced and uninduced states, respectively, and the limited regime of monostability associated with the uninduced state. These results were validated experimentally by showing the presence of bimodality in adhesion assembly and disassembly rates, and demonstrating that Rac activity increases after treating Chinese Hamster Ovary cells with okadaic acid (a paxillin phosphatase inhibitor), followed by a modest recovery after 20 min washout. Spatial gradients of phosphorylated paxillin in a reaction-diffusion model gave rise to distinct regions of Rac and Rho activities, resembling polarization of a cell into front and rear. Perturbing several parameters of the model also revealed important insights into how signalling components upstream and downstream of paxillin phosphorylation affect dynamics.

Cellular migration is crucial in both physiological and pathological functions. Maintenance of proper migration and development of aberrant migration are effectuated by cellular machinery involving protein complexes, called adhesions, that anchor the cell to its environment. Over time, these adhesions assemble at the leading edge, as the cell extends forward, anchoring the front of the cells to its substrate, while those at the cell rear disassemble, allowing detachment and forward movement. Their dynamics are controlled by a number of regulatory factors, occurring on both cell-wide and adhesion-level scales. The coordination of these regulatory factors is complex, but insights about their dynamics can be gained from the use of mathematical modeling techniques which integrate many of these components together. Here, we developed several molecularly explicit models to explore how local regulation of paxillin, an adhesion protein, interacts with the activities of Rac and Rho to produce cell-wide polarization associated with motility and directionality. By altering paxillin phosphorylation/dephosphorylation within such models, we have advanced our understanding of how a shift from a non-motile state to a highly motile state occurs. Deciphering these key processes quantitatively thus helped us gain insight into the subcellular factors underlying polarity and movement.

Rapid valproic acid-induced modulation of the traumatic proteome in a porcine model of traumatic brain injury and hemorrhagic shock

BACKGROUND

Histone deacetylase inhibitors such as valproic acid (VPA) improve survival in lethal models of hemorrhagic shock and polytrauma. Although VPA is known to modulate transcription, its ability to reduce mortality within minutes of administration suggests involvement of a rapid, posttranslational mechanism. We hypothesized that VPA treatment would cause proteomic changes within minutes of treatment including quantitative and/or posttranslational differences in structural and/or effector proteins.

MATERIALS AND METHODS

We used a porcine model of traumatic brain injury (computer-controlled cortical impact, 12 mm depth) and hemorrhagic shock (40% hemorrhage). Animals were kept in shock for 2 h and randomized to two groups (n = 3): normal saline (volume = 3:1 hemorrhage volume) or normal saline + VPA (150 mg/kg, single dose). Peripheral blood mononuclear cells were collected at baseline, postshock, and postresuscitation. Intracellular protein profiles were assessed using 1 dimensional gel electrophoresis, liquid chromatography, mass spectrometry, and analyzed with Ingenuity Pathway Analysis software.

RESULTS

Animals treated with VPA demonstrated significant proteomic changes. Quantitative differences were found in over 200 proteins including effector, regulatory, and structural proteins in critical cell signaling pathways. Posttranslational modification analysis demonstrated differential VPA-induced acetylation of lysine residues in histone and nonhistone proteins. Pathway analysis correlated these changes with significant increases in numerous prosurvival and cytoskeletal intracellular pathways, including Rho GTPase signaling (P = 1.66E-11), integrin signaling (P = 4.19E-21), and a decrease in Rho guanosine nucleotide dissociation inhibitor signaling (P = 4.83E-12).

CONCLUSIONS

In a porcine model of severe injuries, a single dose of VPA is associated with protective changes in the proteome that are measurable within minutes of treatment.

A Boolean network of the crosstalk between IGF and Wnt signaling in aging satellite cells

Aging is a complex biological process, which determines the life span of an organism. Insulin-like growth factor (IGF) and Wnt signaling pathways govern the process of aging. Both pathways share common downstream targets that allow competitive crosstalk between these branches. Of note, a shift from IGF to Wnt signaling has been observed during aging of satellite cells. Biological regulatory networks necessary to recreate aging have not yet been discovered. Here, we established a mathematical in silico model that robustly recapitulates the crosstalk between IGF and Wnt signaling. Strikingly, it predicts critical nodes following a shift from IGF to Wnt signaling. These findings indicate that this shift might cause age-related diseases.

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.

Many roads to symmetry breaking: molecular mechanisms and theoretical models of yeast cell polarity

Mathematical modeling has been instrumental in identifying common principles of cell polarity across diverse systems. These principles include positive feedback loops that are required to destabilize a spatially uniform state of the cell. The conserved small G-protein Cdc42 is a master regulator of eukaryotic cellular polarization. Here we discuss recent developments in studies of Cdc42 polarization in budding and fission yeasts and demonstrate that models describing symmetry-breaking polarization can be classified into six minimal classes based on the structure of positive feedback loops that activate and localize Cdc42. Owing to their generic system-independent nature, these model classes are also likely to be relevant for the G-protein–based symmetry-breaking systems of higher eukaryotes. We review experimental evidence pro et contra different theoretically plausible models and conclude that several parallel and non–mutually exclusive mechanisms are likely involved in cellular polarization of yeasts. This potential redundancy needs to be taken into consideration when interpreting the results of recent cell-rewiring studies.

Stability analysis in spatial modeling of cell signaling

Advances in high-resolution microscopy and other techniques have emphasized the spatio-temporal nature of information transfer through signal transduction pathways. The compartmentalization of signaling molecules and the existence of microdomains are now widely acknowledged as key features in biochemical signaling. To complement experimental observations of spatio-temporal dynamics, mathematical modeling has emerged as a powerful tool. Using modeling, one can not only recapitulate experimentally observed dynamics of signaling molecules, but also gain an understanding of the underlying mechanisms in order to generate experimentally testable predictions. Reaction-diffusion systems are commonly used to this end; however, the analysis of coupled nonlinear systems of partial differential equations, generated by considering large reaction networks is often challenging. Here, we aim to provide an introductory tutorial for the application of reaction-diffusion models to the spatio-temporal dynamics of signaling pathways. In particular, we outline the steps for stability analysis of such models, with a focus on biochemical signal transduction. WIREs Syst Biol Med 2018, 10:e1395. doi: 10.1002/wsbm.1395 This article is categorized under: Biological Mechanisms > Cell Signaling Analytical and Computational Methods > Dynamical Methods Models of Systems Properties and Processes > Mechanistic Models.

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.

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.

Bistable front dynamics in a contractile medium: Travelling wave fronts and cortical advection define stable zones of RhoA signaling at epithelial adherens junctions

Mechanical coherence of cell layers is essential for epithelia to function as tissue barriers and to control active tissue dynamics during morphogenesis. RhoA signaling at adherens junctions plays a key role in this process by coupling cadherin-based cell-cell adhesion together with actomyosin contractility. Here we propose and analyze a mathematical model representing core interactions involved in the spatial localization of junctional RhoA signaling. We demonstrate how the interplay between biochemical signaling through positive feedback, combined with diffusion on the cell membrane and mechanical forces generated in the cortex, can determine the spatial distribution of RhoA signaling at cell-cell junctions. This dynamical mechanism relies on the balance between a propagating bistable signal that is opposed by an advective flow generated by an actomyosin stress gradient. Experimental observations on the behavior of the system when contractility is inhibited are in qualitative agreement with the predictions of the model.

A mathematical model of GTPase pattern formation during single-cell wound repair

Rho GTPases are regulatory proteins whose patterns on the surface of a cell affect cell polarization, cell motility and repair of single-cell wounds. The stereotypical patterns formed by two such proteins, Rho and Cdc42, around laser-injured frog oocytes permit experimental analysis of GTPase activation, inactivation, segregation and crosstalk. Here, we review the development and analysis of a spatial model of GTPase dynamics that describe the formation of concentric zones of Rho and Cdc42 activity around wounds, and describe how this model has provided insights into the roles of the GTPase effector molecules protein kinase C (PKCβ and PKCη) and guanosine nucleotide dissociation inhibitor (GDI) in the wound response. We further demonstrate how the use of a 'sharp switch' model approximation in combination with bifurcation analysis can aid mapping the model behaviour in parameter space (approximate results confirmed with numerical simulation methods). Using these methods in combination with experimental manipulation of PKC activity (PKC overexpression (OE) and dominant negative conditions), we have shown that: (i) PKCβ most probably acts by enhancing existing positive feedbacks (from Rho to itself via the guanosine nucleotide exchange factor domain of Abr, and from Cdc42 to itself), (ii) PKCη most probably increases basal rates of inactivation (or possibly decreases basal rates of activation) of Rho and Cdc42, and (iii) the graded distribution of PKCη and its effect on initial Rho activity accounts for inversion of zones in a fraction (20%) of PKCη OE cells. Finally, we speculate that GDIs (which sequester GTPases) may have a critical role in defining the spatial domain, where the wound response may occur. This paper provides a more thorough exposition of the methods of analysis used in the investigation, whereas previous work on this topic was addressed to biologists and abbreviated such discussion.

Rac1 and RhoA: Networks, loops and bistability

Cell migration requires a precise temporal and spatial coordination of several processes which allow the cell to efficiently move. The extension and retraction of membrane protrusion, as well as adhesion are controlled by the Rho-family small GTPases. Two members of the family, Rac1 and RhoA, can show opposite behaviors and spatial localisations, with RhoA being active toward the rear of the cell and regulating its retraction during migration, whereas Rac1 is active toward the front of the cell. In addition to the spatial segregation, RhoA and Rac1 activity at the leading edge of the cells has an element of temporal segregation, with RhoA and Rac1 activities peaking at separate points during the migratory cycle of protrusion and retraction. Elements of this separation have been explained by the presence of 2 mutually inhibitory feedbacks, where Rac1 inhibits RhoA and RhoA in turn can inhibit Rac1. Recently, it was shown that Rac1 and RhoA activity and downstream signaling respond in a bistable manner to perturbations of this network.