Invasion from a cell aggregate--the roles of active cell motion and mechanical equilibrium

Active phase separation: new phenomenology from non-equilibrium physics

In active systems, whose constituents have non-equilibrium dynamics at local level, fluid-fluid phase separation is widely observed. Examples include the formation of membraneless organelles within cells; the clustering of self-propelled colloidal particles in the absence of attractive forces, and some types of ecological segregation. A schematic understanding of such active phase separation was initially borrowed from what is known for the equilibrium case, in which detailed balance holds at microscopic level. However it has recently become clear that in active systems the absence of detailed balance, although it leave phase separation qualitatively unchanged in some regimes (for example domain growth driven by interfacial tension via Ostwald ripening), can in other regimes radically alter its phenomenology at mechanistic level. For example, microphase separation can be caused by reverse Ostwald ripening, a process that is hard to imagine from an equilibrium perspective. This and other new phenomena arise because, instead of having a single, positive interfacial tension like their equilibrium counterparts, the fluid-fluid interfaces created by active phase separation can have several distinct interfacial tensions governing different properties, some of which can be negative. These phenomena can be broadly understood by studying continuum field theories for a single conserved scalar order parameter (the fluid density), supplemented with a velocity field in cases where momentum conservation is also present. More complex regimes arise in systems described by multiple scalar order parameters (especially with nonreciprocal interactions between these); or when an order parameter undergoes both conserved and non-conserved dynamics (such that the combination breaks detailed balance); or in systems that support orientational long-range order in one or more of the coexisting phases. In this Review, we survey recent progress in understanding the specific role of activity in phase separation, drawing attention to many open questions. We focus primarily on continuum theories, especially those with a single scalar order parameter, reviewing both analytical and numerical work. We compare their predictions with particle-based models, which have mostly been studied numerically although a few have been explicitly coarse-grained to continuum level. We also compare, where possible, with experimental results. In the latter case, qualitative comparisons are broadly encouraging whereas quantitative ones are hindered by the dynamical complexity of most experimental systems relative that of simplified (particle-level or continuum) models of active matter.

A Simple Framework for Agent-Based Modeling with Extracellular Matrix

Extracellular matrix (ECM) is a key component of the cellular microenvironment and critical in multiple disease and developmental processes. Representing ECM and cell-ECM interactions is a challenging multiscale problem as they span molecular-level details to tissue-level dynamics. While several computational frameworks exist for ECM modeling, they often focus on very detailed modeling of individual ECM fibers or represent only a single aspect of the ECM. Using the PhysiCell agent-based modeling platform, we developed a framework of intermediate detail with the ability to capture bidirectional cell-ECM interactions. We represent a small region of ECM, an ECM element, with three variables describing its local microstructure: anisotropy, density, and overall fiber orientation. To spatially model the ECM, we use an array of ECM elements. Cells remodel local ECM microstructure and in turn, local microstructure impacts cellular motility. We demonstrate the utility of this framework and reusability of its core cell-ECM interaction model through examples in cellular invasion, wound healing, basement membrane degradation, and leader-follower collective migration. Despite the relative simplicity of the framework, it is able to capture a broad range of cell-ECM interactions of interest to the modeling community. Furthermore, variables representing the ECM microstructure are accessible through simple programming interfaces. This allows them to impact cell behaviors, such as proliferation and death, without requiring custom code for each interaction, particularly through PhysiCell's modeling grammar, enabling rapid modeling of a diverse range of cell-matrix biology. We make this framework available as a free and open source software package at https://github.com/PhysiCell-Models/collective-invasion .

Minimal cellular automaton model with heterogeneous cell sizes predicts epithelial colony growth

Regulation of cell proliferation is a crucial aspect of tissue development and homeostasis and plays a major role in morphogenesis, wound healing, and tumor invasion. A phenomenon of such regulation is contact inhibition, which describes the dramatic slowing of proliferation, cell migration and individual cell growth when multiple cells are in contact with each other. While many physiological, molecular and genetic factors are known, the mechanism of contact inhibition is still not fully understood. In particular, the relevance of cellular signaling due to interfacial contact for contact inhibition is still debated. Cellular automata (CA) have been employed in the past as numerically efficient mathematical models to study the dynamics of cell ensembles, but they are not suitable to explore the origins of contact inhibition as such agent-based models assume fixed cell sizes. We develop a minimal, data-driven model to simulate the dynamics of planar cell cultures by extending a probabilistic CA to incorporate size changes of individual cells during growth and cell division. We successfully apply this model to previous in-vitro experiments on contact inhibition in epithelial tissue: After a systematic calibration of the model parameters to measurements of single-cell dynamics, our CA model quantitatively reproduces independent measurements of emergent, culture-wide features, like colony size, cell density and collective cell migration. In particular, the dynamics of the CA model also exhibit the transition from a low-density confluent regime to a stationary postconfluent regime with a rapid decrease in cell size and motion. This implies that the volume exclusion principle, a mechanical constraint which is the only inter-cellular interaction incorporated in the model, paired with a size-dependent proliferation rate is sufficient to generate the observed contact inhibition. We discuss how our approach enables the introduction of effective bio-mechanical interactions in a CA framework for future studies.

