How Nucleus Mechanics and ECM Microstructure Influence the Invasion of Single Cells and Multicellular Aggregates

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.

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.

Computational estimates of mechanical constraints on cell migration through the extracellular matrix

Cell migration through a three-dimensional (3D) extracellular matrix (ECM) underlies important physiological phenomena and is based on a variety of mechanical strategies depending on the cell type and the properties of the ECM. By using computer simulations of the cell’s mid-plane, we investigate two such migration mechanisms—‘push-pull’ (forming a finger-like protrusion, adhering to an ECM node, and pulling the cell body forward) and ‘rear-squeezing’ (pushing the cell body through the ECM by contracting the cell cortex and ECM at the cell rear). We present a computational model that accounts for both elastic deformation and forces of the ECM, an active cell cortex and nucleus, and for hydrodynamic forces and flow of the extracellular fluid, cytoplasm, and nucleoplasm. We find that relations between three mechanical parameters—the cortex’s contractile force, nuclear elasticity, and ECM rigidity—determine the effectiveness of cell migration through the dense ECM. The cell can migrate persistently even if its cortical contraction cannot deform a near-rigid ECM, but then the contraction of the cortex has to be able to sufficiently deform the nucleus. The cell can also migrate even if it fails to deform a stiff nucleus, but then it has to be able to sufficiently deform the ECM. Simulation results show that nuclear stiffness limits the cell migration more than the ECM rigidity. Simulations show the rear-squeezing mechanism of motility results in more robust migration with larger cell displacements than those with the push-pull mechanism over a range of parameter values. Additionally, results show that the rear-squeezing mechanism is aided by hydrodynamics through a pressure gradient.

Computational simulations of two different mechanisms of 3D cell migration in an extracellular matrix are presented. One mechanism represents a mesenchymal mode, characterized by finger-like actin protrusions, while the second mode is more amoeboid in that rear contraction of the cortex propels the cell forward. In both mechanisms, the cell generates a thin actin protrusion on the cortex that attaches to an ECM node. The cell is then either pulled (mesenchymal) or pushed (amoeboid) forward. Results show both mechanisms result in successful migration over a range of simulated parameter values as long as the contractile tension of the cortex exceeds either the nuclear stiffness or ECM stiffness, but not necessarily both. However, the distance traveled by the amoeboid migration mode is more robust to changes in parameter values, and is larger than in simulations of the mesenchymal mode. Additionally, cells experience a favorable fluid pressure gradient when migrating in the amoeboid mode, and an adverse fluid pressure gradient in the mesenchymal mode.

A three dimensional model of multicellular aggregate compression

Multicellular aggregates are an excellent model system to explore the role of tissue biomechanics, which has been demonstrated to play a crucial role in many physiological and pathological processes. In this paper, we propose a three-dimensional mechanical model and apply it to the uniaxial compression of a multicellular aggregate in a realistic biological setting. In particular, we consider an aggregate of initially spherical shape and describe both its elastic deformations and the reorganisation of the cells forming the spheroid. The latter phenomenon, understood as remodelling, is accounted for by assuming that the aggregate undergoes plastic-like distortions. The study of the compression of the spheroid, achieved by means of two parallel, compressive plates, needs the formulation of a contact problem between the living spheroid itself and the plates, and is solved with the aid of the augmented Lagrangian method. The results of the performed numerical simulations are in qualitative agreement with the biological observations reported in the literature and can also be used to estimate quantitatively some fundamental aggregate mechanical parameters.

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.

A phenomenological model for cell and nucleus deformation during cancer metastasis

Cell migration plays an essential role in cancer metastasis. In cancer invasion through confined spaces, cells must undergo extensive deformation, which is a capability related to their metastatic potentials. Here, we simulate the deformation of the cell and nucleus during invasion through a dense, physiological microenvironment by developing a phenomenological computational model. In our work, cells are attracted by a generic emitting source (e.g., a chemokine or stiffness signal), which is treated by using Green’s Fundamental solutions. We use an IMEX integration method where the linear parts and the nonlinear parts are treated by using an Euler backward scheme and an Euler forward method, respectively. We develop the numerical model for an obstacle-induced deformation in 2D or/and 3D. Considering the uncertainty in cell mobility, stochastic processes are incorporated and uncertainties in the input variables are evaluated using Monte Carlo simulations. This quantitative study aims at estimating the likelihood for invasion and the length of the time interval in which the cell invades the tissue through an obstacle. Subsequently, the two-dimensional cell deformation model is applied to simplified cancer metastasis processes to serve as a model for in vivo or in vitro biomedical experiments.

On the morphological stability of multicellular tumour spheroids growing in porous media

Multicellular tumour spheroids (MCTSs) are extensively used as in vitro system models for investigating the avascular growth phase of solid tumours. In this work, we propose a continuous growth model of heterogeneous MCTSs within a porous material, taking into account a diffusing nutrient from the surrounding material directing both the proliferation rate and the mobility of tumour cells. At the time scale of interest, the MCTS behaves as an incompressible viscous fluid expanding inside a porous medium. The cell motion and proliferation rate are modelled using a non-convective chemotactic mass flux, driving the cell expansion in the direction of the external nutrients' source. At the early stages, the growth dynamics is derived by solving the quasi-stationary problem, obtaining an initial exponential growth followed by an almost linear regime, in accordance with experimental observations. We also perform a linear-stability analysis of the quasi-static solution in order to investigate the morphological stability of the radially symmetric growth pattern. We show that mechano-biological cues, as well as geometric effects related to the size of the MCTS subdomains with respect to the diffusion length of the nutrient, can drive a morphological transition to fingered structures, thus triggering the formation of complex shapes that might promote tumour invasiveness. The results also point out the formation of a retrograde flow in the MCTS close to the regions where protrusions form, that could describe the initial dynamics of metastasis detachment from the in vivo tumour mass. In conclusion, the results of the proposed model demonstrate that the integration of mathematical tools in biological research could be crucial for better understanding the tumour's ability to invade its host environment.