Modeling angiogenesis: A discrete to continuum description

Computational Insights into the Interplay of Mechanical Forces in Angiogenesis

This study employs a meshless computational model to investigate the impacts of compression and traction on angiogenesis, exploring their effects on vascular endothelial growth factor (VEGF) diffusion and subsequent capillary network formation. Three distinct initial domain geometries were defined to simulate variations in endothelial cell sprouting and VEGF release. Compression and traction were applied, and the ensuing effects on VEGF diffusion coefficients were analysed. Compression promoted angiogenesis, increasing capillary network density. The reduction in the VEGF diffusion coefficient under compression altered VEGF concentration, impacting endothelial cell migration patterns. The findings were consistent across diverse simulation scenarios, demonstrating the robust influence of compression on angiogenesis. This computational study enhances our understanding of the intricate interplay between mechanical forces and angiogenesis. Compression emerges as an effective mediator of angiogenesis, influencing VEGF diffusion and vascular pattern. These insights may contribute to innovative therapeutic strategies for angiogenesis-related disorders, fostering tissue regeneration and addressing diseases where angiogenesis is crucial.

Discrete and continuous mathematical models of sharp-fronted collective cell migration and invasion

Mathematical models describing the spatial spreading and invasion of populations of biological cells are often developed in a continuum modelling framework using reaction-diffusion equations. While continuum models based on linear diffusion are routinely employed and known to capture key experimental observations, linear diffusion fails to predict well-defined sharp fronts that are often observed experimentally. This observation has motivated the use of nonlinear degenerate diffusion; however, these nonlinear models and the associated parameters lack a clear biological motivation and interpretation. Here, we take a different approach by developing a stochastic discrete lattice-based model incorporating biologically inspired mechanisms and then deriving the reaction-diffusion continuum limit. Inspired by experimental observations, agents in the simulation deposit extracellular material, which we call a substrate, locally onto the lattice, and the motility of agents is taken to be proportional to the substrate density. Discrete simulations that mimic a two-dimensional circular barrier assay illustrate how the discrete model supports both smooth and sharp-fronted density profiles depending on the rate of substrate deposition. Coarse-graining the discrete model leads to a novel partial differential equation (PDE) model whose solution accurately approximates averaged data from the discrete model. The new discrete model and PDE approximation provide a simple, biologically motivated framework for modelling the spreading, growth and invasion of cell populations with well-defined sharp fronts. Open-source Julia code to replicate all results in this work is available on GitHub.

Computational modeling of angiogenesis: The importance of cell rearrangements during vascular growth

Angiogenesis is the process wherein endothelial cells (ECs) form sprouts that elongate from the pre-existing vasculature to create new vascular networks. In addition to its essential role in normal development, angiogenesis plays a vital role in pathologies such as cancer, diabetes and atherosclerosis. Mathematical and computational modeling has contributed to unraveling its complexity. Many existing theoretical models of angiogenic sprouting are based on the "snail-trail" hypothesis. This framework assumes that leading ECs positioned at sprout tips migrate toward low-oxygen regions while other ECs in the sprout passively follow the leaders' trails and proliferate to maintain sprout integrity. However, experimental results indicate that, contrary to the snail-trail assumption, ECs exchange positions within developing vessels, and the elongation of sprouts is primarily driven by directed migration of ECs. The functional role of cell rearrangements remains unclear. This review of the theoretical modeling of angiogenesis is the first to focus on the phenomenon of cell mixing during early sprouting. We start by describing the biological processes that occur during early angiogenesis, such as phenotype specification, cell rearrangements and cell interactions with the microenvironment. Next, we provide an overview of various theoretical approaches that have been employed to model angiogenesis, with particular emphasis on recent in silico models that account for the phenomenon of cell mixing. Finally, we discuss when cell mixing should be incorporated into theoretical models and what essential modeling components such models should include in order to investigate its functional role. This article is categorized under: Cardiovascular Diseases > Computational Models Cancer > Computational Models.

