Formation of VEGF isoform-specific spatial distributions governing angiogenesis: computational analysis

Vascular endothelial growth factor (VEGF) delivery approaches in regenerative medicine

The utilization of growth factors in the process of tissue regeneration has garnered significant interest and has been the subject of extensive research. However, despite the fervent efforts invested in recent clinical trials, a considerable number of these studies have produced outcomes that are deemed unsatisfactory. It is noteworthy that the trials that have yielded the most satisfactory outcomes have exhibited a shared characteristic, namely, the existence of a mechanism for the regulated administration of growth factors. Despite the extensive exploration of drug delivery vehicles and their efficacy in delivering certain growth factors, the development of a reliable predictive approach for the delivery of delicate growth factors like Vascular Endothelial Growth Factor (VEGF) remains elusive. VEGF plays a crucial role in promoting angiogenesis; however, the administration of VEGF demands a meticulous approach as it necessitates precise localization and transportation to a specific target tissue. This process requires prolonged and sustained exposure to a low concentration of VEGF. Inaccurate administration of drugs, either through off-target effects or inadequate delivery, may heighten the risk of adverse reactions and potentially result in tumorigenesis. At present, there is a scarcity of technologies available for the accurate encapsulation of VEGF and its subsequent sustained and controlled release. The objective of this review is to present and assess diverse categories of VEGF administration mechanisms. This paper examines various systems, including polymeric, liposomal, hydrogel, inorganic, polyplexes, and microfluidic, and evaluates the appropriate dosage of VEGF for multiple applications.

Monocytes deposit migrasomes to promote embryonic angiogenesis

Pro-angiogenic factors are key regulators of angiogenesis. Here we report that highly migratory cells patrol the area of capillary formation in chick embryo chorioallantoic membrane. These cells deposit migrasomes on their migration tracks, creating migrasome-enriched areas. Single-cell sequencing identified these cells as monocytes. Depletion of monocytes impairs capillary formation. Quantitative mass spectrometry analysis reveals that monocyte migrasomes are enriched with pro-angiogenic factors. Purified migrasomes promote capillary formation and monocyte recruitment in vivo, and endothelial cell tube formation and monocyte chemotaxis in vitro. Knockdown or knockout of TSPAN4 reduces migrasome formation and impairs capillary formation and monocyte recruitment. Mechanistically, monocytes promote angiogenesis via VEGFA and CXCL12 in migrasomes. In summary, monocytes deposit migrasomes enriched in pro-angiogenic factors to promote angiogenesis.

Mechanical regulation of the early stages of angiogenesis

Favouring or thwarting the development of a vascular network is essential in fields as diverse as oncology, cardiovascular disease or tissue engineering. As a result, understanding and controlling angiogenesis has become a major scientific challenge. Mechanical factors play a fundamental role in angiogenesis and can potentially be exploited for optimizing the architecture of the resulting vascular network. Largely focusing on in vitro systems but also supported by some in vivo evidence, the aim of this Highlight Review is dual. First, we describe the current knowledge with particular focus on the effects of fluid and solid mechanical stimuli on the early stages of the angiogenic process, most notably the destabilization of existing vessels and the initiation and elongation of new vessels. Second, we explore inherent difficulties in the field and propose future perspectives on the use of in vitro and physics-based modelling to overcome these difficulties.

A combined experimental and computational framework to evaluate the behavior of therapeutic cells for peripheral nerve regeneration

Recent studies have explored the potential of tissue‐mimetic scaffolds in encouraging nerve regeneration. One of the major determinants of the regenerative success of cellular nerve repair constructs (NRCs) is the local microenvironment, particularly native low oxygen conditions which can affect implanted cell survival and functional performance. In vivo, cells reside in a range of environmental conditions due to the spatial gradients of nutrient concentrations that are established. Here we evaluate in vitro the differences in cellular behavior that such conditions induce, including key biological features such as oxygen metabolism, glucose consumption, cell death, and vascular endothelial growth factor secretion. Experimental measurements are used to devise and parameterize a mathematical model that describes the behavior of the cells. The proposed model effectively describes the interactions between cells and their microenvironment and could in the future be extended, allowing researchers to compare the behavior of different therapeutic cells. Such a combinatorial approach could be used to accelerate the clinical translation of NRCs by identifying which critical design features should be optimized when fabricating engineered nerve repair conduits.

Inverting angiogenesis with interstitial flow and chemokine matrix-binding affinity

The molecular signaling pathways that orchestrate angiogenesis have been widely studied, but the role of biophysical cues has received less attention. Interstitial flow is unavoidable in vivo, and has been shown to dramatically change the neovascular patterns, but the mechanisms by which flow regulates angiogenesis remain poorly understood. Here, we study the complex interactions between interstitial flow and the affinity for matrix binding of different chemokine isoforms. Using a computational model, we find that changing the matrix affinity of the chemokine isoform can invert the effect of interstitial flow on angiogenesis—from preferential growth in the direction of the flow when the chemokine is initially matrix-bound to preferential flow against the flow when it is unbound. Although fluid forces signal endothelial cells directly, our data suggests a mechanism for the inversion based on biotransport arguments only, and offers a potential explanation for experimental results in which interstitial flow produced preferential vessel growth with and against the flow. Our results point to a particularly intricate effect of interstitial flow on angiogenesis in the tumor microenvironment, where the vessel network geometry and the interstitial flow patterns are complex.

Modelling-informed cell-seeded nerve repair construct designs for treating peripheral nerve injuries

Millions of people worldwide are affected by peripheral nerve injuries (PNI), involving billions of dollars in healthcare costs. Common outcomes for patients include paralysis and loss of sensation, often leading to lifelong pain and disability. Engineered Neural Tissue (EngNT) is being developed as an alternative to the current treatments for large-gap PNIs that show underwhelming functional recovery in many cases. EngNT repair constructs are composed of a stabilised hydrogel cylinder, surrounded by a sheath of material, to mimic the properties of nerve tissue. The technology also enables the spatial seeding of therapeutic cells in the hydrogel to promote nerve regeneration. The identification of mechanisms leading to maximal nerve regeneration and to functional recovery is a central challenge in the design of EngNT repair constructs. Using in vivo experiments in isolation is costly and time-consuming, offering a limited insight on the mechanisms underlying the performance of a given repair construct. To bridge this gap, we derive a cell-solute model and apply it to the case of EngNT repair constructs seeded with therapeutic cells which produce vascular endothelial growth factor (VEGF) under low oxygen conditions to promote vascularisation in the construct. The model comprises a set of coupled non-linear diffusion-reaction equations describing the evolving cell population along with its interactions with oxygen and VEGF fields during the first 24h after transplant into the nerve injury site. This model allows us to evaluate a wide range of repair construct designs (e.g. cell-seeding strategy, sheath material, culture conditions), the idea being that designs performing well over a short timescale could be shortlisted for in vivo trials. In particular, our results suggest that seeding cells beyond a certain density threshold is detrimental regardless of the situation considered, opening new avenues for future nerve tissue engineering.