Modeling the extracellular matrix in cell migration and morphogenesis: a guide for the curious biologist

The extracellular matrix (ECM) is a highly complex structure through which biochemical and mechanical signals are transmitted. In processes of cell migration, the ECM also acts as a scaffold, providing structural support to cells as well as points of potential attachment. Although the ECM is a well-studied structure, its role in many biological processes remains difficult to investigate comprehensively due to its complexity and structural variation within an organism. In tandem with experiments, mathematical models are helpful in refining and testing hypotheses, generating predictions, and exploring conditions outside the scope of experiments. Such models can be combined and calibrated with in vivo and in vitro data to identify critical cell-ECM interactions that drive developmental and homeostatic processes, or the progression of diseases. In this review, we focus on mathematical and computational models of the ECM in processes such as cell migration including cancer metastasis, and in tissue structure and morphogenesis. By highlighting the predictive power of these models, we aim to help bridge the gap between experimental and computational approaches to studying the ECM and to provide guidance on selecting an appropriate model framework to complement corresponding experimental studies.

A mesoscale mechanical model of cellular interactions

Computational models of cell mechanics allow the precise interrogation of cell shape change. These morphological changes are required for cells to survive in diverse tissue environments. Here, we present a mesoscale mechanical model of cell-substrate interactions using the level set method based on experimentally measured parameters. By implementing a viscoelastic mechanical equivalent circuit, we accurately model whole-cell deformations that are important for a variety of cellular processes. To effectively model shape changes as a cell interacts with a substrate, we have included receptor-mediated adhesion, which is governed by catch-slip bond behavior. The effect of adhesion was explored by subjecting cells to a variety of different substrates including flat, curved, and deformable surfaces. Finally, we increased the accuracy of our simulations by including a deformable nucleus in our cells. This model sets the foundation for further exploration into computational analyses of multicellular interactions.

Cell and extracellular matrix interaction models in benign mesothelial and malignant pleural mesothelioma cells in 2D and 3D in-vitro

Malignant pleural mesothelioma (MPM) is an aggressive tumour that grows in the pleural cavity. MPM spheroids released in the pleural fluid can form new tumour foci. Cell-cell, cell-extracellular matrix (ECM) interactions in 2D and 3D impact malignant cell behaviour during cell adhesion, migration, proliferation and epithelial-mesenchymal transition (EMT). In this study, epithelioid, biphasic and sarcomatoid MPM cell types as well as benign mesothelial cells were tested with regards to the above phenotypes. Fibronectin (FN) and homologous cell-derived extracellular matrix (hcd-ECM) treated substratum differentially affected the above phenotypes. 3D MPM spheroid invasion was higher in FN-collagen matrices in the epithelioid and biphasic cells, while 3D cell cultures of epithelioid and sarcomatoid MPM cells in FN-collagen showed a higher contractility compared to hcd-ECM-collagen. Cell aggregates demonstrated invasive behaviour in hcd-ECM matrices alone. Our results suggest that ECM and the dimensionality affect malignant cell behaviour during cell culture studies.

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.

Cancer-leukocyte hybrid membrane-cloaked magnetic beads for the ultrasensitive isolation, purification, and non-destructive release of circulating tumor cells

Most of the current circulating tumor cell (CTC) isolation techniques are based on immunomagnetic beads with antibodies or aptamers that specifically target epithelial cell adhesion molecules (EpCAMs). However, these techniques are unsuitable for the isolation and purification of circulating tumor cells because they fail to recognize EpCAM-negative CTCs and thus lead to the non-specific adsorption of background leucocytes and EpCAM-positive circulating epithelial cells. Moreover, releasing the CTCs from the capture platform without disruption is a big challenge. To address these issues, herein, we developed biomimetic magnetic beads (MBs) by cloaking a cancer cell-leukocyte hybrid membrane on the MBs. These biomimetic MBs inherited homologous CTC binding capability from the cancer cell membrane and less affinity for the background cells from the leukocyte membrane, exhibitng a higher CTC capture efficiency and separation purity than EpCAM-based MBs. Importantly, the captured CTCs could be rapidly released by a facile method i.e. co-incubation with a trypsin-EDTA solution. We demonstrated the excellent performance of these MBs for the highly pure separation and non-destructive release of CTCs in metastatic mammary carcinoma models. Our results indicate that the proposed homologous cancer-leukocyte membrane coating strategy may provide a promising method for the ultrahigh-specific and sensitive detection of CTCs.

Collective invasion of glioma cells through OCT1 signalling and interaction with reactive astrocytes after surgery

Glioblastoma multiforme (GBM) is the most aggressive form of brain cancer with a short median survival time. GBM is characterized by the hallmarks of aggressive proliferation and cellular infiltration of normal brain tissue. miR-451 and its downstream molecules are known to play a pivotal role in regulation of the balance of proliferation and aggressive invasion in response to metabolic stress in the tumour microenvironment (TME). Surgery-induced transition in reactive astrocyte populations can play a significant role in tumour dynamics. In this work, we develop a multi-scale mathematical model of miR-451-LKB1-AMPK-OCT1-mTOR pathway signalling and individual cell dynamics of the tumour and reactive astrocytes after surgery. We show how the effects of fluctuating glucose on tumour cells need to be reprogrammed by taking into account the recent history of glucose variations and an AMPK/miR-451 reciprocal feedback loop. The model shows how variations in glucose availability significantly affect the activity of signalling molecules and, in turn, lead to critical cell migration. The model also predicts that microsurgery of a primary tumour induces phenotypical changes in reactive astrocytes and stem cell-like astrocytes promoting tumour cell proliferation and migration by Cxcl5. Finally, we investigated a new anti-tumour strategy by Cxcl5-targeting drugs. This article is part of the theme issue 'Multi-scale analysis and modelling of collective migration in biological systems'.

