Computational modeling of interacting VEGF and soluble VEGF receptor concentration gradients

Angiogenic Microvascular Wall Shear Stress Patterns Revealed Through Three-dimensional Red Blood Cell Resolved Modeling

The wall shear stress (WSS) exerted by blood flowing through microvascular capillaries is an established driver of new blood vessel growth, or angiogenesis. Such adaptations are central to many physiological processes in both health and disease, yet three-dimensional (3D) WSS characteristics in real angiogenic microvascular networks are largely unknown. This marks a major knowledge gap because angiogenesis, naturally, is a 3D process. To advance current understanding, we model 3D red blood cells (RBCs) flowing through rat angiogenic microvascular networks using state-of-the-art simulation. The high-resolution fluid dynamics reveal 3D WSS patterns occurring at sub-endothelial cell (EC) scales that derive from distinct angiogenic morphologies, including microvascular loops and vessel tortuosity. We identify the existence of WSS hot and cold spots caused by angiogenic surface shapes and RBCs, and notably enhancement of low WSS regions by RBCs. Spatiotemporal characteristics further reveal how fluctuations follow timescales of RBC "footprints." Altogether, this work provides a new conceptual framework for understanding how shear stress might regulate EC dynamics in vivo.

The Endothelial Glycocalyx as a Target of Excess Soluble Fms-like Tyrosine Kinase-1

Soluble fms-like tyrosine kinase-1 (sFlt-1) is a secreted protein that binds heparan sulfate expressed on the endothelial glycocalyx (eGC). In this paper we analyze how excess sFlt-1 causes conformational changes in the eGC, leading to monocyte adhesion, a key event triggering vascular dysfunction. In vitro exposure of primary human umbilical vein endothelial cells to excess sFlt-1 decreased eGC height and increased stiffness as determined by atomic force microscopy (AFM). Yet, structural loss of the eGC components was not observed, as indicated by Ulex europaeus agglutinin I and wheat germ agglutinin staining. Moreover, the conformation observed under excess sFlt-1, a collapsed eGC, is flat and stiff with unchanged coverage and sustained content. Functionally, this conformation increased the endothelial adhesiveness to THP-1 monocytes by about 35%. Heparin blocked all these effects, but the vascular endothelial growth factor did not. In vivo administration of sFlt-1 in mice also resulted in the collapse of the eGC in isolated aorta analyzed ex vivo by AFM. Our findings show that excess sFlt-1 causes the collapse of the eGC and favors leukocyte adhesion. This study provides an additional mechanism of action by which sFlt-1 may cause endothelial dysfunction and injury.

The Cardiorespiratory Network in Healthy First-Degree Relatives of Schizophrenic Patients

Impaired heart rate- and respiratory regulatory processes as a sign of an autonomic dysfunction seems to be obviously present in patients suffering from schizophrenia. Since the linear and non-linear couplings within the cardiorespiratory system with respiration as an important homeostatic control mechanism are only partially investigated so far for those subjects, we aimed to characterize instantaneous cardiorespiratory couplings by quantifying the casual interaction between heart rate (HR) and respiration (RESP). Therefore, we investigated causal linear and non-linear cardiorespiratory couplings of 23 patients suffering from schizophrenia (SZO), 20 healthy first-degree relatives (REL) and 23 healthy subjects, who were age-gender matched (CON). From all participants' heart rate (HR) and respirations (respiratory frequency, RESP) were investigated for 30 min under resting conditions. The results revealed highly significant increased HR, reduced HR variability, increased respiration rates and impaired cardiorespiratory couplings in SZO in comparison to CON. SZO were revealed bidirectional couplings, with respiration as the driver (RESP → HR), and with weaker linear and non-linear coupling strengths when RESP influencing HR (RESP → HR) and with stronger linear and non-linear coupling strengths when HR influencing RESP (HR → RESP). For REL we found only significant increased HR and only slightly reduced cardiorespiratory couplings compared to CON. These findings clearly pointing to an underlying disease-inherent genetic component of the cardiac system for SZO and REL, and those respiratory alterations are only clearly present in SZO seem to be connected to their mental emotional states.

