A Review of Mathematical Models for the Formation of\ud Vascular Networks

Mainly two mechanisms are involved in the formation of blood vasculature: vasculogenesis and angiogenesis. The former consists of the formation of a capillary-like network from either a dispersed or a monolayered population of endothelial cells, reproducible also in vitro by specific experimental assays. The latter consists of the sprouting of new vessels from an existing capillary or post-capillary venule. Similar phenomena are also involved in the formation of the lymphatic system through a process generally called lymphangiogenesis.\ud \ud A number of mathematical approaches have analysed these phenomena. This paper reviews the different modelling procedures, with a special emphasis on their ability to reproduce the biological system and to predict measured quantities which describe the overall processes. A comparison between the different methods is also made, highlighting their specific features

Computer simulation of glioma growth and morphology

Despite major advances in the study of glioma, the quantitative links between intra-tumor molecular/cellular properties, clinically observable properties such as morphology, and critical tumor behaviors such as growth and invasiveness remain unclear, hampering more effective coupling of tumor physical characteristics with implications for prognosis and therapy. Although molecular biology, histopathology, and radiological imaging are employed in this endeavor, studies are severely challenged by the multitude of different physical scales involved in tumor growth, i.e., from molecular nanoscale to cell microscale and finally to tissue centimeter scale. Consequently, it is often difficult to determine the underlying dynamics across dimensions. New techniques are needed to tackle these issues. Here, we address this multi-scalar problem by employing a novel predictive three-dimensional mathematical and computational model based on first-principle equations (conservation laws of physics) that describe mathematically the diffusion of cell substrates and other processes determining tumor mass growth and invasion. The model uses conserved variables to represent known determinants of glioma behavior, e.g., cell density and oxygen concentration, as well as biological functional relationships and parameters linking phenomena at different scales whose specific forms and values are hypothesized and calculated based on in vitro and in vivo experiments and from histopathology of tissue specimens from human gliomas. This model enables correlation of glioma morphology to tumor growth by quantifying interdependence of tumor mass on the microenvironment (e.g., hypoxia, tissue disruption) and on the cellular phenotypes (e.g., mitosis and apoptosis rates, cell adhesion strength). Once functional relationships between variables and associated parameter values have been informed, e.g., from histopathology or intra-operative analysis, this model can be used for disease diagnosis/prognosis, hypothesis testing, and to guide surgery and therapy. In particular, this tool identifies and quantifies the effects of vascularization and other cell-scale glioma morphological characteristics as predictors of tumor-scale growth and invasion

Modeling Stem/Progenitor Cell-Induced Neovascularization and\ud Oxygenation around Solid Implants

Tissue engineering constructs and other solid implants with biomedical applications, such as drug delivery devices or bioartificial organs, need oxygen (O2) to function properly. To understand better the vascular integration of such devices, we recently developed a novel model sensor containing O2-sensitive crystals, consisting of a polymeric capsule limited by a nano-porous filter. The sensor was implanted in mice with hydrogel alone (control) or hydrogel embedded with mouse CD117/c-kit+ bone marrow progenitor cells (BMPC) in order to stimulate peri-implant neovascularization. The sensor provided local partial O2 pressure (pO2) using non-invasive electron paramagnetic resonance (EPR) signal measurements. A consistently higher level of per-implant oxygenation was observed in the cell-treatment case as compared to the control over a 10-week period. In order to provide a mechanistic explanation of these experimental observations, we present in this paper a mathematical model, formulated as a system of coupled partial differential equations, that simulates peri-implant vascularization. In the control case, vascularization is considered to be the result of a Foreign Body Reaction (FBR) while in the cell-treatment case, adipogenesis in response to paracrine stimuli produced by the stem cells is assumed to induce neovascularization. The model is validated by fitting numerical predictions of local pO2 to measurements from the implanted sensor. The model is then used to investigate further the potential for using stem cell treatment to enhance the vascular integration of biomedical implants. We thus demonstrate how mathematical modeling combined with experimentation can be used to infer how vasculature develops around biomedical implants in control and stem celltreated cases

Multiphase modelling of vascular tumour growth in two spatial dimensions

