Pattern formation of Rho GTPases in single cell wound healing

Patterning of the cell cortex by Rho GTPases

The Rho GTPases - RHOA, RAC1 and CDC42 - are small GTP binding proteins that regulate basic biological processes such as cell locomotion, cell division and morphogenesis by promoting cytoskeleton-based changes in the cell cortex. This regulation results from active (GTP-bound) Rho GTPases stimulating target proteins that, in turn, promote actin assembly and myosin 2-based contraction to organize the cortex. This basic regulatory scheme, well supported by in vitro studies, led to the natural assumption that Rho GTPases function in vivo in an essentially linear matter, with a given process being initiated by GTPase activation and terminated by GTPase inactivation. However, a growing body of evidence based on live cell imaging, modelling and experimental manipulation indicates that Rho GTPase activation and inactivation are often tightly coupled in space and time via signalling circuits and networks based on positive and negative feedback. In this Review, we present and discuss this evidence, and we address one of the fundamental consequences of coupled activation and inactivation: the ability of the Rho GTPases to self-organize, that is, direct their own transition from states of low order to states of high order. We discuss how Rho GTPase self-organization results in the formation of diverse spatiotemporal cortical patterns such as static clusters, oscillatory pulses, travelling wave trains and ring-like waves. Finally, we discuss the advantages of Rho GTPase self-organization and pattern formation for cell function.

Wound Repair of the Cell Membrane: Lessons from Dictyostelium Cells

The cell membrane is frequently subjected to damage, either through physical or chemical means. The swift restoration of the cell membrane's integrity is crucial to prevent the leakage of intracellular materials and the uncontrolled influx of extracellular ions. Consequently, wound repair plays a vital role in cell survival, akin to the importance of DNA repair. The mechanisms involved in wound repair encompass a series of events, including ion influx, membrane patch formation, endocytosis, exocytosis, recruitment of the actin cytoskeleton, and the elimination of damaged membrane sections. Despite the absence of a universally accepted general model, diverse molecular models have been proposed for wound repair in different organisms. Traditional wound methods not only damage the cell membrane but also impact intracellular structures, including the underlying cortical actin networks, microtubules, and organelles. In contrast, the more recent improved laserporation selectively targets the cell membrane. Studies on Dictyostelium cells utilizing this method have introduced a novel perspective on the wound repair mechanism. This review commences by detailing methods for inducing wounds and subsequently reviews recent developments in the field.

Multilayer omics reveals the molecular mechanism of early infection of Clonorchis sinensis juvenile

BACKGROUND

Clonorchiasis remains a non-negligible global zoonosis, causing serious socioeconomic burdens in endemic areas. Clonorchis sinensis infection typically elicits Th1/Th2 mixed immune responses during the course of biliary injury and periductal fibrosis. However, the molecular mechanism by which C. sinensis juvenile initially infects the host remains poorly understood.

METHODS

The BALB/c mouse model was established to study early infection (within 7 days) with C. sinensis juveniles. Liver pathology staining and observation as well as determination of biochemical enzymes, blood routine and cytokines in blood were conducted. Furthermore, analysis of liver transcriptome, proteome and metabolome changes was performed using multi-omics techniques. Statistical analyses were performed using Student's t-test.

RESULTS

Histopathological analysis revealed that liver injury, characterized by collagen deposition and inflammatory cell infiltration, occurred as early as 24 h of infection. Blood indicators including ALT, AST, WBC, CRP and IL-6 indicated that both liver injury and systemic inflammation worsened as the infection progressed. Proteomic data showed that apoptosis and junction-related pathways were enriched within 3 days of infection, indicating the occurrence of liver injury. Furthermore, proteomic and transcriptomic analysis jointly verified that the detoxification and antioxidant defense system was activated by enrichment of glutathione metabolism and cytochrome P450-related pathways in response to acute liver injury. Proteomic-based GO analysis demonstrated that biological processes such as cell deformation, proliferation, migration and wound healing occurred in the liver during the early infection. Correspondingly, transcriptomic results showed significant enrichment of cell cycle pathway on day 3 and 7. In addition, the KEGG analysis of multi-omics data demonstrated that numerous pathways related to immunity, inflammation, tumorigenesis and metabolism were enriched in the liver. Besides, metabolomic screening identified several metabolites that could promote inflammation and hepatobiliary periductal fibrosis, such as CA7S.