Agent-based models of inflammation in translational systems biology: A decade later

Agent-based modeling is a rule-based, discrete-event, and spatially explicit computational modeling method that employs computational objects that instantiate the rules and interactions among the individual components ("agents") of system. Agent-based modeling is well suited to translating into a computational model the knowledge generated from basic science research, particularly with respect to translating across scales the mechanisms of cellular behavior into aggregated cell population dynamics manifesting at the tissue and organ level. This capacity has made agent-based modeling an integral method in translational systems biology (TSB), an approach that uses multiscale dynamic computational modeling to explicitly represent disease processes in a clinically relevant fashion. The initial work in the early 2000s using agent-based models (ABMs) in TSB focused on examining acute inflammation and its intersection with wound healing; the decade since has seen vast growth in both the application of agent-based modeling to a wide array of disease processes as well as methodological advancements in the use and analysis of ABM. This report presents an update on an earlier review of ABMs in TSB and presents examples of exciting progress in the modeling of various organs and diseases that involve inflammation. This review also describes developments that integrate the use of ABMs with cutting-edge technologies such as high-performance computing, machine learning, and artificial intelligence, with a view toward the future integration of these methodologies. This article is categorized under: Translational, Genomic, and Systems Medicine > Translational Medicine Models of Systems Properties and Processes > Mechanistic Models Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models Models of Systems Properties and Processes > Organismal Models.

Blood Vessel Patterning on Retinal Astrocytes Requires Endothelial Flt-1 (VEGFR-1)

Feedback mechanisms are critical components of many pro-angiogenic signaling pathways that keep vessel growth within a functional range. The Vascular Endothelial Growth Factor-A (VEGF-A) pathway utilizes the decoy VEGF-A receptor Flt-1 to provide negative feedback regulation of VEGF-A signaling. In this study, we investigated how the genetic loss of flt-1 differentially affects the branching complexity of vascular networks in tissues despite similar effects on endothelial sprouting. We selectively ablated flt-1 in the post-natal retina and found that maximum induction of flt-1 loss resulted in alterations in endothelial sprouting and filopodial extension, ultimately yielding hyper-branched networks in the absence of changes in retinal astrocyte architecture. The mosaic deletion of flt-1 revealed that sprouting endothelial cells flanked by flt-1^−/− regions of vasculature more extensively associated with underlying astrocytes and exhibited aberrant sprouting, independent of the tip cell genotype. Overall, our data support a model in which tissue patterning features, such as retinal astrocytes, integrate with flt-1-regulated angiogenic molecular and cellular mechanisms to yield optimal vessel patterning for a given tissue.

Computational Modeling to Quantify the Contributions of VEGFR1, VEGFR2, and Lateral Inhibition in Sprouting Angiogenesis

Sprouting angiogenesis is a necessary process in regeneration and development as well as in tumorigenesis. VEGF-A is the main pro-angiogenic chemoattractant and it can bind to the decoy receptor VEGFR1 or to VEGFR2 to induce sprouting. Active sprout cells express Dll4, which binds to Notch1 on neighboring cells, in turn inhibiting VEGFR2 expression. It is known that the balance between VEGFR2 and VEGFR1 determines tip selection and network architecture, however the quantitative interrelationship of the receptors and their interrelated balances, also with relation to Dll4-Notch1 signaling, remains yet largely unknown. Here, we present an agent-based computer model of sprouting angiogenesis, integrating VEGFR1 and VEGFR2 in a detailed model of cellular signaling. Our model reproduces experimental data on VEGFR1 knockout. We show that soluble VEGFR1 improves the efficiency of angiogenesis by directing sprouts away from existing cells over a wide range of parameters. Our analysis unravels the relevance of the stability of the active notch intracellular domain as a dominating hub in this regulatory network. Our analysis quantitatively dissects the regulatory interactions in sprouting angiogenesis. Because we use a detailed model of intracellular signaling, the results of our analysis are directly linked to biological entities. We provide our computational model and simulation engine for integration in complementary modeling approaches.

