Source author record

Tommaso Lorenzi

Tommaso Lorenzi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

9works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2022arXiv

Individual-based and continuum models of phenotypically heterogeneous growing cell populations

Existing studies comparing individual-based models of growing cell populations and their continuum counterparts have mainly focused on homogeneous populations, in which all cells have the same phenotypic characteristics. However, significant intercellular phenotypic variability is commonly observed in cellular systems. Therefore, we develop here an individual-based model for the growth of phenotypically heterogeneous cell populations. In this model, the phenotypic state of each cell is described by a structuring variable that captures intercellular variability in cell proliferation and migration rates. The model tracks the spatial evolutionary dynamics of single cells, which undergo pressure-dependent proliferation, heritable phenotypic changes and directional movement in response to pressure differentials. We formally show that the continuum limit of this model comprises a non-local partial differential equation for the cell population density, which generalises earlier models of growing cell populations. Results of the individual-based model illustrate how proliferation-migration tradeoffs shaping the evolution of single cells can lead to the formation of travelling waves at the population level where highly-mobile cells locally dominate at the invasive front, while more-proliferative cells are found at the rear. We demonstrate that there is an excellent quantitative agreement between these results and the results of numerical simulations and formal travelling-wave analysis of the continuum model, when sufficiently large cell numbers are considered. We provide numerical evidence of scenarios in which the predictions of the two models may differ due to demographic stochasticity, which cannot be captured by the continuum model. This indicates the importance of integrating individual-based and continuum approaches when modelling the growth of phenotypically heterogeneous cell populations.

preprint2021arXiv

Trade-offs between chemotaxis and proliferation shape the phenotypic structuring of invading waves

Chemotaxis-driven invasions have been proposed across a broad spectrum of biological processes, from cancer to ecology. The influential system of equations introduced by Keller and Segel has proven a popular choice in the modelling of such phenomena, but in its original form restricts to a homogeneous population. To account for the possibility of phenotypic heterogeneity, we extend to the case of a population continuously structured across space, time and phenotype, where the latter determines variation in chemotactic responsiveness, proliferation rate, and the level of chemical environment modulation. The extended model considered here comprises a non-local partial differential equation for the local phenotype distribution of cells which is coupled, through an integral term, with a differential equation for the concentration of an attractant, which is sensed and degraded by the cells. In the framework of this model, we concentrate on a chemotaxis/proliferation trade-off scenario, where the cell phenotypes span a spectrum of states from highly-chemotactic but minimally-proliferative to minimally-chemotactic but highly-proliferative. Using a combination of numerical simulation and formal asymptotic analysis, we explore the properties of travelling-wave solutions. The results of our study demonstrate how incorporating phenotypic heterogeneity may lead to a highly-structured wave profile, where cells in different phenotypic states dominate different spatial positions across the invading wave, and clarify how the phenotypic structuring of the wave can be shaped by trade-offs between chemotaxis and proliferation.

preprint2021arXiv

Travelling-wave analysis of a model of tumour invasion with degenerate, cross-dependent diffusion

In this paper, we carry out a travelling-wave analysis of a model of tumour invasion with degenerate, cross-dependent diffusion. We consider two types of invasive fronts of tumour tissue into extracellular matrix (ECM), which represents healthy tissue. These types differ according to whether the density of ECM far ahead of the wave front is maximal or not. In the former case, we use a shooting argument to prove that there exists a unique travelling wave solution for any positive propagation speed. In the latter case, we further develop this argument to prove that there exists a unique travelling wave solution for any propagation speed greater than or equal to a strictly positive minimal wave speed. Using a combination of analytical and numerical results, we conjecture that the minimal wave speed depends monotonically on the degradation rate of ECM by tumour cells and the ECM density far ahead of the front.

preprint2020arXiv

Evolutionary dynamics in vascularised tumours under chemotherapy

We consider a mathematical model for the evolutionary dynamics of tumour cells in vascularised tumours under chemotherapy. The model comprises a system of coupled partial integro-differential equations for the phenotypic distribution of tumour cells, the concentration of oxygen and the concentration of a chemotherapeutic agent. In order to disentangle the impact of different evolutionary parameters on the emergence of intra-tumour phenotypic heterogeneity and the development of resistance to chemotherapy, we construct explicit solutions to the equation for the phenotypic distribution of tumour cells and provide a detailed quantitative characterisation of the long-time asymptotic behaviour of such solutions. Analytical results are integrated with numerical simulations of a calibrated version of the model based on biologically consistent parameter values. The results obtained provide a theoretical explanation for the observation that the phenotypic properties of tumour cells in vascularised tumours vary with the distance from the blood vessels. Moreover, we demonstrate that lower oxygen levels may correlate with higher levels of phenotypic variability, which suggests that the presence of hypoxic regions supports intra-tumour phenotypic heterogeneity. Finally, the results of our analysis put on a rigorous mathematical basis the idea, previously suggested by formal asymptotic results and numerical simulations, that hypoxia favours the selection for chemoresistant phenotypic variants prior to treatment. Consequently, this facilitates the development of resistance following chemotherapy.

preprint2020arXiv

From a discrete model of chemotaxis with volume-filling to a generalised Patlak-Keller-Segel model

