Multiscale dynamics of a heterotypic cancer cell population within a fibrous extracellular matrix

Dynamics of Fibril Collagen Remodeling by Tumor Cells: A Model of Tumor-Associated Collagen Signatures

Many solid tumors are characterized by a dense extracellular matrix (ECM) composed of various ECM fibril proteins. These proteins provide structural support and a biological context for the residing cells. The reciprocal interactions between growing and migrating tumor cells and the surrounding stroma result in dynamic changes in the ECM architecture and its properties. With the use of advanced imaging techniques, several specific patterns in the collagen surrounding the breast tumor have been identified in both tumor murine models and clinical histology images. These tumor-associated collagen signatures (TACS) include loosely organized fibrils far from the tumor and fibrils aligned either parallel or perpendicular to tumor colonies. They are correlated with tumor behavior, such as benign growth or invasive migration. However, it is not fully understood how one specific fibril pattern can be dynamically remodeled to form another alignment. Here, we present a novel multi-cellular lattice-free (MultiCell-LF) agent-based model of ECM that, in contrast to static histology images, can simulate dynamic changes between TACSs. This model allowed us to identify the rules of cell-ECM physical interplay and feedback that guided the emergence and transition among various TACSs.

Nonlocal multiscale modelling of tumour-oncolytic viruses interactions within a heterogeneous fibrous/non-fibrous extracellular matrix

In this study we investigate computationally tumour-oncolytic virus (OV) interactions that take place within a heterogeneous extracellular matrix (ECM). The ECM is viewed as a mixture of two constitutive phases, namely a fibre phase and a non-fibre phase. The multiscale mathematical model presented here focuses on the nonlocal cell-cell and cell-ECM interactions, and how these interactions might be impacted by the infection of cancer cells with the OV. At macroscale we track the kinetics of cancer cells, virus particles and the ECM. At microscale we track (i) the degradation of ECM by matrix degrading enzymes (MDEs) produced by cancer cells, which further influences the movement of tumour boundary; (ii) the re-arrangement of the microfibres that influences the re-arrangement of macrofibres (i.e., fibres at macroscale). With the help of this new multiscale model, we investigate two questions: (i) whether the infected cancer cell fluxes are the result of local or non-local advection in response to ECM density; and (ii) what is the effect of ECM fibres on the the spatial spread of oncolytic viruses and the outcome of oncolytic virotherapy.

Inverse problem approaches for mutation laws in heterogeneous tumours with local and nonlocal dynamics

Cancer cell mutations occur when cells undergo multiple cell divisions, and these mutations can be spontaneous or environmentally-induced. The mechanisms that promote and sustain these mutations are still not fully understood. This study deals with the identification (or reconstruction) of the usually unknown cancer cell mutation law, which lead to the transformation of a primary tumour cell population into a secondary, more aggressive cell population. We focus on local and nonlocal mathematical models for cell dynamics and movement, and identify these mutation laws from macroscopic tumour snapshot data collected at some later stage in the tumour evolution. In a local cancer invasion model, we first reconstruct the mutation law when we assume that the mutations depend only on the surrounding cancer cells (i.e., the ECM plays no role in mutations). Second, we assume that the mutations depend on the ECM only, and we reconstruct the mutation law in this case. Third, we reconstruct the mutation when we assume that there is no prior knowledge about the mutations. Finally, for the nonlocal cancer invasion model, we reconstruct the mutation law that depends on the cancer cells and on the ECM. For these numerical reconstructions, our approximations are based on the finite difference method combined with the finite elements method. As the inverse problem is ill-posed, we use the Tikhonov regularisation technique in order to regularise the solution. Stability of the solution is examined by adding additive noise into the measurements.

Multiscale modeling in disease

Multiscale computational modeling aims to connect the complex networks of effects at different length and/or time scales. For example, these networks often include intracellular molecular signaling, crosstalk, and other interactions between neighboring cell populations, and higher levels of emergent phenomena across different regions of tissues and among collections of tissues or organs interacting with each other in the whole body. Recent applications of multiscale modeling across intracellular, cellular, and/or tissue levels are highlighted here. These models incorporated the roles of biochemical and biomechanical modulation in processes that are implicated in the mechanisms of several diseases including fibrosis, joint and bone diseases, respiratory infectious diseases, and cancers.

Non-local multiscale approach for the impact of go or grow hypothesis on tumour-viruses interactions