VEGF-A splice variants bind VEGFRs with differential affinities

Vascular endothelial growth factor A (VEGF-A) and its binding to VEGFRs is an important angiogenesis regulator, especially the earliest-known isoform, VEGF-A165a. Yet several additional splice variants play prominent roles in regulating angiogenesis in health and in vascular disease, including VEGF-A121 and an anti-angiogenic variant, VEGF-A165b. Few studies have attempted to distinguish these forms from their angiogenic counterparts, experimentally. Previous studies of VEGF-A:VEGFR binding have measured binding kinetics for VEGFA165 and VEGF-A121, but binding kinetics of the other two pro- and all anti-angiogenic splice variants are not known. We measured the binding kinetics for VEGF-A165, -A165b, and -A121 with VEGFR1 and VEGF-R2 using surface plasmon resonance. We validated our methods by reproducing the known affinities between VEGF-A165a:VEGFR1 and VEGF-A165a:VEGFR2, 1.0 pM and 10 pM respectively, and validated the known affinity VEGF-A121:VEGFR2 as KD = 0.66 nM. We found that VEGF-A121 also binds VEGFR1 with an affinity KD = 3.7 nM. We further demonstrated that the anti-angiogenic variant, VEGF-A165b selectively prefers VEGFR2 binding at an affinity = 0.67 pM while binding VEGFR1 with a weaker affinity-KD = 1.4 nM. These results suggest that the - A165b anti-angiogenic variant would preferentially bind VEGFR2. These discoveries offer a new paradigm for understanding VEGF-A, while further stressing the need to take care in differentiating the splice variants in all future VEGF-A studies.

Small extracellular vesicle-bound vascular endothelial growth factor secreted by carcinoma-associated fibroblasts promotes angiogenesis in a bevacizumab-resistant manner

The blood vessel growth inhibitor bevacizumab targets vascular endothelial growth factor (VEGF), a crucial regulator of angiogenesis. Recently, small extracellular vesicles (sEVs) have been demonstrated to be important vehicles in the transport of growth factors to target cells. In this study, we isolated primary carcinoma-associated fibroblasts (CAFs) from four human oral squamous cell carcinoma (OSCC) specimens. Compared with other non-extracellular vesicle components, CAF-derived sEVs were found to be the main regulators of angiogenesis. The ability of CAF sEVs to activate VEGF receptor 2 (VEGFR2) signaling in human umbilical vein endothelial cells (HUVEC) was dependent on the association between sEVs and VEGF. In addition, sEV-bound VEGF secreted by CAFs further activated VEGFR2 signaling in HUVEC in a bevacizumab-resistant manner. VEGF was found to interact with heparan sulfate proteoglycans on the CAF sEV surface and could be released by heparinase I/III. The bioactivity of the dissociated VEGF was retained in vitro and in vivo and could be neutralized by bevacizumab. These findings suggest that the combined use of heparinase and bevacizumab might inhibit angiogenesis in patients with high levels of sEV-bound VEGF.

Clinical value of TV-CDS combined with serum tumor markers in diagnosis of ovarian cancer

Clinical application value was investigated of transvaginal color Doppler ultrasound (TV-CDS) combined with serum tumor markers carbohydrate antigen 125 (CA125), vascular endothelial growth factor (VEGF) and osteopontin (OPN) in the diagnosis of ovarian cancer (OC). One hundred and six patients with OC [malignant tumor group (MTG)] and fifty patients with benign ovarian diseases [benign control group (BCG)] were selected. Both groups of patients underwent TV-CDS examination. The lesion morphology and internal structure were observed, and the tumor blood flow signal, resistance index (RI) and pulsability index (PI) under ultrasound were determined. Serum CA125 was detected by electrochemiluminescence, and VEGF and OPN levels were detected by enzyme-linked immunosorbent assay. The incidence of irregular lesion morphology, unclear boundary, uneven internal echo, microcalcification and side-acoustic images in OC group (OCG) was significantly higher than that in BCG (P<0.01). As for blood flow grading, most patients in the MTG were in grade II and III, while most patients in the BCG were in grade 0. Compared with BCG, the flow RI and PI in the OCG were significantly reduced (P<0.01). The levels of serum CA125, VEGF and OPN in OCG were significantly higher than those in BCG. The expression levels of serum CA125, VEGF and OPN in OC patients with clinical high stage (stage III and IV), poorly differentiated, ascites, recurrence and metastasis were significantly higher than those in patients with clinical low stage (stage I and II), well differentiated, no ascites and no recurrence and metastasis (P<0.05). With the disappearance of the tumor or the decrease of tumor load, the serum marker levels after treatment were significantly lower than that before treatment (P<0.05). The sensitivity and accuracy of the combined examination in the diagnosis of OC were obviously improved compared with the single and partial combined examinations (P<0.05). In conclusion, combined examination can significantly improve the sensitivity and accuracy of OC, which is conducive to early diagnosis and clinical intervention of OC.

Eyeing the Extracellular Matrix in Vascular Development and Microvascular Diseases and Bridging the Divide between Vascular Mechanics and Function

The extracellular matrix (ECM) is critical in all aspects of vascular development and health: supporting cell anchorage, providing structure, organization and mechanical stability, and serving as a sink for growth factors and sustained survival signals. Abnormal changes in ECM protein expression, organization, and/or properties, and the ensuing changes in vascular compliance affect vasodilator responses, microvascular pressure transmission, and collateral perfusion. The changes in microvascular compliance are independent factors initiating, driving, and/or exacerbating a plethora of microvascular diseases of the eye including diabetic retinopathy (DR) and vitreoretinopathy, retinopathy of prematurity (ROP), wet age-related macular degeneration (AMD), and neovascular glaucoma. Congruently, one of the major challenges with most vascular regenerative therapies utilizing localized growth factor, endothelial progenitor, or genetically engineered cell delivery, is the regeneration of blood vessels with physiological compliance properties. Interestingly, vascular cells sense physical forces, including the stiffness of their ECM, through mechanosensitive integrins, their associated proteins and the actomyosin cytoskeleton, which generates biochemical signals that culminate in a rapid expression of matricellular proteins such as cellular communication network 1 (CCN1) and CCN2 (aka connective tissue growth factor or CTGF). Loss or gain of function of these proteins alters genetic programs of cell growth, ECM biosynthesis, and intercellular signaling, that culminate in changes in cell behavior, polarization, and barrier function. In particular, the function of the matricellular protein CCN2/CTGF is critical during retinal vessel development and regeneration wherein new blood vessels form and invest a preformed avascular neural retina following putative gradients of matrix stiffness. These observations underscore the need for further in-depth characterization of the ECM-derived cues that dictate structural and functional properties of the microvasculature, along with the development of new therapeutic strategies addressing the ECM-dependent regulation of pathophysiological stiffening of blood vessels in ischemic retinopathies.