Can VEGFC Form Turing Patterns in the Zebrafish Embryo?

This paper is concerned with a late stage of lymphangiogenesis in the trunk of the zebrafish embryo. At 48 hours post-fertilisation (HPF), a pool of parachordal lymphangioblasts (PLs) lies in the horizontal myoseptum. Between 48 and 168 HPF, the PLs spread from the horizontal myoseptum to form the thoracic duct, dorsal longitudinal lymphatic vessel, and parachordal lymphatic vessel. This paper deals with the potential of vascular endothelial growth factor C (VEGFC) to guide the differentiation of PLs into the mature lymphatic endothelial cells that form the vessels. We built a mathematical model to describe the biochemical interactions between VEGFC, collagen I, and matrix metalloproteinase 2 (MMP2). We also carried out a linear stability analysis of the model and computer simulations of VEGFC patterning. The results suggest that VEGFC can form Turing patterns due to its relations with MMP2 and collagen I, but the zebrafish embryo needs a separate control mechanism to create the right physiological conditions. Furthermore, this control mechanism must ensure that the VEGFC patterns are useful for lymphangiogenesis: stationary, steep gradients, and reasonably fast forming. Generally, the combination of a patterning species, a matrix protein, and a remodelling species is a new patterning mechanism.

Biology of preeclampsia: Combined actions of angiogenic factors, their receptors and placental proteins

Although massive efforts have been undertaken to elucidate the etiology of the pregnancy syndrome preeclampsia, its developmental origin remains a mystery. Most efforts of the last decade have focused on biomarkers to predict and/or diagnose preeclampsia, including the anti-angiogenic factor sFlt-1 (soluble fms-like tyrosin kinase-1), the angiogenic factor PGF (placental growth factor) and PP13 (placental protein 13). The origins of these marker proteins are still under debate, and so far their actions have only been describe separate from each other. This study will focus on the origins and actions of all three markers during pregnancy and outside pregnancy and will describe a scenario where all three markers act synergistically to rescue the mother from the deleterious effects of the debris that is released from the placenta during preeclampsia. This more holistic approach may open new avenues to think about maternal-fetal interactions and putative therapies.

Assessing the Minimum Time-Period of Normoxic Preincubation for Stable Adipose Stromal Cell-Derived Vascular Networks

Introduction

Pre-vascularization of tissue engineered grafts is a promising strategy to facilitate their improved viability following in vivo implantation. In this process, endothelial cells (ECs) form capillary-like networks that can anastomose with host vasculature. Adipose-derived stromal cells (ASCs) are a commonly used cell population for tissue engineering and contain a subpopulation of ECs capable of assembling into robust vascular networks and anastomosing with the host. However, their initial vascular assembly is significantly impaired in hypoxic conditions (2% O₂). In this study, we explored the minimum period of normoxic (20% O₂) pre-treatment required to enable the formation of stable vascular networks.

Methods

ASC-derived vascular structures were allowed to preassemble in fibrin hydrogels in normoxia for 0, 2, 4, or 6 days and then transplanted into hypoxic environments for 6 days. Total vascular length, pericyte coverage, cell proliferation, apoptosis rates, and ECM production was assessed.

Results

Vascular assembly increased with time over the 6 days of culture. We found that 4 days was the minimum period of time required for stable vascular assembly. We compared the major differences in cell behavior and network structure at Days 2 and 4. Neither proliferation nor apoptosis differed, however, the Day 4 time-point was associated with a significant increase in pericyte coverage (46.1 ± 2.6%) compared to Day 2 (24.3 ± 5.3%).

Conclusions

These data suggest oxygen tension may be a mediator of EC-pericyte interactions during vascular assembly. Pre-vascularization strategies should incorporate a normoxic period of to enable successful vascular formation and development.

Agent-based computational model of retinal angiogenesis simulates microvascular network morphology as a function of pericyte coverage

