A stochastic algorithm for accurately predicting path persistence of cells migrating in 3D matrix environments

Tumor proliferation and invasion are intrinsically coupled and unraveled through tunable spheroid and physics-based models

Proliferation and invasion are two key drivers of tumor growth that are traditionally considered independent multicellular processes. However, these processes are intrinsically coupled through a maximum carrying capacity, i.e., the maximum spatial cell concentration supported by the tumor volume, total cell count, nutrient access, and mechanical properties of the tissue stroma. We explored this coupling of proliferation and invasion through in vitro and in silico methods where we modulated the mechanical properties of the tumor and the surrounding extracellular matrix. E-cadherin expression and stromal collagen concentration were manipulated in a tunable breast cancer spheroid to determine the overall impacts of these tumor variables on net tumor proliferation and continuum invasion. We integrated these results into a mixed-constitutive formulation to computationally delineate the influences of cellular and extracellular adhesion, stiffness, and mechanical properties of the extracellular matrix on net proliferation and continuum invasion. This framework integrates biological in vitro data into concise computational models of invasion and proliferation to provide more detailed physical insights into the coupling of these key tumor processes and tumor growth. STATEMENT OF SIGNIFICANCE: Tumor growth involves expansion into the collagen-rich stroma through intrinsic coupling of proliferation and invasion within the tumor continuum. These processes are regulated by a maximum carrying capacity that is determined by the total cell count, tumor volume, nutrient access, and mechanical properties of the surrounding stroma. The influences of biomechanical parameters (i.e., stiffness, cell elongation, net proliferation rate and cell-ECM friction) on tumor proliferation or invasion cannot be unraveled using experimental methods alone. By pairing a tunable spheroid system with computational modeling, we delineated the interdependencies of each system parameter on tumor proliferation and continuum invasion, and established a concise computational framework for studying tumor mechanobiology.

Thermotaxis in an apolar, non-neuronal animal

Neuronal circuits are hallmarks of complex decision-making processes in the animal world. How animals without neurons process information and respond to environmental cues promises a new window into studying precursors of neuronal control and origin of the nervous system as we know it today. Robust decision making in animals, such as in chemotaxis or thermotaxis, often requires internal symmetry breaking (such as anterior-posterior (AP) axis) provided naturally by a given body plan of an animal. Here we report the discovery of robust thermotaxis behaviour in Trichoplax adhaerens, an early-divergent, enigmatic animal with no anterior-posterior symmetry breaking (apolar) and no known neurons or muscles. We present a quantitative and robust behavioural response assay in Placozoa, which presents an apolar flat geometry. By exposing T. adhaerens to a thermal gradient under a long-term imaging set-up, we observe robust thermotaxis that occurs over timescale of hours, independent of any circadian rhythms. We quantify that T. adhaerens can detect thermal gradients of at least 0.1°C cm-1. Positive thermotaxis is observed for a range of baseline temperatures from 17°C to 22.5°C, and distributions of momentary speeds for both thermotaxis and control conditions are well described by single exponential fits. Interestingly, the organism does not maintain a fixed orientation while performing thermotaxis. Using natural diversity in size of adult organisms (100 µm to a few millimetres), we find no apparent size-dependence in thermotaxis behaviour across an order of magnitude of organism size. Several transient receptor potential (TRP) family homologues have been previously reported to be conserved in metazoans, including in T. adhaerens. We discover naringenin, a known TRPM3 antagonist, inhibits thermotaxis in T. adhaerens. The discovery of robust thermotaxis in T. adhaerens provides a tractable handle to interrogate information processing in a brainless animal. Understanding how divergent marine animals process thermal cues is also critical due to rapid temperature rise in our oceans.

Contact inhibition of locomotion generates collective cell migration without chemoattractants in an open domain

Neural crest cells are embryonic stem cells that migrate throughout embryos and, at different target locations, give rise to the formation of a variety of tissues and organs. The directional migration of the neural crest cells is experimentally described using a process referred to as contact inhibition of locomotion, by which cells redirect their movement upon the cell-cell contacts. However, it is unclear how the migration alignment is affected by the motility properties of the cells. Here, we theoretically model the migration alignment as a function of the motility dynamics and interaction of the cells in an open domain with a channel geometry. The results indicate that by increasing the influx rate of the cells into the domain a transition takes place from random movement to an organized collective migration, where the migration alignment is maximized and the migration time is minimized. This phase transition demonstrates that the cells can migrate efficiently over long distances without any external chemoattractant information about the direction of migration just based on local interactions with each other. The analysis of the dependence of this transition on the characteristic properties of cellular motility shows that the cell density determines the coordination of collective migration whether the migration domain is open or closed. In the open domain, this density is determined by a feedback mechanism between the flux and order parameter, which characterises the alignment of collective migration. The model also demonstrates that the coattraction mechanism proposed earlier is not necessary for collective migration and a constant flux of cells moving into the channel is sufficient to produce directed movement over arbitrary long distances.