Sex differences in cancer mechanisms

We now know that cancer is many different diseases, with great variation even within a single histological subtype. With the current emphasis on developing personalized approaches to cancer treatment, it is astonishing that we have not yet systematically incorporated the biology of sex differences into our paradigms for laboratory and clinical cancer research. While some sex differences in cancer arise through the actions of circulating sex hormones, other sex differences are independent of estrogen, testosterone, or progesterone levels. Instead, these differences are the result of sexual differentiation, a process that involves genetic and epigenetic mechanisms, in addition to acute sex hormone actions. Sexual differentiation begins with fertilization and continues beyond menopause. It affects virtually every body system, resulting in marked sex differences in such areas as growth, lifespan, metabolism, and immunity, all of which can impact on cancer progression, treatment response, and survival. These organismal level differences have correlates at the cellular level, and thus, males and females can fundamentally differ in their protections and vulnerabilities to cancer, from cellular transformation through all stages of progression, spread, and response to treatment. Our goal in this review is to cover some of the robust sex differences that exist in core cancer pathways and to make the case for inclusion of sex as a biological variable in all laboratory and clinical cancer research. We finish with a discussion of lab- and clinic-based experimental design that should be used when testing whether sex matters and the appropriate statistical models to apply in data analysis for rigorous evaluations of potential sex effects. It is our goal to facilitate the evaluation of sex differences in cancer in order to improve outcomes for all patients.

Imaging the Dynamic Interaction Between Sprouting Microvessels and the Extracellular Matrix

Thorough understanding of growth and evolution of tissue vasculature is fundamental to many fields of medicine including cancer therapy, wound healing, and tissue engineering. Angiogenesis, the growth of new vessels from existing ones, is dynamically influenced by a variety of environmental factors, including mechanical and biophysical factors, chemotactic factors, proteolysis, and interaction with stromal cells. Yet, dynamic interactions between neovessels and their environment are difficult to study with traditional fixed time imaging techniques. Advancements in imaging technologies permit time-series and volumetric imaging, affording the ability to visualize microvessel growth over 3D space and time. Time-lapse imaging has led to more informative investigations of angiogenesis. The environmental factors implicated in angiogenesis span a wide range of signals. Neovessels advance through stromal matrices by forming attachments and pulling and pushing on their microenvironment, reorganizing matrix fibers, and inducing large deformations of the surrounding stroma. Concurrently, neovessels secrete proteolytic enzymes to degrade their basement membrane, create space for new vessels to grow, and release matrix-bound cytokines. Growing neovessels also respond to a host of soluble and matrix-bound growth factors, and display preferential growth along a cytokine gradient. Lastly, stromal cells such as macrophages and mesenchymal stem cells (MSCs) interact directly with neovessels and their surrounding matrix to facilitate sprouting, vessel fusion, and tissue remodeling. This review highlights how time-lapse imaging techniques advanced our understanding of the interaction of blood vessels with their environment during sprouting angiogenesis. The technology provides means to characterize the evolution of microvessel behavior, providing new insights and holding great promise for further research on the process of angiogenesis.

Single-Cell Receptor Quantification of an In Vitro Coculture Angiogenesis Model Reveals VEGFR, NRP1, Tie2, and PDGFR Regulation and Endothelial Heterogeneity

Angiogenesis, the formation of new blood vessels from pre-existing ones, is essential for both normal development and numerous pathologies. Systems biology has offered a unique approach to study angiogenesis by profiling tyrosine kinase receptors (RTKs) that regulate angiogenic processes and computationally modeling RTK signaling pathways. Historically, this systems biology approach has been applied on ex vivo angiogenesis assays, however, these assays are difficult to quantify and limited in their potential of temporal analysis. In this study, we adopted a simple two-dimensional angiogenesis assay comprised of human umbilical vein endothelial cells (HUVECs) and human dermal fibroblasts (HDFs) and examined temporal dynamics of a panel of six RTKs and cell heterogeneity up to 17 days. We observed ~2700 VEGFR1 (vascular endothelial growth factor receptor 1) per cell on 24-h-old cocultured HDF plasma membranes, which do not express VEGFR when cultured alone. We observed 4000–8100 VEGFR2 per cell on cocultured HUVEC plasma membranes throughout endothelial tube formation. We showed steady increase of platelet-derived growth factor receptors (PDGFRs) on cocultured HDF plasma membranes, and more interestingly, 1900–2900 PDGFRβ per plasma membrane were found on HUVECs within the first six hours of coculturing. These quantitative findings will offer us insights into molecular regulation during angiogenesis and help assess in vitro tube formation models and their physiological relevance.

Mapping Tyrosine Kinase Receptor Dimerization to Receptor Expression and Ligand Affinities

Tyrosine kinase receptor (RTK) ligation and dimerization is a key mechanism for translating external cell stimuli into internal signaling events. This process is critical to several key cell and physiological processes, such as in angiogenesis and embryogenesis, among others. While modulating RTK activation is a promising therapeutic target, RTK signaling axes have been shown to involve complicated interactions between ligands and receptors both within and across different protein families. In angiogenesis, for example, several signaling protein families, including vascular endothelial growth factors and platelet-derived growth factors, exhibit significant cross-family interactions that can influence pathway activation. Computational approaches can provide key insight to detangle these signaling pathways but have been limited by the sparse knowledge of these cross-family interactions. Here, we present a framework for studying known and potential non-canonical interactions. We constructed generalized models of RTK ligation and dimerization for systems of two, three and four receptor types and different degrees of cross-family ligation. Across each model, we developed parameter-space maps that fully determine relative pathway activation for any set of ligand-receptor binding constants, ligand concentrations and receptor concentrations. Therefore, our generalized models serve as a powerful reference tool for predicting not only known ligand: Receptor axes but also how unknown interactions could alter signaling dimerization patterns. Accordingly, it will drive the exploration of cross-family interactions and help guide therapeutic developments across processes like cancer and cardiovascular diseases, which depend on RTK-mediated signaling.

Soluble VEGFR1 signaling guides vascular patterns into dense branching morphologies