Multiscale modelling of motility wave propagation in cell migration

The collective motion of cell monolayers within a tissue is a fundamental biological process that occurs during tissue formation, wound healing, cancerous invasion, and viral infection. Experiments have shown that at the onset of migration, the motility is self-generated as a polarisation wave starting from the leading edge of the monolayer and progressively propagates into the bulk. However, it is unclear how the propagation of this motility wave is influenced by cellular properties. Here, we investigate this question using a computational model based on the Potts model coupled to the dynamics of intracellular polarisation. The model captures the propagation of the polarisation wave and suggests that the cells cortex can regulate the migration modes: strongly contractile cells may depolarise the monolayer, whereas less contractile cells can form swirling movement. Cortical contractility is further found to limit the cells motility, which (i) decelerates the wave speed and the leading edge progression, and (ii) destabilises the leading edge. Together, our model describes how different mechanical properties of cells can contribute to the regulation of collective cell migration.

MMP Secretion Rate and Inter-invadopodia Spacing Collectively Govern Cancer Invasiveness

Invadopodia are micron-sized invasive structures that mediate extracellular matrix (ECM) degradation through a combination of membrane-bound and soluble matrix metalloproteinases (MMPs). However, how such localized degradation is converted into pores big enough for cancer cells to invade, and the relative contributions of membrane-bound versus soluble MMPs to this process remain unclear. In this article, we address these questions by combining experiments and simulations. We show that in MDA-MB-231 cells, an increase in ECM density enhances invadopodia-mediated ECM degradation and decreases inter-invadopodia spacing. ECM degradation is mostly mediated by soluble MMPs, which are activated by membrane-bound MT1-MMP. We present a computational model of invadopodia-mediated ECM degradation, which recapitulates the above observations and identifies MMP secretion rate as an important regulator of invadopodia stability. Simulations with multiple invadopodia suggest that inter-invadopodia spacing and MMP secretion rate collectively dictate the size of the degraded zones. Taken together, our results suggest that for creating pores conducive for cancer invasion, cells must tune inter-invadopodia spacing and MMP secretion rate in an ECM density-dependent manner, thereby striking a balance between invadopodia penetration and ECM degradation.

Soluble VEGFR1 signaling guides vascular patterns into dense branching morphologies

Vascular patterning is a key process during development and disease. The diffusive decoy receptor sVEGFR1 (sFlt1) is a known regulator of endothelial cell behavior, yet the mechanism by which it controls vascular structure is little understood. We propose computational models to shed light on how vascular patterning is guided by self-organized gradients of the VEGF/sVEGFR1 factors. We demonstrate that a diffusive inhibitor can generate structures with a dense branching morphology in models where the activator elicits directed growth. Inadequate presence of the inhibitor leads to compact growth, while excessive production of the inhibitor blocks expansion and stabilizes existing structures. Model predictions were compared with time-resolved experimental data obtained from endothelial sprout kinetics in fibrin gels. In the presence of inhibitory antibodies against VEGFR1 vascular sprout density increases while the speed of sprout expansion remains unchanged. Thus, the rate of secretion and stability of extracellular sVEGFR1 can modulate vascular sprout density.

Role of extracellular matrix and microenvironment in regulation of tumor growth and LAR-mediated invasion in glioblastoma

The cellular dispersion and therapeutic control of glioblastoma, the most aggressive type of primary brain cancer, depends critically on the migration patterns after surgery and intracellular responses of the individual cancer cells in response to external biochemical cues in the microenvironment. Recent studies have shown that miR-451 regulates downstream molecules including AMPK/CAB39/MARK and mTOR to determine the balance between rapid proliferation and invasion in response to metabolic stress in the harsh tumor microenvironment. Surgical removal of the main tumor is inevitably followed by recurrence of the tumor due to inaccessibility of dispersed tumor cells in normal brain tissue. In order to address this complex process of cell proliferation and invasion and its response to conventional treatment, we propose a mathematical model that analyzes the intracellular dynamics of the miR-451-AMPK- mTOR-cell cycle signaling pathway within a cell. The model identifies a key mechanism underlying the molecular switches between proliferative phase and migratory phase in response to metabolic stress in response to fluctuating glucose levels. We show how up- or down-regulation of components in these pathways affects the key cellular decision to infiltrate or proliferate in a complex microenvironment in the absence and presence of time delays and stochastic noise. Glycosylated chondroitin sulfate proteoglycans (CSPGs), a major component of the extracellular matrix (ECM) in the brain, contribute to the physical structure of the local brain microenvironment but also induce or inhibit glioma invasion by regulating the dynamics of the CSPG receptor LAR as well as the spatiotemporal activation status of resident astrocytes and tumor-associated microglia. Using a multi-scale mathematical model, we investigate a CSPG-induced switch between invasive and non-invasive tumors through the coordination of ECM-cell adhesion and dynamic changes in stromal cells. We show that the CSPG-rich microenvironment is associated with non-invasive tumor lesions through LAR-CSGAG binding while the absence of glycosylated CSPGs induce the critical glioma invasion. We illustrate how high molecular weight CSPGs can regulate the exodus of local reactive astrocytes from the main tumor lesion, leading to encapsulation of non-invasive tumor and inhibition of tumor invasion. These different CSPG conditions also change the spatial profiles of ramified and activated microglia. The complex distribution of CSPGs in the tumor microenvironment can determine the nonlinear invasion behaviors of glioma cells, which suggests the need for careful therapeutic strategies.