CONCLUSIONS

This study revealed that acute inflammatory injury was rapidly triggered by initial infection by C. sinensis juveniles in the host, accompanied by the enrichment of detoxification, inflammation, fibrosis, tumor and metabolism-related pathways in the liver, which provides a new perspective for the early intervention and therapy of clonorchiasis.

Spatiotemporal Coordination of Rac1 and Cdc42 at the Whole Cell Level during Cell Ruffling

Rho-GTPases are central regulators within a complex signaling network that controls cytoskeletal organization and cell movement. The network includes multiple GTPases, such as the most studied Rac1, Cdc42, and RhoA, along with their numerous effectors that provide mutual regulation through feedback loops. Here we investigate the temporal and spatial relationship between Rac1 and Cdc42 during membrane ruffling, using a simulation model that couples GTPase signaling with cell morphodynamics and captures the GTPase behavior observed with FRET-based biosensors. We show that membrane velocity is regulated by the kinetic rate of GTPase activation rather than the concentration of active GTPase. Our model captures both uniform and polarized ruffling. We also show that cell-type specific time delays between Rac1 and Cdc42 activation can be reproduced with a single signaling motif, in which the delay is controlled by feedback from Cdc42 to Rac1. The resolution of our simulation output matches those of time-lapsed recordings of cell dynamics and GTPase activity. Our data-driven modeling approach allows us to validate simulation results with quantitative precision using the same pipeline for the analysis of simulated and experimental data.

Spatiotemporal coordination of Rac1 and Cdc42 at the whole cell level during cell ruffling

Rho-GTPases are central regulators within a complex signaling network that controls the cytoskeletal organization and cell movement. This network includes multiple GTPases, such as the most studied Rac1, Cdc42, and RhoA, and their numerous effectors that provide mutual regulation and feedback loops. Here we investigate the temporal and spatial relationship between Rac1 and Cdc42 during membrane ruffling using a simulation model which couples GTPase signaling with cell morphodynamics to capture the GTPase behavior observed with FRET-based biosensors. We show that membrane velocity is regulated by the kinetic rate of GTPase activation rather than the concentration of active GTPase. Our model captures both uniform and polarized ruffling. We also show that cell-type specific time delays between Rac1 and Cdc42 activation can be reproduced with a single signaling motif, in which the delay is controlled by feedback from Cdc42 to Rac1. The resolution of our simulation output matches those of the time-lapsed recordings of cell dynamics and GTPase activity. This approach allows us to validate simulation results with quantitative precision using the same pipeline for the analysis of simulated and experimental data.

Scar/WAVE has Rac GTPase-independent functions during cell wound repair

Rho family GTPases regulate both linear and branched actin dynamics by activating downstream effectors to facilitate the assembly and function of complex cellular structures such as lamellipodia and contractile actomyosin rings. Wiskott-Aldrich Syndrome (WAS) family proteins are downstream effectors of Rho family GTPases that usually function in a one-to-one correspondence to regulate branched actin nucleation. In particular, the WAS protein Scar/WAVE has been shown to exhibit one-to-one correspondence with Rac GTPase. Here we show that Rac and SCAR are recruited to cell wounds in the Drosophila repair model and are required for the proper formation and maintenance of the dynamic actomyosin ring formed at the wound periphery. Interestingly, we find that SCAR is recruited to wounds earlier than Rac and is still recruited to the wound periphery in the presence of a potent Rac inhibitor. We also show that while Rac is important for actin recruitment to the actomyosin ring, SCAR serves to organize the actomyosin ring and facilitate its anchoring to the overlying plasma membrane. These differing spatiotemporal recruitment patterns and wound repair phenotypes highlight the Rac-independent functions of SCAR and provide an exciting new context in which to investigate these newly uncovered SCAR functions.

EspH interacts with the host active Bcr related (ABR) protein to suppress RhoGTPases

