A novel nonlocal partial differential equation model of endothelial progenitor cell cluster formation during the early stages of vasculogenesis

TL;DR

This paper introduces a new nonlocal PDE model for early endothelial progenitor cell cluster formation during vasculogenesis, combining biological mechanisms and analyzing their roles through stability and numerical studies.

Contribution

The study develops a novel mathematical nonlocal PDE model for EPC cluster formation, incorporating key biological mechanisms and providing insights into their roles.

Findings

Chemotaxis and matrix degradation are crucial for cluster formation.

Previous models of cell-to-cell adhesion are insufficient for stability analysis.

Numerical simulations align qualitatively with in vitro experimental data.

Abstract

Neovascularisation is essential for tissue development and regeneration, in addition to playing a key role in pathological settings such as ischemia and tumour development. Experimental findings in the past two decades have led to the identification of a new mechanism of neovascularisation, cluster-based vasculogenesis, during which endothelial progenitor cells (EPCs) mobilised from the bone marrow are capable of bridging distant vascular beds in a variety of hypoxic settings in vivo. This process is characterised by the formation of EPC clusters during its early stages and, while much progress has been made in identifying various mechanisms underlying cluster formation, we are still far from a comprehensive description of such spatio-temporal dynamics. In order to achieve this, we propose a novel mathematical model of the early stages of cluster-based vasculogenesis, comprising of a system of nonlocal partial differential equations including key mechanisms such as endogenous chemotaxis, matrix degradation, cell proliferation and cell-to-cell adhesion. We conduct a linear stability analysis on the system, solve the equations numerically, conduct a parametric analysis of the numerical solutions of the 1D problem to investigate the role of underlying dynamics on the speed of cluster formation and the size of clusters, and verify the key results of the parametric analysis with simulations of the 2D problem. Our results, which qualitatively compare with data from in vitro experiments, elucidate the complementary role played by endogenous chemotaxis and matrix degradation in the formation of clusters, and they indicate that previous approaches to the nonlocal modelling of cell-to-cell adhesion, while they capture the aggregating effect of cell-to-cell adhesion, are not sufficient to capture its stabilising effect on clusters, and new continuum cell-adhesion modelling strategies are required.