The biology and mathematical modelling of glioma invasion: a review

Adult gliomas are aggressive brain tumours associated with low patient survival rates and limited life expectancy. The most important hallmark of this type of tumour is its invasive behaviour, characterized by a markedly phenotypic plasticity, infiltrative tumour morphologies and the ability of malignant progression from low- to high-grade tumour types. Indeed, the widespread infiltration of healthy brain tissue by glioma cells is largely responsible for poor prognosis and the difficulty of finding curative therapies. Meanwhile, mathematical models have been established to analyse potential mechanisms of glioma invasion. In this review, we start with a brief introduction to current biological knowledge about glioma invasion, and then critically review and highlight future challenges for mathematical models of glioma invasion.

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.

Tumour heterogeneity promotes collective invasion and cancer metastatic dissemination

Heterogeneity within tumour cell populations is commonly observed in most cancers. However, its impact on metastatic dissemination, one of the primary determinants of the disease prognosis, remains poorly understood. Working with a simplified numerical model of tumour spheroids, we investigated the impact of mechanical heterogeneity on the onset of tumour invasion into surrounding tissues. Our work establishes a positive link between tumour heterogeneity and metastatic dissemination, and recapitulates a number of invasion patterns identified in vivo, such as multicellular finger-like protrusions. Two complementary mechanisms are at play in heterogeneous tumours. A small proportion of stronger cells are able to initiate and lead the escape of cells, while collective effects in the bulk of the tumour provide the coordination required to sustain the invasive process through multicellular streaming. This suggests that the multicellular dynamics observed during metastasis is a generic feature of mechanically heterogeneous cell populations and might rely on a limited and generic set of attributes.

The Molecular Basis of Radial Intercalation during Tissue Spreading in Early Development

Radial intercalation is a fundamental process responsible for the thinning of multilayered tissues during large-scale morphogenesis; however, its molecular mechanism has remained elusive. Using amphibian epiboly, the thinning and spreading of the animal hemisphere during gastrulation, here we provide evidence that radial intercalation is driven by chemotaxis of cells toward the external layer of the tissue. This role of chemotaxis in tissue spreading and thinning is unlike its typical role associated with large-distance directional movement of cells. We identify the chemoattractant as the complement component C3a, a factor normally linked with the immune system. The mechanism is explored by computational modeling and tested in vivo, ex vivo, and in vitro. This mechanism is robust against fluctuations of chemoattractant levels and expression patterns and explains expansion during epiboly. This study provides insight into the fundamental process of radial intercalation and could be applied to a wide range of morphogenetic events.

Modeling the two-way feedback between contractility and matrix realignment reveals a nonlinear mode of cancer cell invasion

Cancer cell invasion from primary tumors is mediated by a complex interplay between cellular adhesions, actomyosin-driven contractility, and the physical characteristics of the extracellular matrix (ECM). Here, we incorporate a mechanochemical free-energy-based approach to elucidate how the two-way feedback loop between cell contractility (induced by the activity of chemomechanical interactions such as Ca2+ and Rho signaling pathways) and matrix fiber realignment and strain stiffening enables the cells to polarize and develop contractile forces to break free from the tumor spheroids and invade into the ECM. Interestingly, through this computational model, we are able to identify a critical stiffness that is required by the matrix to break intercellular adhesions and initiate cell invasion. Also, by considering the kinetics of the cell movement, our model predicts a biphasic invasiveness with respect to the stiffness of the matrix. These predictions are validated by analyzing the invasion of melanoma cells in collagen matrices of varying concentration. Our model also predicts a positive correlation between the elongated morphology of the invading cells and the alignment of fibers in the matrix, suggesting that cell polarization is directly proportional to the stiffness and alignment of the matrix. In contrast, cells in nonfibrous matrices are found to be rounded and not polarized, underscoring the key role played by the nonlinear mechanics of fibrous matrices. Importantly, our model shows that mechanical principles mediated by the contractility of the cells and the nonlinearity of the ECM behavior play a crucial role in determining the phenotype of the cell invasion.

Physical models of collective cell motility: from cell to tissue

In this article, we review physics-based models of collective cell motility. We discuss a range of techniques at different scales, ranging from models that represent cells as simple self-propelled particles to phase field models that can represent a cell's shape and dynamics in great detail. We also extensively review the ways in which cells within a tissue choose their direction, the statistics of cell motion, and some simple examples of how cell-cell signaling can interact with collective cell motility. This review also covers in more detail selected recent works on collective cell motion of small numbers of cells on micropatterns, in wound healing, and the chemotaxis of clusters of cells.

