Modeling the influence of nucleus elasticity on cell invasion in fiber networks and microchannels

Learning dynamical models of single and collective cell migration: a review

Single and collective cell migration are fundamental processes critical for physiological phenomena ranging from embryonic development and immune response to wound healing and cancer metastasis. To understand cell migration from a physical perspective, a broad variety of models for the underlying physical mechanisms that govern cell motility have been developed. A key challenge in the development of such models is how to connect them to experimental observations, which often exhibit complex stochastic behaviours. In this review, we discuss recent advances in data-driven theoretical approaches that directly connect with experimental data to infer dynamical models of stochastic cell migration. Leveraging advances in nanofabrication, image analysis, and tracking technology, experimental studies now provide unprecedented large datasets on cellular dynamics. In parallel, theoretical efforts have been directed towards integrating such datasets into physical models from the single cell to the tissue scale with the aim of conceptualising the emergent behaviour of cells. We first review how this inference problem has been addressed in both freely migrating and confined cells. Next, we discuss why these dynamics typically take the form of underdamped stochastic equations of motion, and how such equations can be inferred from data. We then review applications of data-driven inference and machine learning approaches to heterogeneity in cell behaviour, subcellular degrees of freedom, and to the collective dynamics of multicellular systems. Across these applications, we emphasise how data-driven methods can be integrated with physical active matter models of migrating cells, and help reveal how underlying molecular mechanisms control cell behaviour. Together, these data-driven approaches are a promising avenue for building physical models of cell migration directly from experimental data, and for providing conceptual links between different length-scales of description.

The Influence of Nucleus Mechanics in Modelling Adhesion-independent Cell Migration in Structured and Confined Environments

Recent biological experiments (Lämmermann et al. in Nature 453(7191):51-55, 2008; Reversat et al. in Nature 7813:582-585, 2020; Balzer et al. in ASEB J Off Publ Fed Am Soc Exp Biol 26(10):4045-4056, 2012) have shown that certain types of cells are able to move in structured and confined environments even without the activation of focal adhesion. Focusing on this particular phenomenon and based on previous works (Jankowiak et al. in Math Models Methods Appl Sci 30(03):513-537, 2020), we derive a novel two-dimensional mechanical model, which relies on the following physical ingredients: the asymmetrical renewal of the actin cortex supporting the membrane, resulting in a backward flow of material; the mechanical description of the nuclear membrane and the inner nuclear material; the microtubule network guiding nucleus location; the contact interactions between the cell and the external environment. The resulting fourth order system of partial differential equations is then solved numerically to conduct a study of the qualitative effects of the model parameters, mainly those governing the mechanical properties of the nucleus and the geometry of the confining structure. Coherently with biological observations, we find that cells characterized by a stiff nucleus are unable to migrate in channels that can be crossed by cells with a softer nucleus. Regarding the geometry, cell velocity and ability to migrate are influenced by the width of the channel and the wavelength of the external structure. Even though still preliminary, these results may be potentially useful in determining the physical limit of cell migration in confined environments and in designing scaffolds for tissue engineering.

Calculation of the force field required for nucleus deformation during cell migration through constrictions

During cell migration in confinement, the nucleus has to deform for a cell to pass through small constrictions. Such nuclear deformations require significant forces. A direct experimental measure of the deformation force field is extremely challenging. However, experimental images of nuclear shape are relatively easy to obtain. Therefore, here we present a method to calculate predictions of the deformation force field based purely on analysis of experimental images of nuclei before and after deformation. Such an inverse calculation is technically non-trivial and relies on a mechanical model for the nucleus. Here we compare two simple continuum elastic models of a cell nucleus undergoing deformation. In the first, we treat the nucleus as a homogeneous elastic solid and, in the second, as an elastic shell. For each of these models we calculate the force field required to produce the deformation given by experimental images of nuclei in dendritic cells migrating in microchannels with constrictions of controlled dimensions. These microfabricated channels provide a simplified confined environment mimicking that experienced by cells in tissues. Our calculations predict the forces felt by a deforming nucleus as a migrating cell encounters a constriction. Since a direct experimental measure of the deformation force field is very challenging and has not yet been achieved, our numerical approaches can make important predictions motivating further experiments, even though all the parameters are not yet available. We demonstrate the power of our method by showing how it predicts lateral forces corresponding to actin polymerisation around the nucleus, providing evidence for actin generated forces squeezing the sides of the nucleus as it enters a constriction. In addition, the algorithm we have developed could be adapted to analyse experimental images of deformation in other situations.