We present a discrete model of chemotaxis whereby cells responding to a chemoattractant are seen as individual agents whose movement is described through a set of rules that result in a biased random walk. In order to take into account possible alterations in cellular motility observed at high cell densities (i.e. volume-filling), we let the probabilities of cell movement be modulated by a decaying function of the cell density. We formally show that a general form of the celebrated Patlak-Keller-Segel (PKS) model of chemotaxis can be formally derived as the appropriate continuum limit of this discrete model. The family of steady-state solutions of such a generalised PKS model are characterised and the conditions for the emergence of spatial patterns are studied via linear stability analysis. Moreover, we carry out a systematic quantitative comparison between numerical simulations of the discrete model and numerical solutions of the corresponding PKS model, both in one and in two spatial dimensions. The results obtained indicate that there is excellent quantitative agreement between the spatial patterns produced by the two models. Finally, we numerically show that the outcomes of the two models faithfully replicate those of the classical PKS model in a suitable asymptotic regime.

preprint2020arXiv

From individual-based mechanical models of multicellular systems to free-boundary problems

In this paper we present an individual-based mechanical model that describes the dynamics of two contiguous cell populations with different proliferative and mechanical characteristics. An off-lattice modelling approach is considered whereby: (i) every cell is identified by the position of its centre; (ii) mechanical interactions between cells are described via generic nonlinear force laws; and (iii) cell proliferation is contact inhibited. We formally show that the continuum counterpart of this discrete model is given by a free-boundary problem for the cell densities. The results of the derivation demonstrate how the parameters of continuum mechanical models of multicellular systems can be related to biophysical cell properties. We prove an existence result for the free-boundary problem and construct travelling-wave solutions. Numerical simulations are performed in the case where the cellular interaction forces are described by the celebrated Johnson-Kendall-Roberts model of elastic contact, which has been previously used to model cell-cell interactions. The results obtained indicate excellent agreement between the simulation results for the individual-based model, the numerical solutions of the corresponding free-boundary problem and the travelling-wave analysis.

preprint2019arXiv

Evolutionary dynamics of competing phenotype-structured populations in periodically fluctuating environments

Living species, ranging from bacteria to animals, exist in environmental conditions that exhibit spatial and temporal heterogeneity which requires them to adapt. Risk-spreading through spontaneous phenotypic variations is a known concept in ecology, which is used to explain how species may survive when faced with the evolutionary risks associated with temporally varying environments. In order to support a deeper understanding of the adaptive role of spontaneous phenotypic variations in fluctuating environments, we consider a system of non-local partial differential equations modelling the evolutionary dynamics of two competing phenotype-structured populations in the presence of periodically oscillating nutrient levels. The two populations undergo spontaneous phenotypic variations at different rates. The phenotypic state of each individual is represented by a continuous variable, and the phenotypic landscape of the populations evolves in time due to variations in the nutrient level. Exploiting the analytical tractability of our model, we study the long-time behaviour of the solutions to obtain a detailed mathematical depiction of evolutionary dynamics. The results suggest that when nutrient levels undergo small and slow oscillations, it is evolutionarily more convenient to rarely undergo spontaneous phenotypic variations. Conversely, under relatively large and fast periodic oscillations in the nutrient levels, which bring about alternating cycles of starvation and nutrient abundance, higher rates of spontaneous phenotypic variations confer a competitive advantage. We discuss the implications of our results in the context of cancer metabolism.

preprint2013arXiv

Effects of space structure and combination therapies on phenotypic heterogeneity and drug resistance in solid tumors

Histopathological evidence supports the idea that the emergence of phenotypic heterogeneity and resistance to cytotoxic drugs can be considered as a process of adaptation, or evolution, in tumor cell populations. In this framework, can we explain intra-tumor heterogeneity in terms of cell adaptation to local conditions? How do anti-cancer therapies affect the outcome of cell competition for nutrients within solid tumors? Can we overcome the emergence of resistance and favor the eradication of cancer cells by using combination therapies? Bearing these questions in mind, we develop a model describing cell dynamics inside a tumor spheroid under the effects of cytotoxic and cytostatic drugs. Cancer cells are assumed to be structured as a population by two real variables standing for space position and the expression level of a cytotoxic resistant phenotype. The model takes explicitly into account the dynamics of resources and anti-cancer drugs as well as their interactions with the cell population under treatment. We analyze the effects of space structure and combination therapies on phenotypic heterogeneity and chemotherapeutic resistance. Furthermore, we study the efficacy of combined therapy protocols based on constant infusion and/or bang-bang delivery of cytotoxic and cytostatic drugs.

preprint2012arXiv

Populational adaptive evolution, chemotherapeutic resistance and multiple anti-cancer therapies

Resistance to chemotherapies, particularly to anticancer treatments, is an increasing medical concern. Among the many mechanisms at work in cancers, one of the most important is the selection of tumor cells expressing resistance genes or phenotypes. Motivated by the theory of mutation-selection in adaptive evolution, we propose a model based on a continuous variable that represents the expression level of a resistance gene (or genes, yielding a phenotype) influencing in healthy and tumor cells birth/death rates, effects of chemotherapies (both cytotoxic and cytostatic) and mutations. We extend previous work by demonstrating how qualitatively different actions of chemotherapeutic and cytostatic treatments may induce different levels of resistance. The mathematical interest of our study is in the formalism of constrained Hamilton-Jacobi equations in the framework of viscosity solutions. We derive the long-term temporal dynamics of the fittest traits in the regime of small mutations. In the context of adaptive cancer management, we also analyse whether an optimal drug level is better than the maximal tolerated dose.