A multiscale model of complex endothelial cell dynamics in early angiogenesis

Modeling the extracellular matrix in cell migration and morphogenesis: a guide for the curious biologist

The extracellular matrix (ECM) is a highly complex structure through which biochemical and mechanical signals are transmitted. In processes of cell migration, the ECM also acts as a scaffold, providing structural support to cells as well as points of potential attachment. Although the ECM is a well-studied structure, its role in many biological processes remains difficult to investigate comprehensively due to its complexity and structural variation within an organism. In tandem with experiments, mathematical models are helpful in refining and testing hypotheses, generating predictions, and exploring conditions outside the scope of experiments. Such models can be combined and calibrated with in vivo and in vitro data to identify critical cell-ECM interactions that drive developmental and homeostatic processes, or the progression of diseases. In this review, we focus on mathematical and computational models of the ECM in processes such as cell migration including cancer metastasis, and in tissue structure and morphogenesis. By highlighting the predictive power of these models, we aim to help bridge the gap between experimental and computational approaches to studying the ECM and to provide guidance on selecting an appropriate model framework to complement corresponding experimental studies.

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.

Endothelial cells signaling and patterning under hypoxia: a mechanistic integrative computational model including the Notch-Dll4 pathway

Introduction: Several signaling pathways are activated during hypoxia to promote angiogenesis, leading to endothelial cell patterning, interaction, and downstream signaling. Understanding the mechanistic signaling differences between endothelial cells under normoxia and hypoxia and their response to different stimuli can guide therapies to modulate angiogenesis. We present a novel mechanistic model of interacting endothelial cells, including the main pathways involved in angiogenesis. Methods: We calibrate and fit the model parameters based on well-established modeling techniques that include structural and practical parameter identifiability, uncertainty quantification, and global sensitivity. Results: Our results indicate that the main pathways involved in patterning tip and stalk endothelial cells under hypoxia differ, and the time under hypoxia interferes with how different stimuli affect patterning. Additionally, our simulations indicate that Notch signaling might regulate vascular permeability and establish different Nitric Oxide release patterns for tip/stalk cells. Following simulations with various stimuli, our model suggests that factors such as time under hypoxia and oxygen availability must be considered for EC pattern control. Discussion: This project provides insights into the signaling and patterning of endothelial cells under various oxygen levels and stimulation by VEGFA and is our first integrative approach toward achieving EC control as a method for improving angiogenesis. Overall, our model provides a computational framework that can be built on to test angiogenesis-related therapies by modulation of different pathways, such as the Notch pathway.

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.

Enhanced perfusion following exposure to radiotherapy: A theoretical investigation

Tumour angiogenesis leads to the formation of blood vessels that are structurally and spatially heterogeneous. Poor blood perfusion, in conjunction with increased hypoxia and oxygen heterogeneity, impairs a tumour's response to radiotherapy. The optimal strategy for enhancing tumour perfusion remains unclear, preventing its regular deployment in combination therapies. In this work, we first identify vascular architectural features that correlate with enhanced perfusion following radiotherapy, using in vivo imaging data from vascular tumours. Then, we present a novel computational model to determine the relationship between these architectural features and blood perfusion in silico. If perfusion is defined to be the proportion of vessels that support blood flow, we find that vascular networks with small mean diameters and large numbers of angiogenic sprouts show the largest increases in perfusion post-irradiation for both biological and synthetic tumours. We also identify cases where perfusion increases due to the pruning of hypoperfused vessels, rather than blood being rerouted. These results indicate the importance of considering network composition when determining the optimal irradiation strategy. In the future, we aim to use our findings to identify tumours that are good candidates for perfusion enhancement and to improve the efficacy of combination therapies.

Optimal mechanical interactions direct multicellular network formation on elastic substrates

Cells self-organize into functional, ordered structures during tissue morphogenesis, a process that is evocative of colloidal self-assembly into engineered soft materials. Understanding how intercellular mechanical interactions may drive the formation of ordered and functional multicellular structures is important in developmental biology and tissue engineering. Here, by combining an agent-based model for contractile cells on elastic substrates with endothelial cell culture experiments, we show that substrate deformation-mediated mechanical interactions between cells can cluster and align them into branched networks. Motivated by the structure and function of vasculogenic networks, we predict how measures of network connectivity like percolation probability and fractal dimension as well as local morphological features including junctions, branches, and rings depend on cell contractility and density and on substrate elastic properties including stiffness and compressibility. We predict and confirm with experiments that cell network formation is substrate stiffness dependent, being optimal at intermediate stiffness. We also show the agreement between experimental data and predicted cell cluster types by mapping a combined phase diagram in cell density substrate stiffness. Overall, we show that long-range, mechanical interactions provide an optimal and general strategy for multicellular self-organization, leading to more robust and efficient realizations of space-spanning networks than through just local intercellular interactions.

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.

