Computational model of VEGFR2 pathway to ERK activation and modulation through receptor trafficking

CD133-containing microvesicles promote colorectal cancer progression by inducing tumor angiogenesis

Angiogenesis is an indispensable mechanism in cancer progression, as cancer cells need to establish blood vessels to supply oxygen and nutrients. Extracellular vesicles (EVs) derived from cancer cells act as messengers in the tumor microenvironment and induce resistance to anti-angiogenic cancer treatment. EVs can be classified into two categories: exosomes and microvesicles (MVs). Although exosomes are involved in angiogenesis, the role of MVs in angiogenesis and cancer progression remains unclear. CD133 plays a key role in MV formation and oncoprotein trafficking. In this study, we investigated the role of CD133-containing MVs derived from colorectal cancer (CRC) in angiogenesis and cancer progression. CRC-derived MVs were incorporated into endothelial cells and increased the mesh area and tube length of endothelial cells. CD133-containing MVs also stimulate vessel sprouting in endothelial cell spheroids and mouse thoracic aortas. However, MVs derived from CD133-knockdown CRC cells exerted a limited effect on tube formation and vessel sprouting. CD133-containing MVs induced angiogenesis through p38 activation and angiogenesis induced by CD133-containing MVs was insensitive to the anti-vascular endothelial growth factor antibody bevacizumab. Survival analysis revealed that high expression level of CD133 correlated with poor prognosis in patients with metastatic CRC. These findings suggest that CD133-containing MVs act as key regulators of angiogenesis and are related to the prognosis of CRC patients.

Trafficking dynamics of VEGFR1, VEGFR2, and NRP1 in human endothelial cells

The vascular endothelial growth factor (VEGF) family of cytokines are key drivers of blood vessel growth and remodeling. These ligands act via multiple VEGF receptors (VEGFR) and co-receptors such as Neuropilin (NRP) expressed on endothelial cells. These membrane-associated receptors are not solely expressed on the cell surface, they move between the surface and intracellular locations, where they can function differently. The location of the receptor alters its ability to 'see' (access and bind to) its ligands, which regulates receptor activation; location also alters receptor exposure to subcellularly localized phosphatases, which regulates its deactivation. Thus, receptors in different subcellular locations initiate different signaling, both in terms of quantity and quality. Similarly, the local levels of co-expression of other receptors alters competition for ligands. Subcellular localization is controlled by intracellular trafficking processes, which thus control VEGFR activity; therefore, to understand VEGFR activity, we must understand receptor trafficking. Here, for the first time, we simultaneously quantify the trafficking of VEGFR1, VEGFR2, and NRP1 on the same cells-specifically human umbilical vein endothelial cells (HUVECs). We build a computational model describing the expression, interaction, and trafficking of these receptors, and use it to simulate cell culture experiments. We use new quantitative experimental data to parameterize the model, which then provides mechanistic insight into the trafficking and localization of this receptor network. We show that VEGFR2 and NRP1 trafficking is not the same on HUVECs as on non-human ECs; and we show that VEGFR1 trafficking is not the same as VEGFR2 trafficking, but rather is faster in both internalization and recycling. As a consequence, the VEGF receptors are not evenly distributed between the cell surface and intracellular locations, with a very low percentage of VEGFR1 being on the cell surface, and high levels of NRP1 on the cell surface. Our findings have implications both for the sensing of extracellular ligands and for the composition of signaling complexes at the cell surface versus inside the cell.

Regulation of cells of the arterial wall by hypoxia and its role in the development of atherosclerosis

The cell's response to hypoxia depends on stabilization of the hypoxia-inducible factor 1 complex and transactivation of nuclear factor kappa-B (NF-κB). HIF target gene transcription in cells resident to atherosclerotic lesions adjoins a complex interplay of cytokines and mediators of inflammation affecting cholesterol uptake, migration, and inflammation. Maladaptive activation of the HIF-pathway and transactivation of nuclear factor kappa-B causes monocytes to invade early atherosclerotic lesions, maintaining inflammation and aggravating a low-oxygen environment. Meanwhile HIF-dependent upregulation of the ATP-binding cassette transporter ABCA1 causes attenuation of cholesterol efflux and ultimately macrophages becoming foam cells. Hypoxia facilitates neovascularization by upregulation of vascular endothelial growth factor (VEGF) secreted by endothelial cells and vascular smooth muscle cells lining the arterial wall destabilizing the plaque. HIF-knockout animal models and inhibitor studies were able to show beneficial effects on atherogenesis by counteracting the HIF-pathway in the cell wall. In this review the authors elaborate on the up-to-date literature on regulation of cells of the arterial wall through activation of HIF-1α and its effect on atherosclerotic plaque formation.

Brain arteriovenous malformation in hereditary hemorrhagic telangiectasia: Recent advances in cellular and molecular mechanisms

