A Cellular Potts Model of single cell migration in presence of durotaxis

Piezo1-induced durotaxis of pancreatic stellate cells depends on TRPC1 and TRPV4 channels

Pancreatic stellate cells (PSCs) are primarily responsible for producing the stiff tumor tissue in pancreatic ductal adenocarcinoma (PDAC). Thereby, PSCs generate a stiffness gradient between the healthy pancreas and the tumor. This gradient induces durotaxis, a form of directional cell migration driven by differential stiffness. The molecular sensors behind durotaxis are still unclear. To investigate the role of mechanosensitive ion channels in PSC durotaxis, we established a two-dimensional stiffness gradient mimicking PDAC. Using pharmacological and genetic methods, we investigated the role of the ion channels Piezo1, TRPC1, and TRPV4 in PSC durotaxis. We found that PSC migration towards a stiffer substrate is diminished by altering Piezo1 activity. Moreover, disrupting TRPC1 along with TRPV4 abolishes PSC durotaxis even when Piezo1 is functional. Hence, PSC durotaxis is optimal with an intermediary level of mechanosensitive channel activity, which we simulated using a numerically discretized mathematical model. Our findings suggest that mechanosensitive ion channels, particularly Piezo1, detect the mechanical microenvironment to guide PSC migration.

A non local model for cell migration in response to mechanical stimuli

Cell migration is one of the most studied phenomena in biology since it plays a fundamental role in many physiological and pathological processes such as morphogenesis, wound healing and tumorigenesis. In recent years, researchers have performed experiments showing that cells can migrate in response to mechanical stimuli of the substrate they adhere to. Motion towards regions of the substrate with higher stiffness is called durotaxis, while motion guided by the stress or the deformation of the substrate itself is called tensotaxis. Unlike chemotaxis (i.e. the motion in response to a chemical stimulus), these migratory processes are not yet fully understood from a biological point of view. In this respect, we present a mathematical model of single-cell migration in response to mechanical stimuli, in order to simulate these two processes. Specifically, the cell moves by changing its direction of polarization and its motility according to material properties of the substrate (e.g., stiffness) or in response to proper scalar measures of the substrate strain or stress. The equations of motion of the cell are non-local integro-differential equations, with the addition of a stochastic term to account for random Brownian motion. The mechanical stimulus to be integrated in the equations of motion is defined according to experimental measurements found in literature, in the case of durotaxis. Conversely, in the case of tensotaxis, substrate strain and stress are given by the solution of the mechanical problem, assuming that the extracellular matrix behaves as a hyperelastic Yeoh's solid. In both cases, the proposed model is validated through numerical simulations that qualitatively reproduce different experimental scenarios.

Mathematical model of mechano-sensing and mechanically induced collective motility of cells on planar elastic substrates

Cells mechanically interact with their environment to sense, for example, topography, elasticity and mechanical cues from other cells. Mechano-sensing has profound effects on cellular behaviour, including motility. The current study aims to develop a mathematical model of cellular mechano-sensing on planar elastic substrates and demonstrate the model's predictive capabilities for the motility of individual cells in a colony. In the model, a cell is assumed to transmit an adhesion force, derived from a dynamic focal adhesion integrin density, that locally deforms a substrate, and to sense substrate deformation originating from neighbouring cells. The substrate deformation from multiple cells is expressed as total strain energy density with a spatially varying gradient. The magnitude and direction of the gradient at the cell location define the cell motion. Cell-substrate friction, partial motion randomness, and cell death and division are included. The substrate deformation by a single cell and the motility of two cells are presented for several substrate elasticities and thicknesses. The collective motility of 25 cells on a uniform substrate mimicking the closure of a circular wound of 200 µm is predicted for deterministic and random motion. Cell motility on substrates with varying elasticity and thickness is explored for four cells and 15 cells, the latter again mimicking wound closure. Wound closure by 45 cells is used to demonstrate the simulation of cell death and division during migration. The mathematical model can adequately simulate the mechanically induced collective cell motility on planar elastic substrates. The model is suitable for extension to other cell and substrates shapes and the inclusion of chemotactic cues, offering the potential to complement in vitro and in vivo studies.

Matrix adhesion and remodeling diversifies modes of cancer invasion across spatial scales