In this paper we present a continuum mathematical model of vascular tumour growth which is based on a multiphase framework in which the tissue is decomposed into four distinct phases and the principles of conservation of mass and momentum are applied to the normal/healthy cells, tumour cells, blood vessels and extracellular material. The inclusion of a diffusible nutrient, supplied by the blood vessels, allows the vasculature to have a nonlocal influence on the other phases. Two-dimensional computational simulations are carried out on unstructured, triangular meshes to allow a natural treatment of irregular geometries, and the tumour boundary is captured as a diffuse interface on this mesh, thereby obviating the need to explicitly track the (potentially highly irregular and ill-defined) tumour boundary. A hybrid finite volume/finite element algorithm is used to discretise the continuum model: the application of a conservative, upwind, finite volume scheme to the hyperbolic mass balance equations and a finite element scheme with a stable element pair to the generalised Stokes equations derived from momentum balance, leads to a robust algorithm which does not use any form of artificial stabilisation. The use of a matrix-free Newton iteration with a finite element scheme for the nutrient reaction-diffusion equations allows full nonlinearity in the source terms of the mathematical model. Numerical simulations reveal that this four-phase model reproduces the characteristic pattern of tumour growth in which a necrotic core forms behind an expanding rim of well-vascularised proliferating tumour cells. The simulations consistently predict linear tumour growth rates. The dependence of both the speed with which the tumour grows and the irregularity of the invading tumour front on the model parameters are investigated

A review of mathematical models for the formation of vascular networks

Two major mechanisms are involved in the formation of blood vasculature: vasculogenesis and angiogenesis. The former term describes the formation of a capillary-like network from either a dispersed or a monolayered population of endothelial cells, reproducible also in vitro by specific experimental assays. The latter term describes the sprouting of new vessels from an existing capillary or post-capillary venule. Similar mechanisms are also involved in the formation of the lymphatic system through a process generally called lymphangiogenesis. A number of mathematical approaches have been used to analyse these phenomena. In this article, we review the different types of models, with special emphasis on their ability to reproduce different biological systems and to predict measurable quantities which describe the overall processes. Finally, we highlight the advantages specific to each of the different modelling approaches. The research that led to the present paper was partially supported by a grant of the group GNFM of INdA

Multiscale modeling of palisade formation in gliobastoma multiforme

Palisades are characteristic tissue aberrations that arise in glioblastomas. Observation of palisades is considered as a clinical indicator of the transition from a noninvasive to an invasive tumour. In this article we propose a computational model to study the influence of genotypic and phenotypic heterogeneity in palisade formation. For this we produced three dimensional realistic simulations, based on a multiscale hybrid model, coupling the evolution of tumour cells and the oxygen diffusion in tissue, that depict the shape of palisades during its formation. Our results can be summarized as the following: (1) we show that cell heterogeneity is a crucial factor in palisade formation and tumour growth; (2) we present results that can explain the observed fact that recursive tumours are more malignant than primary tumours; and (3) the presented simulations can provide to clinicians and biologists for a better understanding of palisades 3D structure as well as glioblastomas growth dynamics

The role of malignant tissue on the thermal distribution of cancerous breast

The present work focuses on the integration of analytical and numerical strategies to investigate the thermal distribution of cancerous breasts. Coupled stationary bioheat transfer equations are considered for the glandular and heterogeneous tumor regions, which are characterized by different thermophysical properties. The cross-section of the cancerous breast is identified by a homogeneous glandular tissue that surrounds the heterogeneous tumor tissue, which is assumed to be a two-phase periodic composite with non-overlapping circular inclusions and a square lattice distribution, wherein the constituents exhibit isotropic thermal conductivity behavior. Asymptotic periodic homogenization method is used to find the effective properties in the heterogeneous region. The tissue effective thermal conductivities are computed analytically and then used in the homogenized model, which is solved numerically. Results are compared with appropriate experimental data reported in the literature. In particular, the tissue scale temperature profile agrees with experimental observations. Moreover, as a novelty result we find that the tumor volume fraction in the heterogeneous zone influences the breast surface temperature

Multiscale modeling of palisade formation in gliobastoma multiforme

Palisades are characteristic tissue aberrations that arise in glioblastomas. Observation of palisades is considered as a clinical indicator of the transition from a noninvasive to an invasive tumour. In this article we propose a computational model to study the influence of genotypic and phenotypic heterogeneity in palisade formation. For this we produced three dimensional realistic simulations, based on a multiscale hybrid model, coupling the evolution of tumour cells and the oxygen diffusion in tissue, that depict the shape of palisades during its formation. Our results can be summarized as the following: (1) we show that cell heterogeneity is a crucial factor in palisade formation and tumour growth; (2) we present results that can explain the observed fact that recursive tumours are more malignant than primary tumours; and (3) the presented simulations can provide to clinicians and biologists for a better understanding of palisades 3D structure as well as glioblastomas growth dynamic