Many cell types are able to migrate and squeeze through constrictions that are narrower than the cell’s resting radius. For example, both immune cells and metastatic cancer cells change their shape to migrate through small holes in the complex tissue media they move in. During migration the cell nucleus is more difficult to deform than the cell cytoplasm and therefore significant forces are required for a cell to pass through spaces that are smaller than the resting size of the nucleus. Experimental measurements of these forces are extremely challenging but experimental images of nuclear deformation are regularly obtained in many labs. Therefore we present a computational method to analyse experimental images of nuclear deformation to deduce the forces required to produce such deformations. A mechanical model of the nucleus is necessary for this analysis and here we present two different models. Our computational tool enables us to obtain detailed information about forces causing deformation from microscopy images and consequently provide evidence for actin generated forces squeezing the sides of the nucleus as it enters a constriction.

A Cellular Potts Model for Analyzing Cell Migration across Constraining Pillar Arrays

Cell migration in highly constrained environments is fundamental in a wide variety of physiological and pathological phenomena. In particular, it has been experimentally shown that the migratory capacity of most cell lines depends on their ability to transmigrate through narrow constrictions, which in turn relies on their deformation capacity. In this respect, the nucleus, which occupies a large fraction of the cell volume and is substantially stiffer than the surrounding cytoplasm, imposes a major obstacle. This aspect has also been investigated with the use of microfluidic devices formed by dozens of arrays of aligned polymeric pillars that limit the available space for cell movement. Such experimental systems, in particular, in the designs developed by the groups of Denais and of Davidson, were here reproduced with a tailored version of the Cellular Potts model, a grid-based stochastic approach where cell dynamics are established by a Metropolis algorithm for energy minimization. The proposed model allowed quantitatively analyzing selected cell migratory determinants (e.g., the cell and nuclear speed and deformation, and forces acting at the nuclear membrane) in the case of different experimental setups. Most of the numerical results show a remarkable agreement with the corresponding empirical data.

Evaluation of Cell's Passability in the ECM Network

The microstructure of the extracellular matrix (ECM) plays a key role in affecting cell migration, especially nonproteolytic migration. It is difficult, however, to measure some properties of the ECM, such as stiffness and the passability for cell migration. On the basis of a network model of collagen fiber in the ECM, which has been well applied to simulate mechanical behaviors such as the stress-strain relationship, damage, and failure, we proposed a series of methods to study the microstructural properties containing pore size and pore stiffness and to search for the possible migration paths for cells. Finally, with a given criterion, we quantitatively evaluated the passability of the ECM network for cell migration. The fiber network model with a microstructure and the analysis method presented in this study further our understanding of and ability to evaluate the properties of an ECM network.

Kinetic models with non-local sensing determining cell polarization and speed according to independent cues

Cells move by run and tumble, a kind of dynamics in which the cell alternates runs over straight lines and re-orientations. This erratic motion may be influenced by external factors, like chemicals, nutrients, the extra-cellular matrix, in the sense that the cell measures the external field and elaborates the signal eventually adapting its dynamics. We propose a kinetic transport equation implementing a velocity-jump process in which the transition probability takes into account a double bias, which acts, respectively, on the choice of the direction of motion and of the speed. The double bias depends on two different non-local sensing cues coming from the external environment. We analyze how the size of the cell and the way of sensing the environment with respect to the variation of the external fields affect the cell population dynamics by recovering an appropriate macroscopic limit and directly integrating the kinetic transport equation. A comparison between the solutions of the transport equation and of the proper macroscopic limit is also performed.

Modelling actin polymerization: the effect on confined cell migration

The aim of this work is to model cell motility under conditions of mechanical confinement. This cell migration mode may occur in extravasation of tumour and neutrophil-like cells. Cell migration is the result of the complex action of different forces exerted by the interplay between myosin contractility forces and actin processes. Here, we propose and implement a finite element model of the confined migration of a single cell. In this model, we consider the effects of actin and myosin in cell motility. Both filament and globular actin are modelled. We model the cell considering cytoplasm and nucleus with different mechanical properties. The migration speed in the simulation is around 0.1 μm/min, which is in agreement with existing literature. From our simulation, we observe that the nucleus size has an important role in cell migration inside the channel. In the simulation the cell moves further when the nucleus is smaller. However, this speed is less sensitive to nucleus stiffness. The results show that the cell displacement is lower when the nucleus is stiffer. The degree of adhesion between the channel walls and the cell is also very important in confined migration. We observe an increment of cell velocity when the friction coefficient is higher.