Vascular patterning is a key process during development and disease. The diffusive decoy receptor sVEGFR1 (sFlt1) is a known regulator of endothelial cell behavior, yet the mechanism by which it controls vascular structure is little understood. We propose computational models to shed light on how vascular patterning is guided by self-organized gradients of the VEGF/sVEGFR1 factors. We demonstrate that a diffusive inhibitor can generate structures with a dense branching morphology in models where the activator elicits directed growth. Inadequate presence of the inhibitor leads to compact growth, while excessive production of the inhibitor blocks expansion and stabilizes existing structures. Model predictions were compared with time-resolved experimental data obtained from endothelial sprout kinetics in fibrin gels. In the presence of inhibitory antibodies against VEGFR1 vascular sprout density increases while the speed of sprout expansion remains unchanged. Thus, the rate of secretion and stability of extracellular sVEGFR1 can modulate vascular sprout density.

A sulfonated reversible thermal gel for the spatiotemporal control of VEGF delivery to promote therapeutic angiogenesis

Despite medical and surgical advancements for the treatment of cardiovascular disease, mortality and morbidity remain high. Therapeutic angiogenesis has been one approach to address the major clinical need for a more effective treatment to restoring blood flow in ischemic organs and tissues, but current progress in angiogenic drug delivery is inadequate at providing sufficient bioavailability without causing safety concerns. An injectable sulfonated reversible thermal gel composed of a polyurea conjugated with poly(N-isopropylacrylamide) and sulfonate groups has been developed for the delivery of angiogenic factors. The thermal gel allowed for the spatiotemporal control of vascular endothelial growth factor release with a decreased initial burst release and reduced release rate in vitro. A subcutaneous injection mouse model was used to evaluate efficacious vascularization and assess the inflammatory response due to a foreign body. Thermal gel injections showed substantial vascularization properties by inducing vessel formation, recruitment and differentiation of vascular endothelial cells, and vessel stabilization by perivascular cells, while infiltrating macrophages due to the thermal gel injections decreased over time. These results demonstrated effective localization and delivery of angiogenic factors for therapeutic angiogenesis. © 2018 Wiley Periodicals, Inc. J Biomed Mater Res Part A: 106A: 3053-3064, 2018.

The Role of Biomaterials as Angiogenic Modulators of Spinal Cord Injury: Mimetics of the Spinal Cord, Cell and Angiogenic Factor Delivery Agents

Spinal cord injury (SCI) represents an extremely debilitating condition for which no efficacious treatment is available. One of the main contributors to the inhospitable environment found in SCI is the vascular disruption that happens at the moment of injury that compromises the blood-spinal cord barrier (BSCB) and triggers a cascade of events that includes infiltration of inflammatory cells, ischemia and intraparenchymal hemorrhage. Due to the unsatisfactory nature of revascularization following SCI, restoring vascular perfusion and the BSCB seems an interesting way of modulating the lesion environment into a regenerative phenotype, with a potential increase in functional recovery. Certain biomaterials possess interesting features to enhance SCI therapies, and in fact have been applied as angiogenic promoters in other pathologies. The present mini-review intends to highlight the contribution that biomaterials could make in the development of novel therapeutic solutions able to restore proper vascularization and the BSCB.

A Poroelasticity Theory Approach to Study the Mechanisms Leading to Elevated Interstitial Fluid Pressure in Solid Tumours

Although the mechanisms responsible for elevated interstitial fluid pressure (IFP) in tumours remain obscure, it seems clear that high IFP represents a barrier to drug delivery (since the resulting adverse pressure gradient implies a reduction in the driving force for transvascular exchange of both fluid and macromolecules). R. Jain and co-workers studied this problem, and although the conclusions drawn from their idealized mathematical models offered useful insights into the causes of elevated IFP, they by no means gave a definitive explanation for this phenomenon. In this paper, we use poroelasticity theory to also develop a macroscopic mathematical model to describe the time evolution of a solid tumour, but focus our attention on the mechanisms responsible for the rise of the IFP, from that for a healthy interstitium to that measured in malignant tumours. In particular, we discuss a number of possible time scales suggested by our mathematical model and propose a tumour-dependent time scale that leads to results in agreement with experimental observations. We apply our mathematical model to simulate the effect of "vascular normalization" (as proposed by Jain in Nat Med 7:987-989, 2001) on the IFP profile and discuss and contrast our conclusions with those of previous work in the literature.

Retinal vasculature development in health and disease

Development of the retinal vasculature is based on highly coordinated signalling between different cell types of the retina, integrating internal metabolic requirements with external influences such as the supply of oxygen and nutrients. The developing mouse retinal vasculature is a useful model system to study these interactions because it is experimentally accessible for intra ocular injections and genetic manipulations, can be easily imaged and develops in a similar fashion to that of humans. Research using this model has provided insights about general principles of angiogenesis as well as pathologies that affect the developing retinal vasculature. In this review, we discuss recent advances in our understanding of the molecular and cellular mechanisms that govern the interactions between neurons, glial and vascular cells in the developing retina. This includes a review of mechanisms that shape the retinal vasculature, such as sprouting angiogenesis, vascular network remodelling and vessel maturation. We also explore how the disruption of these processes in mice can lead to pathology - such as oxygen induced retinopathy - and how this translates to human retinopathy of prematurity.

A computational analysis of in vivo VEGFR activation by multiple co-expressed ligands

The splice isoforms of vascular endothelial growth A (VEGF) each have different affinities for the extracellular matrix (ECM) and the coreceptor NRP1, which leads to distinct vascular phenotypes in model systems expressing only a single VEGF isoform. ECM-immobilized VEGF can bind to and activate VEGF receptor 2 (VEGFR2) directly, with a different pattern of site-specific phosphorylation than diffusible VEGF. To date, the way in which ECM binding alters the distribution of isoforms of VEGF and of the related placental growth factor (PlGF) in the body and resulting angiogenic signaling is not well-understood. Here, we extend our previous validated cell-level computational model of VEGFR2 ligation, intracellular trafficking, and site-specific phosphorylation, which captured differences in signaling by soluble and immobilized VEGF, to a multi-scale whole-body framework. This computational systems pharmacology model captures the ability of the ECM to regulate isoform-specific growth factor distribution distinctly for VEGF and PlGF, and to buffer free VEGF and PlGF levels in tissue. We show that binding of immobilized growth factor to VEGF receptors, both on endothelial cells and soluble VEGFR1, is likely important to signaling in vivo. Additionally, our model predicts that VEGF isoform-specific properties lead to distinct profiles of VEGFR1 and VEGFR2 binding and VEGFR2 site-specific phosphorylation in vivo, mediated by Neuropilin-1. These predicted signaling changes mirror those observed in murine systems expressing single VEGF isoforms. Simulations predict that, contrary to the 'ligand-shifting hypothesis,' VEGF and PlGF do not compete for receptor binding at physiological concentrations, though PlGF is predicted to slightly increase VEGFR2 phosphorylation when over-expressed by 10-fold. These results are critical to design of appropriate therapeutic strategies to control VEGF availability and signaling in regenerative medicine applications.