Hereditary hemorrhagic telangiectasia (HHT) is a genetic disorder characterized by vessel dilatation, such as telangiectasia in skin and mucosa and arteriovenous malformations (AVM) in internal organs such as the gastrointestinal tract, lungs, and brain. AVMs are fragile and tortuous vascular anomalies that directly connect arteries and veins, bypassing healthy capillaries. Mutations in transforming growth factor β (TGFβ) signaling pathway components, such as ENG (ENDOGLIN), ACVRL1 (ALK1), and SMAD4 (SMAD4) genes, account for most of HHT cases. 10-20% of HHT patients develop brain AVMs (bAVMs), which can lead to vessel wall rupture and intracranial hemorrhages. Though the main mutations are known, mechanisms leading to AVM formation are unclear, partially due to lack of animal models. Recent mouse models allowed significant advances in our understanding of AVMs. Endothelial-specific deletion of either Acvrl1, Eng or Smad4 is sufficient to induce AVMs, identifying endothelial cells (ECs) as primary targets of BMP signaling to promote vascular integrity. Loss of ALK1/ENG/SMAD4 signaling is associated with NOTCH signaling defects and abnormal arteriovenous EC differentiation. Moreover, cumulative evidence suggests that AVMs originate from venous ECs with defective flow-migration coupling and excessive proliferation. Mutant ECs show an increase of PI3K/AKT signaling and inhibitors of this signaling pathway rescue AVMs in HHT mouse models, revealing new therapeutic avenues. In this review, we will summarize recent advances and current knowledge of mechanisms controlling the pathogenesis of bAVMs, and discuss unresolved questions.

Non-Muscle MLCK Contributes to Endothelial Cell Hyper-Proliferation through the ERK Pathway as a Mechanism for Vascular Remodeling in Pulmonary Hypertension

Pulmonary arterial hypertension (PAH) is characterized by endothelial dysfunction, uncontrolled proliferation and migration of pulmonary arterial endothelial cells leading to increased pulmonary vascular resistance resulting in great morbidity and poor survival. Bone morphogenetic protein receptor II (BMPR2) plays an important role in the pathogenesis of PAH as the most common genetic mutation. Non-muscle myosin light chain kinase (nmMLCK) is an essential component of the cellular cytoskeleton and recent studies have shown that increased nmMLCK activity regulates biological processes in various pulmonary diseases such as asthma and acute lung injury. In this study, we aimed to discover the role of nmMLCK in the proliferation and migration of pulmonary arterial endothelial cells (HPAECs) in the pathogenesis of PAH. We used two cellular models relevant to the pathobiology of PAH including BMPR2 silenced and vascular endothelial growth factor (VEGF) stimulated HPAECs. Both models demonstrated an increase in nmMLCK activity along with a robust increase in cellular proliferation, inflammation, and cellular migration. The upregulated nmMLCK activity was also associated with increased ERK expression pointing towards a potential integral cytoplasmic interaction. Mechanistically, we confirmed that when nmMLCK is inhibited by MLCK selective inhibitor (ML-7), proliferation and migration are attenuated. In conclusion, our results demonstrate that nmMLCK upregulation in association with increased ERK expression may contribute to the pathogenesis of PAHby stimulating cellular proliferation and migration.

Safranal Inhibits Angiogenesis via Targeting HIF-1α/VEGF Machinery: In Vitro and Ex Vivo Insights

Nature has a nearly infinite inventory of unexplored phytochemicals and biomolecules that have the potential to treat a variety of diseases. Safranal exhibits anti-cancer property and the present study explores its antiangiogenic property. Hepatocellular carcinoma (HCC) ranks as the sixth deadliest among all cancer types. Targeting the non-tumor vasculature supporting system is very promising as it has less plasticity, unlike malignant cells that are often associated with issues like drug resistance, poor prognosis, and relapse. In this study, we successfully inhibited the proliferation of primary human umbilical vein endothelial cells (HUVEC) with an IC50 of 300μM and blocked VEGF secretion in HepG2 cells. Furthermore, safranal inhibited VEGF-induced angiogenesis in vitro and ex vivo via scratch wound assay, tube formation assay, transmembrane assay, and aortic ring assay. In addition, safranal downregulated the in vitro expression of HIF-1α, VEGF, VEGFR2, p-AKT, p-ERK1/2, MMP9, p-FAK, and p-STAT3. The present study is the first to reveal the antiangiogenic potential of safranal and propose its possible underlying mechanism in HCC.

DHA-enriched phosphatidylcholine suppressed angiogenesis by activating PPARγ and modulating the VEGFR2/Ras/ERK pathway in human umbilical vein endothelial cells

Docosahexaenoic acid-enriched phosphatidylcholine (DHA-PC) is a new generation of omega-3 lipids, which contains an ester bond linking DHA at the sn-2 position of phospholipid. DHA-PC has become the interest recently as its better bioavailability and anti-oxidation capacity. In this study, the anti-angiogenic effect of DHA-PC was evaluated. The capacities of proliferation, migration, tube formation of human umbilical vein endothelial cells were significantly declined after DHA-PC treatment. Furthermore, DHA-PC inhibited the neovascularization of the chick chorioallantoic membrane in vivo. Mechanism results indicated that DHA-PC enhances the expression of peroxisome proliferator-activated receptor γ (PPARγ) at transcriptional and translational level, subsequently down-regulates the VEGFR2 expression and VEGFR2-mediated downstream Ras/ERK pathway, resulting in significant reduction in proliferation and differentiation. Additionally, PPARγ-specific antagonist GW9662 partly reversed the inhibition effects of DHA-PC on tube formation and neovascularization, suggesting that DHA-PC exerts anti-angiogenesis effect through activating PPARγ. These findings indicated that DHA-PC has a great prospect of anti-tumor angiogenesis therapy.