A multiscale model of the role of microenvironmental factors in cell segregation and heterogeneity in breast cancer development

We analyzed a quantitative multiscale model that describes the epigenetic dynamics during the growth and evolution of an avascular tumor. A gene regulatory network (GRN) formed by a set of ten genes that are believed to play an important role in breast cancer development was kinetically coupled to the microenvironmental agents: glucose, estrogens, and oxygen. The dynamics of spontaneous mutations was described by a Yule-Furry master equation whose solution represents the probability that a given cell in the tissue undergoes a certain number of mutations at a given time. We assumed that the mutation rate is modified by a spatial gradient of nutrients. The tumor mass was simulated by means of cellular automata supplemented with a set of reaction diffusion equations that described the transport of microenvironmental agents. By analyzing the epigenetic state space described by the GRN dynamics, we found three attractors that were identified with cellular epigenetic states: normal, precancer and cancer. For two-dimensional (2D) and three-dimensional (3D) tumors we calculated the spatial distribution of the following quantities: (i) number of mutations, (ii) mutation of each gene and, (iii) phenotypes. Using estrogen as the principal microenvironmental agent that regulates cell proliferation process, we obtained tumor shapes for different values of estrogen consumption and supply rates. It was found that he majority of mutations occurred in cells that were located close to the 2D tumor perimeter or close to the 3D tumor surface. Also, it was found that the occurrence of different phenotypes in the tumor are controlled by estrogen concentration levels since they can change the individual cell threshold and gene expression levels. All results were consistently observed for 2D and 3D tumors.

Self-organized and directed branching results in optimal coverage in developing dermal lymphatic networks

Branching morphogenesis is a ubiquitous process that gives rise to high exchange surfaces in the vasculature and epithelial organs. Lymphatic capillaries form branched networks, which play a key role in the circulation of tissue fluid and immune cells. Although mouse models and correlative patient data indicate that the lymphatic capillary density directly correlates with functional output, i.e., tissue fluid drainage and trafficking efficiency of dendritic cells, the mechanisms ensuring efficient tissue coverage remain poorly understood. Here, we use the mouse ear pinna lymphatic vessel network as a model system and combine lineage-tracing, genetic perturbations, whole-organ reconstructions and theoretical modeling to show that the dermal lymphatic capillaries tile space in an optimal, space-filling manner. This coverage is achieved by two complementary mechanisms: initial tissue invasion provides a non-optimal global scaffold via self-organized branching morphogenesis, while VEGF-C dependent side-branching from existing capillaries rapidly optimizes local coverage by directionally targeting low-density regions. With these two ingredients, we show that a minimal biophysical model can reproduce quantitatively whole-network reconstructions, across development and perturbations. Our results show that lymphatic capillary networks can exploit local self-organizing mechanisms to achieve tissue-scale optimization.

Combined computational modeling and experimental study of the biomechanical mechanisms of platelet-driven contraction of fibrin clots

While blood clot formation has been relatively well studied, little is known about the mechanisms underlying the subsequent structural and mechanical clot remodeling called contraction or retraction. Impairment of the clot contraction process is associated with both life-threatening bleeding and thrombotic conditions, such as ischemic stroke, venous thromboembolism, and others. Recently, blood clot contraction was observed to be hindered in patients with COVID-19. A three-dimensional multiscale computational model is developed and used to quantify biomechanical mechanisms of the kinetics of clot contraction driven by platelet-fibrin pulling interactions. These results provide important biological insights into contraction of platelet filopodia, the mechanically active thin protrusions of the plasma membrane, described previously as performing mostly a sensory function. The biomechanical mechanisms and modeling approach described can potentially apply to studying other systems in which cells are embedded in a filamentous network and exert forces on the extracellular matrix modulated by the substrate stiffness.