Enteropathogenic Escherichia coli are bacterial pathogens that colonize the gut and cause severe diarrhea in humans. Upon intimate attachment to the intestinal epithelium, these pathogens translocate via a type III secretion system virulent proteins, termed effectors, into the host cells. These effectors manipulate diverse host cell organelles and functions for the pathogen's benefit. However, the precise mechanisms underlying their activities are not fully understood despite intensive research. EspH, a critical effector protein, has been previously reported to disrupt the host cell actin cytoskeleton by suppressing RhoGTPase guanine exchange factors. However, native host proteins targeted by EspH to mediate these activities remained unknown. Here, we identified the active Bcr related (ABR), a protein previously characterized to possess dual Rho guanine nucleotide exchange factor and GTPase activating protein (GAP) domains, as a native EspH interacting partner. These interactions are mediated by the effector protein's C-terminal 38 amino acid segment. The effector primarily targets the GAP domain of ABR to suppress Rac1 and Cdc42, host cell cytotoxicity, bacterial invasion, and filopodium formation at infection sites. Knockdown of ABR expression abolished the ability of EspH to suppress Rac1, Cdc42. Our studies unravel a novel mechanism by which host RhoGTPases are hijacked by bacterial effectors.

Spatial models of pattern formation during phagocytosis

Phagocytosis, the biological process in which cells ingest large particles such as bacteria, is a key component of the innate immune response. Fcγ receptor (FcγR)-mediated phagocytosis is initiated when these receptors are activated after binding immunoglobulin G (IgG). Receptor activation initiates a signaling cascade that leads to the formation of the phagocytic cup and culminates with ingestion of the foreign particle. In the experimental system termed "frustrated phagocytosis", cells attempt to internalize micropatterned disks of IgG. Cells that engage in frustrated phagocytosis form "rosettes" of actin-enriched structures called podosomes around the IgG disk. The mechanism that generates the rosette pattern is unknown. We present data that supports the involvement of Cdc42, a member of the Rho family of GTPases, in pattern formation. Cdc42 acts downstream of receptor activation, upstream of actin polymerization, and is known to play a role in polarity establishment. Reaction-diffusion models for GTPase spatiotemporal dynamics exist. We demonstrate how the addition of negative feedback and minor changes to these models can generate the experimentally observed rosette pattern of podosomes. We show that this pattern formation can occur through two general mechanisms. In the first mechanism, an intermediate species forms a ring of high activity around the IgG disk, which then promotes rosette organization. The second mechanism does not require initial ring formation but relies on spatial gradients of intermediate chemical species that are selectively activated over the IgG patch. Finally, we analyze the models to suggest experiments to test their validity.

A computational model of mutual antagonism in the mechano-signaling network of RhoA and nitric oxide

BACKGROUND

RhoA is a master regulator of cytoskeletal contractility, while nitric oxide (NO) is a master regulator of relaxation, e.g., vasodilation. There are multiple forms of cross-talk between the RhoA/ROCK pathway and the eNOS/NO/cGMP pathway, but previous work has not studied their interplay at a systems level. Literature review suggests that the majority of their cross-talk interactions are antagonistic, which motivates us to ask whether the RhoA and NO pathways exhibit mutual antagonism in vitro, and if so, to seek the theoretical implications of their mutual antagonism.

RESULTS

Experiments found mutual antagonism between RhoA and NO in epithelial cells. Since mutual antagonism is a common motif for bistability, we sought to explore through theoretical simulations whether the RhoA-NO network is capable of bistability. Qualitative modeling showed that there are parameters that can cause bistable switching in the RhoA-NO network, and that the robustness of the bistability would be increased by positive feedback between RhoA and mechanical tension.

CONCLUSIONS

We conclude that the RhoA-NO bistability is robust enough in silico to warrant the investment of further experimental testing. Tension-dependent bistability has the potential to create sharp concentration gradients, which could contribute to the localization and self-organization of signaling domains during cytoskeletal remodeling and cell migration.

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.

Tension modulation of actomyosin ring assembly and RhoGTPases activity: Perspectives from the Xenopus oocyte wound healing model

Cells are remarkably resilient structures; they are able to recover from injuries to their plasma membrane (PM) and cytoskeleton that would normally constitute existential threats. This capacity is exemplified by Xenopus laevis oocytes which can recover from very large PM defects through exocytotic and endocytic events and can repair damaged cortical cytoskeleton structures through the formation of a contractile actomyosin ring (AMR). Formation of the AMR involves the localized Ca²⁺ -dependent activation of RhoA and Cdc42, and the pre-patterning of guanine nucleotide exchange factors (GEFs) and GTPase-activating proteins (GAPs). However, this model fails to account for observations that suggest a link between cytoskeletal dynamics, intracellular tension, and AMR formation. It also does not explain why the formation of an AMR is not involved in the cytoskeletal repair program of adherent cells. We show here evidence for the support of tension as an essential regulatory signal for the formation of AMR. Indeed, oocytes in which global tension has been experimentally reduced were unable to form a functional AMR following injury, showing severely diminished RhoA activity at the wound site. These new insights place the cytoskeleton at the center of events involving changes in cell shape such as cytokinesis which also involves the formation and closure of an AMR.