Effects of Fruquintinib on the Pluripotency Maintenance and Differentiation Potential of Mouse Embryonic Stem Cells

Mouse embryonic stem cells (mESCs) can maintain self-renewal and differentiate into any cell type of the three primary germ layers. The vascular endothelial growth factor (VEGF) is involved in the regulation of mESC differentiation and induces the activation of a series of kinase responses and several cell signaling pathways by binding to its respective transmembrane receptors, vascular endothelial growth factor receptor VEGFR1, and VEGFR2. Fruquintinib is a selective inhibitor of VEGFRs, and we used it to investigate the effects on the maintenance of pluripotency and differentiation potential of mESCs in this study. Our results showed that fruquintinib-treated cells expressed higher levels of pluripotent markers, including Oct4, Nanog, Sox2, and Esrrb under serum and leukemia inhibitory factor (LIF) condition, whereas the expression of phosphorylated Erk1/2 was restricted. Mitogen-activated protein kinase (MAPK)/extracellular signal-regulated kinase (MEK) signaling inhibitor (PD0325901) and glycogen synthase kinase 3 (GSK3) signaling inhibitor (CHIR99021) (also known as 2i) enable cells to maintain naive pluripotency with LIF, and fruquintinib can also promote cells to maintain naive pluripotent state even under serum/LIF condition, whereas VEGF addition limits the pluripotency characteristics in serum/LIF mESCs. Furthermore, fruquintinib could inhibit the three-germ layer establishment in embryoid body formation and maintain the undifferentiated characteristics of mESCs, indicating that fruquintinib could promote the maintenance of naive pluripotency and inhibit early differentiation programs.

Anti-angiogenic effects of testis-specific gene antigen 10 on primary endothelial cells