Deterministic and Stochastic Study for a Microscopic Angiogenesis Model: Applications to the Lewis Lung Carcinoma

Angiogenesis modelling is an important tool to understand the underlying mechanisms yielding tumour growth. Nevertheless, there is usually a gap between models and experimental data. We propose a model based on the intrinsic microscopic reactions defining the angiogenesis process to link experimental data with previous macroscopic models. The microscopic characterisation can describe the macroscopic behaviour of the tumour, which stability analysis reveals a set of predicted tumour states involving different morphologies. Additionally, the microscopic description also gives a framework to study the intrinsic stochasticity of the reactive system through the resulting Langevin equation. To follow the goal of the paper, we use available experimental information on the Lewis lung carcinoma to infer meaningful parameters for the model that are able to describe the different stages of the tumour growth. Finally we explore the predictive capabilities of the fitted model by showing that fluctuations are determinant for the survival of the tumour during the first week and that available treatments can give raise to new stable tumour dormant states with a reduced vascular network.

Agent-based model of angiogenesis simulates capillary sprout initiation in multicellular networks

Many biological processes are controlled by both deterministic and stochastic influences. However, efforts to model these systems often rely on either purely stochastic or purely rule-based methods. To better understand the balance between stochasticity and determinism in biological processes a computational approach that incorporates both influences may afford additional insight into underlying biological mechanisms that give rise to emergent system properties. We apply a combined approach to the simulation and study of angiogenesis, the growth of new blood vessels from existing networks. This complex multicellular process begins with selection of an initiating endothelial cell, or tip cell, which sprouts from the parent vessels in response to stimulation by exogenous cues. We have constructed an agent-based model of sprouting angiogenesis to evaluate endothelial cell sprout initiation frequency and location, and we have experimentally validated it using high-resolution time-lapse confocal microscopy. ABM simulations were then compared to a Monte Carlo model, revealing that purely stochastic simulations could not generate sprout locations as accurately as the rule-informed agent-based model. These findings support the use of rule-based approaches for modeling the complex mechanisms underlying sprouting angiogenesis over purely stochastic methods.

The Force at the Tip--Modelling Tension and Proliferation in Sprouting Angiogenesis

Sprouting angiogenesis, where new blood vessels grow from pre-existing ones, is a complex process where biochemical and mechanical signals regulate endothelial cell proliferation and movement. Therefore, a mathematical description of sprouting angiogenesis has to take into consideration biological signals as well as relevant physical processes, in particular the mechanical interplay between adjacent endothelial cells and the extracellular microenvironment. In this work, we introduce the first phase-field continuous model of sprouting angiogenesis capable of predicting sprout morphology as a function of the elastic properties of the tissues and the traction forces exerted by the cells. The model is very compact, only consisting of three coupled partial differential equations, and has the clear advantage of a reduced number of parameters. This model allows us to describe sprout growth as a function of the cell-cell adhesion forces and the traction force exerted by the sprout tip cell. In the absence of proliferation, we observe that the sprout either achieves a maximum length or, when the traction and adhesion are very large, it breaks. Endothelial cell proliferation alters significantly sprout morphology, and we explore how different types of endothelial cell proliferation regulation are able to determine the shape of the growing sprout. The largest region in parameter space with well formed long and straight sprouts is obtained always when the proliferation is triggered by endothelial cell strain and its rate grows with angiogenic factor concentration. We conclude that in this scenario the tip cell has the role of creating a tension in the cells that follow its lead. On those first stalk cells, this tension produces strain and/or empty spaces, inevitably triggering cell proliferation. The new cells occupy the space behind the tip, the tension decreases, and the process restarts. Our results highlight the ability of mathematical models to suggest relevant hypotheses with respect to the role of forces in sprouting, hence underlining the necessary collaboration between modelling and molecular biology techniques to improve the current state-of-the-art.

Sprouting angiogenesis—a process by which new blood vessels grow from existing ones—is an ubiquitous phenomenon in health and disease of higher organisms, playing a crucial role in organogenesis, wound healing, inflammation, as well as on the onset and progression of over 50 different diseases such as cancer, rheumatoid arthritis and diabetes. Mathematical models have the ability to suggest relevant hypotheses with respect to the mechanisms of cell movement and rearrangement within growing vessel sprouts. The inclusion of both biochemical and mechanical processes in a mathematical model of sprouting angiogenesis permits to describe sprout extension as a function of the forces exerted by the cells in the tissue. It also allows to question the regulation of biochemical processes by mechanical forces and vice-versa. In this work we present a compact model of sprouting angiogenesis that includes the mechanical characteristics of the vessel and the tissue. We use this model to suggest the mechanism for the regulation of proliferation within sprout formation. We conclude that the tip cell has the role of creating a tension in the cells that follow its lead. On those first cells of the stalk, this tension produces strain and/or empty spaces, inevitably triggering cell proliferation. The new cells occupy the space behind the tip, the tension decreases, and the process restarts. The modelling strategy used, deemed phase-field, permits to describe the evolution of the shape of different domains in complex systems. It is focused on the movement of the interfaces between the domains, and not on an exhaustive description of the transport properties within each domain. For this reason, it requires a reduced number of parameters, and has been used extensively in modelling other biological phenomena such as tumor growth. The coupling of mechanical and biochemical processes in a compact mathematical model of angiogenesis will enable the study of lumen formation and aneurisms in the near future. Also, this framework will allow the study of the action of flow in vessel remodelling, since local forces can readily be coupled with cell movement to obtain the final vessel morphology.

Site-Specific Phosphorylation of VEGFR2 Is Mediated by Receptor Trafficking: Insights from a Computational Model

Matrix-binding isoforms and non-matrix-binding isoforms of vascular endothelial growth factor (VEGF) are both capable of stimulating vascular remodeling, but the resulting blood vessel networks are structurally and functionally different. Here, we develop and validate a computational model of the binding of soluble and immobilized ligands to VEGF receptor 2 (VEGFR2), the endosomal trafficking of VEGFR2, and site-specific VEGFR2 tyrosine phosphorylation to study differences in induced signaling between these VEGF isoforms. In capturing essential features of VEGFR2 signaling and trafficking, our model suggests that VEGFR2 trafficking parameters are largely consistent across multiple endothelial cell lines. Simulations demonstrate distinct localization of VEGFR2 phosphorylated on Y1175 and Y1214. This is the first model to clearly show that differences in site-specific VEGFR2 activation when stimulated with immobilized VEGF compared to soluble VEGF can be accounted for by altered trafficking of VEGFR2 without an intrinsic difference in receptor activation. The model predicts that Neuropilin-1 can induce differences in the surface-to-internal distribution of VEGFR2. Simulations also show that ligated VEGFR2 and phosphorylated VEGFR2 levels diverge over time following stimulation. Using this model, we identify multiple key levers that alter how VEGF binding to VEGFR2 results in different coordinated patterns of multiple downstream signaling pathways. Specifically, simulations predict that VEGF immobilization, interactions with Neuropilin-1, perturbations of VEGFR2 trafficking, and changes in expression or activity of phosphatases acting on VEGFR2 all affect the magnitude, duration, and relative strength of VEGFR2 phosphorylation on tyrosines 1175 and 1214, and they do so predictably within our single consistent model framework.