Cross-talk-dependent cortical patterning of Rho GTPases during cell repair

Rho GTPases such as Rho, Rac, and Cdc42 are important regulators of the cortical cytoskeleton in processes including cell division, locomotion, and repair. In these processes, Rho GTPases assume characteristic patterns wherein the active GTPases occupy mutually exclusive "zones" in the cell cortex. During cell wound repair, for example, a Rho zone encircles the wound edge and is in turn encircled by a Cdc42 zone. Here we evaluated the contributions of cross-talk between Rho and Cdc42 to the patterning of their respective zones in wounded Xenopus oocytes using experimental manipulations in combination with mathematical modeling. The results show that the position of the Cdc42 zone relative to the Rho zone and relative to the wound edge is controlled by the level of Rho activity. In contrast, the outer boundary of the Rho zone is limited by the level of Cdc42 activity. Models based on positive feedback within zones and negative feedback from Rho to the GEF-GAP Abr to Cdc42 capture some, but not all, of the observed behaviors. We conclude that GTPase zone positioning is controlled at the level of Rho activity and we speculate that the Cdc42 zone or something associated with it limits the spread of Rho activity.

Transplantation of Wnt5a-modified NSCs promotes tissue repair and locomotor functional recovery after spinal cord injury

Traditional therapeutic strategies for spinal cord injury (SCI) are insufficient to repair locomotor function because of the failure of axonal reconnection and neuronal regeneration in the injured central nervous system (CNS). Neural stem cell (NSC) transplantation has been considered a potential strategy and is generally feasible for repairing the neural circuit after SCI; however, the most formidable problem is that the neuronal differentiation rate of NSCs is quite limited. Therefore, it is essential to induce the neuronal differentiation of NSCs and improve the differentiation rate of NSCs in spinal cord repair. Our results demonstrate that both Wnt5a and miRNA200b-3p could promote NSC differentiation into neurons and that Wnt5a upregulated miRNA200b-3p expression through MAPK/JNK signaling to promote NSC differentiation into neurons. Wnt5a could reduce RhoA expression by upregulating miRNA200b-3p expression to inhibit activation of the RhoA/Rock signaling pathway, which has been reported to suppress neuronal differentiation. Overexpression of RhoA abolished the neurogenic capacity of Wnt5a and miRNA200b-3p. In vivo, miRNA200b-3p was critical for Wnt5a-induced NSC differentiation into neurons to promote motor functional and histological recovery after SCI by suppressing RhoA/Rock signaling. These findings provide more insight into SCI and help with the identification of novel treatment strategies.

Incorporating key molecules into neural stem cells enhances their ability to differentiate correctly and promote repair following spinal cord injury. Spinal cord injuries can have a debilitating effect on patients’ lives, yet there are no therapies that fully restore movement and sensation. Therapies based on neural stem cells (NSCs) show promise, but initial studies show many NSCs differentiate into astrocytes, supportive cells that do not conduct nerve impulses, instead of neurons, leading to treatment failure. Yong Wan and Le Wang at Sun Yat-sen University in Guangzhou, China, and co-workers demonstrated that adding a protein called Wnt5a and a specific microRNA molecule to NSCs significantly increases differentiation into neurons. Wnt5a suppresses a signalling pathway that otherwise interferes with NSC differentiation. Experiments on rat models showed that the therapy improved locomotor function and tissue repair after injury.

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.

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.

Dynamics of actin polymerisation during the mammalian single-cell wound healing response

OBJECTIVE

The contribution of actomyosin contractile rings in the wound healing program of somatic cells as never been directly assessed. This contrast with the events characterising the wound healing response of in wounded Xenopus oocytes, in which formation and contraction of an actomyosin ring provides a platform for cytoskeletal repair and drives the restoration of proper plasma membrane composition at the site of injury. As such, we aimed to characterize, using high-resolution live-cell confocal microscopy, the cytoskeletal repair dynamics of HeLa cells.