Multicompartment cell-based modeling of confined migration: regulation by cell intrinsic and extrinsic factors

Though cell and nuclear deformability are expected to influence efficiency of confined migration, their individual and collective influence on migration efficiency remains incompletely understood. In addition to cell intrinsic properties, the relevance of cell extrinsic factors on confined migration, if any, has not been adequately explored. Here we address these questions using a statistical mechanics-based stochastic modeling approach where cell/nuclear dimensions and their deformability are explicitly taken into consideration. In addition to demonstrating the importance of cell softness in sustaining confined migration, our results suggest that dynamic tuning of cell and nuclear properties at different stages of migration is essential for maximizing migration efficiency. Our simulations also implicate confinement shape and confinement history as two important cell extrinsic regulators of cell invasiveness. Together, our findings illustrate the strength of a multicompartment model in dissecting the contributions of multiple factors that collectively influence confined cell migration.

Role of the Nucleus as a Sensor of Cell Environment Topography

The proper integration of biophysical cues from the cell vicinity is crucial for cells to maintain homeostasis, cooperate with other cells within the tissues, and properly fulfill their biological function. It is therefore crucial to fully understand how cells integrate these extracellular signals for tissue engineering and regenerative medicine. Topography has emerged as a prominent component of the cellular microenvironment that has pleiotropic effects on cell behavior. This progress report focuses on the recent advances in the understanding of the topography sensing mechanism with a special emphasis on the role of the nucleus. Here, recent techniques developed for monitoring the nuclear mechanics are reviewed and the impact of various topographies and their consequences on nuclear organization, gene regulation, and stem cell fate is summarized. The role of the cell nucleus as a sensor of cell-scale topography is further discussed.

TRPM8 inhibits endothelial cell migration via a non-channel function by trapping the small GTPase Rap1

Endothelial cell adhesion and migration are critical steps of the angiogenic process, whose dysfunction is associated with tumor growth and metastasis. The TRPM8 channel has recently been proposed to play a protective role in prostate cancer by impairing cell motility. However, the mechanisms by which it could influence vascular behavior are unknown. Here, we reveal a novel non-channel function for TRPM8 that unexpectedly acts as a Rap1 GTPase inhibitor, thereby inhibiting endothelial cell motility, independently of pore function. TRPM8 retains Rap1 intracellularly through direct protein-protein interaction, thus preventing its cytoplasm-plasma membrane trafficking. In turn, this mechanism impairs the activation of a major inside-out signaling pathway that triggers the conformational activation of integrin and, consequently, cell adhesion, migration, in vitro endothelial tube formation, and spheroid sprouting. Our results bring to light a novel, pore-independent molecular mechanism by which endogenous TRPM8 expression inhibits Rap1 GTPase and thus plays a critical role in the behavior of vascular endothelial cells by inhibiting migration.

Physical role for the nucleus in cell migration

Cell migration is important for the function of many eukaryotic cells. Recently the nucleus has been shown to play an important role in cell motility. After giving an overview of cell motility mechanisms we review what is currently known about the mechanical properties of the nucleus and the connections between it and the cytoskeleton. We also discuss connections to the extracellular matrix and mechanotransduction. We identify key physical roles of the nucleus in cell migration.

A node-based version of the cellular Potts model

The cellular Potts model (CPM) is a lattice-based Monte Carlo method that uses an energetic formalism to describe the phenomenological mechanisms underlying the biophysical problem of interest. We here propose a CPM-derived framework that relies on a node-based representation of cell-scale elements. This feature has relevant consequences on the overall simulation environment. First, our model can be implemented on any given domain, provided a proper discretization (which can be regular or irregular, fixed or time evolving). Then, it allowed an explicit representation of cell membranes, whose displacements realistically result in cell movement. Finally, our node-based approach can be easily interfaced with continuous mechanics or fluid dynamics models. The proposed computational environment is here applied to some simple biological phenomena, such as cell sorting and chemotactic migration, also in order to achieve an analysis of the performance of the underlying algorithm. This work is finally equipped with a critical comparison between the advantages and disadvantages of our model with respect to the traditional CPM and to some similar vertex-based approaches.