Vascular endothelial growth factor (VEGF) is an important regulator of blood vessel growth. To date, therapies attempting to harness the VEGF system to promote blood vessel growth (e.g. for wound healing or ischemic disease) have achieved only limited success. To improve VEGF-based therapies, we need to better understand how VEGF promotes development of functional blood vessels. We have developed a computational model of VEGF binding to the receptor VEGFR2, trafficking of VEGFR2 through endosomal compartments in the cell, and activation of VEGFR2 on several tyrosine residues. The pattern of tyrosines activated on VEGFR2 influences cell behavior, promoting cell survival, proliferation, or migration. The combination of these cues influences the diameter of vessels, degree of branching, and leakiness of the resultant vessel network. Our model shows that changes in VEGFR2 trafficking as a result of VEGF immobilization to the extracellular matrix are sufficient to describe observed changes in the pattern of VEGFR2 activation compared to stimulation with purely soluble VEGF. This model can be used to predict how VEGF immobilization, interactions with co-receptors or proteins that deactivate VEGFR2, and changes to VEGFR2 trafficking can be tuned to promote development of functional blood vessel networks for tissue engineering applications.

Systems biology of the microvasculature

The vascular network carries blood throughout the body, delivering oxygen to tissues and providing a pathway for communication between distant organs. The network is hierarchical and structured, but also dynamic, especially at the smaller scales. Remodeling of the microvasculature occurs in response to local changes in oxygen, gene expression, cell-cell communication, and chemical and mechanical stimuli from the microenvironment. These local changes occur as a result of physiological processes such as growth and exercise, as well as acute and chronic diseases including stroke, cancer, and diabetes, and pharmacological intervention. While the vasculature is an important therapeutic target in many diseases, drugs designed to inhibit vascular growth have achieved only limited success, and no drug has yet been approved to promote therapeutic vascular remodeling. This highlights the challenges involved in identifying appropriate therapeutic targets in a system as complex as the vasculature. Systems biology approaches provide a means to bridge current understanding of the vascular system, from detailed signaling dynamics measured in vitro and pre-clinical animal models of vascular disease, to a more complete picture of vascular regulation in vivo. This will translate to an improved ability to identify multi-component biomarkers for diagnosis, prognosis, and monitoring of therapy that are easy to measure in vivo, as well as better drug targets for specific disease states. In this review, we summarize systems biology approaches that have advanced our understanding of vascular function and dysfunction in vivo, with a focus on computational modeling.

Applications of computational models to better understand microvascular remodelling: a focus on biomechanical integration across scales

Microvascular network remodelling is a common denominator for multiple pathologies and involves both angiogenesis, defined as the sprouting of new capillaries, and network patterning associated with the organization and connectivity of existing vessels. Much of what we know about microvascular remodelling at the network, cellular and molecular scales has been derived from reductionist biological experiments, yet what happens when the experiments provide incomplete (or only qualitative) information? This review will emphasize the value of applying computational approaches to advance our understanding of the underlying mechanisms and effects of microvascular remodelling. Examples of individual computational models applied to each of the scales will highlight the potential of answering specific questions that cannot be answered using typical biological experimentation alone. Looking into the future, we will also identify the needs and challenges associated with integrating computational models across scales.

Galectins in angiogenesis: consequences for gestation

Members of the galectin family have been shown to exert several roles in the context of reproduction. They contribute to placentation, maternal immune regulation and facilitate angiogenesis encompassing decidualisation and placenta formation during pregnancy. In the context of neo-vascularisation, galectins have been shown to augment signalling pathways that lead to endothelial cell activation, cell proliferation, migration and tube formation in vitro in addition to angiogenesis in vivo. Angiogenesis during gestation ensures not only proper foetal growth and development, but also maternal health. Consequently, restriction of placental blood flow has major consequences for both foetus and mother, leading to pregnancy diseases. In this review we summarise both the established and the emerging roles of galectin in angiogenesis and discuss the possible implications during healthy and pathological gestation.

Pericellular proteolysis in cancer

Pericellular proteases have long been associated with cancer invasion and metastasis due to their ability to degrade extracellular matrix components. Recent studies demonstrate that proteases also modulate tumor progression and metastasis through highly regulated and complex processes involving cleavage, processing, or shedding of cell adhesion molecules, growth factors, cytokines, and kinases. In this review, we address how cancer cells, together with their surrounding microenvironment, regulate pericellular proteolysis. We dissect the multitude of mechanisms by which pericellular proteases contribute to cancer progression and discuss how this knowledge can be integrated into therapeutic opportunities.

Design of growth factor sequestering biomaterials

Growth factors (GFs) are major regulatory proteins that can govern cell fate, migration, and organization. Numerous aspects of the cell milieu can modulate cell responses to GFs, and GF regulation is often achieved by the native extracellular matrix (ECM). For example, the ECM can sequester GFs and thereby control GF bioavailability. In addition, GFs can exert distinct effects depending on whether they are sequestered in solution, at two-dimensional interfaces, or within three-dimensional matrices. Understanding how the context of GF sequestering impacts cell function in the native ECM can instruct the design of soluble or insoluble GF sequestering moieties, which can then be used in a variety of bioengineering applications. This Feature Article provides an overview of the natural mechanisms of GF sequestering in the cell milieu, and reviews the recent bioengineering approaches that have sequestered GFs to modulate cell function. Results to date demonstrate that the cell response to GF sequestering depends on the affinity of the sequestering interaction, the spatial proximity of sequestering in relation to cells, the source of the GF (supplemented or endogenous), and the phase of the sequestering moiety (soluble or insoluble). We highlight the importance of context for the future design of biomaterials that can leverage endogenous molecules in the cell milieu and mitigate the need for supplemented factors.

Computational systems biology approaches to anti-angiogenic cancer therapeutics

Angiogenesis is an exquisitely regulated process that is required for physiological processes and is also important in numerous diseases. Tumors utilize angiogenesis to generate the vascular network needed to supply the cancer cells with nutrients and oxygen, and many cancer drugs aim to inhibit tumor angiogenesis. Anti-angiogenic therapy involves inhibiting multiple cell types, molecular targets, and intracellular signaling pathways. Computational tools are useful in guiding treatment strategies, predicting the response to treatment, and identifying new targets of interest. Here, we describe progress that has been made in applying mathematical modeling and bioinformatics approaches to study anti-angiogenic therapeutics in cancer.

