Phenomenological modeling of durotaxis

Durotaxis and extracellular matrix degradation promote the clustering of cancer cells

Early stages of metastasis depend on the collective behavior of cancer cells and their interaction with the extracellular matrix (ECM). Cancer cell clusters are known to exhibit higher metastatic potential than single cells. To explore clustering dynamics, we developed a calibrated computational model describing how motile cancer cells biochemically and biomechanically interact with the ECM during the initial invasion phase, including ECM degradation and mechanical remodeling. The model reveals that cluster formation time, size, and shape are influenced by ECM degradation rates and cellular compliance to external stresses (durotaxis). The results align with experimental observations, demonstrating distinct cell trajectories and cluster morphologies shaped by biomechanical parameters. The simulations provide valuable insights into cancer invasion dynamics and may suggest potential therapeutic strategies targeting early-stage invasive cells.

Elastic interactions compete with persistent cell motility to drive durotaxis

Many animal cells that crawl on extracellular substrates exhibit durotaxis, i.e., directed migration toward stiffer substrate regions. This has implications in several biological processes including tissue development and tumor progression. Here, we introduce a phenomenological model for single-cell durotaxis that incorporates both elastic deformation-mediated cell-substrate interactions and the stochasticity of cell migration. Our model is motivated by a key observation in an early demonstration of durotaxis: a single, contractile cell at a sharp interface between a softer and a stiffer region of an elastic substrate reorients and migrates toward the stiffer region. We model migrating cells as self-propelling, persistently motile agents that exert contractile traction forces on their elastic substrate. The resulting substrate deformations induce elastic interactions with mechanical boundaries, captured by an elastic potential. The dynamics is determined by two crucial parameters: the strength of the cellular traction-induced boundary elastic interaction (A), and the persistence of cell motility (Pe). Elastic forces and torques resulting from the potential orient cells perpendicular (parallel) to the boundary and accumulate (deplete) them at the clamped (free) boundary. Thus, a clamped boundary induces an attractive potential that drives durotaxis, while a free boundary induces a repulsive potential that prevents antidurotaxis. By quantifying the steady-state position and orientation probability densities, we show how the extent of accumulation (depletion) depends on the strength of the elastic potential and motility. We compare and contrast crawling cells with biological microswimmers and other synthetic active particles, where accumulation at confining boundaries is well known. We define metrics quantifying boundary accumulation and durotaxis, and present a phase diagram that identifies three possible regimes: durotaxis, and adurotaxis with and without motility-induced accumulation at the boundary. Overall, our model predicts how durotaxis depends on cell contractility and motility, successfully explains some previous observations, and provides testable predictions to guide future experiments.

Positive and negative durotaxis - mechanisms and emerging concepts

Cell migration is controlled by the coordinated action of cell adhesion, cytoskeletal dynamics, contractility and cell extrinsic cues. Integrins are the main adhesion receptors to ligands of the extracellular matrix (ECM), linking the actin cytoskeleton to the ECM and enabling cells to sense matrix rigidity and mount a directional cell migration response to stiffness gradients. Most models studied show preferred migration of single cells or cell clusters towards increasing rigidity. This is referred to as durotaxis, and since its initial discovery in 2000, technical advances and elegant computational models have provided molecular level details of stiffness sensing in cell migration. However, modeling has long predicted that, depending on cell intrinsic factors, such as the balance of cell adhesion molecules (clutches) and the motor proteins pulling on them, cells might also prefer adhesion to intermediate rigidity. Recently, experimental evidence has supported this notion and demonstrated the ability of cells to migrate towards lower rigidity, in a process called negative durotaxis. In this Review, we discuss the significant conceptual advances that have been made in our appreciation of cell plasticity and context dependency in stiffness-guided directional cell 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.

Three-dimensional cancer cell migration directed by dual mechanochemical guidance

Directed cell migration guided by external cues plays a central role in many physiological and pathophysiological processes. The microenvironment of cells often simultaneously contains various cues and the motility response of cells to multiplexed guidance is poorly understood. Here we combine experiments and mathematical models to study the three-dimensional migration of breast cancer cells in the presence of both contact guidance and a chemoattractant gradient. We find that the chemotaxis of cells is complicated by the presence of contact guidance as the microstructure of extracellular matrix (ECM) vary spatially. In the presence of dual guidance, the impact of ECM alignment is determined externally by the coherence of ECM fibers and internally by cell mechanosensing Rho/Rock pathways. When contact guidance is parallel to the chemical gradient, coherent ECM fibers significantly increase the efficiency of chemotaxis. When contact guidance is perpendicular to the chemical gradient, cells exploit the ECM disorder to locate paths for chemotaxis. Our results underscore the importance of fully characterizing the cancer cell microenvironment in order to better understand invasion and metastasis.

Computational modelling of cell motility modes emerging from cell-matrix adhesion dynamics