Adhesion strength and contractility enable metastatic cells to become adurotactic

Significant changes in cell stiffness, contractility, and adhesion, i.e., mechanotype, are observed during a variety of biological processes. Whether cell mechanics merely change as a side effect of or driver for biological processes is still unclear. Here, we sort genotypically similar metastatic cancer cells into strongly adherent (SA) versus weakly adherent (WA) phenotypes to study how contractility and adhesion differences alter the ability of cells to sense and respond to gradients in material stiffness. We observe that SA cells migrate up a stiffness gradient, or durotax, while WA cells largely ignore the gradient, i.e., adurotax. Biophysical modeling and experimental validation suggest that differences in cell migration and durotaxis between weakly and strongly adherent cells are driven by differences in intra-cellular actomyosin activity. These results provide a direct relationship between cell phenotype and durotaxis and suggest how, unlike other senescent cells, metastatic cancer cells navigate against stiffness gradients.

Computational models of migration modes improve our understanding of metastasis

Tumor cells migrate through changing microenvironments of diseased and healthy tissue, making their migration particularly challenging to describe. To better understand this process, computational models have been developed for both the ameboid and mesenchymal modes of cell migration. Here, we review various approaches that have been used to account for the physical environment's effect on cell migration in computational models, with a focus on their application to understanding cancer metastasis and the related phenomenon of durotaxis. We then discuss how mesenchymal migration models typically simulate complex cell–extracellular matrix (ECM) interactions, while ameboid migration models use a cell-focused approach that largely ignores ECM when not acting as a physical barrier. This approach greatly simplifies or ignores the mechanosensing ability of ameboid migrating cells and should be reevaluated in future models. We conclude by describing future model elements that have not been included to date but would enhance model accuracy.

Matrix degradation regulates osteoblast protrusion dynamics and individual migration

Protrusions are one of the structures that cells use to sense their surrounding environment in a probing and exploratory manner as well as to communicate with other cells. In particular, osteoblasts embedded within a 3D matrix tend to originate a large number of protrusions compared to other type of cells. In this work, we study the role that mechanochemical properties of the extracellular matrix (ECM) play on the dynamics of these protrusions, namely, the regulation of the size and number of emanating structures. In addition, we also determine how the dynamics of the protrusions may lead the 3D movement of the osteoblasts. Significant differences were found in protrusion size and cell velocity, when degradation activity due to metalloproteases was blocked by means of an artificial broad-spectrum matrix metalloproteinase inhibitor, whereas stiffening of the matrix by introducing transglutaminase crosslinking, only induced slight changes in both protrusion size and cell velocity, suggesting that the ability of cells to create a path through the matrix is more critical than the matrix mechanical properties themselves. To confirm this, we developed a cell migration computational model in 3D including both the mechanical and chemical properties of the ECM as well as the protrusion mechanics, obtaining good agreement with experimental results.

To lead or to herd: optimal strategies for 3D collective migration of cell clusters

Cells migrating in clusters play a significant role in a number of biological processes such as embryogenesis, wound healing, and tumor metastasis during cancer progression. A variety of environmental and biochemical factors can influence the collective migration of cells with differing degrees of cell autonomy and inter-cellular coupling strength. For example, weakly coupled cells can move collectively under the influence of contact guidance from neighboring cells or the environment. Alternatively strongly coupled cells might follow one or more leader cells to move as a single cohesive unit. Additionally, chemical and mechanical signaling between these cells may alter the degree of coupling and determine effective cluster sizes. Being able to understand this collective cell migration process is critical in the prediction and manipulation of outcomes of key biological processes. Here we focus on understanding how various environmental and cellular factors influence small clusters of cells migrating collectively within a 3D fibrous matrix. We combine existing knowledge of single-cell migration in 2D and 3D environments, prior experimental observations of cell-cell interactions and collective migration, and a newly developed stochastic model of cell migration in 3D matrices, to simulate the migration of cell clusters in different physiologically relevant environments. Our results show that based on the extracellular environment and the strength of cell-cell mechanical coupling, two distinct optimal approaches to driving collective cell migration emerge. The ability to effectively employ these two distinct migration strategies might be critical for cells to collectively migrate through the heterogeneous tissue environments within the body.