Analyzing the effect of cell rearrangement on Delta-Notch pattern formation

The Delta-Notch system plays a vital role in many areas of biology and typically forms a salt and pepper pattern in which cells strongly expressing Delta and cells strongly expressing Notch are alternately aligned via lateral inhibition. In this study, we consider cell rearrangement events, such as cell mixing and proliferation, that alter the spatial structure itself and affect the pattern dynamics. We model cell rearrangement events by a Poisson process and analyze the model while preserving the discrete properties of the spatial structure. We investigate the effects of the intermittent perturbations arising from these cell rearrangement events on the discrete spatial structure itself in the context of pattern formation and by using an analytical approach, coupled with numerical simulation. We find that the homogeneous expression pattern is stabilized if the frequency of cell rearrangement events is sufficiently large. We analytically obtain the balanced frequencies of the cell rearrangement events where the decrease of the pattern amplitude, as a result of cell rearrangement, is balanced by the increase in amplitude due to the Delta-Notch interaction dynamics. Our framework, while applied here to the specific case of the Delta-Notch system, is applicable more widely to other pattern formation mechanisms.

Endothelial cells signaling and patterning under hypoxia: a mechanistic integrative computational model including the Notch-Dll4 pathway

Several signaling pathways are activated during hypoxia to promote angiogenesis, leading to endothelial cell patterning, interaction, and downstream signaling. Understanding the mechanistic signaling differences between normoxia and hypoxia can guide therapies to modulate angiogenesis. We present a novel mechanistic model of interacting endothelial cells, including the main pathways involved in angiogenesis. We calibrate and fit the model parameters based on well-established modeling techniques. Our results indicate that the main pathways involved in the patterning of tip and stalk endothelial cells under hypoxia differ, and the time under hypoxia affects how a reaction affects patterning. Interestingly, the interaction of receptors with Neuropilin1 is also relevant for cell patterning. Our simulations under different oxygen concentrations indicate time- and oxygen-availability-dependent responses for the two cells. Following simulations with various stimuli, our model suggests that factors such as period under hypoxia and oxygen availability must be considered for pattern control. This project provides insights into the signaling and patterning of endothelial cells under hypoxia, contributing to studies in the field.

A new model of Notch signalling: Control of Notch receptor cis-inhibition via Notch ligand dimers

All tissue development and replenishment relies upon the breaking of symmetries leading to the morphological and operational differentiation of progenitor cells into more specialized cells. One of the main engines driving this process is the Notch signal transduction pathway, a ubiquitous signalling system found in the vast majority of metazoan cell types characterized to date. Broadly speaking, Notch receptor activity is governed by a balance between two processes: 1) intercellular Notch transactivation triggered via interactions between receptors and ligands expressed in neighbouring cells; 2) intracellular cis inhibition caused by ligands binding to receptors within the same cell. Additionally, recent reports have also unveiled evidence of cis activation. Whilst context-dependent Notch receptor clustering has been hypothesized, to date, Notch signalling has been assumed to involve an interplay between receptor and ligand monomers. In this study, we demonstrate biochemically, through a mutational analysis of DLL4, both in vitro and in tissue culture cells, that Notch ligands can efficiently self-associate. We found that the membrane proximal EGF-like repeat of DLL4 was necessary and sufficient to promote oligomerization/dimerization. Mechanistically, our experimental evidence supports the view that DLL4 ligand dimerization is specifically required for cis-inhibition of Notch receptor activity. To further substantiate these findings, we have adapted and extended existing ordinary differential equation-based models of Notch signalling to take account of the ligand dimerization-dependent cis-inhibition reported here. Our new model faithfully recapitulates our experimental data and improves predictions based upon published data. Collectively, our work favours a model in which net output following Notch receptor/ligand binding results from ligand monomer-driven Notch receptor transactivation (and cis activation) counterposed by ligand dimer-mediated cis-inhibition.

Matrix stiffness regulates Notch signaling activity in endothelial cells