3D hybrid modelling of vascular network formation

We develop an off-lattice, agent-based model to describe vasculogenesis, the de novo formation of blood vessels from endothelial progenitor cells during development. The endothelial cells that comprise our vessel network are viewed as linearly elastic spheres that move in response to the forces they experience. We distinguish two types of endothelial cells: vessel elements are contained within the network and tip cells are located at the ends of vessels. Tip cells move in response to mechanical forces caused by interactions with neighbouring vessel elements and the local tissue environment, chemotactic forces and a persistence force which accounts for their tendency to continue moving in the same direction. Vessel elements are subject to similar mechanical forces but are insensitive to chemotaxis. An angular persistence force representing interactions with the local tissue is introduced to stabilise buckling instabilities caused by cell proliferation. Only vessel elements proliferate, at rates which depend on their degree of stretch: elongated elements have increased rates of proliferation, and compressed elements have reduced rates. Following division, the fate of the new cell depends on the local mechanical environment: the probability of forming a new sprout is increased if the parent vessel is highly compressed and the probability of being incorporated into the parent vessel increased if the parent is stretched. Simulation results reveal that our hybrid model can reproduce the key qualitative features of vasculogenesis. Extensive parameter sensitivity analyses show that significant changes in network size and morphology are induced by varying the chemotactic sensitivity of tip cells, and the sensitivities of the proliferation rate and the sprouting probability to mechanical stretch. Varying the chemotactic sensitivity directly influences the directionality of the networks. The degree of branching, and thereby the density of the networks, is influenced by the sprouting probability. Glyphs that simultaneously depict several network properties are introduced to show how these and other network quantities change over time and also as model parameters vary. We also show how equivalent glyphs constructed from in vivo data could be used to discriminate between normal and tumour vasculature and, in the longer term, for model validation. We conclude that our biomechanical hybrid model can generate vascular networks that are qualitatively similar to those generated from in vitro and in vivo experiments.

Collective cancer cell invasion induced by coordinated contractile stresses

The physical underpinnings of fibrosarcoma cell dissemination from a tumor in a surrounding collagen-rich matrix are poorly understood. Here we show that a tumor spheroid embedded in a 3D collagen matrix exerts large contractile forces on the matrix before invasion. Cell invasion is accompanied by complex spatially and temporally dependent patterns of cell migration within and at the surface of the spheroids that are fundamentally different from migratory patterns of individual fibrosarcoma cells homogeneously distributed in the same type of matrix. Cells display a continuous transition from a round morphology at the spheroid core, to highly aligned elongated morphology at the spheroid periphery, which depends on both β₁-integrin-based cell-matrix adhesion and myosin II/ROCK-based cell contractility. This isotropic-to-anisotropic transition corresponds to a shift in migration, from a slow and unpolarized movement at the core, to a fast, polarized and persistent one at the periphery. Our results also show that the ensuing collective invasion of fibrosarcoma cells is induced by anisotropic contractile stresses exerted on the surrounding matrix.

Modeling cell shape and dynamics on micropatterns

Adhesive micropatterns have become a standard tool to study cells under defined conditions. Applications range from controlling the differentiation and fate of single cells to guiding the collective migration of cell sheets. In long-term experiments, single cell normalization is challenged by cell division. For all of these setups, mathematical models predicting cell shape and dynamics can guide pattern design. Here we review recent advances in predicting and explaining cell shape, traction forces and dynamics on micropatterns. Starting with contour models as the simplest approach to explain concave cell shapes, we move on to network and continuum descriptions as examples for static models. To describe dynamic processes, cellular Potts, vertex and phase field models can be used. Different types of model are appropriate to address different biological questions and together, they provide a versatile tool box to predict cell behavior on micropatterns.

Collective Signal Processing in Cluster Chemotaxis: Roles of Adaptation, Amplification, and Co-attraction in Collective Guidance

Single eukaryotic cells commonly sense and follow chemical gradients, performing chemotaxis. Recent experiments and theories, however, show that even when single cells do not chemotax, clusters of cells may, if their interactions are regulated by the chemoattractant. We study this general mechanism of "collective guidance" computationally with models that integrate stochastic dynamics for individual cells with biochemical reactions within the cells, and diffusion of chemical signals between the cells. We show that if clusters of cells use the well-known local excitation, global inhibition (LEGI) mechanism to sense chemoattractant gradients, the speed of the cell cluster becomes non-monotonic in the cluster's size-clusters either larger or smaller than an optimal size will have lower speed. We argue that the cell cluster speed is a crucial readout of how the cluster processes chemotactic signals; both amplification and adaptation will alter the behavior of cluster speed as a function of size. We also show that, contrary to the assumptions of earlier theories, collective guidance does not require persistent cell-cell contacts and strong short range adhesion. If cell-cell adhesion is absent, and the cluster cohesion is instead provided by a co-attraction mechanism, e.g. chemotaxis toward a secreted molecule, collective guidance may still function. However, new behaviors, such as cluster rotation, may also appear in this case. Co-attraction and adaptation allow for collective guidance that is robust to varying chemoattractant concentrations while not requiring strong cell-cell adhesion.