Plasticity of Cell Migration In Vivo and In Silico

Cell migration results from stepwise mechanical and chemical interactions between cells and their extracellular environment. Mechanistic principles that determine single-cell and collective migration modes and their interconversions depend upon the polarization, adhesion, deformability, contractility, and proteolytic ability of cells. Cellular determinants of cell migration respond to extracellular cues, including tissue composition, topography, alignment, and tissue-associated growth factors and cytokines. Both cellular determinants and tissue determinants are interdependent; undergo reciprocal adjustment; and jointly impact cell decision making, navigation, and migration outcome in complex environments. We here review the variability, decision making, and adaptation of cell migration approached by live-cell, in vivo, and in silico strategies, with a focus on cell movements in morphogenesis, repair, immune surveillance, and cancer metastasis.

Migration in Confined 3D Environments Is Determined by a Combination of Adhesiveness, Nuclear Volume, Contractility, and Cell Stiffness

In cancer metastasis and other physiological processes, cells migrate through the three-dimensional (3D) extracellular matrix of connective tissue and must overcome the steric hindrance posed by pores that are smaller than the cells. It is currently assumed that low cell stiffness promotes cell migration through confined spaces, but other factors such as adhesion and traction forces may be equally important. To study 3D migration under confinement in a stiff (1.77 MPa) environment, we use soft lithography to fabricate polydimethylsiloxane (PDMS) devices consisting of linear channel segments with 20 μm length, 3.7 μm height, and a decreasing width from 11.2 to 1.7 μm. To study 3D migration in a soft (550 Pa) environment, we use self-assembled collagen networks with an average pore size of 3 μm. We then measure the ability of four different cancer cell lines to migrate through these 3D matrices, and correlate the results with cell physical properties including contractility, adhesiveness, cell stiffness, and nuclear volume. Furthermore, we alter cell adhesion by coating the channel walls with different amounts of adhesion proteins, and we increase cell stiffness by overexpression of the nuclear envelope protein lamin A. Although all cell lines are able to migrate through the smallest 1.7 μm channels, we find significant differences in the migration velocity. Cell migration is impeded in cell lines with larger nuclei, lower adhesiveness, and to a lesser degree also in cells with lower contractility and higher stiffness. Our data show that the ability to overcome the steric hindrance of the matrix cannot be attributed to a single cell property but instead arises from a combination of adhesiveness, nuclear volume, contractility, and cell stiffness.

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

Cell migration is a fundamental biological phenomenon during which cells sense their surroundings and respond to different types of signals. In presence of durotaxis, cells preferentially crawl from soft to stiff substrates by reorganizing their cytoskeleton from an isotropic to an anisotropic distribution of actin filaments. In the present paper, we propose a Cellular Potts Model to simulate single cell migration over flat substrates with variable stiffness. We have tested five configurations: (i) a substrate including a soft and a stiff region, (ii) a soft substrate including two parallel stiff stripes, (iii) a substrate made of successive stripes with increasing stiffness to create a gradient and (iv) a stiff substrate with four embedded soft squares. For each simulation, we have evaluated the morphology of the cell, the distance covered, the spreading area and the migration speed. We have then compared the numerical results to specific experimental observations showing a consistent agreement.

Cell Invasion Dynamics into a Three Dimensional Extracellular Matrix Fibre Network

The dynamics of filopodia interacting with the surrounding extracellular matrix (ECM) play a key role in various cell-ECM interactions, but their mechanisms of interaction with the ECM in 3D environment remain poorly understood. Based on first principles, here we construct an individual-based, force-based computational model integrating four modules of 1) filopodia penetration dynamics; 2) intracellular mechanics of cellular and nuclear membranes, contractile actin stress fibers, and focal adhesion dynamics; 3) structural mechanics of ECM fiber networks; and 4) reaction-diffusion mass transfers of seven biochemical concentrations in related with chemotaxis, proteolysis, haptotaxis, and degradation in ECM to predict dynamic behaviors of filopodia that penetrate into a 3D ECM fiber network. The tip of each filopodium crawls along ECM fibers, tugs the surrounding fibers, and contracts or retracts depending on the strength of the binding and the ECM stiffness and pore size. This filopodium-ECM interaction is modeled as a stochastic process based on binding kinetics between integrins along the filopodial shaft and the ligands on the surrounding ECM fibers. This filopodia stochastic model is integrated into migratory dynamics of a whole cell in order to predict the cell invasion into 3D ECM in response to chemotaxis, haptotaxis, and durotaxis cues. Predicted average filopodia speed and that of the cell membrane advance agreed with experiments of 3D HUVEC migration at r(2) > 0.95 for diverse ECMs with different pore sizes and stiffness.