OBJECTIVE

Define a role for perivascular cells during developmental retinal angiogenesis in the context of EC Notch1-DLL4 signaling at the multicellular network level.

METHODS

The retinal vasculature is highly sensitive to growth factor-mediated intercellular signaling. Although EC signaling has been explored in detail, it remains unclear how PC function to modulate these signals that lead to a diverse set of vascular network patterns in health and disease. We have developed an ABM of retinal angiogenesis that incorporates both ECs and PCs to investigate the formation of vascular network patterns as a function of pericyte coverage. We use our model to test the hypothesis that PC modulate Notch1-DLL4 signaling in endothelial cell-endothelial cell interactions.

RESULTS

Agent-based model (ABM) simulations that include PCs more accurately predict experimentally observed vascular network morphologies than simulations that lack PCs, suggesting that PCs may influence sprouting behaviors through physical blockade of endothelial intercellular connections.

CONCLUSIONS

This study supports a role for PCs as a physical buffer to signal propagation during vascular network formation-a barrier that may be important for generating healthy microvascular network patterns.

A Mathematical Model of Lymphangiogenesis in a Zebrafish Embryo

The lymphatic system of a vertebrate is important in health and diseases. We propose a novel mathematical model to elucidate the lymphangiogenic processes in zebrafish embryos. Specifically, we are interested in the period when lymphatic endothelial cells (LECs) exit the posterior cardinal vein and migrate to the horizontal myoseptum of a zebrafish embryo. We wonder whether vascular endothelial growth factor C (VEGFC) is a morphogen and a chemotactic factor for these LECs. The model considers the interstitial flow driving convection, the reactive transport of VEGFC, and the changing dynamics of the extracellular matrix in the embryo. Simulations of the model illustrate that VEGFC behaves very differently in diffusion and convection-dominant scenarios. In the former case, it must bind to the matrix to establish a functional morphogen gradient. In the latter case, the opposite is true and the pressure field is the key determinant of what VEGFC may do to the LECs. Degradation of collagen I, a matrix component, by matrix metallopeptidase 2 controls the spatiotemporal dynamics of VEGFC. It controls whether diffusion or convection is dominant in the embryo; it can create channels of abundant VEGFC and scarce collagen I to facilitate lymphangiogenesis; when collagen I is insufficient, VEGFC cannot influence the LECs at all. We predict that VEGFC is a morphogen for the migrating LECs, but it is not a chemotactic factor for them.

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.

Computational Modeling of Muscle Regeneration and Adaptation to Advance Muscle Tissue Regeneration Strategies

Skeletal muscle has an exceptional ability to regenerate and adapt following injury. Tissue engineering approaches (e.g. cell therapy, scaffolds, and pharmaceutics) aimed at enhancing or promoting muscle regeneration from severe injuries are a promising and active field of research. Computational models are beginning to advance the field by providing insight into regeneration mechanisms and therapies. In this paper, we summarize the contributions computational models have made to understanding muscle remodeling and the functional implications thereof. Next, we describe a new agent-based computational model of skeletal muscle inflammation and regeneration following acute muscle injury. Our computational model simulates the recruitment and cellular behaviors of key inflammatory cells (e.g. neutrophils and M1 and M2 macrophages) and their interactions with native muscle cells (muscle fibers, satellite stem cells, and fibroblasts) that result in the clearance of necrotic tissue and muscle fiber regeneration. We demonstrate the ability of the model to track key regeneration metrics during both unencumbered regeneration and in the case of impaired macrophage function. We also use the model to simulate regeneration enhancement when muscle is primed with inflammatory cells prior to injury, which is a putative therapeutic intervention that has not yet been investigated experimentally. Computational modeling of muscle regeneration, pursued in combination with experimental analyses, provides a quantitative framework for evaluating and predicting muscle regeneration and enables the rational design of therapeutic strategies for muscle recovery.

Flt-1 (VEGFR-1) coordinates discrete stages of blood vessel formation

AIMS

