Reseach Interests
Reseach Interests
My main focus lies in the development and the numerical solution of mathematical models in biological and medical applications. In particular, my works include:
Multiscale models for cancer invasion and metastasis
Hybrid models of SDEs-PDEs for collective cellular motion
Mean field and Macroscopic limits in cancer models
Hyperbolic blood flow models on networks
Publications
Cell-cell adhesion is a key regulator of cancer invasion. In this work, we extend a pre-existing individual based cancer invasion model by introducing a stochastic representation of N-cadherin-mediated adhesion, where the lifetime of a cell-cell bond depends on the pulling force acting on the bond. Using experimental data, we derive expressions for the mean and standard deviation of N-cadherin bond lifetimes and fit them to Gamma distributions, enabling their treatment as force-dependent random variables. These distributions are then used to modify the diffusion coefficient of mesenchymal cancer cells. The model predicts reduced random motility with increasing adhesion and incorporates a dynamic transition between catch- and slip-bond behaviour. Along with this model for cell motility, we propose a preliminary physical framework, that can be used to model pattern formation as a result of the new adhesion mechanic.
In predicting metastatic potential and improving treatment outcomes in cancer research, it is crucial that we understand the dynamics of cancer cell dormancy and reactivation. In this paper we propose, study, and evaluate a cancer growth model that incorporates cell death, dormancy, reactivation, and proliferation in the secondary sites. Using experimental data from murine models, we test various statistical distributions and identify models that represent the asymmetry and variability observed in dormancy durations. Notably, the estimated cancer cell death rate remained consistent across all tested distributions, supporting its biological relevance as a robust parameter for modelling dormancy survival dynamics. When post-reactivation cell death is present, the most suitable among the distributions we studied exhibit heavy-tails and asymmetric skewness; this aligns with the prolonged and rare dormancy periods expected of cancer cells. In contrast, analysis of complementary data in which no cancer cell death is observed indicates that dormancy can be equally well represented by distributions with finite variance and comparatively light tails. Our findings highlight a close interplay between reactivation/dormancy-time and cancer cell death mechanisms. Furthermore, they stress the importance of selecting appropriate statistical models for dormancy, both in predicting cancer cell reactivation, and in informing therapeutic strategies that focus on dormancy-driven metastasis.
2. D. Katsaounis, N. Harbour, T. Williams, M.A.J. Chaplain & N. Sfakianakis, A genuinely hybrid, multiscale 3D cancer invasion and metastasis modelling framework. Bull Math Biol 86, 64 (2024). https://doi.org/10.1007/s11538-024-01286-0
We introduce in this paper substantial enhancements to a previously proposed hybrid multiscale cancer invasion modelling framework to better reflect the biological reality and dynamics of cancer. These model updates contribute to a more accurate representation of cancer dynamics, they provide deeper insights and enhance our predictive capabilities. Key updates include the integration of porous medium-like diffusion for the evolution of Epithelial-like Cancer Cells and other essential cellular constituents of the system, more realistic modelling of Epithelial-Mesenchymal Transition and Mesenchymal-Epithelial Transition models with the inclusion of Transforming Growth Factor beta within the tumour microenvironment, and the introduction of Compound Poisson Process in the Stochastic Differential Equations that describe the migration behaviour of the Mesenchymal-like Cancer Cells. Another innovative feature of the model is its extension into a multi-organ metastatic framework. This framework connects various organs through a circulatory network, enabling the study of how cancer cells spread to secondary sites.
Invasion of the surrounding tissue is a key aspect of cancer growth and spread involving a coordinated effort between cell migration and matrix degradation, and has been the subject of mathematical modelling for almost 30 years. In this current paper we address a long-standing question in the field of cancer cell migration modelling. Namely, identify the migratory pattern and spread of individual cancer cells, or small clusters of cancer cells, when the macroscopic evolution of the cancer cell colony is dictated by a specific partial differential equation (PDE). We show that the usual heuristic understanding of the diffusion and advection terms of the PDE being one-to-one responsible for the random and biased motion of the solitary cancer cells, respectively, is not precise. On the contrary, we show that the drift term of the correct stochastic differential equation (SDE) scheme that dictates the individual cancer cell migration, should account also for the divergence of the diffusion of the PDE. We support our claims with a number of numerical experiments and computational simulations.