The metastasis of malignant epithelial tumors begins with the egress of transformed cells from the confines of their basement membrane (BM) to their surrounding collagen-rich stroma. Invasion can be morphologically diverse: when breast cancer cells are separately cultured within BM-like matrix, collagen I (Coll I), or a combination of both, they exhibit collective-, dispersed mesenchymal-, and a mixed collective-dispersed (multimodal)- invasion, respectively. In this paper, we asked how distinct these invasive modes are with respect to the cellular and microenvironmental cues that drive them. A rigorous computational exploration of invasion was performed within an experimentally motivated Cellular Potts-based modeling environment. The model comprised of adhesive interactions between cancer cells, BM- and Coll I-like extracellular matrix (ECM), and reaction-diffusion-based remodeling of ECM. The model outputs were parameters cognate to dispersed- and collective- invasion. A clustering analysis of the output distribution curated through a careful examination of subsumed phenotypes suggested at least four distinct invasive states: dispersed, papillary-collective, bulk-collective, and multimodal, in addition to an indolent/non-invasive state. Mapping input values to specific output clusters suggested that each of these invasive states are specified by distinct input signatures of proliferation, adhesion and ECM remodeling. In addition, specific input perturbations allowed transitions between the clusters and revealed the variation in the robustness between the invasive states. Our systems-level approach proffers quantitative insights into how the diversity in ECM microenvironments may steer invasion into diverse phenotypic modes during early dissemination of breast cancer and contributes to tumor heterogeneity.

Enhanced persistence and collective migration in cooperatively aligning cell clusters

Most cells possess the capacity to locomote. Alone or collectively, this allows them to adapt, to rearrange, and to explore their surroundings. The biophysical characterization of such motile processes, in health and in disease, has so far focused mostly on two limiting cases: single-cell motility on the one hand and the dynamics of confluent tissues such as the epithelium on the other. The in-between regime of clusters, composed of relatively few cells moving as a coherent unit, has received less attention. Such small clusters are, however, deeply relevant in development but also in cancer metastasis. In this work, we use cellular Potts models and analytical active matter theory to understand how the motility of small cell clusters changes with N, the number of cells in the cluster. Modeling and theory reveal our two main findings: cluster persistence time increases with N, whereas the intrinsic diffusivity decreases with N. We discuss a number of settings in which the motile properties of more complex clusters can be analytically understood, revealing that the focusing effects of small-scale cooperation and cell-cell alignment can overcome the increased bulkiness and internal disorder of multicellular clusters to enhance overall migrational efficacy. We demonstrate this enhancement for small-cluster collective durotaxis, which is shown to proceed more effectively than for single cells. Our results may provide some novel, to our knowledge, insights into the connection between single-cell and large-scale collective motion and may point the way to the biophysical origins of the enhanced metastatic potential of small tumor cell clusters.

Cell Shape and Durotaxis Explained from Cell-Extracellular Matrix Forces and Focal Adhesion Dynamics

Many cells are small and rounded on soft extracellular matrices (ECM), elongated on stiffer ECMs, and flattened on hard ECMs. Cells also migrate up stiffness gradients (durotaxis). Using a hybrid cellular Potts and finite-element model extended with ODE-based models of focal adhesion (FA) turnover, we show that the full range of cell shape and durotaxis can be explained in unison from dynamics of FAs, in contrast to previous mathematical models. In our 2D cell-shape model, FAs grow due to cell traction forces. Forces develop faster on stiff ECMs, causing FAs to stabilize and, consequently, cells to spread on stiff ECMs. If ECM stress further stabilizes FAs, cells elongate on substrates of intermediate stiffness. We show that durotaxis follows from the same set of assumptions. Our model contributes to the understanding of the basic responses of cells to ECM stiffness, paving the way for future modeling of more complex cell-ECM interactions.

A Cellular Potts energy-based approach to analyse the influence of the surface topography on single cell motility

The surface shape is an important aspect to take into account to ensure the success of an implant. At the cellular scale level, the cell behaviour, especially its migration, is affected by the specificities of the surface of the substrate, such as the stiffness of the surface and its roughness topography. The latter has been shown to have a great impact on various cell mechanisms, such as the cell adhesion, migration, or proliferation. In fact, the mere presence of micro roughness leads to an improvement of those mechanisms, with a better integration of the implants. However, the phenomena behind those improvements are still not clear. In this paper, we propose a three-dimensional (3D) model of a single cell migration using a Cellular Potts (CP) model to study the influence of the surface topography on cell motility. To do so, various configurations were tested, such as: (i) a substrate with a random roughness, (ii) a substrate with a rectangular groove pattern (parallel and perpendicular to the direction of motion), (ii) a substrate with a sinusoidal groove pattern. To evaluate the influence of the surface topography on cell motility, for each configuration, the cell speed and shape as well as the contact surface between the cell and the substrate have been quantified. Our numerical results demonstrate that, in agreement with the experimental observations of the literature, the substrate topography has an influence on the cell efficiency (i.e. cell speed), orientation and shape. Besides, we also show that the increase of the contact surface alone in presence of roughness is not enough to explain the improvement of cell migration on the various rough surfaces. Finally, we highlight the importance of the roughness dimension on cell motility. This could be a critical aspect to consider for further analyses and applications, such as surface treatments for medical applications.

The Impact of Elastic Deformations of the Extracellular Matrix on Cell Migration