Growing evidence indicates the antitumor and antiangiogenesis activities of testis-specific gene antigen 10 (TSGA10). However, the underlying mechanisms and precise role of TSGA10 in angiogenesis are still elusive. In this study, we isolated human umbilical cord vein endothelial cells (HUVECs) and stably transfected with pcDNA3.1 carrying TSGA10 coding sequence. We demonstrated that TSGA10 over-expression significantly decreases HUVEC tubulogenesis and interconnected capillary network formation. HUVECs over-expressing TSGA10 exhibited a significant decrease in migration and proliferation rates. TSGA10 over-expression markedly decreased expression of angiogenesis-related genes, including VEGF-A, VEGFR-2, Ang-1, Ang-2, and Tie-2. Our ELISA results showed the decrease in VEGF-A mRNA expression level is associated with a significant decrease in its protein secretion. Additionally, over-expressing TSGA10 decreased expression levels of marker genes of cell migration (MMP-2, MMP-9, and SDF-1a) and proliferation (PCNA and Ki-67. Furthermore, ERK-1 and AKT phosphorylation significantly reduced in HUVECs over-expressing TSGA10. Our findings suggest a potent anti-angiogenesis activity of TSGA10 in HUVECs through down-regulation of ERK and AKT signalling pathways, and may provide therapeutic benefits for the management of different pathological angiogenesis.

ERK and Akt exhibit distinct signaling responses following stimulation by pro-angiogenic factors

Background

Angiogenesis plays an important role in the survival of tissues, as blood vessels provide oxygen and nutrients required by the resident cells. Thus, targeting angiogenesis is a prominent strategy in many different settings, including both tissue engineering and cancer treatment. However, not all of the approaches that modulate angiogenesis lead to successful outcomes. Angiogenesis-based therapies primarily target pro-angiogenic factors such as vascular endothelial growth factor-A (VEGF) or fibroblast growth factor (FGF) in isolation, and there is a limited understanding of how these promoters combine together to stimulate angiogenesis. Targeting one pathway could be insufficient, as alternative pathways may compensate, diminishing the overall effect of the treatment strategy.

Methods

To gain mechanistic insight and identify novel therapeutic strategies, we have developed a detailed mathematical model to quantitatively characterize the crosstalk of FGF and VEGF intracellular signaling. The model focuses on FGF- and VEGF-induced mitogen-activated protein kinase (MAPK) signaling to promote cell proliferation and the phosphatidylinositol 3-kinase/protein kinase B (PI3K/Akt) pathway, which promotes cell survival and migration. We fit the model to published experimental datasets that measure phosphorylated extracellular regulated kinase (pERK) and Akt (pAkt) upon FGF or VEGF stimulation. We validate the model with separate sets of data.

Results

We apply the trained and validated mathematical model to characterize the dynamics of pERK and pAkt in response to the mono- and co-stimulation by FGF and VEGF. The model predicts that for certain ranges of ligand concentrations, the maximum pERK level is more responsive to changes in ligand concentration compared to the maximum pAkt level. Also, the combination of FGF and VEGF indicates a greater effect in increasing the maximum pERK compared to the summation of individual effects, which is not seen for maximum pAkt levels. In addition, our model identifies the influential species and kinetic parameters that specifically modulate the pERK and pAkt responses, which represent potential targets for angiogenesis-based therapies.

Conclusions

Overall, the model predicts the combination effects of FGF and VEGF stimulation on ERK and Akt quantitatively and provides a framework to mechanistically explain experimental results and guide experimental design. Thus, this model can be utilized to study the effects of pro- and anti-angiogenic therapies that particularly target ERK and/or Akt activation upon stimulation with FGF and VEGF.

Mathematical Model Predicts Effective Strategies to Inhibit VEGF-eNOS Signaling

The endothelial nitric oxide synthase (eNOS) signaling pathway in endothelial cells has multiple physiological significances. It produces nitric oxide (NO), an important vasodilator, and enables a long-term proliferative response, contributing to angiogenesis. This signaling pathway is mediated by vascular endothelial growth factor (VEGF), a pro-angiogenic species that is often targeted to inhibit tumor angiogenesis. However, inhibiting VEGF-mediated eNOS signaling can lead to complications such as hypertension. Therefore, it is important to understand the dynamics of eNOS signaling in the context of angiogenesis inhibitors. Thrombospondin-1 (TSP1) is an important angiogenic inhibitor that, through interaction with its receptor CD47, has been shown to redundantly inhibit eNOS signaling. However, the exact mechanisms of TSP1's inhibitory effects on this pathway remain unclear. To address this knowledge gap, we established a molecular-detailed mechanistic model to describe VEGF-mediated eNOS signaling, and we used the model to identify the potential intracellular targets of TSP1. In addition, we applied the predictive model to investigate the effects of several approaches to selectively target eNOS signaling in cells experiencing high VEGF levels present in the tumor microenvironment. This work generates insights for pharmacologic targets and therapeutic strategies to inhibit tumor angiogenesis signaling while avoiding potential side effects in normal vasoregulation.

The impact of tumor receptor heterogeneity on the response to anti-angiogenic cancer treatment

Multiple promoters and inhibitors mediate angiogenesis, the formation of new blood vessels, and these factors represent potential targets for impeding vessel growth in tumors. Vascular endothelial growth factor (VEGF) is a potent angiogenic factor targeted in anti-angiogenic cancer therapies. In addition, thrombospondin-1 (TSP1) is a major endogenous inhibitor of angiogenesis, and TSP1 mimetics are being developed as an alternative type of anti-angiogenic agent. The combination of bevacizumab, an anti-VEGF agent, and ABT-510, a TSP1 mimetic, has been tested in clinical trials to treat advanced solid tumors. However, the patients' responses are highly variable and show disappointing outcomes. To obtain mechanistic insight into the effects of this combination anti-angiogenic therapy, we have constructed a novel whole-body systems biology model including the VEGF and TSP1 reaction networks. Using this molecular-detailed model, we investigated how the combination anti-angiogenic therapy changes the amounts of pro-angiogenic and anti-angiogenic complexes in cancer patients. We particularly focus on answering the question of how the effect of the combination therapy is influenced by tumor receptor expression, one aspect of patient-to-patient variability. Overall, this model complements the clinical administration of combination anti-angiogenic therapy, highlights the role of tumor receptor variability in the heterogeneous responses to anti-angiogenic therapy, and identifies the tumor receptor profiles that correlate with a high likelihood of a positive response to the combination therapy. Our model provides novel understanding of the VEGF-TSP1 balance in cancer patients at the systems-level and could be further used to optimize combination anti-angiogenic therapy.

Mechanistic insight into activation of MAPK signaling by pro-angiogenic factors

Background

Angiogenesis is important in physiological and pathological conditions, as blood vessels provide nutrients and oxygen needed for tissue growth and survival. Therefore, targeting angiogenesis is a prominent strategy in both tissue engineering and cancer treatment. However, not all of the approaches to promote or inhibit angiogenesis lead to successful outcomes. Angiogenesis-based therapies primarily target pro-angiogenic factors such as vascular endothelial growth factor-A (VEGF) or fibroblast growth factor (FGF) in isolation. However, pre-clinical and clinical evidence shows these therapies often have limited effects. To improve therapeutic strategies, including targeting FGF and VEGF in combination, we need a quantitative understanding of the how the promoters combine to stimulate angiogenesis.

Results

In this study, we trained and validated a detailed mathematical model to quantitatively characterize the crosstalk of FGF and VEGF intracellular signaling. This signaling is initiated by FGF binding to the FGF receptor 1 (FGFR1) and heparan sulfate glycosaminoglycans (HSGAGs) or VEGF binding to VEGF receptor 2 (VEGFR2) to promote downstream signaling. The model focuses on FGF- and VEGF-induced mitogen-activated protein kinase (MAPK) signaling and phosphorylation of extracellular regulated kinase (ERK), which promotes cell proliferation. We apply the model to predict the dynamics of phosphorylated ERK (pERK) in response to the stimulation by FGF and VEGF individually and in combination. The model predicts that FGF and VEGF have differential effects on pERK. Additionally, since VEGFR2 upregulation has been observed in pathological conditions, we apply the model to investigate the effects of VEGFR2 density and trafficking parameters. The model predictions show that these parameters significantly influence the response to VEGF stimulation.

Conclusions

The model agrees with experimental data and is a framework to synthesize and quantitatively explain experimental studies. Ultimately, the model provides mechanistic insight into FGF and VEGF interactions needed to identify potential targets for pro- or anti-angiogenic therapies.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0668-5) contains supplementary material, which is available to authorized users.