A mechanical modeling framework to study endothelial permeability

The inner lining of blood vessels, the endothelium, is made up of endothelial cells. Vascular endothelial (VE)-cadherin protein forms a bond with VE-cadherin from neighboring cells to determine the size of gaps between the cells and thereby regulate the size of particles that can cross the endothelium. Chemical cues such as thrombin, along with mechanical properties of the cell and extracellular matrix are known to affect the permeability of endothelial cells. Abnormal permeability is found in patients suffering from diseases including cardiovascular diseases, cancer, and COVID-19. Even though some of the regulatory mechanisms affecting endothelial permeability are well studied, details of how several mechanical and chemical stimuli acting simultaneously affect endothelial permeability are not yet understood. In this article, we present a continuum-level mechanical modeling framework to study the highly dynamic nature of the VE-cadherin bonds. Taking inspiration from the catch-slip behavior that VE-cadherin complexes are known to exhibit, we model the VE-cadherin homophilic bond as cohesive contact with damage following a traction-separation law. We explicitly model the actin cytoskeleton and substrate to study their role in permeability. Our studies show that mechanochemical coupling is necessary to simulate the influence of the mechanical properties of the substrate on permeability. Simulations show that shear between cells is responsible for the variation in permeability between bicellular and tricellular junctions, explaining the phenotypic differences observed in experiments. An increase in the magnitude of traction force due to disturbed flow that endothelial cells experience results in increased permeability, and it is found that the effect is higher on stiffer extracellular matrix. Finally, we show that the cylindrical monolayer exhibits higher permeability than the planar monolayer under unconstrained cases. Thus, we present a contact mechanics-based mechanochemical model to investigate the variation in the permeability of endothelial monolayer due to multiple loads acting simultaneously.

Soliton approximation in continuum models of leader-follower behavior

Complex biological processes involve collective behavior of entities (bacteria, cells, animals) over many length and time scales and can be described by discrete models that track individuals or by continuum models involving densities and fields. We consider hybrid stochastic agent-based models of branching morphogenesis and angiogenesis (new blood vessel creation from preexisting vasculature), which treat cells as individuals that are guided by underlying continuous chemical and/or mechanical fields. In these descriptions, leader (tip) cells emerge from existing branches and follower (stalk) cells build the new sprout in their wake. Vessel branching and fusion (anastomosis) occur as a result of tip and stalk cell dynamics. Coarse graining these hybrid models in appropriate limits produces continuum partial differential equations (PDEs) for endothelial cell densities that are more analytically tractable. While these models differ in nonlinearity, they produce similar equations at leading order when chemotaxis is dominant. We analyze this leading order system in a simple quasi-one-dimensional geometry and show that the numerical solution of the leading order PDE is well described by a soliton wave that evolves from vessel to source. This wave is an attractor for intermediate times until it arrives at the hypoxic region releasing the growth factor. The mathematical techniques used here thus identify common features of discrete and continuum approaches and provide insight into general biological mechanisms governing their collective dynamics.

Single-cell characterization of deformation and dynamics of mesenchymal stem cells in microfluidic systems: A computational study

Understanding the homing dynamics of individual mesenchymal stem cells (MSCs) in physiologically relevant microenvironments is crucial for improving the efficacy of MSC-based therapies for therapeutic and targeting purposes. This study investigates the passive homing behavior of individual MSCs in micropores that mimic interendothelial clefts through predictive computational simulations informed by previous microfluidic experiments. Initially, we quantified the size-dependent behavior of MSCs in micropores and elucidated the underlying mechanisms. Subsequently, we analyzed the shape deformation and traversal dynamics of each MSC. In addition, we conducted a systematic investigation to understand how the mechanical properties of MSCs impact their traversal process. We considered geometric and mechanical parameters, such as reduced cell volume, cell-to-nucleus diameter ratio, and cytoskeletal prestress states. Furthermore, we quantified the changes in the MSC traversal process and identified the quantitative limits in their response to variations in micropore length. Taken together, the computational results indicate the complex dynamic behavior of individual MSCs in the confined microflow. This finding offers an objective way to evaluate the homing ability of MSCs in an interendothelial-slit-like microenvironment.