The mechanical properties of the extracellular matrix, in particular its stiffness, are known to impact cell migration. In this paper, we develop a mathematical model of a single cell migrating on an elastic matrix, which accounts for the deformation of the matrix induced by forces exerted by the cell, and investigate how the stiffness impacts the direction and speed of migration. We model a cell in 1D as a nucleus connected to a number of adhesion sites through elastic springs. The cell migrates by randomly updating the position of its adhesion sites. We start by investigating the case where the cell springs are constant, and then go on to assuming that they depend on the matrix stiffness, on matrices of both uniform stiffness as well as those with a stiffness gradient. We find that the assumption that cell springs depend on the substrate stiffness is necessary and sufficient for an efficient durotactic response. We compare simulations to recent experimental observations of human cancer cells exhibiting durotaxis, which show good qualitative agreement.

Mechanisms of endothelial cell coverage by pericytes: computational modelling of cell wrapping and in vitro experiments

Pericytes (PCs) wrap around endothelial cells (ECs) and perform diverse functions in physiological and pathological processes. Although molecular interactions between ECs and PCs have been extensively studied, the morphological processes at the cellular level and their underlying mechanisms have remained elusive. In this study, using a simple cellular Potts model, we explored the mechanisms for EC wrapping by PCs. Based on the observed in vitro cell wrapping in three-dimensional PC-EC coculture, the model identified four putative contributing factors: preferential adhesion of PCs to the extracellular matrix (ECM), strong cell-cell adhesion, PC surface softness and larger PC size. While cell-cell adhesion can contribute to the prevention of cell segregation and the degree of cell wrapping, it cannot determine the orientation of cell wrapping alone. While atomic force microscopy revealed that PCs have a larger Young's modulus than ECs, the experimental analyses supported preferential ECM adhesion and size asymmetry. We also formulated the corresponding energy minimization problem and numerically solved this problem for specific cases. These results give biological insights into the role of PC-ECM adhesion in PC coverage. The modelling framework presented here should also be applicable to other cell wrapping phenomena observed in vivo.

Cell migration driven by substrate deformation gradients

Identifying the cues followed by cells is key to understand processes as embryonic development, tissue homeostasis, or several pathological conditions. Based on a durotaxis model, it is shown that cells moving on predeformed thin elastic membrane follow the direction of increasing strain of the substrate. This mechanism, straintaxis, does not distinguish the origin of the strain, but the active stresses produce large strains on cells or tissues being used as substrates. Hence, straintaxis is the natural realization of duratoaxis in vivo. Considering a circular geometry for the substrate cells, it is shown that if the annular component of the active stress component increases with the radial distance, cells migrate toward the substrate cell borders. With appropriate estimation for the different parameters, the migration speeds are similar to those obtained in recent experiments (Reig et al 2017 Nat. Commun. 8 15431). In these, during the annual killifish epiboly, deep cells that move in contact with the epithelial enveloping cell layer (EVL), migrate toward the EVL cell borders with speeds of microns per minute.

Substrate-rigidity dependent migration of an idealized twitching bacterium

Mechanical properties of the extracellular matrix are important determinants of cellular migration in diverse processes, such as immune response, wound healing, and cancer metastasis. Moreover, recent studies indicate that even bacterial surface colonization can depend on the mechanics of the substrate. Here, we focus on physical mechanisms that can give rise to substrate-rigidity dependent migration. We study a "twitcher", a cell driven by extension-retraction cycles, to idealize bacteria and perhaps eukaryotic cells that employ a slip-stick mode of motion. The twitcher is asymmetric and always pulls itself forward at its front. Analytical calculations show that the migration speed of a twitcher depends non-linearly on substrate rigidity. For soft substrates, deformations do not lead to build-up of significant force and the migration speed is therefore determined by stochastic adhesion unbinding. For rigid substrates, forced adhesion rupture determines the migration speed. Depending on the force-sensitivity of front and rear adhesions, forced bond rupture implies an increase or a decrease of the migration speed. A requirement for the occurrence of rigidity-dependent stick-slip migration is a "sticky" substrate, with binding rates being an order of magnitude larger than unbinding rates in absence of force. Computer simulations show that small stall forces of the driving machinery lead to a reduced movement on high rigidities, regardless of force-sensitivities of bonds. The simulations also confirm the occurrence of rigidity-dependent migration speed in a generic model for slip-stick migration of cells on a sticky substrate.

Mechanical Model for Durotactic Cell Migration