Quantifying the phosphorylation timescales of receptor-ligand complexes: a Markovian matrix-analytic approach

Cells interact with the extracellular environment by means of receptor molecules on their surface. Receptors can bind different ligands, leading to the formation of receptor-ligand complexes. For a subset of receptors, called receptor tyrosine kinases, binding to ligand enables sequential phosphorylation of intra-cellular residues, which initiates a signalling cascade that regulates cellular function and fate. Most mathematical modelling approaches employed to analyse receptor signalling are deterministic, especially when studying scenarios of high ligand concentration or large receptor numbers. There exist, however, biological scenarios where low copy numbers of ligands and/or receptors need to be considered, or where signalling by a few bound receptor-ligand complexes is enough to initiate a cellular response. Under these conditions stochastic approaches are appropriate, and in fact, different attempts have been made in the literature to measure the timescales of receptor signalling initiation in receptor-ligand systems. However, these approaches have made use of numerical simulations or approximations, such as moment-closure techniques. In this paper, we study, from an analytical perspective, the stochastic times to reach a given signalling threshold for two receptor-ligand models. We identify this time as an extinction time for a conveniently defined auxiliary absorbing continuous time Markov process, since receptor-ligand association/dissociation events can be analysed in terms of quasi-birth-and-death processes. We implement algorithmic techniques to compute the different order moments of this time, as well as the steady-state probability distribution of the system. A novel feature of the approach introduced here is that it allows one to quantify the role played by each kinetic rate in the timescales of signal initiation, and in the steady-state probability distribution of the system. Finally, we illustrate our approach by carrying out numerical studies for the vascular endothelial growth factor and one of its receptors, the vascular endothelial growth factor receptor of human endothelial cells.

Computational modeling of synergistic interaction between αVβ3 integrin and VEGFR2 in endothelial cells: Implications for the mechanism of action of angiogenesis-modulating integrin-binding peptides

Cooperation between VEGFR2 and integrin αVβ3 is critical for neovascularization in wound healing, cardiovascular ischemic diseases, ocular diseases, and tumor angiogenesis. In the present study, we developed a rule-based computational model to investigate the potential mechanism by which the Src-induced integrin association with VEGFR2 enhances VEGFR2 activation. Simulations demonstrated that the main function of integrin is to reduce the degradation of VEGFR2 and hence stabilize the activation signal. In addition, receptor synthesis rate and recruitment from internal compartment were found to be sensitive determinants of the activation state of VEGFR2. The model was then applied to simulate the effect of integrin-binding peptides such as tumstatin and cilengitide on VEGFR2 signaling. Further, computational modeling proposed potential molecular mechanisms for the angiogenesis-modulating activity of other integrin-binding peptides. The model highlights the complexity of the crosstalk between αVβ3 integrin and VEGFR2 and the necessity of utilizing models to elucidate potential mechanisms in angiogenesis-modulating peptide therapy.

Computer Simulation of TSP1 Inhibition of VEGF-Akt-eNOS: An Angiogenesis Triple Threat

The matricellular protein thrombospondin-1 (TSP1) is a potent inhibitor of angiogenesis. Specifically, TSP1 has been experimentally shown to inhibit signaling downstream of vascular endothelial growth factor (VEGF). The molecular mechanism of this inhibition is not entirely clear. We developed a detailed computational model of VEGF signaling to Akt-endothelial nitric oxide synthase (eNOS) to investigate the quantitative molecular mechanism of TSP1 inhibition. The model demonstrated that TSP1 acceleration of VEGFR2 degradation is sufficient to explain the inhibition of VEGFR2 and eNOS phosphorylation. However, Akt inhibition requires TSP1-induced phosphatase recruitment to VEGFR2. The model was then utilized to test various strategies for the rescue of VEGF signaling to Akt and eNOS. Inhibiting TSP1 was predicted to be not as effective as CD47 depletion in rescuing signaling to Akt. The model further predicts that combination strategy involving depletion of CD47 and inhibition of TSP1 binding to CD47 is necessary for effective recovery of signaling to eNOS. In all, computational modeling offers insight to molecular mechanisms involving TSP1 interaction with VEGF signaling and provides strategies for rescuing angiogenesis by targeting TSP1-CD47 axis.

Integrative meta-modeling identifies endocytic vesicles, late endosome and the nucleus as the cellular compartments primarily directing RTK signaling