Ring-Shaped Microlanes and Chemical Barriers as a Platform for Probing Single-Cell Migration

Quantification and discrimination of pharmaceutical and disease-related effects on cell migration requires detailed characterization of single-cell motility. In this context, micropatterned substrates that constrain cells within defined geometries facilitate quantitative readout of locomotion. Here, we study quasi-one-dimensional cell migration in ring-shaped microlanes. We observe bimodal behavior in form of alternating states of directional migration (run state) and reorientation (rest state). Both states show exponential lifetime distributions with characteristic persistence times, which, together with the cell velocity in the run state, provide a set of parameters that succinctly describe cell motion. By introducing PEGylated barriers of different widths into the lane, we extend this description by quantifying the effects of abrupt changes in substrate chemistry on migrating cells. The transit probability decreases exponentially as a function of barrier width, thus specifying a characteristic penetration depth of the leading lamellipodia. Applying this fingerprint-like characterization of cell motion, we compare different cell lines, and demonstrate that the cancer drug candidate salinomycin affects transit probability and resting time, but not run time or run velocity. Hence, the presented assay allows to assess multiple migration-related parameters, permits detailed characterization of cell motility, and has potential applications in cell biology and advanced drug screening.

In vivo confinement promotes collective migration of neural crest cells

Collective cell migration is fundamental throughout development and in many diseases. Spatial confinement using micropatterns has been shown to promote collective cell migration in vitro, but its effect in vivo remains unclear. Combining computational and experimental approaches, we show that the in vivo collective migration of neural crest cells (NCCs) depends on such confinement. We demonstrate that confinement may be imposed by the spatiotemporal distribution of a nonpermissive substrate provided by versican, an extracellular matrix molecule previously proposed to have contrasting roles: barrier or promoter of NCC migration. We resolve the controversy by demonstrating that versican works as an inhibitor of NCC migration and also acts as a guiding cue by forming exclusionary boundaries. Our model predicts an optimal number of cells in a given confinement width to allow for directional migration. This optimum coincides with the width of neural crest migratory streams analyzed across different species, proposing an explanation for the highly conserved nature of NCC streams during development.

Dynamics of Cell Ensembles on Adhesive Micropatterns: Bridging the Gap between Single Cell Spreading and Collective Cell Migration

The collective dynamics of multicellular systems arise from the interplay of a few fundamental elements: growth, division and apoptosis of single cells; their mechanical and adhesive interactions with neighboring cells and the extracellular matrix; and the tendency of polarized cells to move. Micropatterned substrates are increasingly used to dissect the relative roles of these fundamental processes and to control the resulting dynamics. Here we show that a unifying computational framework based on the cellular Potts model can describe the experimentally observed cell dynamics over all relevant length scales. For single cells, the model correctly predicts the statistical distribution of the orientation of the cell division axis as well as the final organisation of the two daughters on a large range of micropatterns, including those situations in which a stable configuration is not achieved and rotation ensues. Large ensembles migrating in heterogeneous environments form non-adhesive regions of inward-curved arcs like in epithelial bridge formation. Collective migration leads to swirl formation with variations in cell area as observed experimentally. In each case, we also use our model to predict cell dynamics on patterns that have not been studied before.

Proteolytic and non-proteolytic regulation of collective cell invasion: tuning by ECM density and organization

Cancer cells manoeuvre through extracellular matrices (ECMs) using different invasion modes, including single cell and collective cell invasion. These modes rely on MMP-driven ECM proteolysis to make space for cells to move. How cancer-associated alterations in ECM influence the mode of invasion remains unclear. Further, the sensitivity of the two invasion modes to MMP dynamics remains unexplored. In this paper, we address these open questions using a multiscale hybrid computational model combining ECM density-dependent MMP secretion, MMP diffusion, ECM degradation by MMP and active cell motility. Our results demonstrate that in randomly aligned matrices, collective cell invasion is more efficient than single cell invasion. Although increase in MMP secretion rate enhances invasiveness independent of cell-cell adhesion, sustenance of collective invasion in dense matrices requires high MMP secretion rates. However, matrix alignment can sustain both single cell and collective cell invasion even without ECM proteolysis. Similar to our in-silico observations, increase in ECM density and MMP inhibition reduced migration of MCF-7 cells embedded in sandwich gels. Together, our results indicate that apart from cell intrinsic factors (i.e., high cell-cell adhesion and MMP secretion rates), ECM density and organization represent two important extrinsic parameters that govern collective cell invasion and invasion plasticity.

Cell-based computational modeling of vascular morphogenesis using Tissue Simulation Toolkit

Computational modeling has become a widely used tool for unraveling the mechanisms of higher level cooperative cell behavior during vascular morphogenesis. However, experimenting with published simulation models or adding new assumptions to those models can be daunting for novice and even for experienced computational scientists. Here, we present a step-by-step, practical tutorial for building cell-based simulations of vascular morphogenesis using the Tissue Simulation Toolkit (TST). The TST is a freely available, open-source C++ library for developing simulations with the two-dimensional cellular Potts model, a stochastic, agent-based framework to simulate collective cell behavior. We will show the basic use of the TST to simulate and experiment with published simulations of vascular network formation. Then, we will present step-by-step instructions and explanations for building a recent simulation model of tumor angiogenesis. Demonstrated mechanisms include cell-cell adhesion, chemotaxis, cell elongation, haptotaxis, and haptokinesis.