Cell migration is a fundamental process in biological systems, playing an important role for diverse physiological processes. Cells often exhibit directed migration in a specific direction in response to various types of cues. In particular, cells are able to sense the rigidity of surrounding environments and then migrate toward stiffer regions. To understand this mechanosensitive behavior called durotaxis, several computational models have been developed. However, most of the models employed cell decision making to recapitulate durotactic behaviors, significantly limiting insights provided from these studies. In this study, we developed a computational biomechanical model without any cell decision making to illuminate intrinsic mechanisms of durotactic behaviors of cells migrating on a two-dimensional substrate. The model consists of a simplified cell generating contractile forces and a deformable substrate coarse-grained into an irregular triangulated mesh. Using the model, we demonstrated that durotactic behaviors emerge from purely mechanical interactions between the cell and the underlying substrate. We investigated how durotactic migration is regulated by biophysical properties of the substrate, including elasticity, viscosity, and stiffness profile.

A Biophysical Model for Curvature-Guided Cell Migration

The latest experiments have shown that adherent cells can migrate according to cell-scale curvature variations via a process called curvotaxis. Despite identification of key cellular factors, a clear understanding of the mechanism is lacking. We employ a mechanical model featuring a detailed description of the cytoskeleton filament networks, the viscous cytosol, the cell adhesion dynamics, and the nucleus. We simulate cell adhesion and migration on sinusoidal substrates. We show that cell adhesion on three-dimensional curvatures induces a gradient of pressure inside the cell that triggers the internal motion of the nucleus. We propose that the resulting out-of-equilibrium position of the nucleus alters cell migration directionality, leading to cell motility toward concave regions of the substrate, resulting in lower potential energy states. Altogether, we propose a simple mechanism explaining how intracellular mechanics enable the cells to react to substratum curvature, induce a deterministic cell polarization, and break down cells basic persistent random walk, which correlates with latest experimental evidences.

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.

An investigation of the distribution and location of mast cells affected by the stiffness of substrates as a mechanical niche

The distribution and location of mast cells are closely related to their physiological and pathological functions, such as allergic responses, immunity, and fibrosis, and are used in acupuncture. In this study, the distribution of mast cells in vivo was observed, and mechanical clues for understanding their distribution based on mechanical niches were explored. By toluidine blue staining and immunohistochemical staining, we examined the distribution and location of mast cells in rat skin and found that mast cells are distributed in a spatially nonuniform manner, preferring to locate at regions in the tissue and extracellular matrix with stiffness changes. In vitro experiments for studying the distribution of rat basophilic leukemia (RBL-2H3) mast cell line on poly-di-methyl-siloxane (PDMS) substrates with stiffness variations were performed. It was found that RBL-2H3 cells migrate and tend to remain in the areas with stiffness variations. The present research suggests that changing the stiffness of local tissues may stimulate mast cell recruitment, which may be the method by which some traditional Chinese medicine treatments, such as acupuncture. On the basis of the origin of mast cells and our experimental results, we predict that mast cells exist in tissues that contain permeable capillaries and prefer regions with stiffness changes. We discussed this prediction using examples of specific tissues from some cases.

A minimal mechanics model for mechanosensing of substrate rigidity gradient in durotaxis

Durotaxis refers to the phenomenon in which cells can sense the spatial gradient of the substrate rigidity in the process of cell migration. A conceptual two-part theory consisting of the focal adhesion force generation and mechanotransduction has been proposed previously by Lo et al. to explain the mechanism underlying durotaxis. In the present work, we are concerned with the first part of the theory: how exactly is the larger focal adhesion force generated in the part of the cell adhering to the stiffer region of the substrate? Using a simple elasticity model and by assuming the cell adheres to the substrate continuously underneath the whole cell body, we show that the mechanics principle of static equilibrium alone is sufficient to account for the generation of the larger traction stress on the stiffer region of the substrate. We believe that our model presents a simple mechanistic understanding of mechanosensing of substrate stiffness gradient at the cellular scale, which can be incorporated in more sophisticated mechanobiochemical models to address complex problems in mechanobiology and bioengineering.

Cell migration in microengineered tumor environments

Recent advances in microengineered cell migration platforms are discussed critically with a focus on how cell migration is influenced by engineered tumor microenvironments, the medical relevance being to understand how tumor microenvironments may promote or suppress the progression of cancer. We first introduce key findings in cancer cell migration under the influence of the physical environment, which is systematically controlled by microengineering technology, followed by multi-cues of physico-chemical factors, which represent the complexity of the tumor environment. Recognizing that cancer cells constantly communicate not only with each other but also with tumor-associated cells such as vascular, fibroblast, and immune cells, and also with non-cellular components, it follows that cell motility in tumor microenvironments, especially metastasis via the invasion of cancer cells into the extracellular matrix and other tissues, is closely related to the malignancy of cancer-related mortality. Medical relevance of forefront research realized in microfabricated devices, such as single cell sorting based on the analysis of cell migration behavior, may assist personalized theragnostics based on the cell migration phenotype. Furthermore, we urge development of theory and numerical understanding of single or collective cell migration in microengineered platforms to gain new insights in cancer metastasis and in therapeutic strategies.