Recently, intracellular receptor signaling has been identified as a key component mediating cell responses for various receptor tyrosine kinases (RTKs). However, the extent each endocytic compartment (endocytic vesicle, early endosome, recycling endosome, late endosome, lysosome and nucleus) contributes to receptor signaling has not been quantified. Furthermore, our understanding of endocytosis and receptor signaling is complicated by cell- or receptor-specific endocytosis mechanisms. Therefore, towards understanding the differential endocytic compartment signaling roles, and identifying how to achieve signal transduction control for RTKs, we delineate how endocytosis regulates RTK signaling. We achieve this via a meta-analysis across eight RTKs, integrating computational modeling with experimentally derived cell (compartment volume, trafficking kinetics and pH) and ligand-receptor (ligand/receptor concentration and interaction kinetics) physiology. Our simulations predict the abundance of signaling from eight RTKs, identifying the following hierarchy in RTK signaling: PDGFRβ > IGFR1 > EGFR > PDGFRα > VEGFR1 > VEGFR2 > Tie2 > FGFR1. We find that endocytic vesicles are the primary cell signaling compartment; over 43% of total receptor signaling occurs within the endocytic vesicle compartment for these eight RTKs. Mechanistically, we found that high RTK signaling within endocytic vesicles may be attributed to their low volume (5.3 × 10^-19 L) which facilitates an enriched ligand concentration (3.2 μM per ligand molecule within the endocytic vesicle). Under the analyzed physiological conditions, we identified extracellular ligand concentration as the most sensitive parameter to change; hence the most significant one to modify when regulating absolute compartment signaling. We also found that the late endosome and nucleus compartments are important contributors to receptor signaling, where 26% and 18%, respectively, of average receptor signaling occurs across the eight RTKs. Conversely, we found very low membrane-based receptor signaling, exhibiting <1% of the total receptor signaling for these eight RTKs. Moreover, we found that nuclear translocation, mechanistically, requires late endosomal transport; when we blocked receptor trafficking from late endosomes to the nucleus we found a 57% reduction in nuclear translocation. In summary, our research has elucidated the significance of endocytic vesicles, late endosomes and the nucleus in RTK signal propagation.

VEGFR1 promotes cell migration and proliferation through PLCγ and PI3K pathways

The ability to control vascular endothelial growth factor (VEGF) signaling offers promising therapeutic potential for vascular diseases and cancer. Despite this promise, VEGF-targeted therapies are not clinically effective for many pathologies, such as breast cancer. VEGFR1 has recently emerged as a predictive biomarker for anti-VEGF efficacy, implying a functional VEGFR1 role beyond its classically defined decoy receptor status. Here we introduce a computational approach that accurately predicts cellular responses elicited via VEGFR1 signaling. Aligned with our model prediction, we show empirically that VEGFR1 promotes macrophage migration through PLCγ and PI3K pathways and promotes macrophage proliferation through a PLCγ pathway. These results provide new insight into the basic function of VEGFR1 signaling while offering a computational platform to quantify signaling of any receptor.

Netrin-1 Prevents Rat Primary Cortical Neurons from Apoptosis via the DCC/ERK Pathway

In the nervous system, Netrin-1 serves as a neural guide, mediating the neuronal development. However, it remains blurred whether Netrin-1 can protect neurons from apoptosis induced by cerebral stroke. In the current study, the cultured rat primary cortical neurons were transfected with Netrin-1-encoding lentivirus before the oxygen-glucose-deprivation (OGD) treatment. Cell death and apoptosis were evaluated by lactate dehydrogenase (LDH) release and flow cytometry. We found that Netrin-1 attenuated OGD-induced cell death and neuronal apoptosis at 24 h after OGD treatment, and that the overexpression of Netrin-1 activated the ERK signaling pathway. These effects were partly abolished by blocking its receptor deleted in colorectal cancer (DCC) or U0126, an inhibitor of the ERK signaling pathway. Netrin-1 overexpression in neurons elevated the expression of DCC, on mRNA level and protein level. Netrin-1 also reduced DNA damage. Taken together, our findings suggest that Netrin-1 attenuates cell death and neuronal apoptosis via the DCC/ERK signaling pathway in the cultured primary cortical neurons after OGD injury, which may involve the mediation of DNA damage in the neurons.

Computational investigation of sphingosine kinase 1 (SphK1) and calcium dependent ERK1/2 activation downstream of VEGFR2 in endothelial cells

Vascular endothelial growth factor (VEGF) is a powerful regulator of neovascularization. VEGF binding to its cognate receptor, VEGFR2, activates a number of signaling pathways including ERK1/2. Activation of ERK1/2 is experimentally shown to involve sphingosine kinase 1 (SphK1) activation and its calcium-dependent translocation downstream of ERK1/2. Here we construct a rule-based computational model of signaling downstream of VEGFR2, by including SphK1 and calcium positive feedback mechanisms, and investigate their consequences on ERK1/2 activation. The model predicts the existence of VEGF threshold in ERK1/2 activation that can be continuously tuned by cellular concentrations of SphK1 and sphingosine 1 phosphate (S1P). The computer model also predicts powerful effects of perturbations in plasma and ER calcium pump rates and the current through the CRAC channels on ERK1/2 activation dynamics, highlighting the critical role of intracellular calcium in shaping the pERK1/2 signal. The model is then utilized to simulate anti-angiogenic therapeutic interventions targeting VEGFR2-ERK1/2 axis. Simulations indicate that monotherapies that exclusively target VEGFR2 phosphorylation, VEGF, or VEGFR2 are ineffective in shutting down signaling to ERK1/2. By simulating therapeutic strategies that target multiple nodes of the pathway such as Raf and SphK1, we conclude that combination therapy should be much more effective in blocking VEGF signaling to EKR1/2. The model has important implications for interventions that target signaling pathways in angiogenesis relevant to cancer, vascular diseases, and wound healing.