Serum-dependence of affinity-mediated VEGF release from biomimetic microspheres

Vascular endothelial growth factor (VEGF) activity is highly regulated via sequestering within the ECM and cell-demanded proteolysis to release the sequestered VEGF. Numerous studies have demonstrated that VEGF activity mediates cellular events leading to angiogenesis and capillary formation in vivo. This has motivated the study of biomaterials to sustain VEGF release, and in many cases, the materials are inspired by the structure and function of the native ECM. However, there remains a need for materials that can bind to VEGF with high specificity, as the in vivo environment is rich in a variety of growth factors (GFs) and GF-binding moieties. Here we describe a strategy to control VEGF release using hydrogel microspheres with tethered peptides derived from VEGF receptor 2 (VEGFR2). Using biomaterials covalently modified with varying concentrations of two distinct VEGFR2-derived peptides with varying serum stability, we analyzed both biomaterial and environmental variables that influence VEGF release and activity. The presence of tethered VEGF-binding peptides (VBPs) resulted in significantly extended VEGF release relative to control conditions, and the resulting released VEGF significantly increased the expansion of human umbilical vein endothelial cells in culture. VEGF release rates were also strongly influenced by the concentration of serum. The presence of Feline McDonough Sarcoma-like tyrosine kinase 1 (sFlt-1), a serum-borne receptor fragment derived from VEGF receptor 1, increased VEGF release rates, although sFlt-1 was not sufficient to recapitulate the release profile of VEGF in serum. Further, the influence of serum on VEGF release was not due to protease activity or nonspecific VEGF interactions in the presence of serum-borne heparin. VEGF release kinetics correlated well with a generalizable mathematical model describing affinity-mediated release of VEGF from hydrogel microspheres in defined conditions. Modeling results suggest a potential mechanism whereby competition between VEGF and multiple VEGF-binding serum proteins including sFlt-1, soluble kinase insert domain receptor (sKDR), and α2-macroglobulin (α2-M) likely influenced VEGF release from microspheres. The materials and mathematical model described in this approach may be useful in a range of applications in which sustained, biologically active GF release of a specific GF is desirable.

Analysing calcium dependent and independent regulation of eNOS in endothelium triggered by extracellular signalling events

The vascular endothelium, the intima of blood vessels, coordinately interacts with several biochemical factors expressing endothelial nitric oxide synthase (eNOS) to produce nitric oxide (NO), a potent endogenous vasodilator. The present study investigated the regulation of eNOS by multiple molecular signal transduction pathways, namely vascular endothelial growth factor (VEGF-A) and shear stress which are implicated in the process of angiogenesis and vascular remodelling respectively. In response to signal transduction upstream by VEGF-A and shear stress, different signalling pathways mediated by kinases and intracellular calcium potentiates eNOS activation leading to nitric oxide release. Our study revealed a distinct pattern of eNOS activation driven by VEGF-A and shear stress, maintaining the signalling specificity of the respective pathways. A transient response to eNOS activation was observed under VEGF-A and shear stress stimulus when mediated by calcium dependent cascades, whereas a sustained response was produced by calcium independent vascular signalling kinases. Furthermore, we found that the basal arterial shear stress enhanced eNOS activity when stimulated synergistically even at low VEGF-A levels which might be utilized to facilitate specific endothelial cell functions. Moreover, our study revealed that the presence of PI3K imparted transient behaviour to PLCγ1 supporting the hypothesis that regression and formation of tube structures are mediated by PLCγ1 and PI3K respectively in endothelial cells. This fact is corroborated by the absence of transient behaviour when PI3K is inhibited. We therefore obtained subtle insights into the control mechanism governing the role of specific signalling proteins which are obligate for the regulation of endothelial cell function and the consequent modulation of the nitric oxide release pattern.

Extracellular regulation of VEGF: isoforms, proteolysis, and vascular patterning

The regulation of vascular endothelial growth factor A (VEGF) is critical to neovascularization in numerous tissues under physiological and pathological conditions. VEGF has multiple isoforms, created by alternative splicing or proteolytic cleavage, and characterized by different receptor-binding and matrix-binding properties. These isoforms are known to give rise to a spectrum of angiogenesis patterns marked by differences in branching, which has functional implications for tissues. In this review, we detail the extensive extracellular regulation of VEGF and the ability of VEGF to dictate the vascular phenotype. We explore the role of VEGF-releasing proteases and soluble carrier molecules on VEGF activity. While proteases such as MMP9 can 'release' matrix-bound VEGF and promote angiogenesis, for example as a key step in carcinogenesis, proteases can also suppress VEGF's angiogenic effects. We explore what dictates pro- or anti-angiogenic behavior. We also seek to understand the phenomenon of VEGF gradient formation. Strong VEGF gradients are thought to be due to decreased rates of diffusion from reversible matrix binding, however theoretical studies show that this scenario cannot give rise to lasting VEGF gradients in vivo. We propose that gradients are formed through degradation of sequestered VEGF. Finally, we review how different aspects of the VEGF signal, such as its concentration, gradient, matrix-binding, and NRP1-binding can differentially affect angiogenesis. We explore how this allows VEGF to regulate the formation of vascular networks across a spectrum of high to low branching densities, and from normal to pathological angiogenesis. A better understanding of the control of angiogenesis is necessary to improve upon limitations of current angiogenic therapies.

A systems biology view of blood vessel growth and remodelling

Blood travels throughout the body in an extensive network of vessels – arteries, veins and capillaries. This vascular network is not static, but instead dynamically remodels in response to stimuli from cells in the nearby tissue. In particular, the smallest vessels – arterioles, venules and capillaries – can be extended, expanded or pruned, in response to exercise, ischaemic events, pharmacological interventions, or other physiological and pathophysiological events. In this review, we describe the multi-step morphogenic process of angiogenesis – the sprouting of new blood vessels – and the stability of vascular networks in vivo. In particular, we review the known interactions between endothelial cells and the various blood cells and plasma components they convey. We describe progress that has been made in applying computational modelling, quantitative biology and high-throughput experimentation to the angiogenesis process.

VEGF-A mRNA processing, stability and translation: a paradigm for intricate regulation of gene expression at the post-transcriptional level