RESULTS

We confirm here that the F-actin enrichment that characterizes the late repair program of laser-wounded cells is mostly uniform and is not associated with co-enrichment of myosin-II or the formation of concentric zones of RhoA and Cdc42 activity.

Spatially Adaptive Colocalization Analysis in Dual-Color Fluorescence Microscopy

Colocalization analysis aims to study complex spatial associations between bio-molecules via optical imaging techniques. However, existing colocalization analysis workflows only assess an average degree of colocalization within a certain region of interest and ignore the unique and valuable spatial information offered by microscopy. In the current work, we introduce a new framework for colocalization analysis that allows us to quantify colocalization levels at each individual location and automatically identify pixels or regions where colocalization occurs. The framework, referred to as spatially adaptive colocalization analysis (SACA), integrates a pixel-wise local kernel model for colocalization quantification and a multi-scale adaptive propagation-separation strategy for utilizing spatial information to detect colocalization in a spatially adaptive fashion. Applications to simulated and real biological datasets demonstrate the practical merits of SACA in what we hope to be an easily applicable and robust colocalization analysis method. In addition, theoretical properties of SACA are investigated to provide rigorous statistical justification.

Regulation and Assembly of Actomyosin Contractile Rings in Cytokinesis and Cell Repair

Cytokinesis and single-cell wound repair both involve contractile assemblies of filamentous actin (F-actin) and myosin II organized into characteristic ring-like arrays. The assembly of these actomyosin contractile rings (CRs) is specified spatially and temporally by small Rho GTPases, which trigger local actin polymerization and myosin II contractility via a variety of downstream effectors. We now have a much clearer view of the Rho GTPase signaling cascade that leads to the formation of CRs, but some factors involved in CR positioning, assembly, and function remain poorly understood. Recent studies show that this regulation is multifactorial and goes beyond the long-established Ca²⁺ -dependent processes. There is substantial evidence that the Ca²⁺ -independent changes in cell shape, tension, and plasma membrane composition that characterize cytokinesis and single-cell wound repair also regulate CR formation. Elucidating the regulation and mechanistic properties of CRs is important to our understanding of basic cell biology and holds potential for therapeutic applications in human disease. In this review, we present a primer on the factors influencing and regulating CR positioning, assembly, and contraction as they occur in a variety of cytokinetic and single-cell wound repair models. Anat Rec, 301:2051-2066, 2018. © 2018 Wiley Periodicals, Inc.

Cellular mechanisms and signals that coordinate plasma membrane repair

Plasma membrane forms the barrier between the cytoplasm and the environment. Cells constantly and selectively transport molecules across their plasma membrane without disrupting it. Any disruption in the plasma membrane compromises its selective permeability and is lethal, if not rapidly repaired. There is a growing understanding of the organelles, proteins, lipids, and small molecules that help cells signal and efficiently coordinate plasma membrane repair. This review aims to summarize how these subcellular responses are coordinated and how cellular signals generated due to plasma membrane injury interact with each other to spatially and temporally coordinate repair. With the involvement of calcium and redox signaling in single cell and tissue repair, we will discuss how these and other related signals extend from single cell repair to tissue level repair. These signals link repair processes that are activated immediately after plasma membrane injury with longer term processes regulating repair and regeneration of the damaged tissue. We propose that investigating cell and tissue repair as part of a continuum of wound repair mechanisms would be of value in treating degenerative diseases.

A Rho-based reaction-diffusion system governs cell wall patterning in metaxylem vessels

Rho GTPases play crucial roles in cell polarity and pattern formation. ROPs, Rho of plant GTPases, are widely involved in cell wall patterning in plants, yet the molecular mechanism underlying their action remains unknown. Arabidopsis ROP11 is locally activated to form plasma membrane domains, which direct formation of cell wall pits in metaxylem vessel cells through interaction with cortical microtubules. Here, we show that the pattern formation of cell wall pits is governed by ROP activation via a reaction-diffusion mechanism. Genetic analysis and reconstructive assays revealed that ROPGEF4/7 and ROPGAP3/4, which encode ROP activators and inactivators, respectively, regulated the formation of ROP-activated domains; these in turn determined the pattern of cell wall pits. Mathematical modelling showed that ROP-activation cycle generated ROP domains by reaction-diffusion mechanism. The model predicted that a positive feedback and slow diffusion of ROP11-ROPGEF4 complex were required to generate ROP-activated domains. ROPGEF4 formed a dimer that interacted with activated ROP11 in vivo, which could provide positive feedback for ROP activation. ROPGEF4 was highly stable on the plasma membrane and inhibited ROP11 diffusion. Our study indicated that ROP-based reaction-diffusion system self-organizes ROP-activated domains, thereby determines the pit pattern of metaxylem vessels.

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.