In developing blood vessel networks, the overall level of vessel branching often correlates with angiogenic sprout initiations, but in some pathological situations, increased sprout initiations paradoxically lead to reduced vessel branching and impaired vascular function. We examine the hypothesis that defects in the discrete stages of angiogenesis can uniquely contribute to vessel branching outcomes.

METHODS AND RESULTS

Time-lapse movies of mammalian blood vessel development were used to define and quantify the dynamics of angiogenic sprouting. We characterized the formation of new functional conduits by classifying discrete sequential stages-sprout initiation, extension, connection, and stability-that are differentially affected by manipulation of vascular endothelial growth factor-A (VEGF-A) signalling via genetic loss of the receptor flt-1 (vegfr1). In mouse embryonic stem cell-derived vessels genetically lacking flt-1, overall branching is significantly decreased while sprout initiations are significantly increased. Flt-1(-/-) mutant sprouts are less likely to retract, and they form increased numbers of connections with other vessels. However, loss of flt-1 also leads to vessel collapse, which reduces the number of new stable conduits. Computational simulations predict that loss of flt-1 results in ectopic Flk-1 signalling in connecting sprouts post-fusion, causing protrusion of cell processes into avascular gaps and collapse of branches. Thus, defects in stabilization of new vessel connections offset increased sprout initiations and connectivity in flt-1(-/-) vascular networks, with an overall outcome of reduced numbers of new conduits.

CONCLUSIONS

These results show that VEGF-A signalling has stage-specific effects on vascular morphogenesis, and that understanding these effects on dynamic stages of angiogenesis and how they integrate to expand a vessel network may suggest new therapeutic strategies.

Multi-scale Modeling of the Cardiovascular System: Disease Development, Progression, and Clinical Intervention

Cardiovascular diseases (CVDs) are the leading cause of death in the western world. With the current development of clinical diagnostics to more accurately measure the extent and specifics of CVDs, a laudable goal is a better understanding of the structure-function relation in the cardiovascular system. Much of this fundamental understanding comes from the development and study of models that integrate biology, medicine, imaging, and biomechanics. Information from these models provides guidance for developing diagnostics, and implementation of these diagnostics to the clinical setting, in turn, provides data for refining the models. In this review, we introduce multi-scale and multi-physical models for understanding disease development, progression, and designing clinical interventions. We begin with multi-scale models of cardiac electrophysiology and mechanics for diagnosis, clinical decision support, personalized and precision medicine in cardiology with examples in arrhythmia and heart failure. We then introduce computational models of vasculature mechanics and associated mechanical forces for understanding vascular disease progression, designing clinical interventions, and elucidating mechanisms that underlie diverse vascular conditions. We conclude with a discussion of barriers that must be overcome to provide enhanced insights, predictions, and decisions in pre-clinical and clinical applications.

Elevated Soluble Fms-Like Tyrosine Kinase-1 and Placental-Like Growth Factor Levels Are Associated With Development and Mortality Risk in Heart Failure

BACKGROUND

Vascular endothelial dysfunction may play an important role in the progression of heart failure (HF). We hypothesize that elevated levels of vascular markers, placental-like growth factor, and soluble Fms-like tyrosine kinase-1 (sFlt-1) are associated with adverse outcomes in patients with HF. We also assessed possible triggers of sFlt-1 elevation in animal HF models.

METHODS AND RESULTS

We measured plasma placental-like growth factor and sFlt-1 in 791 HF patients undergoing elective coronary angiogram. Median (interquartile range) placental-like growth factor and sFlt-1 levels were 24 (20-29) and 382 (277-953) pg/mL, respectively. After 5 years of follow-up, and after using receiver operator characteristic curves to determine optimal cutoffs, high levels of sFlt-1 (≥ 280 pg/mL; adjusted hazard ratio, 1.47; 95% confidence interval, 1.03-2.09; P=0.035) but not placental-like growth factor (≥ 25 pg/mL; adjusted hazard ratio, 1.26; 95% confidence interval, 0.94-1.71, P=0.12) were associated with adverse cardiovascular outcomes. In addition, significant elevation of sFlt-1 levels was observed in left anterior descending artery ligation and transverse aortic constriction HF mouse models after 4 and 8 weeks of follow-up, suggesting vascular stress and ischemia as triggers for sFlt-1 elevation in HF.