Approaches to vascular network, blood flow, and metabolite distribution modeling in brain tissue

The cardiovascular system plays a key role in the transport of nutrients, ensuring a continuous supply of all cells of the body with the metabolites necessary for life. The blood supply to the brain is carried out by the large arteries located on its surface, which branch into smaller arterioles that penetrate the cerebral cortex and feed the capillary bed, thereby forming an extensive branching network. The formation of blood vessels is carried out via vasculogenesis and angiogenesis, which play an important role in both embryo and adult life. The review presents approaches to modeling various aspects of both the formation of vascular networks and the construction of the formed arterial tree. In addition, a brief description of models that allows one to study the blood flow in various parts of the circulatory system and the spatiotemporal metabolite distribution in brain tissues is given. Experimental study of these issues is not always possible due to both the complexity of the cardiovascular system and the mechanisms through which the perfusion of all body cells is carried out. In this regard, mathematical models are a good tool for studying hemodynamics and can be used in clinical practice to diagnose vascular diseases and assess the need for treatment.

Development of a coupled modeling for tumor growth, angiogenesis, oxygen delivery, and phenotypic heterogeneity

Analysis of the evolution and growth dynamics of tumors is crucial for understanding cancer and the development of individually optimized therapies. During tumor growth, a hypoxic microenvironment around cancer cells caused by excessive non-vascular tumor growth induces tumor angiogenesis that plays a key role in the ensuing tumor growth and its progression into higher stages. Various mathematical simulation models have been introduced to simulate these biologically and physically complex hallmarks of cancer. Here, we developed a hybrid two-dimensional computational model that integrates spatiotemporally different components of the tumor system to investigate both angiogenesis and tumor growth/proliferation. This spatiotemporal evolution is based on partial diffusion equations, the cellular automation method, transition and probabilistic rules, and biological assumptions. The new vascular network provided by angiogenesis affects tumor microenvironmental conditions and drives individual cells to adapt themselves to spatiotemporal conditions. Furthermore, some stochastic rules are involved besides microenvironmental conditions. Overall, the conditions promote some commonly observed cellular states, i.e., proliferative, migrative, quiescent, and cell death, depending on the condition of each cell. Altogether, our results offer a theoretical basis for the biological evidence that regions of the tumor tissue near blood vessels are densely populated by proliferative phenotypic variants, while poorly oxygenated regions are sparsely populated by hypoxic phenotypic variants.

Data-driven spatio-temporal modelling of glioblastoma

Mathematical oncology provides unique and invaluable insights into tumour growth on both the microscopic and macroscopic levels. This review presents state-of-the-art modelling techniques and focuses on their role in understanding glioblastoma, a malignant form of brain cancer. For each approach, we summarize the scope, drawbacks and assets. We highlight the potential clinical applications of each modelling technique and discuss the connections between the mathematical models and the molecular and imaging data used to inform them. By doing so, we aim to prime cancer researchers with current and emerging computational tools for understanding tumour progression. By providing an in-depth picture of the different modelling techniques, we also aim to assist researchers who seek to build and develop their own models and the associated inference frameworks. Our article thus strikes a unique balance. On the one hand, we provide a comprehensive overview of the available modelling techniques and their applications, including key mathematical expressions. On the other hand, the content is accessible to mathematicians and biomedical scientists alike to accommodate the interdisciplinary nature of cancer research.

Computational models for generating microvascular structures: Investigations beyond medical imaging resolution