A computational mechanics approach to assess the link between cell morphology and forces during confined migration

Confined migration plays a fundamental role during several biological phenomena such as embryogenesis, immunity and tumorogenesis. Here, we propose a two-dimensional mechanical model to simulate the migration of a HeLa cell through a micro-channel. As in our previous works, the cell is modelled as a continuum and a standard Maxwell model is used to describe the mechanical behaviour of both the cytoplasm (including active strains) and the nucleus. The cell cyclically protrudes and contracts and develops viscous forces to adhere to the substrate. The micro-channel is represented by two rigid walls, and it exerts an additional viscous force on the cell boundaries. We test four channels whose dimensions in terms of width are i) larger than the cell diameter, ii) sub-cellular, ii) sub-nuclear and iv) much smaller than the nucleus diameter. The main objective of the work is to assess the necessary conditions for the cell to enter into the channel and migrate through it. Therefore, we evaluate both the evolution of the cell morphology and the cell-channel and cell-substrate surface forces, and we show that there exists a link between the two, which is the essential parameter determining whether the cell is permeative, invasive or penetrating.

The effects of cell compressibility, motility and contact inhibition on the growth of tumor cell clusters using the Cellular Potts Model

There are numerous biological examples where genes associated with migratory ability of cells also confer the cells with an increased fitness even though these genes may not have any known effect on the cell mitosis rates. Here, we provide insight into these observations by analyzing the effects of cell migration, compression, and contact inhibition on the growth of tumor cell clusters using the Cellular Potts Model (CPM) in a monolayer geometry. This is a follow-up of a previous study (Thalhauser et al. 2010) in which a Moran-type model was used to study the interaction of cell proliferation, migratory potential and death on the emergence of invasive phenotypes. Here, we extend the study to include the effects of cell size and shape. In particular, we investigate the interplay between cell motility and compressibility within the CPM and find that the CPM predicts that increased cell motility leads to smaller cells. This is an artifact in the CPM. An analysis of the CPM reveals an explicit inverse-relationship between the cell stiffness and motility parameters. We use this relationship to compensate for motility-induced changes in cell size in the CPM so that in the corrected CPM, cell size is independent of the cell motility. We find that subject to comparable levels of compression, clusters of motile cells grow faster than clusters of less motile cells, in qualitative agreement with biological observations and our previous study. Increasing compression tends to reduce growth rates. Contact inhibition penalizes clumped cells by halting their growth and gives motile cells an even greater advantage. Finally, our model predicts cell size distributions that are consistent with those observed in clusters of neuroblastoma cells cultured in low and high density conditions.

Influence of nucleus deformability on cell entry into cylindrical structures

The mechanical properties of cell nuclei have been demonstrated to play a fundamental role in cell movement across extracellular networks and micro-channels. In this work, we focus on a mathematical description of a cell entering a cylindrical channel composed of extracellular matrix. An energetic approach is derived in order to obtain a necessary condition for which cells enter cylindrical structures. The nucleus of the cell is treated either (i) as an elastic membrane surrounding a liquid droplet or (ii) as an incompressible elastic material with Neo-Hookean constitutive equation. The results obtained highlight the importance of the interplay between mechanical deformability of the nucleus and the capability of the cell to establish adhesive bonds and generate active forces in the cytoskeleton due to myosin action.

Modeling the mechanics of cancer: effect of changes in cellular and extra-cellular mechanical properties

Malignant transformation, though primarily driven by genetic mutations in cells, is also accompanied by specific changes in cellular and extra-cellular mechanical properties such as stiffness and adhesivity. As the transformed cells grow into tumors, they interact with their surroundings via physical contacts and the application of forces. These forces can lead to changes in the mechanical regulation of cell fate based on the mechanical properties of the cells and their surrounding environment. A comprehensive understanding of cancer progression requires the study of how specific changes in mechanical properties influences collective cell behavior during tumor growth and metastasis. Here we review some key results from computational models describing the effect of changes in cellular and extra-cellular mechanical properties and identify mechanistic pathways for cancer progression that can be targeted for the prediction, treatment, and prevention of cancer.