CONCLUSIONS

Circulating sFlt-1 is generated as a result of myocardial injury and subsequent HF development. Elevated levels of sFlt-1 are associated with adverse outcomes in stable patients with HF.

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.

Collective motion of cells: from experiments to models

Swarming or collective motion of living entities is one of the most common and spectacular manifestations of living systems that have been extensively studied in recent years. A number of general principles have been established. The interactions at the level of cells are quite different from those among individual animals, therefore the study of collective motion of cells is likely to reveal some specific important features which we plan to overview in this paper. In addition to presenting the most appealing results from the quickly growing related literature we also deliver a critical discussion of the emerging picture and summarize our present understanding of collective motion at the cellular level. Collective motion of cells plays an essential role in a number of experimental and real-life situations. In most cases the coordinated motion is a helpful aspect of the given phenomenon and results in making a related process more efficient (e.g., embryogenesis or wound healing), while in the case of tumor cell invasion it appears to speed up the progression of the disease. In these mechanisms cells both have to be motile and adhere to one another, the adherence feature being the most specific to this sort of collective behavior. One of the central aims of this review is to present the related experimental observations and treat them in light of a few basic computational models so as to make an interpretation of the phenomena at a quantitative level as well.

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.

A double-edged sword: the role of VEGF in wound repair and chemoattraction of opportunist pathogens

Wound healing is a complex process essential to repairing damaged tissues and preventing infection. Skin is the first line of defense, a chief physical barrier to microbe entry. Wound healing is a physical rebuilding process, but at the same time it is an inflammatory event. In turn, molecules for wound repair are secreted by fibroblasts and others present at the wound site. Vascular endothelial growth factor (VEGF) is a critical cytokine that exhibits chemoattractant properties, recruiting other immune cells to the site. Although generally beneficial, VEGF may also act as a chemoattractant for invading microorganisms, such as Pseudomonas aeruginosa.P. aeruginosa is problematic during wound infection due to its propensity to form biofilms and exhibit heightened antimicrobial resistance. Here, we explored the influence of VEGF gradients (in a microfluidic device wound model) on the motility and chemotactic properties of P. aeruginosa. At lower concentrations, VEGF had little effect on motility, but as the maximal concentration within the gradient increased, P. aeruginosa cells exhibited directed movement along the gradient. Our data provide evidence that while beneficial, VEGF, in excess, may aid colonization by P. aeruginosa. This highlights the necessity for the efficient resolution of inflammation. Understanding the dynamics of wound colonization may lead to new/enhanced therapeutics to hasten recovery.

Do endothelial cells dream of eclectic shape?

Endothelial cells (ECs) exhibit dramatic plasticity of form at the single- and collective-cell level during new vessel growth, adult vascular homeostasis, and pathology. Understanding how, when, and why individual ECs coordinate decisions to change shape, in relation to the myriad of dynamic environmental signals, is key to understanding normal and pathological blood vessel behavior. However, this is a complex spatial and temporal problem. In this review we show that the multidisciplinary field of Adaptive Systems offers a refreshing perspective, common biological language, and straightforward toolkit that cell biologists can use to untangle the complexity of dynamic, morphogenetic systems.

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.

Identical effects of VEGF and serum-deprivation on phenotype and function of adipose-derived stromal cells from healthy donors and patients with ischemic heart disease

BACKGROUND

Adipose-derived stromal cells (ASCs) stimulated with vascular endothelial growth factor (VEGF) and serum-deprived, are applied in the first in-man double-blind placebo-controlled MyStromalCell Trial, as a novel therapeutic option for treatment of ischemic heart disease (IHD). This in vitro study explored the effect of VEGF and serum deprivation on endothelial differentiation capacity of ASCs from healthy donors and IHD patients.