Angiogenesis, arteriogenesis, and pruning are revascularization processes essential to our natural vascular development and adaptation, as well as central players in the onset and development of pathologies such as tumoral growth and stroke recovery. Computational modeling allows for repeatable experimentation and exploration of these complex biological processes. In this review, we provide an introduction to the biological understanding of the vascular adaptation processes of sprouting angiogenesis, intussusceptive angiogenesis, anastomosis, pruning, and arteriogenesis, discussing some of the more significant contributions made to the computational modeling of these processes. Each computational model represents a theoretical framework for how biology functions, and with rises in computing power and study of the problem these frameworks become more accurate and complete. We highlight physiological, pathological, and technological applications that can be benefit from the advances performed by these models, and we also identify which elements of the biology are underexplored in the current state-of-the-art computational models. This article is categorized under: Cancer > Computational Models Cardiovascular Diseases > Computational Models.

A Model for Cell Proliferation in a Developing Organism

In mathematical biology, there is a great deal of interest in producing continuum models by scaling discrete agent-based models governed by local stochastic rules. We discuss a particular example of this approach: a model for the proliferation of neural crest cells that can help us understand the development of Hirschprung's disease, a potentially-fatal condition in which the enteric nervous system of a new-born child does not extend all the way through the intestine and colon. Our starting point is a discrete-state, continuous-time Markov chain model proposed by Hywood et al. (2013a) for the location of the neural crest cells that make up the enteric nervous system. Hywood et al. (2013a) scaled their model to derive an approximate second order partial differential equation describing how the limiting expected number of neural crest cells evolve in space and time. In contrast, we exploit the relationship between the above-mentioned Markov chain model and the well-known Yule-Furry process to derive the exact form of the scaled version of the process. Furthermore, we provide expressions for other features of the domain agent occupancy process, such as the variance of the marginal occupancy at a particular site, the distribution of the number of agents that are yet to reach a given site and a stochastic description of the process itself.

Mathematical models of neuronal growth

The establishment of a functioning neuronal network is a crucial step in neural development. During this process, neurons extend neurites-axons and dendrites-to meet other neurons and interconnect. Therefore, these neurites need to migrate, grow, branch and find the correct path to their target by processing sensory cues from their environment. These processes rely on many coupled biophysical effects including elasticity, viscosity, growth, active forces, chemical signaling, adhesion and cellular transport. Mathematical models offer a direct way to test hypotheses and understand the underlying mechanisms responsible for neuron development. Here, we critically review the main models of neurite growth and morphogenesis from a mathematical viewpoint. We present different models for growth, guidance and morphogenesis, with a particular emphasis on mechanics and mechanisms, and on simple mathematical models that can be partially treated analytically.

Comparative analysis of continuum angiogenesis models