Notch signaling is critical for many developmental and disease-related processes. It is widely accepted that Notch has a mechanotransduction module that regulates receptor cleavage. However, the role of biomechanical properties of the cellular environment in Notch signaling in general is still poorly understood. During angiogenesis, differentiation of endothelial cells into tip and stalk cells is regulated by Notch signaling, and remodeling of the extracellular matrix occurs. We investigated the influence of substrate stiffness on the Notch signaling pathway in endothelial cells. Using stiffness-tuned polydimethylsiloxane (PDMS) substrates, we show that activity of the Notch signaling pathway inversely correlates with a physiologically relevant range of substrate stiffness (i.e. increased Notch signaling activity on softer substrates). Trans-endocytosis of the Notch extracellular domain, but not the overall endocytosis, is regulated by substrate stiffness, and integrin cell-matrix connections are both stiffness dependent and influenced by Notch signaling. We conclude that mechanotransduction of Notch activation is modulated by substrate stiffness, highlighting the role of substrate rigidity as an important cue for signaling. This might have implications in pathological situations associated with stiffening of the extracellular matrix, such as tumor growth.

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.

Towards integration of time-resolved confocal microscopy of a 3D in vitro microfluidic platform with a hybrid multiscale model of tumor angiogenesis

The goal of this study is to calibrate a multiscale model of tumor angiogenesis with time-resolved data to allow for systematic testing of mathematical predictions of vascular sprouting. The multi-scale model consists of an agent-based description of tumor and endothelial cell dynamics coupled to a continuum model of vascular endothelial growth factor concentration. First, we calibrate ordinary differential equation models to time-resolved protein concentration data to estimate the rates of secretion and consumption of vascular endothelial growth factor by endothelial and tumor cells, respectively. These parameters are then input into the multiscale tumor angiogenesis model, and the remaining model parameters are then calibrated to time resolved confocal microscopy images obtained within a 3D vascularized microfluidic platform. The microfluidic platform mimics a functional blood vessel with a surrounding collagen matrix seeded with inflammatory breast cancer cells, which induce tumor angiogenesis. Once the multi-scale model is fully parameterized, we forecast the spatiotemporal distribution of vascular sprouts at future time points and directly compare the predictions to experimentally measured data. We assess the ability of our model to globally recapitulate angiogenic vasculature density, resulting in an average relative calibration error of 17.7% ± 6.3% and an average prediction error of 20.2% ± 4% and 21.7% ± 3.6% using one and four calibrated parameters, respectively. We then assess the model's ability to predict local vessel morphology (individualized vessel structure as opposed to global vascular density), initialized with the first time point and calibrated with two intermediate time points. In this study, we have rigorously calibrated a mechanism-based, multiscale, mathematical model of angiogenic sprouting to multimodal experimental data to make specific, testable predictions.

On the role of mechanical signals on sprouting angiogenesis through computer modeling approaches

Sprouting angiogenesis, the formation of new vessels from preexisting vasculature, is an essential process in the regeneration of new tissues as well as in the development of some diseases like cancer. Although early studies identified chemical signaling as the main driver of this process, many recent studies have shown a strong role of mechanical signals in the formation of new capillaries. Different types of mechanical signals (e.g., external forces, cell traction forces, and blood flow-induced shear forces) have been shown to play distinct roles in the process; however, their interplay remains still largely unknown. During the last decades, mathematical and computational modeling approaches have been developed to investigate and better understand the mechanisms behind mechanically driven angiogenesis. In this manuscript, we review computational models of angiogenesis with a focus on models investigating the role of mechanics on the process. Our aim is not to provide a detailed review on model methodology but to describe what we have learnt from these models. We classify models according to the mechanical signals being investigated and describe how models have looked into their role on the angiogenic process. We show that a better understanding of the mechanobiology of the angiogenic process will require the development of computer models that incorporate the interactions between the multiple mechanical signals and their effect on cellular responses, since they all seem to play a key in sprout patterning. In the end, we describe some of the remaining challenges of computational modeling of angiogenesis and discuss potential avenues for future research.

Computational model of brain endothelial cell signaling pathways predicts therapeutic targets for cerebral pathologies