Automated and Robust Quantification of Colocalization in Dual-Color Fluorescence Microscopy: A Nonparametric Statistical Approach

Colocalization is a powerful tool to study the interactions between fluorescently labeled molecules in biological fluorescence microscopy. However, existing techniques for colocalization analysis have not undergone continued development especially in regards to robust statistical support. In this paper, we examine two of the most popular quantification techniques for colocalization and argue that they could be improved upon using ideas from nonparametric statistics and scan statistics. In particular, we propose a new colocalization metric that is robust, easily implementable, and optimal in a rigorous statistical testing framework. Application to several benchmark data sets, as well as biological examples, further demonstrates the usefulness of the proposed technique.

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.

Host cell perforation by listeriolysin O (LLO) activates a Ca2+-dependent cPKC/Rac1/Arp2/3 signaling pathway that promotes Listeria monocytogenes internalization independently of membrane resealing

Host cell invasion is an indispensable step for a successful infection by intracellular pathogens. Recent studies identified pathogen-induced host cell plasma membrane perforation as a novel mechanism used by diverse pathogens (Trypanosoma cruzi, Listeria monocytogenes, and adenovirus) to promote their internalization into target cells. It was concluded that T. cruzi and adenovirus damage the host cell plasma membrane to hijack the endocytic-dependent membrane resealing machinery, thereby invading the host cell. We studied L. monocytogenes and its secreted pore-forming toxin listeriolysin O (LLO) to identify key signaling events activated upon plasma membrane perforation that lead to bacterial internalization. Using various approaches, including fluorescence resonance energy transfer imaging, we found that the influx of extracellular Ca²⁺ subsequent to LLO-mediated plasma membrane perforation is required for the activation of a conventional protein kinase C (cPKC). cPKC is positioned upstream of Rac1 and the Arp2/3 complex, which activation leads to F-actin--dependent bacterial internalization. Inhibition of this pathway did not prevent membrane resealing, revealing that perforation-dependent L. monocytogenes endocytosis is distinct from the resealing machinery. These studies identified the LLO-dependent endocytic pathway of L. monocytogenes and support a novel model for pathogen uptake promoted by plasma membrane injury that is independent of membrane resealing.

Antagonistic Behaviors of NMY-1 and NMY-2 Maintain Ring Channels in the C. elegans Gonad

Contractile rings play critical roles in a number of biological processes, including oogenesis, wound healing, and cytokinesis. In many cases, the activity of motor proteins such as nonmuscle myosins is required for appropriate constriction of these contractile rings. In the gonad of the nematode worm Caenorhabditis elegans, ring channels are a specialized form of contractile ring that are maintained at a constant diameter before oogenesis. We propose a model of ring channel maintenance that explicitly incorporates force generation by motor proteins that can act normally or tangentially to the ring channel opening. We find that both modes of force generation are needed to maintain the ring channels. We demonstrate experimentally that the type II myosins NMY-1 and NMY-2 antagonize each other in the ring channels by producing force in perpendicular directions: the experimental depletion of NMY-1/theoretical decrease in orthogonal force allows premature ring constriction and cellularization, whereas the experimental depletion of NMY-2/theoretical decrease in tangential force opens the ring channels and prevents cellularization. Together, our experimental and theoretical results show that both forces, mediated by NMY-1 and NMY-2, are crucial for maintaining the appropriate ring channel diameter and dynamics throughout the gonad.

Molecular perturbation strategies to examine spatiotemporal features of Rho GEF and Rho GTPase activity in living cells

Much of our current knowledge of Rho GTPase networks and the regulation by Rho guanine exchange factors (Rho GEFs) and Rho GTPase activating proteins (Rho GAPs) is based on population-based techniques. Over the last decades, technologies that enable single cell analysis with high spatial and temporal resolution have revealed that Rho GTPase activity in cells is regulated on second timescales and at submicrometer length scales. Therefore, perturbation methods with matching spatial and temporal resolution are crucial to further our understanding of Rho GTPase signaling. Here, we give a brief overview of the components of Rho GTPase signaling networks and review a range of existing perturbation strategies that target a specific component of the Rho GTPase signaling module. The advantages and limitations of each perturbation method are discussed. Several recommendations are formulated to guide future studies aimed at addressing spatiotemporal aspects of Rho GEF and Rho GTPase signaling.