METHODS

ASCs stimulated with rhVEGF(A165) in serum-deprived medium for one to three weeks were compared with ASCs in serum-deprived (2% fetal bovine serum) or complete medium (10% fetal bovine serum). Expression of VEGF receptors, endothelial and stem cell markers was measured using qPCR, flow cytometry and immunocytochemistry. In vitro tube formation and proliferation was also measured.

RESULTS

ASCs from VEGF-stimulated and serum-deprived medium significantly increased transcription of transcription factor FOXF1, endothelial marker vWF and receptor VEGFR1 compared with ASCs from complete medium. ASCs maintained stem cell characteristics in all conditions. Tube formation of ASCs occurred in VEGF-stimulated and serum-deprived medium. The only difference between healthy and patient ASCs was a variation in proliferation rate.

CONCLUSIONS

ASCs from IHD patients and healthy donors proved equally inclined to differentiate in endothelial direction by serum-deprivation, however with no visible additive effect of VEGF stimulation. The treatment did not result in complete endothelial differentiation, but priming towards endothelial lineage.

Endothelial cell motility, coordination and pattern formation during vasculogenesis

How vascular networks assemble is a fundamental problem of developmental biology that also has medical importance. To explain the organizational principles behind vascular patterning, we must understand how can tissue level structures be controlled through cell behavior patterns like motility and adhesion that, in turn, are determined by biochemical signal transduction processes? We discuss the various ideas that have been proposed as mechanisms for vascular network assembly: cell motility guided by extracellular matrix alignment (contact guidance), chemotaxis guided by paracrine and autocrine morphogens, and multicellular sprouting guided by cell-cell contacts. All of these processes yield emergent patterns, thus endothelial cells can form an interconnected structure autonomously, without guidance from an external pre-pattern.

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.

Overcoming resistance to antiangiogenic therapies

The concept of targeting new blood vessel formation, or angiogenesis, in tumors is an important advancement in cancer therapy, resulting, in part, from the development of such biologic agents as bevacizumab, a monoclonal antibody directed against vascular endothelial growth factor (VEGF)-A. The rationale for antiangiogenic therapy is based on the hypothesis that if tumors are limited in their capacity to obtain a new blood supply, so too is their capacity for growth and metastasis. Additional evidence suggests that pruning and/or "normalization" of irregular tumor vasculature and reduction of hypoxia may facilitate greater access of cytotoxic chemotherapy (CT) to the tumor. Indeed, for metastatic colorectal cancer, bevacizumab in combination with established CT regimens has efficacy superior to that of CT alone. Despite ~2-month longer progression-free and overall survival times than with CT alone, patients still progress, possibly because of alternative angiogenic "escape" pathways that emerge independent of VEGF-A, or are driven by hypoxic stress on the tumor. Other VEGF family members may contribute to resistance, and many factors that contribute to the regulation of tumor angiogenesis function as part of a complex network, existing in different concentrations and spatiotemporal gradients and producing a wide range of biologic responses. Integrating these concepts into the design and evaluation of new antiangiogenic therapies may help overcome resistance mechanisms and allow for greater efficacy over longer treatment periods.

Pattern formation during vasculogenesis

Vasculogenesis, the assembly of the first vascular network, is an intriguing developmental process that yields the first functional organ system of the embryo. In addition to being a fundamental part of embryonic development, vasculogenic processes also have medical importance. To explain the organizational principles behind vascular patterning, we must understand how morphogenesis of tissue level structures can be controlled through cell behavior patterns that, in turn, are determined by biochemical signal transduction processes. Mathematical analyses and computer simulations can help conceptualize how to bridge organizational levels and thus help in evaluating hypotheses regarding the formation of vascular networks. Here, we discuss the ideas that have been proposed to explain the formation of the first vascular pattern: cell motility guided by extracellular matrix alignment (contact guidance), chemotaxis guided by paracrine and autocrine morphogens, and sprouting guided by cell-cell contacts.

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.