Mechanical cell-matrix feedback explains pairwise and collective endothelial cell behavior in vitro

In vitro cultures of endothelial cells are a widely used model system of the collective behavior of endothelial cells during vasculogenesis and angiogenesis. When seeded in an extracellular matrix, endothelial cells can form blood vessel-like structures, including vascular networks and sprouts. Endothelial morphogenesis depends on a large number of chemical and mechanical factors, including the compliancy of the extracellular matrix, the available growth factors, the adhesion of cells to the extracellular matrix, cell-cell signaling, etc. Although various computational models have been proposed to explain the role of each of these biochemical and biomechanical effects, the understanding of the mechanisms underlying in vitro angiogenesis is still incomplete. Most explanations focus on predicting the whole vascular network or sprout from the underlying cell behavior, and do not check if the same model also correctly captures the intermediate scale: the pairwise cell-cell interactions or single cell responses to ECM mechanics. Here we show, using a hybrid cellular Potts and finite element computational model, that a single set of biologically plausible rules describing (a) the contractile forces that endothelial cells exert on the ECM, (b) the resulting strains in the extracellular matrix, and (c) the cellular response to the strains, suffices for reproducing the behavior of individual endothelial cells and the interactions of endothelial cell pairs in compliant matrices. With the same set of rules, the model also reproduces network formation from scattered cells, and sprouting from endothelial spheroids. Combining the present mechanical model with aspects of previously proposed mechanical and chemical models may lead to a more complete understanding of in vitro angiogenesis.

Velocity alignment leads to high persistence in confined cells

Many cell types display random motility on two-dimensional substrates but crawl persistently in a single direction when confined in a microchannel or on an adhesive micropattern. Does this imply that the motility mechanism of confined cells is fundamentally different from that of unconfined cells? We argue that both free- and confined-cell migration may be described by a generic model of cells as "velocity-aligning" active Brownian particles previously proposed to solve a completely separate problem in collective cell migration. Our model can be mapped to a diffusive escape over a barrier and analytically solved to determine the cell's orientation distribution and repolarization rate. In quasi-one-dimensional confinement, velocity-aligning cells maintain their direction for times that can be exponentially larger than their persistence time in the absence of confinement. Our results suggest an important connection between single- and collective-cell migration: high persistence in confined cells corresponds with fast alignment of velocity to cell-cell forces.

Synergy of cell-cell repulsion and vacuolation in a computational model of lumen formation

A key step in blood vessel development (angiogenesis) is lumen formation: the hollowing of vessels for blood perfusion. Two alternative lumen formation mechanisms are suggested to function in different types of blood vessels. The vacuolation mechanism is suggested for lumen formation in small vessels by coalescence of intracellular vacuoles, a view that was extended to extracellular lumen formation by exocytosis of vacuoles. The cell–cell repulsion mechanism is suggested to initiate extracellular lumen formation in large vessels by active repulsion of adjacent cells, and active cell shape changes extend the lumen. We used an agent-based computer model, based on the cellular Potts model, to compare and study both mechanisms separately and combined. An extensive sensitivity analysis shows that each of the mechanisms on its own can produce lumens in a narrow region of parameter space. However, combining both mechanisms makes lumen formation much more robust to the values of the parameters, suggesting that the mechanisms may work synergistically and operate in parallel, rather than in different vessel types.

Endothelial cell motility, coordination and pattern formation during vasculogenesis

How vascular networks assemble is a fundamental problem of developmental biology that also has medical importance. To explain the organizational principles behind vascular patterning, we must understand how can tissue level structures be controlled through cell behavior patterns like motility and adhesion that, in turn, are determined by biochemical signal transduction processes? We discuss the various ideas that have been proposed as mechanisms for vascular network assembly: cell motility guided by extracellular matrix alignment (contact guidance), chemotaxis guided by paracrine and autocrine morphogens, and multicellular sprouting guided by cell-cell contacts. All of these processes yield emergent patterns, thus endothelial cells can form an interconnected structure autonomously, without guidance from an external pre-pattern.

Collective cell streams in epithelial monolayers depend on cell adhesion

We report a spontaneously emerging, randomly oriented, collective streaming behavior within a monolayer culture of a human keratinocyte cell line, and explore the effect of modulating cell adhesions by perturbing the function of calcium-dependent cell adhesion molecules. We demonstrate that decreasing cell adhesion induces narrower and more anisotropic cell streams, reminiscent of decreasing the Taylor scale of turbulent liquids. To explain our empirical findings, we propose a cell-based model that represents the dual nature of cell-cell adhesions. Spring-like connections provide mechanical stability, while a cellular Potts model formalism represents surface-tension driven attachment. By changing the relevance and persistence of mechanical links between cells, we are able to explain the experimentally observed changes in emergent flow patterns.

Geometrical constraints and physical crowding direct collective migration of fibroblasts