Although discrete approaches are increasingly employed to model biological phenomena, it remains unclear how complex, population-level behaviours in such frameworks arise from the rules used to represent interactions between individuals. Discrete-to-continuum approaches, which are used to derive systems of coarse-grained equations describing the mean-field dynamics of a microscopic model, can provide insight into such emergent behaviour. Coarse-grained models often contain nonlinear terms that depend on the microscopic rules of the discrete framework, however, and such nonlinearities can make a model difficult to mathematically analyse. By contrast, models developed using phenomenological approaches are typically easier to investigate but have a more obscure connection to the underlying microscopic system. To our knowledge, there has been little work done to compare solutions of phenomenological and coarse-grained models. Here we address this problem in the context of angiogenesis (the creation of new blood vessels from existing vasculature). We compare asymptotic solutions of a classical, phenomenological “snail-trail” model for angiogenesis to solutions of a nonlinear system of partial differential equations (PDEs) derived via a systematic coarse-graining procedure (Pillay et al. in Phys Rev E 95(1):012410, 2017. https://doi.org/10.1103/PhysRevE.95.012410). For distinguished parameter regimes corresponding to chemotaxis-dominated cell movement and low branching rates, both continuum models reduce at leading order to identical PDEs within the domain interior. Numerical and analytical results confirm that pointwise differences between solutions to the two continuum models are small if these conditions hold, and demonstrate how perturbation methods can be used to determine when a phenomenological model provides a good approximation to a more detailed coarse-grained system for the same biological process.

A multiscale model of complex endothelial cell dynamics in early angiogenesis

We introduce a hybrid two-dimensional multiscale model of angiogenesis, the process by which endothelial cells (ECs) migrate from a pre-existing vascular bed in response to local environmental cues and cell-cell interactions, to create a new vascular network. Recent experimental studies have highlighted a central role of cell rearrangements in the formation of angiogenic networks. Our model accounts for this phenomenon via the heterogeneous response of ECs to their microenvironment. These cell rearrangements, in turn, dynamically remodel the local environment. The model reproduces characteristic features of angiogenic sprouting that include branching, chemotactic sensitivity, the brush border effect, and cell mixing. These properties, rather than being hardwired into the model, emerge naturally from the gene expression patterns of individual cells. After calibrating and validating our model against experimental data, we use it to predict how the structure of the vascular network changes as the baseline gene expression levels of the VEGF-Delta-Notch pathway, and the composition of the extracellular environment, vary. In order to investigate the impact of cell rearrangements on the vascular network structure, we introduce the mixing measure, a scalar metric that quantifies cell mixing as the vascular network grows. We calculate the mixing measure for the simulated vascular networks generated by ECs of different lineages (wild type cells and mutant cells with impaired expression of a specific receptor). Our results show that the time evolution of the mixing measure is directly correlated to the generic features of the vascular branching pattern, thus, supporting the hypothesis that cell rearrangements play an essential role in sprouting angiogenesis. Furthermore, we predict that lower cell rearrangement leads to an imbalance between branching and sprout elongation. Since the computation of this statistic requires only individual cell trajectories, it can be computed for networks generated in biological experiments, making it a potential biomarker for pathological angiogenesis.

Angiogenesis, the process by which new blood vessels are formed by sprouting from the pre-existing vascular bed, plays a key role in both physiological and pathological processes, including tumour growth. The structure of a growing vascular network is determined by the coordinated behaviour of endothelial cells in response to various signalling cues. Recent experimental studies have highlighted the importance of cell rearrangements as a driver for sprout elongation. However, the functional role of this phenomenon remains unclear. We formulate a new multiscale model of angiogenesis which, by accounting explicitly for the complex dynamics of endothelial cells within growing angiogenic sprouts, is able to reproduce generic features of angiogenic structures (branching, chemotactic sensitivity, cell mixing, etc.) as emergent properties of its dynamics. We validate our model against experimental data and then use it to quantify the phenomenon of cell mixing in vascular networks generated by endothelial cells of different lineages. Our results show that there is a direct correlation between the time evolution of cell mixing in a growing vascular network and its branching structure, thus paving the way for understanding the functional role of cell rearrangements in angiogenesis.

Evaluating snail-trail frameworks for leader-follower behavior with agent-based modeling

Branched networks constitute a ubiquitous structure in biology, arising in plants, lungs, and the circulatory system; however, the mechanisms behind their creation are not well understood. A commonly used model for network morphogenesis proposes that sprouts develop through interactions between leader (tip) cells and follower (stalk) cells. In this description, tip cells emerge from existing structures, travel up chemoattractant gradients, and form new networks by guiding the movement of stalk cells. Such dynamics have been mathematically represented by continuum "snail-trail" models in which the tip cell flux contributes to the stalk cell proliferation rate. Although snail-trail models constitute a classical depiction of leader-follower behavior, their accuracy has yet to be evaluated in a rigorous quantitative setting. Here, we extend the snail-trail modeling framework to two spatial dimensions by introducing a novel multiplicative factor to the stalk cell rate equation, which corrects for neglected network creation in directions other than that of the migrating front. Our derivation of this factor demonstrates that snail-trail models are valid descriptions of cell dynamics when chemotaxis dominates cell movement. We confirm that our snail-trail model accurately predicts the dynamics of tip and stalk cells in an existing agent-based model (ABM) for network formation [Pillay et al., Phys. Rev. E 95, 012410 (2017)10.1103/PhysRevE.95.012410]. We also derive conditions for which it is appropriate to use a reduced, one-dimensional snail-trail model to analyze ABM results. Our analysis identifies key metrics for cell migration that may be used to anticipate when simple snail-trail models will accurately describe experimentally observed cell dynamics in network formation.

Therapeutic Potential of Brassinosteroids in Biomedical and Clinical Research

Steroids are a pivotal class of hormones with a key role in growth modulation and signal transduction in multicellular organisms. Synthetic steroids are widely used to cure large array of viral, fungal, bacterial, and cancerous infections. Brassinosteroids (BRs) are a natural collection of phytosterols, which have structural similarity with animal steroids. BRs are dispersed universally throughout the plant kingdom. These plant steroids are well known to modulate a plethora of physiological responses in plants leading to improvement in quality as well as yield of food crops. Moreover, they have been found to play imperative role in stress-fortification against various stresses in plants. Over a decade, BRs have conquered worldwide interest due to their diverse biological activities in animal systems. Recent studies have indicated anticancerous, antiangiogenic, antiviral, antigenotoxic, antifungal, and antibacterial bioactivities of BRs in the animal test systems. BRs inhibit replication of viruses and induce cytotoxic effects on cancerous cell lines. Keeping in view the biological activities of BRs, this review is an attempt to update the information about prospects of BRs in biomedical and clinical application.

Modelling skin wound healing angiogenesis: A review

The occurrence of wounds is a main health concern in Western society due to their high frequency and treatment cost. During wound healing, the formation of a functional blood vessel network through angiogenesis is an essential process. Angiogenesis allows the reestablishment of the normal blood flow, the sufficient exchange of oxygen and nutrients and the removal of metabolic waste, necessary for cell proliferation and viability. Mathematical and computational models provide new tools to improve the healing process. In fact, over the last thirty years, in silico models have been continuously formulated to describe the effect of several biological and mechanical factors in angiogenesis during wound healing. Additionally, with different levels of complexity, these models allow coupling the human skin structure, to distinct cell types and growth factors, to study extracellular matrix composition and to understand its deformation. This paper discusses how in silico models, which are more economical and less time-consuming comparatively to laboratory methodologies, can help test new strategies to promote/optimize angiogenesis. The continuum, cell-based and hybrid mathematical models of wound healing angiogenesis are reviewed in the present paper, in order to identify possible improvements. Accordingly, the development of higher dimension models incorporating multiscale analysis at molecular, cellular and tissue level remains a challenge that future models should consider.

The impact of exclusion processes on angiogenesis models

Angiogenesis is the process by which new blood vessels form from existing vessels. During angiogenesis, tip cells migrate via diffusion and chemotaxis, new tip cells are introduced through branching, loops form via tip-to-tip and tip-to-sprout anastomosis, and a vessel network forms as endothelial cells, known as stalk cells, follow the paths of tip cells (a process known as the snail-trail). Using a mean-field approximation, we systematically derive one-dimensional non-linear continuum models from a lattice-based cellular automaton model of angiogenesis in the corneal assay, explicitly accounting for cell volume. We compare our continuum models and a well-known phenomenological snail-trail model that is linear in the diffusive, chemotactic and branching terms, with averaged cellular automaton simulation results to distinguish macroscale volume exclusion effects and determine whether linear models can capture them. We conclude that, in general, both linear and non-linear models can be used at low cell densities when single or multi-species exclusion effects are negligible at the macroscale. When cell densities increase, our non-linear model should be used to capture non-linear tip cell behavior that occurs when single-species exclusion effects are pronounced, and alternative models should be derived for non-negligible multi-species exclusion effects.