Brain endothelial cells serve many critical homeostatic functions. In addition to sensing and regulating blood flow, they maintain blood-brain barrier function, including precise control of nutrient exchange and efflux of xenobiotics. Many signaling pathways in brain endothelial cells have been implicated in both health and disease; however, our understanding of how these signaling pathways functionally integrate is limited. A model capable of integrating these signaling pathways could both advance our understanding of brain endothelial cell signaling networks and potentially identify promising molecular targets for endothelial cell-based drug or gene therapies. To this end, we developed a large-scale computational model, wherein brain endothelial cell signaling pathways were reconstructed from the literature and converted into a network of logic-based differential equations. The model integrates 63 nodes (including proteins, mRNA, small molecules, and cell phenotypes) and 82 reactions connecting these nodes. Specifically, our model combines signaling pathways relating to VEGF-A, BDNF, NGF, and Wnt signaling, in addition to incorporating pathways relating to focused ultrasound as a therapeutic delivery tool. To validate the model, independently established relationships between selected inputs and outputs were simulated, with the model yielding correct predictions 73% of the time. We identified influential and sensitive nodes under different physiological or pathological contexts, including altered brain endothelial cell conditions during glioma, Alzheimer's disease, and ischemic stroke. Nodes with the greatest influence over combinations of desired model outputs were identified as potential druggable targets for these disease conditions. For example, the model predicts therapeutic benefits from inhibiting AKT, Hif-1α, or cathepsin D in the context of glioma - each of which are currently being studied in clinical or pre-clinical trials. Notably, the model also permits testing multiple combinations of node alterations for their effects on the network and the desired outputs (such as inhibiting AKT and overexpressing the P75 neurotrophin receptor simultaneously in the context of glioma), allowing for the prediction of optimal combination therapies. In all, our approach integrates results from over 100 past studies into a coherent and powerful model, capable of both revealing network interactions unapparent from studying any one pathway in isolation and predicting therapeutic targets for treating devastating brain pathologies.

A multiscale cell-based model of tumor growth for chemotherapy assessment and tumor-targeted therapy through a 3D computational approach

OBJECTIVES

Computational modeling of biological systems is a powerful tool to clarify diverse processes contributing to cancer. The aim is to clarify the complex biochemical and mechanical interactions between cells, the relevance of intracellular signaling pathways in tumor progression and related events to the cancer treatments, which are largely ignored in previous studies.

MATERIALS AND METHODS

A three-dimensional multiscale cell-based model is developed, covering multiple time and spatial scales, including intracellular, cellular, and extracellular processes. The model generates a realistic representation of the processes involved from an implementation of the signaling transduction network.

RESULTS

Considering a benign tumor development, results are in good agreement with the experimental ones, which identify three different phases in tumor growth. Simulating tumor vascular growth, results predict a highly vascularized tumor morphology in a lobulated form, a consequence of cells' motile behavior. A novel systematic study of chemotherapy intervention, in combination with targeted therapy, is presented to address the capability of the model to evaluate typical clinical protocols. The model also performs a dose comparison study in order to optimize treatment efficacy and surveys the effect of chemotherapy initiation delays and different regimens.

CONCLUSIONS

Results not only provide detailed insights into tumor progression, but also support suggestions for clinical implementation. This is a major step toward the goal of predicting the effects of not only traditional chemotherapy but also tumor-targeted therapies.

High content Image Analysis to study phenotypic heterogeneity in endothelial cell monolayers

Endothelial cells (EC) are heterogeneous across and within tissues, reflecting distinct, specialised functions. EC heterogeneity has been proposed to underpin EC plasticity independently from vessel microenvironments. However, heterogeneity driven by contact-dependent or short-range cell-cell crosstalk cannot be evaluated with single cell transcriptomic approaches as spatial and contextual information is lost. Nonetheless, quantification of EC heterogeneity and understanding of its molecular drivers is key to developing novel therapeutics for cancer, cardiovascular diseases and for revascularisation in regenerative medicine. Here, we developed an EC profiling tool (ECPT) to examine individual cells within intact monolayers. We used ECPT to characterise different phenotypes in arterial, venous and microvascular EC populations. In line with other studies, we measured heterogeneity in terms of cell cycle, proliferation, and junction organisation. ECPT uncovered a previously under-appreciated single-cell heterogeneity in NOTCH activation. We correlated cell proliferation with different NOTCH activation states at the single cell and population levels. The positional and relational information extracted with our novel approach is key to elucidating the molecular mechanisms underpinning EC heterogeneity.

Systems biology of angiogenesis signaling: Computational models and omics