Vascular endothelial growth factor (VEGF) signaling is a potent regulator of angiogenesis, the growth and development of new vessels out of a preexisting vascular network. Angiogenesis requires enhanced survival, proliferation, and motility of the vascular endothelial cells. Crucial signaling endpoints in VEGF-mediated angiogenic response include elevation in intracellular calcium and the activation of the proteins ERK1 and 2 (ERK1/2). In this study, we have developed a novel computer model for the activation of ERK1/2 and calcium downstream of VEGF receptor type 2 (VEGFR2). Our model is the first of its kind to incorporate and investigate the consequences of calcium elevation and the role of a cellular lipid modifier known as sphingosine kinase 1 (SphK1). We also utilize the model to simulate therapeutic strategies targeting VEGF signaling to ERK1/2 indicating inefficiency of single therapies known as tyrosine kinase inhibitors (TKI) that target receptor phosphorylation. Computer simulations indicate that combination therapy is essential for effective blockade of this important pathway. Our results have important implications for human diseases such as cancer where plethora of anti-VEGF therapies are currently employed. Overall, our computer model sheds new light on a complex feedback involving SphK1 and calcium that radically alters the response of cells to VEGF.

Fibronectin fibrils regulate TGF-β1-induced Epithelial-Mesenchymal Transition

Epithelial-Mesenchymal Transition (EMT) is a dynamic process through which epithelial cells transdifferentiate from an epithelial phenotype into a mesenchymal phenotype. Previous studies have demonstrated that both mechanical signaling and soluble growth factor signaling facilitate this process. One possible point of integration for mechanical and growth factor signaling is the extracellular matrix. Here we investigate the role of the extracellular matrix (ECM) protein fibronectin (FN) in this process. We demonstrate that inhibition of FN fibrillogenesis blocks activation of the Transforming Growth Factor-Beta (TGF-β) signaling pathway via Smad2 signaling, decreases cell migration and ultimately leads to inhibition of EMT. Results show that soluble FN, FN fibrils, or increased contractile forces are insufficient to independently induce EMT. We further demonstrate that inhibition of latent TGF-β1 binding to FN fibrils via either a monoclonal blocking antibody against the growth factor binding domain of FN or through use of a FN deletion mutant that lacks the growth factor binding domains of FN blocks EMT progression, indicating a novel role for FN in EMT in which the assembly of FN fibrils serves to localize TGF-β1 signaling to drive EMT.

VEGF regulates local inhibitory complement proteins in the eye and kidney

Outer retinal and renal glomerular functions rely on specialized vasculature maintained by VEGF that is produced by neighboring epithelial cells, the retinal pigment epithelium (RPE) and podocytes, respectively. Dysregulation of RPE- and podocyte-derived VEGF is associated with neovascularization in wet age-related macular degeneration (ARMD), choriocapillaris degeneration, and glomerular thrombotic microangiopathy (TMA). Since complement activation and genetic variants in inhibitory complement factor H (CFH) are also features of both ARMD and TMA, we hypothesized that VEGF and CFH interact. Here, we demonstrated that VEGF inhibition decreases local CFH and other complement regulators in the eye and kidney through reduced VEGFR2/PKC-α/CREB signaling. Patient podocytes and RPE cells carrying disease-associated CFH genetic variants had more alternative complement pathway deposits than controls. These deposits were increased by VEGF antagonism, a common wet ARMD treatment, suggesting that VEGF inhibition could reduce cellular complement regulatory capacity. VEGF antagonism also increased markers of endothelial cell activation, which was partially reduced by genetic complement inhibition. Together, these results suggest that VEGF protects the retinal and glomerular microvasculature, not only through VEGFR2-mediated vasculotrophism, but also through modulation of local complement proteins that could protect against complement-mediated damage. Though further study is warranted, these findings could be relevant for patients receiving VEGF antagonists.

Selective Targeting of a Novel Epsin-VEGFR2 Interaction Promotes VEGF-Mediated Angiogenesis

RATIONALE

We previously reported that vascular endothelial growth factor (VEGF)-induced binding of VEGF receptor 2 (VEGFR2) to epsins 1 and 2 triggers VEGFR2 degradation and attenuates VEGF signaling. The epsin ubiquitin interacting motif (UIM) was shown to be required for the interaction with VEGFR2. However, the molecular determinants that govern how epsin specifically interacts with and regulates VEGFR2 were unknown.

OBJECTIVE

The goals for the present study were as follows: (1) to identify critical molecular determinants that drive the specificity of the epsin and VEGFR2 interaction and (2) to ascertain whether such determinants were critical for physiological angiogenesis in vivo.

METHODS AND RESULTS