The cellular and molecular mechanisms of tissue repair and regeneration as revealed by studies in Xenopus

Survival of any living organism critically depends on its ability to repair and regenerate damaged tissues and/or organs during its lifetime following injury, disease, or aging. Various animal models from invertebrates to vertebrates have been used to investigate the molecular and cellular mechanisms of wound healing and tissue regeneration. It is hoped that such studies will form the framework for identifying novel clinical treatments that will improve the healing and regenerative capacity of humans. Amongst these models, Xenopus stands out as a particularly versatile and powerful system. This review summarizes recent findings using this model, which have provided fundamental knowledge of the mechanisms responsible for efficient and perfect tissue repair and regeneration.

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.

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.

Paradoxical signaling regulates structural plasticity in dendritic spines

Transient spine enlargement (3- to 5-min timescale) is an important event associated with the structural plasticity of dendritic spines. Many of the molecular mechanisms associated with transient spine enlargement have been identified experimentally. Here, we use a systems biology approach to construct a mathematical model of biochemical signaling and actin-mediated transient spine expansion in response to calcium influx caused by NMDA receptor activation. We have identified that a key feature of this signaling network is the paradoxical signaling loop. Paradoxical components act bifunctionally in signaling networks, and their role is to control both the activation and the inhibition of a desired response function (protein activity or spine volume). Using ordinary differential equation (ODE)-based modeling, we show that the dynamics of different regulators of transient spine expansion, including calmodulin-dependent protein kinase II (CaMKII), RhoA, and Cdc42, and the spine volume can be described using paradoxical signaling loops. Our model is able to capture the experimentally observed dynamics of transient spine volume. Furthermore, we show that actin remodeling events provide a robustness to spine volume dynamics. We also generate experimentally testable predictions about the role of different components and parameters of the network on spine dynamics.

The dynamics of spatio-temporal Rho GTPase signaling: formation of signaling patterns

Rho GTPases are crucial signaling molecules that regulate a plethora of biological functions. Traditional biochemical, cell biological, and genetic approaches have founded the basis of Rho GTPase biology. The development of biosensors then allowed measuring Rho GTPase activity with unprecedented spatio-temporal resolution. This revealed that Rho GTPase activity fluctuates on time and length scales of tens of seconds and micrometers, respectively. In this review, we describe Rho GTPase activity patterns observed in different cell systems. We then discuss the growing body of evidence that upstream regulators such as guanine nucleotide exchange factors and GTPase-activating proteins shape these patterns by precisely controlling the spatio-temporal flux of Rho GTPase activity. Finally, we comment on additional mechanisms that might feed into the regulation of these signaling patterns and on novel technologies required to dissect this spatio-temporal complexity.

Modeling the roles of protein kinase Cβ and η in single-cell wound repair

Wounded cells such as Xenopus oocytes respond to damage by assembly and closure of an array of actin filaments and myosin-2 controlled by Rho GTPases, including Rho and Cdc42. Rho and Cdc42 are patterned around wounds in a characteristic manner, with active Rho concentrating in a ring-like zone inside a larger, ring-like zone of active Cdc42. How this patterning is achieved is unknown, but Rho and Cdc42 at wounds are subject to regulation by other proteins, including the protein kinases C. Specifically, Cdc42 and Rho activity are enhanced by PKCβ and inhibited by PKCη. We adapt a mathematical model of Simon and coworkers to probe the possible roles of these kinases. We show that PKCβ likely affects the magnitude of positive Rho-Abr feedback, whereas PKCη acts on Cdc42 inactivation. The model explains both qualitative and some overall quantitative features of PKC-Rho GTPase regulation. It also accounts for the previous, peculiar observation that ∼ 20% of cells overexpressing PKCη display zone inversions--that is, displacement of active Rho to the outside of the active Cdc42.

Measurement of extracellular ion fluxes using the ion-selective self-referencing microelectrode technique