Angiogenesis is a highly regulated multiscale process that involves a plethora of cells, their cellular signal transduction, activation, proliferation, differentiation, as well as their intercellular communication. The coordinated execution and integration of such complex signaling programs is critical for physiological angiogenesis to take place in normal growth, development, exercise, and wound healing, while its dysregulation is critically linked to many major human diseases such as cancer, cardiovascular diseases, and ocular disorders; it is also crucial in regenerative medicine. Although huge efforts have been devoted to drug development for these diseases by investigation of angiogenesis‐targeted therapies, only a few therapeutics and targets have proved effective in humans due to the innate multiscale complexity and nonlinearity in the process of angiogenic signaling. As a promising approach that can help better address this challenge, systems biology modeling allows the integration of knowledge across studies and scales and provides a powerful means to mechanistically elucidate and connect the individual molecular and cellular signaling components that function in concert to regulate angiogenesis. In this review, we summarize and discuss how systems biology modeling studies, at the pathway‐, cell‐, tissue‐, and whole body‐levels, have advanced our understanding of signaling in angiogenesis and thereby delivered new translational insights for human diseases.

This article is categorized under:

Cardiovascular Diseases > Computational Models

Cancer > Computational Models

The role of calcium oscillations in the phenotype selection in endothelial cells

Angiogenesis is an important process in the formation and maintenance of tissues which is driven by a complex system of intracellular and intercellular signaling mechanisms. Endothelial cells taking part in early angiogenesis must select their phenotype as either a tip cells (leading, migratory) or a stalk cells (following). Recent experiments have demonstrated that rapid calcium oscillations within active cells characterize this phenotype selection process and that these oscillations play a necessary role in governing phenotype selection and eventual vessel architecture. In this work, we develop a mathematical model capable of describing these oscillations and their role in phenotype selection then use it to improve our understanding of the biological mechanisms at play. We developed a model based on two previously published and experimentally validated mathematical models of calcium and angiogenesis then use our resulting model to simulate various multi-cell scenarios. We are able to capture essential calcium oscillation dynamics and intercellular communication between neighboring cells. The results of our model show that although the late DLL4 (a transmembrane protein that activates Notch pathway) levels of a cell are connected with its initial IP3 (Inositol 1,4,5-trisphosphate) level, cell-to-cell communication determines its eventual phenotype.

Biologically-Based Mathematical Modeling of Tumor Vasculature and Angiogenesis via Time-Resolved Imaging Data

Tumor-associated vasculature is responsible for the delivery of nutrients, removal of waste, and allowing growth beyond 2-3 mm3. Additionally, the vascular network, which is changing in both space and time, fundamentally influences tumor response to both systemic and radiation therapy. Thus, a robust understanding of vascular dynamics is necessary to accurately predict tumor growth, as well as establish optimal treatment protocols to achieve optimal tumor control. Such a goal requires the intimate integration of both theory and experiment. Quantitative and time-resolved imaging methods have emerged as technologies able to visualize and characterize tumor vascular properties before and during therapy at the tissue and cell scale. Parallel to, but separate from those developments, mathematical modeling techniques have been developed to enable in silico investigations into theoretical tumor and vascular dynamics. In particular, recent efforts have sought to integrate both theory and experiment to enable data-driven mathematical modeling. Such mathematical models are calibrated by data obtained from individual tumor-vascular systems to predict future vascular growth, delivery of systemic agents, and response to radiotherapy. In this review, we discuss experimental techniques for visualizing and quantifying vascular dynamics including magnetic resonance imaging, microfluidic devices, and confocal microscopy. We then focus on the integration of these experimental measures with biologically based mathematical models to generate testable predictions.

Topological data analysis distinguishes parameter regimes in the Anderson-Chaplain model of angiogenesis

Angiogenesis is the process by which blood vessels form from pre-existing vessels. It plays a key role in many biological processes, including embryonic development and wound healing, and contributes to many diseases including cancer and rheumatoid arthritis. The structure of the resulting vessel networks determines their ability to deliver nutrients and remove waste products from biological tissues. Here we simulate the Anderson-Chaplain model of angiogenesis at different parameter values and quantify the vessel architectures of the resulting synthetic data. Specifically, we propose a topological data analysis (TDA) pipeline for systematic analysis of the model. TDA is a vibrant and relatively new field of computational mathematics for studying the shape of data. We compute topological and standard descriptors of model simulations generated by different parameter values. We show that TDA of model simulation data stratifies parameter space into regions with similar vessel morphology. The methodologies proposed here are widely applicable to other synthetic and experimental data including wound healing, development, and plant biology.