We propose and study computationally a novel non-local multiscale moving boundary mathematical model for tumour and oncolytic virus (OV) interactions when we consider the go or grow hypothesis for cancer dynamics. This spatio-temporal model focuses on two cancer cell phenotypes that can be infected with the OV or remain uninfected, and which can either move in response to the extracellular-matrix (ECM) density or proliferate. The interactions between cancer cells, those among cancer cells and ECM, and those among cells and OV occur at the macroscale. At the micro-scale, we focus on the interactions between cells and matrix degrading enzymes (MDEs) that impact the movement of tumour boundary. With the help of this multiscale model we explore the impact on tumour invasion patterns of two different assumptions that we consider in regard to cell-cell and cell-matrix interactions. In particular we investigate model dynamics when we assume that cancer cell fluxes are the result of local advection in response to the density of extracellular matrix (ECM), or of non-local advection in response to cell-ECM adhesion. We also investigate the role of the transition rates between mainly-moving and mainly-growing cancer cell sub-populations, as well as the role of virus infection rate and virus replication rate on the overall tumour dynamics.

Collective Cell Migration in a Fibrous Environment: A Hybrid Multiscale Modelling Approach

The specific structure of the extracellular matrix (ECM), and in particular the density and orientation of collagen fibres, plays an important role in the evolution of solid cancers. While many experimental studies discussed the role of ECM in individual and collective cell migration, there are still unanswered questions about the impact of nonlocal cell sensing of other cells on the overall shape of tumour aggregation and its migration type. There are also unanswered questions about the migration and spread of tumour that arises at the boundary between different tissues with different collagen fibre orientations. To address these questions, in this study we develop a hybrid multi-scale model that considers the cells as individual entities and ECM as a continuous field. The numerical simulations obtained through this model match experimental observations, confirming that tumour aggregations are not moving if the ECM fibres are distributed randomly, and they only move when the ECM fibres are highly aligned. Moreover, the stationary tumour aggregations can have circular shapes or irregular shapes (with finger-like protrusions), while the moving tumour aggregations have elongate shapes (resembling to clusters, strands or files). We also show that the cell sensing radius impacts tumour shape only when there is a low ratio of fibre to non-fibre ECM components. Finally, we investigate the impact of different ECM fibre orientations corresponding to different tissues, on the overall tumour invasion of these neighbouring tissues.

Directionality of Macrophages Movement in Tumour Invasion: A Multiscale Moving-Boundary Approach

Invasion of the surrounding tissue is one of the recognised hallmarks of cancer (Hanahan and Weinberg in Cell 100: 57–70, 2000. 10.1016/S0092-8674(00)81683-9), which is accomplished through a complex heterotypic multiscale dynamics involving tissue-scale random and directed movement of the population of both cancer cells and other accompanying cells (including here, the family of tumour-associated macrophages) as well as the emerging cell-scale activity of both the matrix-degrading enzymes and the rearrangement of the cell-scale constituents of the extracellular matrix (ECM) fibres. The involved processes include not only the presence of cell proliferation and cell adhesion (to other cells and to the extracellular matrix), but also the secretion of matrix-degrading enzymes. This is as a result of cancer cells as well as macrophages, which are one of the most abundant types of immune cells in the tumour micro-environment. In large tumours, these tumour-associated macrophages (TAMs) have a tumour-promoting phenotype, contributing to tumour proliferation and spread. In this paper, we extend a previous multiscale moving-boundary mathematical model for cancer invasion, by considering also the multiscale effects of TAMs, with special focus on the influence that their directional movement exerts on the overall tumour progression. Numerical investigation of this new model shows the importance of the interactions between pro-tumour TAMs and the fibrous ECM, highlighting the impact of the fibres on the spatial structure of solid tumour.

Cell-Scale Degradation of Peritumoural Extracellular Matrix Fibre Network and Its Role Within Tissue-Scale Cancer Invasion

Local cancer invasion of tissue is a complex, multiscale process which plays an essential role in tumour progression. During the complex interaction between cancer cell population and the extracellular matrix (ECM), of key importance is the role played by both bulk two-scale dynamics of ECM fibres within collective movement of the tumour cells and the multiscale leading edge dynamics driven by proteolytic activity of the matrix-degrading enzymes (MDEs) that are secreted by the cancer cells. As these two multiscale subsystems share and contribute to the same tumour macro-dynamics, in this work we develop further the model introduced in Shuttleworth and Trucu (Bull Math Biol 81:2176–2219, 2019. 10.1007/s11538-019-00598-w) by exploring a new aspect of their interaction that occurs at the cell scale. Specifically, here we will focus on understanding the cell-scale cross talk between the micro-scale parts of these two multiscale subsystems which get to interact directly in the peritumoural region, with immediate consequences both for MDE micro-dynamics occurring at the leading edge of the tumour and for the cell-scale rearrangement of the naturally oriented ECM fibres in the peritumoural region, ultimately influencing the way tumour progresses in the surrounding tissue. To that end, we will propose a new modelling that captures the ECM fibres degradation not only at macro-scale in the bulk of the tumour but also explicitly in the micro-scale neighbourhood of the tumour interface as a consequence of the interactions with molecular fluxes of MDEs that exercise their spatial dynamics at the invasive edge of the tumour.