Cells from animals, plants and single cells are enclosed by a barrier called the cell membrane that separates the cytoplasm from the outside. Cell layers such as epithelia also form a barrier that separates the inside from the outside or different compartments of multicellular organisms. A key feature of these barriers is the differential distribution of ions across cell membranes or cell layers. Two properties allow this distribution: 1) membranes and epithelia display selective permeability to specific ions; 2) ions are transported through pumps across cell membranes and cell layers. These properties play crucial roles in maintaining tissue physiology and act as signaling cues after damage, during repair, or under pathological condition. The ion-selective self-referencing microelectrode allows measurements of specific fluxes of ions such as calcium, potassium or sodium at single cell and tissue levels. The microelectrode contains an ionophore cocktail which is selectively permeable to a specific ion. The internal filling solution contains a set concentration of the ion of interest. The electric potential of the microelectrode is determined by the outside concentration of the ion. As the ion concentration varies, the potential of the microelectrode changes as a function of the log of the ion activity. When moved back and forth near a source or sink of the ion (i.e. in a concentration gradient due to ion flux) the microelectrode potential fluctuates at an amplitude proportional to the ion flux/gradient. The amplifier amplifies the microelectrode signal and the output is recorded on computer. The ion flux can then be calculated by Fick's law of diffusion using the electrode potential fluctuation, the excursion of microelectrode, and other parameters such as the specific ion mobility. In this paper, we describe in detail the methodology to measure extracellular ion fluxes using the ion-selective self-referencing microelectrode and present some representative results.

Lipid domain-dependent regulation of single-cell wound repair

After damage, cells reseal their plasma membrane and repair the underlying cortical cytoskeleton. Although many different proteins have been implicated in cell repair, the potential role of specific lipids has not been explored. Here we report that cell damage elicits rapid formation of spatially organized lipid domains around the damage site, with different lipids concentrated in different domains as a result of both de novo synthesis and transport. One of these lipids-diacylglycerol (DAG)-rapidly accumulates in a broad domain that overlaps the zones of active Rho and Cdc42, GTPases that regulate repair of the cortical cytoskeleton. Formation of the DAG domain is required for Cdc42 and Rho activation and healing. Two DAG targets, protein kinase C (PKC) β and η, are recruited to cell wounds and play mutually antagonistic roles in the healing process: PKCβ participates in Rho and Cdc42 activation, whereas PKCη inhibits Rho and Cdc42 activation. The results reveal an unexpected diversity in subcellular lipid domains and the importance of such domains for a basic cellular process.

Single cell wound generates electric current circuit and cell membrane potential variations that requires calcium influx

Breaching of the cell membrane is one of the earliest and most common causes of cell injury, tissue damage, and disease. If the compromise in cell membrane is not repaired quickly, irreversible cell damage, cell death and defective organ functions will result. It is therefore fundamentally important to efficiently repair damage to the cell membrane. While the molecular aspects of single cell wound healing are starting to be deciphered, its bio-physical counterpart has been poorly investigated. Using Xenopus laevis oocytes as a model for single cell wound healing, we describe the temporal and spatial dynamics of the wound electric current circuitry and the temporal dynamics of cell membrane potential variation. In addition, we show the role of calcium influx in controlling electric current circuitry and cell membrane potential variations. (i) Upon wounding a single cell: an inward electric current appears at the wound center while an outward electric current is observed at its sides, illustrating the wound electric current circuitry; the cell membrane is depolarized; calcium flows into the cell. (ii) During cell membrane re-sealing: the wound center current density is maintained for a few minutes before decreasing; the cell membrane gradually re-polarizes; calcium flow into the cell drops. (iii) In conclusion, calcium influx is required for the formation and maintenance of the wound electric current circuitry, for cell membrane re-polarization and for wound healing.

Single cell pattern formation and transient cytoskeletal arrays

A major goal of developmental biology is to explain the emergence of pattern in cell layers, tissues and organs. Developmental biologists now accept that reaction diffusion-based mechanisms are broadly employed in developing organisms to direct pattern formation. Here we briefly consider these mechanisms and then apply some of the concepts derived from them to several processes that occur in single cells: wound repair, yeast budding, and cytokinesis. Two conclusions emerge from this analysis: first, there is considerable overlap at the level of general mechanisms between developmental and single cell pattern formation; second, dynamic structures based on the actin cytoskeleton may be far more ordered than is generally recognized.