Structural modeling uncovered 2 novel binding surfaces within VEGFR2 that mediate specific interactions with epsin UIM. Three glutamic acid residues in epsin UIM were found to interact with residues in VEGFR2. Furthermore, we found that the VEGF-induced VEGFR2-epsin interaction promoted casitas B-lineage lymphoma-mediated ubiquitination of epsin, and uncovered a previously unappreciated ubiquitin-binding surface within VEGFR2. Mutational analysis revealed that the VEGFR2-epsin interaction is supported by VEGFR2 interacting specifically with the UIM and with ubiquitinated epsin. An epsin UIM peptide, but not a mutant UIM peptide, potentiated endothelial cell proliferation, migration and angiogenic properties in vitro, increased postnatal retinal angiogenesis, and enhanced VEGF-induced physiological angiogenesis and wound healing.

CONCLUSIONS

Distinct residues in the epsin UIM and VEGFR2 mediate specific interactions between epsin and VEGFR2, in addition to UIM recognition of ubiquitin moieties on VEGFR2. These novel interactions are critical for pathophysiological angiogenesis, suggesting that these sites could be selectively targeted by therapeutics to modulate angiogenesis.

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.

Systems pharmacology approaches for optimization of antiangiogenic therapies: challenges and opportunities

Targeted therapies have become an important therapeutic paradigm for multiple malignancies. The rapid development of resistance to these therapies impedes the successful management of advanced cancer. Due to the redundancy in angiogenic signaling, alternative proangiogenic factors are activated upon treatment with anti-VEGF agents. Higher doses of the agents lead to greater stimulation of compensatory proangiogenic pathways that limit the therapeutic efficacy of VEGF-targeted drugs and produce escape mechanisms for tumor. Evidence suggests that dose intensity and schedules affect the dynamics of the development of this resistance. Thus, an optimal dosing regimen is crucial to maximizing the therapeutic benefit of antiangiogenic agents and limiting treatment resistance. A systems pharmacology approach using multiscale computational modeling can facilitate a mechanistic understanding of these dynamics of angiogenic biomarkers and their impacts on tumor reduction and resistance. Herein, we discuss a systems pharmacology approach integrating the biology of VEGF-targeted therapy resistance, including circulating biomarkers, and pharmacodynamics to enable the optimization of antiangiogenic therapy for therapeutic gains.

Quantitative insight in utilizing circulating angiogenic factors as biomarkers for antiangiogenic therapy: systems pharmacology approach

Circulating angiogenic factors (CAF) like vascular endothelial growth factor (VEGF), placental growth factor (PlGF), and sVEGFR2 have potential as biomarkers for antiangiogenic therapy. The interpretation of changes in CAF is complicated by the dynamic nature of the tumor and host cells emanating CAF in response to VEGF pathway inhibition. We developed a systems pharmacology model of anti-VEGF agents to investigate CAF modulation by tumor and host cells, and the relationship between overall CAF changes in response to sunitinib and antitumor efficacy. This model distinguishes between the tumor cells' contributions from tumor-independent response to therapy and total plasma CAF correlating with antitumor activity. Altered VEGF is more likely to serve as a useful biomarker reflecting tumor responses in cancer patients whose pretreatment VEGF is higher than baseline VEGF in healthy subjects. Our findings provide a mechanistic insight into tumor modulation of angiogenic molecules, and may explain the inconsistent results found in previous biomarker studies.

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.

Quantitative characterization of cellular membrane-receptor heterogeneity through statistical and computational modeling

Cell population heterogeneity can affect cellular response and is a major factor in drug resistance. However, there are few techniques available to represent and explore how heterogeneity is linked to population response. Recent high-throughput genomic, proteomic, and cellomic approaches offer opportunities for profiling heterogeneity on several scales. We have recently examined heterogeneity in vascular endothelial growth factor receptor (VEGFR) membrane localization in endothelial cells. We and others processed the heterogeneous data through ensemble averaging and integrated the data into computational models of anti-angiogenic drug effects in breast cancer. Here we show that additional modeling insight can be gained when cellular heterogeneity is considered. We present comprehensive statistical and computational methods for analyzing cellomic data sets and integrating them into deterministic models. We present a novel method for optimizing the fit of statistical distributions to heterogeneous data sets to preserve important data and exclude outliers. We compare methods of representing heterogeneous data and show methodology can affect model predictions up to 3.9-fold. We find that VEGF levels, a target for tuning angiogenesis, are more sensitive to VEGFR1 cell surface levels than VEGFR2; updating VEGFR1 levels in the tumor model gave a 64% change in free VEGF levels in the blood compartment, whereas updating VEGFR2 levels gave a 17% change. Furthermore, we find that subpopulations of tumor cells and tumor endothelial cells (tEC) expressing high levels of VEGFR (>35,000 VEGFR/cell) negate anti-VEGF treatments. We show that lowering the VEGFR membrane insertion rate for these subpopulations recovers the anti-angiogenic effect of anti-VEGF treatment, revealing new treatment targets for specific tumor cell subpopulations. This novel method of characterizing heterogeneous distributions shows for the first time how different representations of the same data set lead to different predictions of drug efficacy.