Vascular Endothelial Growth Factor A (VEGF-A) is a potent secreted mitogen crucial for physiological and pathological angiogenesis. Post-transcriptional regulation of VEGF-A occurs at multiple levels. Firstly, alternative splicing gives rise to different transcript variants encoding diverse isoforms that exhibit distinct biological properties with regard to receptor binding and extra-cellular localization. Secondly, VEGF-A mRNA stability is regulated by effectors such as hypoxia or growth factors through the binding of stabilizing and destabilizing proteins at AU-rich elements located in the 3′-untranslated region. Thirdly, translation of VEGF-A mRNA is a controlled process involving alternative initiation codons, internal ribosome entry sites (IRESs), an upstream open reading frame (uORF), miRNA targeting and a riboswitch in the 3′ untranslated region. These different levels of regulation cooperate for the crucial fine-tuning of the expression of VEGF-A variants. This review will be focused on our current knowledge of the complex post-transcriptional regulatory switches that modulate the cellular VEGF-A level, a paradigmatic model of post-transcriptional regulation.

Multiscale computational models of complex biological systems

Integration of data across spatial, temporal, and functional scales is a primary focus of biomedical engineering efforts. The advent of powerful computing platforms, coupled with quantitative data from high-throughput experimental methodologies, has allowed multiscale modeling to expand as a means to more comprehensively investigate biological phenomena in experimentally relevant ways. This review aims to highlight recently published multiscale models of biological systems, using their successes to propose the best practices for future model development. We demonstrate that coupling continuous and discrete systems best captures biological information across spatial scales by selecting modeling techniques that are suited to the task. Further, we suggest how to leverage these multiscale models to gain insight into biological systems using quantitative biomedical engineering methods to analyze data in nonintuitive ways. These topics are discussed with a focus on the future of the field, current challenges encountered, and opportunities yet to be realized.

MOSAIC: a multiscale model of osteogenesis and sprouting angiogenesis with lateral inhibition of endothelial cells

The healing of a fracture depends largely on the development of a new blood vessel network (angiogenesis) in the callus. During angiogenesis tip cells lead the developing sprout in response to extracellular signals, amongst which vascular endothelial growth factor (VEGF) is critical. In order to ensure a correct development of the vasculature, the balance between stalk and tip cell phenotypes must be tightly controlled, which is primarily achieved by the Dll4-Notch1 signaling pathway. This study presents a novel multiscale model of osteogenesis and sprouting angiogenesis, incorporating lateral inhibition of endothelial cells (further denoted MOSAIC model) through Dll4-Notch1 signaling, and applies it to fracture healing. The MOSAIC model correctly predicted the bone regeneration process and recapitulated many experimentally observed aspects of tip cell selection: the salt and pepper pattern seen for cell fates, an increased tip cell density due to the loss of Dll4 and an excessive number of tip cells in high VEGF environments. When VEGF concentration was even further increased, the MOSAIC model predicted the absence of a vascular network and fracture healing, thereby leading to a non-union, which is a direct consequence of the mutual inhibition of neighboring cells through Dll4-Notch1 signaling. This result was not retrieved for a more phenomenological model that only considers extracellular signals for tip cell migration, which illustrates the importance of implementing the actual signaling pathway rather than phenomenological rules. Finally, the MOSAIC model demonstrated the importance of a proper criterion for tip cell selection and the need for experimental data to further explore this. In conclusion, this study demonstrates that the MOSAIC model creates enhanced capabilities for investigating the influence of molecular mechanisms on angiogenesis and its relation to bone formation in a more mechanistic way and across different time and spatial scales.

Vascular endothelial growth factor-delivery systems for cardiac repair: an overview

Since the discovery of the Vascular Endothelial Growth Factor (VEGF) and its leading role in the angiogenic process, this has been seen as a promising molecule for promoting neovascularization in the infarcted heart. However, even though several clinical trials were initiated, no therapeutic effects were observed, due in part to the short half life of this factor when administered directly to the tissue. In this context, drug delivery systems appear to offer a promising strategy to overcome limitations in clinical trials of VEGF.The aim of this paper is to review the principal drug delivery systems that have been developed to administer VEGF in cardiovascular disease. Studies published in the last 5 years are reviewed and the main features of these systems are explained. The tissue engineering concept is introduced as a therapeutic alternative that holds promise for the near future.

Integration of experimental and computational approaches to sprouting angiogenesis

PURPOSE OF REVIEW

We summarize recent experimental and computational studies that investigate molecular and cellular mechanisms of sprouting angiogenesis. We discuss how experimental tools have unveiled new opportunities for computational modeling by providing detailed phenomenological descriptions and conceptual models of cell-level behaviors underpinned by high-quality molecular data. Using recent examples, we show how new understanding results from bridging computational and experimental approaches.

RECENT FINDINGS

Experimental data extends beyond the tip cell vs. stalk cell paradigm, and involves numerous molecular inputs such as vascular endothelial growth factor and Notch. This data is being used to generate and validate computational models, which can then be used to predict the results of hypothetical experiments that are difficult to perform in the laboratory, and to generate new hypotheses that account for system-wide interactions. As a result of this integration, descriptions of critical gradients of growth factor-receptor complexes have been generated, and new modulators of cell behavior have been described.

SUMMARY

We suggest that the recent emphasis on the different stages of sprouting angiogenesis, and integration of experimental and computational approaches, should provide a way to manage the complexity of this process and help identify new regulatory paradigms and therapeutic targets.

A two-compartment model of VEGF distribution in the mouse

Vascular endothelial growth factor (VEGF) is a key regulator of angiogenesis – the growth of new microvessels from existing microvasculature. Angiogenesis is a complex process involving numerous molecular species, and to better understand it, a systems biology approach is necessary. In vivo preclinical experiments in the area of angiogenesis are typically performed in mouse models; this includes drug development targeting VEGF. Thus, to quantitatively interpret such experimental results, a computational model of VEGF distribution in the mouse can be beneficial. In this paper, we present an in silico model of VEGF distribution in mice, determine model parameters from existing experimental data, conduct sensitivity analysis, and test the validity of the model.

The multiscale model is comprised of two compartments: blood and tissue. The model accounts for interactions between two major VEGF isoforms (VEGF₁₂₀ and VEGF₁₆₄) and their endothelial cell receptors VEGFR-1, VEGFR-2, and co-receptor neuropilin-1. Neuropilin-1 is also expressed on the surface of parenchymal cells. The model includes transcapillary macromolecular permeability, lymphatic transport, and macromolecular plasma clearance. Simulations predict that the concentration of unbound VEGF in the tissue is approximately 50-fold greater than in the blood. These concentrations are highly dependent on the VEGF secretion rate. Parameter estimation was performed to fit the simulation results to available experimental data, and permitted the estimation of VEGF secretion rate in healthy tissue, which is difficult to measure experimentally. The model can provide quantitative interpretation of preclinical animal data and may be used in conjunction with experimental studies in the development of pro- and anti-angiogenic agents. The model approximates the normal tissue as skeletal muscle and includes endothelial cells to represent the vasculature. As the VEGF system becomes better characterized in other tissues and cell types, the model can be expanded to include additional compartments and vascular elements.