Lymphocytes have been described to perform different motility patterns such as Brownian random walks, persistent random walks, and Lévy walks. Depending on the conditions, such as confinement or the distribution of target cells, either Brownian or Lévy walks lead to more efficient interaction with the targets. The diversity of these motility patterns may be explained by an adaptive response to the surrounding extracellular matrix (ECM). Indeed, depending on the ECM composition, lymphocytes either display a floating motility without attaching to the ECM, or sliding and stepping motility with respectively continuous or discontinuous attachment to the ECM, or pivoting behaviour with sustained attachment to the ECM. Moreover, on the long term, lymphocytes either perform a persistent random walk or a Brownian-like movement depending on the ECM composition. How the ECM affects cell motility is still incompletely understood. Here, we integrate essential mechanistic details of the lymphocyte-matrix adhesions and lymphocyte intrinsic cytoskeletal induced cell propulsion into a Cellular Potts model (CPM). We show that the combination of de novo cell-matrix adhesion formation, adhesion growth and shrinkage, adhesion rupture, and feedback of adhesions onto cell propulsion recapitulates multiple lymphocyte behaviours, for different lymphocyte subsets and various substrates. With an increasing attachment area and increased adhesion strength, the cells’ speed and persistence decreases. Additionally, the model predicts random walks with short-term persistent but long-term subdiffusive properties resulting in a pivoting type of motility. For small adhesion areas, the spatial distribution of adhesions emerges as a key factor influencing cell motility. Small adhesions at the front allow for more persistent motility than larger clusters at the back, despite a similar total adhesion area. In conclusion, we present an integrated framework to simulate the effects of ECM proteins on cell-matrix adhesion dynamics. The model reveals a sufficient set of principles explaining the plasticity of lymphocyte motility.

During immunosurveillance, lymphocytes patrol through tissues to interact with cancer cells, other immune cells, and pathogens. The efficiency of this process depends on the kinds of trajectories taken, ranging from simple Brownian walks to Lévy walks. The composition of the extracellular matrix (ECM), a network of macromolecules, affects the formation of cell-matrix adhesions, thus strongly influencing the way lymphocytes move. Here, we present a model of lymphocyte motility driven by adhesions that grow, shrink and rupture in response to the ECM and cellular forces. Compared to other models, our model is computationally light making it suitable for generating long term cell track data, while still capturing actin dynamics and adhesion turnover. Our model suggests that cell motility is affected by the force required to break adhesions and the rate at which new adhesions form. Adhesions can promote cell protrusion by inhibiting retrograde actin flow. After introducing this effect into the model, we found that it reduces the cellular diffusivity and that it promotes stick-slip behaviour. Furthermore, location and size of adhesion clusters determined cell persistence. Overall, our model explains the plasticity of lymphocyte behaviour in response to the ECM.

Durotaxis: The Hard Path from In Vitro to In Vivo

Durotaxis, the process by which cells follow gradients of extracellular mechanical stiffness, has been proposed as a mechanism driving directed migration. Despite the lack of evidence for its existence in vivo, durotaxis has become an active field of research, focusing on the mechanism by which cells respond to mechanical stimuli from the environment. In this review, we describe the technical and conceptual advances in the study of durotaxis in vitro, discuss to what extent the evidence suggests durotaxis may occur in vivo, and emphasize the urgent need for in vivo demonstration of durotaxis.

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.

All Roads Lead to Directional Cell Migration

Directional cell migration normally relies on a variety of external signals, such as chemical, mechanical, or electrical, which instruct cells in which direction to move. Many of the major molecular and physical effects derived from these cues are now understood, leading to questions about whether directional cell migration is alike or distinct under these different signals, and how cells might be directed by multiple simultaneous cues, which would be expected in complex in vivo environments. In this review, we compare how different stimuli are spatially distributed, often as gradients, to direct cell movement and the mechanisms by which they steer cells. A comparison of the downstream effectors of directional cues suggests that different external signals regulate a common set of components: small GTPases and the actin cytoskeleton, which implies that the mechanisms downstream of different signals are likely to be closely related and underlies the idea that cell migration operates by a common set of physical principles, irrespective of the input.

Topotaxis of active Brownian particles

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

Force chains in cell-cell mechanical communication

Force chains (FCs) are a key determinant of the micromechanical properties and behaviour of heterogeneous materials, such as granular systems. However, less is known about FCs in fibrous materials, such as the networks composing the extracellular matrix (ECM) of biological systems. Using a finite-element computational model, we simulated the contraction of a single cell and two nearby cells embedded in two-dimensional fibrous elastic networks and analysed the tensile FCs that developed in the ECM. The role of ECM nonlinear elasticity on FC formation was evaluated by considering linear and nonlinear, i.e. exhibiting 'buckling' and/or 'strain-stiffening', stress-strain curves. The effect of the degree of cell contraction and network coordination value was assessed. We found that nonlinear elasticity of the ECM fibres influenced the structure of the FCs, facilitating the transition towards more distinct chains that were less branched and more radially oriented than the chains formed in linear elastic networks. When two neighbouring cells contract, a larger number of FCs bridged between the cells in nonlinear networks, and these chains had a larger effective rigidity than the chains that did not reach a neighbouring cell. These results suggest that FCs function as a route for mechanical communication between distant cells and highlight the contribution of ECM fibre nonlinear elasticity to the formation of FCs.

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 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.