Computational Model for Migration of a Cell Cluster in Three-Dimensional Matrices

Mechanochemical Coupling and Junctional Forces during Collective Cell Migration

Cell migration, a fundamental physiological process in which cells sense and move through their surrounding physical environment, plays a critical role in development and tissue formation, as well as pathological processes, such as cancer metastasis and wound healing. During cell migration, dynamics are governed by the bidirectional interplay between cell-generated mechanical forces and the activity of Rho GTPases, a family of small GTP-binding proteins that regulate actin cytoskeleton assembly and cellular contractility. These interactions are inherently more complex during the collective migration of mechanically coupled cells because of the additional regulation of cell-cell junctional forces. In this study, we adapted a recent minimal modeling framework to simulate the interactions between mechanochemical signaling in individual cells and interactions with cell-cell junctional forces during collective cell migration. We find that migration of individual cells depends on the feedback between mechanical tension and Rho GTPase activity in a biphasic manner. During collective cell migration, waves of Rho GTPase activity mediate mechanical contraction/extension and thus synchronization throughout the tissue. Further, cell-cell junctional forces exhibit distinct spatial patterns during collective cell migration, with larger forces near the leading edge. Larger junctional force magnitudes are associated with faster collective cell migration and larger tissue size. Simulations of heterogeneous tissue migration exhibit a complex dependence on the properties of both leading and trailing cells. Computational predictions demonstrate that collective cell migration depends on both the emergent dynamics and interactions between cellular-level Rho GTPase activity and contractility and multicellular-level junctional forces.

Phenomenological modelling and simulation of cell clusters in 3D cultures

Cell clustering and aggregation are fundamental processes in the development of several tissues and the progression of many diseases. The formation of these aggregates also has a direct impact on the oxygen concentration in their surroundings due to cellular respiration and poor oxygen diffusion through clusters. In this work, we propose a mathematical model that is capable of simulating cell cluster formation in 3D cultures through combining a particle-based and a finite element approach to recreate complex experimental conditions. Cells are modelled considering cell proliferation, cell death and cell-cell mechanical interactions. Additionally, the oxygen concentration profile is calculated through finite element analysis using a reaction-diffusion model that considers cell oxygen consumption and diffusion through the extracellular matrix and the cell clusters. In our model, the local oxygen concentration in the medium determines both cell proliferation and cell death. Numerical predictions are also compared with experimental data from the literature. The simulation results indicate that our model can predict cell clustering, cluster growth and oxygen distribution in 3D cultures. We conclude that the initial cell distribution, cell death and cell proliferation dynamics determine the size and density of clusters. Moreover, these phenomena are directly affected by the oxygen transport in the 3D culture.

A mathematical model of collective cell migration in a three-dimensional, heterogeneous environment

Cell migration is essential in animal development, homeostasis, and disease progression, but many questions remain unanswered about how this process is controlled. While many kinds of individual cell movements have been characterized, less effort has been directed towards understanding how clusters of cells migrate collectively through heterogeneous, cellular environments. To explore this, we have focused on the migration of the border cells during Drosophila egg development. In this case, a cluster of different cell types coalesce and traverse as a group between large cells, called nurse cells, in the center of the egg chamber. We have developed a new model for this collective cell migration based on the forces of adhesion, repulsion, migration and stochastic fluctuation to generate the movement of discrete cells. We implement the model using Identical Math Cells, or IMCs. IMCs can each represent one biological cell of the system, or can be aggregated using increased adhesion forces to model the dynamics of larger biological cells. The domain of interest is filled with IMCs, each assigned specific biophysical properties to mimic a diversity of cell types. Using this system, we have successfully simulated the migration of the border cell cluster through an environment filled with larger cells, which represent nurse cells. Interestingly, our simulations suggest that the forces utilized in this model are sufficient to produce behaviors of the cluster that are observed in vivo, such as rotation. Our framework was developed to capture a heterogeneous cell population, and our implementation strategy allows for diverse, but precise, initial position specification over a three- dimensional domain. Therefore, we believe that this model will be useful for not only examining aspects of Drosophila oogenesis, but also for modeling other two or three-dimensional systems that have multiple cell types and where investigating the forces between cells is of interest.

Computational model to probe cellular mechanics during epithelial-mesenchymal transition

During the epithelial to mesenchymal transition (EMT), polarized cells in the epithelium can undergo a transformation characterized by the loss of cell-cell junctions and increased migratory activity into nonpolarized invasive cells. These cells adopt a mesenchymal shape and migrate into the basal lamina. Such transitions have been observed in developmental processes and have been linked to cancer cell metastasis. Most experimental studies on EMT search for molecular markers indicating an epithelial or mesenchymal conformation, focussing on afferent signaling pathways received by cells undergoing this transformation; however, these approaches are unable to track mechanical changes in the cell and the possible role this plays in EMT. In order to address this gap in our understanding, we have used a quantitative approach to study population level effects of single cell changes typically occurring during EMT. We have developed a computational model making use of the advantages of both single cell migratory models and agent-based cell population models to study the effect of cellular molecular processes in EMT. The disruption of a cell sheet representing the epithelium over a dense extracellular matrix (ECM) is simulated using interaction forces between different cells and between cells and discrete fibers representing the ECM. In our study, two different parameters were varied: protrusion force magnitude and E-cadherin (cell-cell junction) concentration. The cell population was tracked for 3 days and the number of cells that leave the layer, the depth of invasion, and the percentage of initial number of cells that remain in the layer (a measure of epithelium disruption) were monitored. Our studies suggest that having a high protrusion force or a reduction in cell-cell attachments is enough to cause EMT. Our results also demonstrate that the morphological progression in membrane disruption has an effect on the number of cells becoming invasive, with epithelial layers broken into clusters hindering the further exodus of cells. The results reveal the quantitative interplay between two key parameters involved in EMT and suggest potential avenues for further exploration of a systems level understanding of EMT.

A computational model for collective cellular motion in three dimensions: general framework and case study for cell pair dynamics

Cell migration in healthy and diseased systems is a combination of single and collective cell motion. While single cell motion has received considerable attention, our understanding of collective cell motion remains elusive. A new computational framework for the migration of groups of cells in three dimensions is presented, which focuses on the forces acting at the microscopic scale and the interactions between cells and their extracellular matrix (ECM) environment. Cell-cell adhesion, resistance due to the ECM and the factors regulating the propulsion of each cell through the matrix are considered. In particular, our approach emphasizes the role of receptors that mediate cell-cell and cell-matrix interactions, and examines how variation in their properties induces changes in cellular motion. As an important case study, we analyze two interacting cells. Our results show that the dynamics of cell pairs depends on the magnitude and the stochastic nature of the forces. Stronger intercellular stability is generally promoted by surface receptors that move. We also demonstrate that matrix resistance, cellular stiffness and intensity of adhesion contribute to migration behaviors in different ways, with memory effects present that can alter pair motility. If adhesion weakens with time, our findings show that cell pair break-up depends strongly on the way cells interact with the matrix. Finally, the motility for cells in a larger cluster (size 50 cells) is examined to illustrate the full capabilities of the model and to stress the role of cellular pairs in complex cellular structures. Overall, our framework shows how properties of cells and their environment influence the stability and motility of cellular assemblies. This is an important step in the advancement of the understanding of collective motility, and can contribute to knowledge of complex biological processes involving migration, aggregation and detachment of cells in healthy and diseased systems.