Migrating cells constantly interact with their immediate microenvironment and neighbors. Although studies on single cell migration offer us insights into the molecular and biochemical signaling pathways, they cannot predict the influence of cell crowding and geometrical cues. Using microfabrication techniques, we examine the influence of cell density and geometrical constraints on migrating fibroblasts. Fibroblasts were allowed to migrate on fibronectin strips of different widths. Under such conditions, cells experience various physical guidance cues including boundary effect, confinement and contact inhibition from neighboring cells. Fibroblasts migrating along the edge of the fibronectin pattern exhibit spindle-like morphology, reminiscent of migrating cells within confined space and high cell density are associated with increased alignment and higher speed in migrating fibroblasts.

Cellular potts modeling of tumor growth, tumor invasion, and tumor evolution

Despite a growing wealth of available molecular data, the growth of tumors, invasion of tumors into healthy tissue, and response of tumors to therapies are still poorly understood. Although genetic mutations are in general the first step in the development of a cancer, for the mutated cell to persist in a tissue, it must compete against the other, healthy or diseased cells, for example by becoming more motile, adhesive, or multiplying faster. Thus, the cellular phenotype determines the success of a cancer cell in competition with its neighbors, irrespective of the genetic mutations or physiological alterations that gave rise to the altered phenotype. What phenotypes can make a cell "successful" in an environment of healthy and cancerous cells, and how? A widely used tool for getting more insight into that question is cell-based modeling. Cell-based models constitute a class of computational, agent-based models that mimic biophysical and molecular interactions between cells. One of the most widely used cell-based modeling formalisms is the cellular Potts model (CPM), a lattice-based, multi particle cell-based modeling approach. The CPM has become a popular and accessible method for modeling mechanisms of multicellular processes including cell sorting, gastrulation, or angiogenesis. The CPM accounts for biophysical cellular properties, including cell proliferation, cell motility, and cell adhesion, which play a key role in cancer. Multiscale models are constructed by extending the agents with intracellular processes including metabolism, growth, and signaling. Here we review the use of the CPM for modeling tumor growth, tumor invasion, and tumor progression. We argue that the accessibility and flexibility of the CPM, and its accurate, yet coarse-grained and computationally efficient representation of cell and tissue biophysics, make the CPM the method of choice for modeling cellular processes in tumor development.

Quantifying stretching and rearrangement in epithelial sheet migration

Although understanding the collective migration of cells, such as that seen in epithelial sheets, is essential for understanding diseases such as metastatic cancer, this motion is not yet as well characterized as individual cell migration. Here we adapt quantitative metrics used to characterize the flow and deformation of soft matter to contrast different types of motion within a migrating sheet of cells. Using a Finite-Time Lyapunov Exponent (FTLE) analysis, we find that - in spite of large fluctuations - the flow field of an epithelial cell sheet is not chaotic. Stretching of a sheet of cells (i.e., positive FTLE) is localized at the leading edge of migration and increases when the cells are more highly stimulated. By decomposing the motion of the cells into affine and non-affine components using the metric D2min , we quantify local plastic rearrangements and describe the motion of a group of cells in a novel way. We find an increase in plastic rearrangements with increasing cell densities, whereas inanimate systems tend to exhibit less non-affine rearrangements with increasing density.

The mechanics of metastasis: insights from a computational model

Although it may seem obvious that mechanical forces are required to drive metastatic cell movements, understanding of the mechanical aspects of metastasis has lagged far behind genetic and biochemical knowledge. The goal of this study is to learn about the mechanics of metastasis using a cell-based finite element model that proved useful for advancing knowledge about the forces that drive embryonic cell and tissue movements. Metastasis, the predominant cause of cancer-related deaths, involves a series of mechanical events in which one or more cells dissociate from a primary tumour, migrate through normal tissue, traverse in and out of a multi-layer circulatory system vessel and resettle. The present work focuses on the dissemination steps, from dissociation to circulation. The model shows that certain surface tension relationships must be satisfied for cancerous cells to dissociate from a primary tumour and that these equations are analogous to those that govern dissociation of embryonic cells. For a dissociated cell to then migrate by invadopodium extension and contraction and exhibit the shapes seen in experiments, the invadopodium must generate a contraction equal to approximately twice that produced by the interfacial tension associated with surrounding cells. Intravasation through the wall of a vessel is governed by relationships akin to those in the previous two steps, while release from the vessel wall is governed by equations that involve surface and interfacial tensions. The model raises a number of potential research questions. It also identifies how specific mechanical properties and the sub-cellular structural components that give rise to them might be changed so as to thwart particular metastatic steps and thereby block the spread of cancer.

Pattern formation during vasculogenesis

Vasculogenesis, the assembly of the first vascular network, is an intriguing developmental process that yields the first functional organ system of the embryo. In addition to being a fundamental part of embryonic development, vasculogenic processes also have medical importance. To explain the organizational principles behind vascular patterning, we must understand how morphogenesis of tissue level structures can be controlled through cell behavior patterns that, in turn, are determined by biochemical signal transduction processes. Mathematical analyses and computer simulations can help conceptualize how to bridge organizational levels and thus help in evaluating hypotheses regarding the formation of vascular networks. Here, we discuss the ideas that have been proposed to explain the formation of the first vascular pattern: cell motility guided by extracellular matrix alignment (contact guidance), chemotaxis guided by paracrine and autocrine morphogens, and sprouting guided by cell-cell contacts.