PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of npjsbaLink to Publisher's site
 
NPJ Syst Biol Appl. 2018; 4: 1.
Published online 2017 December 19. doi:  10.1038/s41540-017-0037-9
PMCID: PMC5736688

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

Abstract

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.

Introduction

Vascular endothelial growth factor (VEGF) plays a critical role in many pathologies, including vascular disease and cancer.15 Despite this role, VEGF-targeted therapies are not clinically effective for many patients.6,7 As such, there is an urgent need to develop a better understanding of how VEGF-promoted pathologies can be controlled, mechanistically, to improve the efficiency and specificity of current VEGF treatments.

VEGF receptor-1 (VEGFR1) has emerged as a predictive biomarker for anti-VEGF therapeutics in cancer,810 but its signaling mechanisms and function remain incompletely defined. VEGFR1 is conventionally described as a decoy receptor that does not produce intracellular signals, due to its high VEGF affinity but low phosphorylation compared to VEGFR2.11 However, emerging evidence suggests an active VEGFR1 signaling role: membrane VEGFR1 is upregulated during vascular reperfusion stages in ischemic tissue,12 in both hypoxic tumor cells and tumor endothelial cells,13 and mice lacking VEGFR1 signaling exhibit reduced tumor vascularization.14 Furthermore, VEGFR1 demonstrates tumor activity via placental growth factor (PlGF);15,16 PlGF inhibition has shown promise to prevent tumor growth and metastasis.17 Given such emerging evidence, we believe that VEGFR1 must have an important signaling role, and we aim to delineate it.

VEGFR1 signaling can be characterized by systems biology to mathematically define receptor signaling mechanisms. The power of this mechanistic approach is its faithfulness to the biological structure of the receptor. Toward this end, the two key signaling mechanisms post-VEGFR1 ligation include: (1) carboxy-terminal VEGFR1 phosphorylation at specific tyrosine sites and (2) adapter binding at these sites. We define these as the key steps because they structurally facilitate the second messenger signaling that directs the angiogenic hallmarks of cell proliferation and migration;1820 as such, these steps may together predict those hallmarks. Indeed, there is evidence that tyrosine site phosphorylation is linked to cell response: cell proliferation results from phosphorylation at the VEGFR2 Tyr1175 site, whereas phosphorylation at the VEGFR2 Tyr1214 site is linked to cell migration.18 Cell responses are similarly associated with adapter binding and phosphorylation at receptor phosphor-tyrosine sites:2125 receptor-induced phosphoinositide-3 kinase (PI3K)-p85α/γ regulatory subunit phosphorylation (hence simply called PI3K phosphorylation) is known to result in cell migration.25 While these tyrosine site- and adapter-based approaches are useful to predict cell response, they are often analyzed separately, which does not enable a unified understanding of how RTK structure directs cell function.26,27 Therefore, computational models that integrate adapter binding and phosphorylation at specific receptor tyrosine sites would advance structure-based predictions of VEGFR1 signaling.

Here we predict how VEGFR1 directs cell responses by developing, comparing, and validating a structure-based model of carboxy-terminal VEGFR activation and a general VEGFR activation model. We validate our modeling approach experimentally by quantifying adapter phosphorylation and cell migration and proliferation stemming from both VEGFR1 and VEGFR2 signaling and computationally parse out VEGFR1 signaling alone to map the VEGFR1 function. The models quantitatively rank adapter protein contributions to VEGFR1-mediated cell migration and cell proliferation. Model comparison reveals how degrees of model complexity affect predictions of receptor activation and cell response. Computational predictions of cell response to drug treatment are validated via functional assays. Together, our modeling approach provides a new, validated tool for structure-based prediction of cell signaling, applied to grant the exigent mapping of VEGFR1.

Results

VEGFR1 primarily induces cell migration

Following VEGF binding, the initial intracellular VEGFR signal transduction steps include: receptor dimerization; autophosphorylation, a post-translational modification (PTM) of carboxy-terminal tyrosines; adapter binding to phospho-tyrosine residues; and adapter phosphorylation. Here we model these receptor signaling mechanisms using mass-action kinetics: specifically, we computationally model VEGF-induced VEGFR phosphorylation, specific adapters binding the VEGFRs, and adapter phosphorylation (Fig. 1). To identify the importance of individual receptor sites in directing aggregated cell responses, we model adapter binding and PTMs occurring non-specifically (nonspecific model) at a single tyrosine site (a common receptor modeling approach) compared to adapter binding and PTMs occurring at specific receptor tyrosine sites (specific model, representing complete receptor physiology) (Fig. 2a). VEGFR-induced cell migration and proliferation were modeled by calibrating adapter phosphorylation to each specific cell response, which interested readers can find a full mathematical description of in the Supplementary Information. To understand the VEGFR1 function, we computationally predict cell signaling stemming from VEGFR1 alone (Figs. 23). Both the nonspecific and specific models predict that VEGFR1 primarily induces cell migration (Fig. 2b). This is evidenced by migration exhibiting both the highest integrated cell response (Fig. 2c) and the highest phosphorylation amplitude (Fig. 2d). The specific model reveals mechanistic insight into the migratory cell response: the VEGFR1 tyrosine sites specify cell migration signaling. This is evidenced by the specific model exhibiting a greater contribution to migration signaling; the integrated migration response, relative to proliferation and degradation, increases 16% in the specific model, relative to the nonspecific model (Fig. 2c). Furthermore, the migration phosphorylation amplitude increases 23% in the specific model, relative to the nonspecific model (Fig. 2d). Therefore, we predict that VEGFR1 tyrosine sites are structured to specify cell migration signaling and aim to identify which adapters result in cell migration.

Fig. 1
VEGFR–adapter interaction schematics. This schematic depicts the known biology of VEGFR1, VEGFR2, and their related adapters. Adapters bind specific tyrosine (Tyr) sites on a VEGFR1 and b VEGFR2 (SI Table S4). VEGFR1 and VEGFR2 kinase ...
Fig. 2
The VEGFR1 structure preferentially activates PLCγ and PI3K. a Schematics for the VEGFR–adapter interaction models: (left) adapters bind a single nonspecific VEGFR tyrosine site vs. (right) adapters binding specific VEGFR tyrosine sites. ...
Fig. 3
VEGFR1-induced cell responses are primarily directed by PLCγ, PI3K, and Src concentrations. HUVEC ac Integrated responses (area under the activation–time curve) and df phosphorylation amplitudes for each cell response, ...

VEGFR1 tyrosine sites specify PLCγ, and PI3K activation through adapter binding competition

VEGFR1 tyrosine sites specify cell migration signaling through phospholipase C-gamma (PLCγ) and PI3K phosphorylation (Fig. 2e). PLCγ and PI3K are the only adapters with increased integrated responses (Fig. 2f) and phosphorylation amplitudes (Fig. 2g) between nonspecific and specific models. This unique increase in PLCγ and PI3K activation is due to their binding preference with the VEGFR1 phospho-tyrosine sites (Fig. 1a); only two adapters bind VEGFR1 simultaneously (SI Tables S4 and S5): one adapter at Tyr794 and a second adapter at another tyrosine site. PI3K and PLCγ are the only adapters that bind Tyr794, thus experiencing less VEGFR1-binding competition than the other adapters, resulting in greater activation. This is evidenced by PLCγ and PI3K activation preferentially occurring at Tyr794 (SI Fig S1).

VEGFR1-promoted cell responses are regulated by coordinated PLCγ, PI3K, and Src activation

To predict which adapters primarily direct VEGFR1 cell responses, we perform sensitivity analyses between adapter concentrations and cell responses with the specific site model. We predict that cell proliferation and migration are primarily mediated by PLCγ, PI3K, and Src concentrations, in that order (Fig. 3a, b, d, e). Conversely, degradation signaling is not highly altered by adapter concentration (Fig. 3c, f). These three adapters direct VEGFR1 signaling in a coordinated fashion: increasing the PLCγ (Fig. 3g, j), PI3K (Fig. 3h, k), or Src (Fig. 3i, l) concentration to ~2×104 molecules/cell increases phosphorylation of the other two adapters. Increasing PI3K (Fig. 3h) and Src (Fig. 3i) concentrations above ~2×104 molecules/cell increases the PLCγ integrated response, indicating that PI3K and Src promote PLCγ phosphorylation. Together with our result that VEGFR1 is structured to preferentially activate PLCγ and PI3K, we predict that PLCγ and PI3K mediate VEGFR1 cell responses through coordinated activation involving Src.

Specific tyrosine site modeling captures adapter phosphorylation dynamics

The specific model accurately predicts PI3K phosphorylation dynamics and magnitude in VEGF-treated RAW 264.7 macrophages, evidenced by the Χ2 goodness-of-fit test (Fig. 4a).28 All model validations were performed by modeling both VEGFR1 and VEGFR2 signaling (Figs. 4 and and5),5), as RAWs express both these receptors (SI Fig S3). However, we identify that VEGFR1 signaling dominates VEGF signaling in RAWs computationally (SI Fig S4), and focus on the VEGFR1 signaling contribution. The specific model accurately predicts that PI3K phosphorylation is abrogated by the PI3K-specific inhibitor Wortmannin, while relatively unaffected by inhibiting other adapters (Fig. 4a). Conversely, the nonspecific model accurately predicts relative phosphorylation trends (SI Fig S2) but not phosphorylation magnitudes; the nonspecific model underestimates PI3K phosphorylation by 81% and fails the Χ2 goodness-of-fit test (Fig. 4a). Model-predicted PLCγ phosphorylation shows the same trend: the site-specific model accurately predicts PLCγ phosphorylation given VEGF and inhibitor treatments, whereas the nonspecific model fails validation (Fig. 4b). The specific model also accurately identifies which VEGFR1-associated adapters are not critical to VEGFR1 signaling: Abl phosphorylation is not detected as predicted (Fig. 4c). This validation highlights that modeling-specific receptor tyrosine sites is essential to capturing adapter phosphorylation magnitudes and is translatable across cell lines, whereas the conventional approach to model a nonspecific receptor tyrosine site fails physiological validation.

Fig. 4
VEGFR1 phosphorylates PI3K and PLCγ with model predicted dynamics. a PI3K, b PLCγ, and c Abl phosphorylation in RAWs were quantified with ELISAs at multiple time points given treatment with VEGF-A164 (50 ng/mL), 100 nM Wortmannin ...
Fig. 5
PLCγ and PI3K regulate VEGFR1-induced cell responses in vitro. a RAW migration was measured in wound-healing assays at 0 and 24 h post scratch. Scale bars represent 50 µm. b Analyzed wound-healing assays show that inhibiting ...

PI3K and PLCγ are critical to VEGFR1-induced cell migration

We validate our prediction that VEGFR1 promotes cell migration, which is primarily regulated by PLCγ followed by PI3K. We find that VEGFR1 does promote cell migration: VEGF induces significant RAW migration in vitro (Fig. 5a, b). Furthermore, VEGFR1-induced migration is primarily regulated by PLCγ, followed by PI3K (Fig. 5a, b). Our specific VEGFR1 tyrosine site model accurately quantifies adapter contributions to RAW migration; RAW migration decreases 79% in vitro with PLCγ inhibition (72% predicted) and 64% with PI3K inhibition (64% predicted) (Fig. 5b). Additionally, our model accurately identifies that Abl is insignificant to VEGFR1-induced migration (Fig. 5b).

VEGFR1-induced cell proliferation is primarily mediated via PLCγ

We validate our prediction that VEGFR1 promotes cell proliferation, primarily through PLCγ activation. VEGFR1 promotes cell proliferation: VEGF induces significant RAW proliferation in vitro (Fig. 5c). We validate our prediction that VEGFR1-induced migration is only significantly regulated by PLCγ; RAW proliferation decreases 50% in vitro with PLCγ inhibition (Fig. 5c). Conversely, PI3K and Abl inhibition do not significantly affect cell proliferation, accurately predicted by the specific VEGFR1 site model.

Discussion

The VEGFR1 status as a decoy receptor may not fully capture its signaling role;8 however, few studies have probed VEGFR1 signaling,11 which is due to the low phosphorylation levels VEGFR1 exhibits. As tyrosine kinase receptor family is known to signal through coupling with the SH2 domain of adapters,29 examining VEGFR1-adapter binding can offer new insight into VEGFR1 signal propagation. To this end, we developed and validated a receptor-adapter interaction modeling approach, which accurately predicts cell responses from adapter phosphorylation, and is translatable across receptor and cell types. Combining this modeling approach with experimental validation identified that VEGFR1 induces macrophage migration via PLCγ and PI3K pathways and induces proliferation via a PLCγ pathway.

Modeling techniques allow prediction of receptor signaling roles

Our modeling approach quantifies adapter phosphorylation and cell responses simultaneously to map unknown receptor signaling pathways. Our modeling approach integrates the pioneered approaches that accurately predict select adapter–receptor interactions3033 and cell responses3436 from external stimuli. We additionally advance receptor signaling models by providing the ability to map unknown receptor pathways. Furthermore, we show that our approach to model specific receptor tyrosine sites offers physiological relevancy; both nonspecific and specific VEGFR1 tyrosine site models are validated when only the shape of adapter phosphorylation over time is considered (SI Fig S2), but only the specific tyrosine site model accurately predicts adapter phosphorylation magnitudes (Fig. 4). Additionally, our receptor–adapter modeling approach can be easily integrated into pharmacokinetic/pharmacodynamic models that quantify extracellular VEGF dynamics in response to anti-VEGF drugs.3739 One of the major challenges for developing personalized, clinically relevant computational platforms is the difficulty of capturing all relevant physiological processes at multiple scales.40,41 By integrating extracellular VEGF dynamics, VEGF–VEGFR interactions, and subsequent intracellular VEGFR signaling, we can provide a clinically relevant platform to explore how anti-VEGF drugs mediate VEGFR signaling simultaneously at the tissue macroscale and intracellular microscale.

Mapping the native VEGFR1 function requires VEGFR2 signaling

Our modeling approach accurately quantified adapter phosphorylation and cell proliferation and migration in macrophages, which was then extended to map the VEGFR1 function. While macrophages express high VEGFR1, they also lowly express VEGFR2 (SI Fig S3). While our model captures signaling stemming from both VEGFR1 and VEGFR2, the question may arise as to why these experiments were not conducted with VEGFR2 inhibition to measure VEGFR1 signaling alone. The reason is that we seek to map the native function of VEGF-stimulated VEGFR1, which cannot be identified if VEGFR2 signaling is abolished. Indeed, it is well established that VEGFR1 and VEGFR2 expression are inversely related;12,4244 abolishing VEGFR2 signaling would increase VEGFR1 expression, which may affect VEGFR1-adapter binding kinetics, altering adapter phosphorylation and cell response dynamics. Furthermore, VEGFR2 knockdown has been shown to alter basal cell physiology and function, altering cell proliferation potential of unstimulated cells.42,45 VEGFR2 knockdown also alters basal levels of total intracellular kinase expression and baseline phosphorylation: VEGFR2 knockdown increases baseline c-Jun phosphorylation in unstimulated cells,42 a kinase downstream PI3K that regulates cell proliferation and migration.46,47 By contrast, VEGFR2 knockdown decreases baseline phosphorylation of extracellular signal–regulated kinase and Akt,45 kinases of downstream adapters such as Src and PI3K that promote cell migration/proliferation.48 We expect such adapter and downstream kinase expression alterations from VEGFR2 knockdown to alter VEGFR1 signal propagation and cell response dynamics. Rather than studying VEGFR2 inhibition for these reasons, we computationally parse out the RAW signaling contribution from VEGFR1 compared to VEGFR2, identifying that VEGFR1 signaling dominates VEGF signaling in RAWS; VEGFR1 exhibits a 2.4-fold higher RAW migration, and 2.6-fold higher RAW proliferation, integrated response than VEGFR2. Thus VEGFR2 knockdown presents dynamic changes in cell physiology that may alter the VEGFR1 function, which should be explored in a future study.

Translating receptor signaling across cell lines advances translational research

Our modeling approach accurately quantified VEGFR1-induced macrophage migration and proliferation through adapter phosphorylation. Our modeling approach required a calibration step where adapters VEGFR1 signals through were identified in endothelial cells (SI Table S3), due to the rich VEGF–VEGFR data available in endothelial cells,49,50 and then translated to simulate VEGFR1 signaling in macrophages, which are data poorer in regards to VEGF signaling.51 The rationale for this approach is our hypothesis that VEGFR1 is the exact same protein in endothelial cells and macrophages, and thus differential roles VEGFR1 has in these cells must be dependent on the intracellular (adapter expression) and extracellular (VEGFR2 surface expression) environments, as opposed to an inherent difference in VEGFR1 function. Indeed, we show that there is a sizeable difference between VEGFR1-interacting adapter expression (SI Tables S1 and S6) and VEGFR2 expression (SI Table S1, Fig S4) for human umblical vein endothelial cells (HUVECs) and RAW macrophages, agreeing with this hypothesis. We performed this translation across cell lines for two reasons: (1) it allows higher confidence in our model validations, since we are seeding and validating our model from two independent data sources, and (2) it highlights the translatable feature of this modeling approach. Indeed, the ability to translate across cell types has the potential to reduce experimental costs and increase model development time, as the richest data source available can be used even when modeling a new condition. This increased model development time and reduced experimental costs advances translational research by allowing model-informed decisions to be made quicker and with higher confidence.

qFlow cytometry accurately quantifies membrane receptors

Our ability to accurately quantify VEGFR1 signaling highlights the power of integrating experiment and computation to provide new biology insight: empirical evidence defined VEGFR1–adapter reactions, kinetics, and concentrations for our model, which in turn provided testable VEGFR1 signaling predictions that we confirmed experimentally. This first step, model parameterization, is essential to develop physiologically relevant models, as previously described.5254 We achieved VEGFR concentration parameterization with quantitative flow (qFlow) cytometry,12,55,56 a recently established high-throughput approach that detects receptor expression with a fluorescent affinity probe and quantifies absolute receptor concentrations using fluorescent calibration standards.55 While qFlow cytometry is becoming an essential tool for parameterizing receptor concentrations in computational models,8,32,37,39,5760 analogous methods for quantifying other receptor signaling parameters, such as adapter phosphorylation rates, are not well established. As such, most computational models contain parameters that are estimated or generalized across multiple species or interactions;61 Bose and Janes recently developed one such method for high-throughput characterization of signal molecule dephosphorylation kinetics via phosphatase activity.62 Development of such high-throughput methods to completely parameterize receptor signaling models, from species concentrations to specific kinetics for every interaction, would unlock additional options for tuning receptor signaling, such as by targeting specific phosphatases, while maintaining high physiological relevancy.

VEGFR1 preferentially activates PLCγ in burst activation to induce cell migration and proliferation, possibly through Ca2+ signaling

We show that VEGFR1-induced PLCγ activation is required for macrophage migration and proliferation and hypothesize that this VEGFR1-PLCγ-mediated migration involves Ca2+ signaling. PLCγ phosphorylation is known to activate Ca2+ influx63,64 in oscillatory bursts.6568 Furthermore, directed cell migration requires Ca2+ pulses near the leading edge of the cell.6971 From this prior knowledge, combined with the delta function-like PLCγ activation we observe here, we hypothesize that VEGFR1 phosphorylates PLCγ in quick bursts to induce Ca2+ pulses and direct cell migration. This burst PLCγ activation could explain how cells migrate toward a VEGF gradient, with a possible mechanism being as follows: (1) VEGF binds plasma membrane VEGFR1 on the cell facing the gradient; (2) VEGFR1 recruits and phosphorylates PLCγ; and (3) phosphorylated PLCγ causes Ca2+ pulses by activating Ca2+ channels, a well-established mechanism7274 reviewed by Mikoshiba,75 initiating migration toward the VEGF gradient. This mechanism is further supported by experimental data showing that Ca2+ pulse following VEGF simulation is required for HUVEC migration.73 As the extent of directed cell migration is dependent on growth factor gradient patterns,76 we hypothesize that VEGFR1-PLCγ activation acts as a VEGF gradient sensor to determine both cell migration direction and magnitude. Additionally, PLCγ-Ca2+ signaling promotes cell proliferation through downstream activation of protein kinase C44. These studies showcase that cell migration and proliferation mechanisms through PLCγ-Ca2+ signaling have been established. Future work experimentally probing PLCγ-Ca2+ signaling through VEGFR1 to mediate cell migration and proliferation is necessary to validate our VEGFR1-PLCγ-Ca2+ signaling hypothesis.

Ca2+ signaling may indirectly regulate PI3K activation by VEGFR1

We identified PI3K as a primary adapter directing VEGFR1-mediated macrophage migration. Primarily, PI3K is known to promote cell migration through Akt activation,77,78 which also involves Ca2+ signaling; PI3K/Akt activation translocates Ca2+ channels to the cell membrane, inducing Ca2+ entry into cells, and subsequent cell migration.79 However, PI3K activation does not induce Ca2+ signaling in HUVECs;80 rather, PI3K is activated by Ca2+ to promote HUVEC migration.81 Since VEGF–VEGFR–adapter phosphorylation data from HUVECs was used to calibrate our computational model, and validated in macrophages, this indicates that VEGFR1-PI3K may play an important role in indirectly activating Ca2+ signaling in macrophages.

The PLCγ-, PI3K-, and Src-dependent relationship may form a Ca2+ signaling regulatory loop

We observed a dependent relationship between VEGFR1-induced PLCγ, PI3K, and Src phosphorylation. As PI3K and PLCγ cooperate to initiate Ca2+ signaling,82 we hypothesize that PI3K, PLCγ, and Ca2+ have a dependent relationship to robustly mediate VEGFR1-induced cell migration and proliferation. Furthermore, PLCγ-induced Ca2+ signaling phosphorylates Src,83 and Src phosphorylates PLCγ 80,83,84 and PI3K.8587 Thus we hypothesize from these studies and our results that VEGFR1 is structured to preferentially activate a PLCγ, PI3K, and Src regulatory loop mediating Ca2+ signaling (Fig. 6) and subsequent cell migration and proliferation.

Fig. 6
VEGFR1 preferentially activates PLCγ, PI3K, and Src, possibly to form a Ca2+ signaling regulatory loop. Our simulations predict that VEGFR1 tyrosine sites are structured to preferentially associate with PLCγ or PI3K at Tyr794 and Src at ...

VEGFR1-promoted hematopoietic progenitor cell migration may be required for tumor cell metastasis

The strong VEGFR1 migratory signal we identify here indicates that VEGFR1 signaling may be required for hematopoietic progenitor cell (HPC) migration to form premetastatic niche clusters. Metastasis from the primary tumor site requires circulating tumor cells to extravaste into secondary sites.88 Prior to this process, the tumor primes premetastatic niches, sites receptive to recruiting circulating tumor cells, to direct at which secondary sites metastasis occurs.89 These premetastatic niches are characterized by clustering of VEGFR1-positive HPCs; inhibiting VEGFR1 on HPCs prevents premetastatic niche formation and tumor cell metastasis.90 This effect of premetastatic niche formation being prevented with VEGFR1 inhibition may be explained by HPC migration requiring VEGFR1 signaling; thus inhibiting VEGFR1 signaling would prevent HPC migration, HPC clustering, and subsequent tumor cell metastasis. Furthermore, Akt activation has been implicated in macrophage-assisted cancer cell invasion,91 supporting our hypothesis that VEGFR1-PI3K-Ca2+ signaling (Fig. 6) promotes macrophage migration. Therefore, targeting VEGFR1-induced HPC migration may be a therapeutic option to prevent tumor cell metastasis.

Receptor signaling can be comprehensively quantified by modeling adapter–adapter interactions and specific phosphatases

Our modeling approach accurately and quantitatively predicted adapter phosphorylation and cell responses through complex formation between specific VEGFR1 tyrosine sites and single adapters, with adapter dephosphorylation occurring through a generalized phosphatase. Building upon this validated model to include adapter–adapter interactions and specific phosphatases would comprehensively represent VEGFR1 signaling. Modeling adapter–adapter interactions would identify how VEGFR1 signaling is directed through adapter cooperativity; adapter–adapter interactions occur via adapter SH3 domains92 to form larger signaling complexes that direct differential cell outcomes.93,94 Our ability to accurately model multi-adapter complex formation with VEGFR1 is currently limited, as no known experimental or computational studies have mapped the adapter–adapter interactions downstream of VEGFR1. This limitation may be overcome by identifying VEGFR1-associated adapter–adapter interactions from VEGF-induced protein phosphorylation dynamics, a predictive approach validated with the epidermal growth factor receptor signaling axis.95

Modeling specific phosphatases would identify additional VEGFR1-targeting therapeutics; since different phosphatases bind specific adapters to dynamically regulate receptor signaling,96 VEGFR1-induced adapter phosphorylation and cell responses could be directed by targeting specific phosphatases. Our ability to model specific phosphatases is currently limited however, as the specific phosphatases involved in VEGFR1 signaling, and their adapter interaction kinetics, have not been determined. This limitation may be overcome using the high-throughput assay for identifying phosphoprotein-specific phosphatases and kinetics developed by the Janes' laboratory.62 Overall, incorporating adapter–adapter interactions and phosphatase specificity into our VEGFR1 model would provide further insight into how VEGFR1 signaling is directed systemically and identify additional proteins or interactions that can be targeted to tune VEGFR1 signaling.

Conclusions

Our modeling approach has identified that VEGFR1 actively promotes macrophage migration and proliferation primarily via the PLCγ and PI3K pathways and has posited a new hypothesis that adapter coordination and Ca2+ signaling may regulate this VEGFR1-mediated migratory response. These findings critically advance our understanding of VEGF signaling by providing a structurally based mechanism for VEGFR1 function. Our findings and our modeling platform also offer mechanistic guidance for developing therapeutics targeting VEGFR1 signaling. This also represents a paradigm shift, since VEGF, generally, and VEGFR2 are primary targets for drug discovery. This modeling approach provides a foundation to fully understand signaling mechanisms for any receptor, an essential step to develop effective therapeutics for a wealth of pathologies.

Methods

Here we provide a brief overview of the computational and experimental methods. Interested readers can find complete description of the computational model, including parameters, and experimental procedures, including vendor information, in the Supplementary Information.

Computational models

VEGFR–adapter interaction models are defined by mass-action kinetics using ordinary differential equations and solved with the SimBiology toolbox in MATLAB. In general, the VEGFR–adapter scheme interaction scheme follows:

VEGF + VEGFRkoffVEGFkonVEGFpVEGFRpVEGFR + AkoffAkonApVEGFR : ApVEGFR : AkppVEGFR : pApVEGFR : pA+PTPNkoffPTPNkonPTPNpVEGFR : pA : PTPNpVEGFR : pA : PTPNkdPTPNpVEGFR + [A : PTPN]

for each adapter A and both VEGFRs, where PTPN are phosphatases. Model-predicted adapter phosphorylation in HUVECs shows good agreement to previous experimental data (SI Fig S2). VEGFR1 and VEGFR2 are both modeled for this validation (Fig. 1), as HUVECs express both receptors. Following this validation, we examine adapter–VEGFR1 interactions specifically to determine the VEGFR1 function. See SI Materials and Methods for details.

Code availability

Model code is available upon request.

Model assumptions

The following assumptions are used for model development—each assumption is described in SI Materials and Methods

Protein concentrations

(1) HUVEC protein concentrations are determined by western blot intensity, relative to a known protein concentration, assuming a linear relationship between protein band intensities (SI Appendix, Table S1). (2) PTPN acts as an “infinite reservoir”; the PTPN concentration is sufficiently high to not be a limiting species in any reaction.

Kinetics parameters

(1) Each adapter has the same interaction kinetics (on-rate and off-rate) for both VEGFR1 and VEGFR2 and is the same for all tyrosine sites (SI Table S2). (2) Adapter–VEGFR interaction kinetics are identical to adapter–EGFR interaction kinetics. (3) If adapter–VEGFR or adapter–EGFR interaction rates are unavailable, we assume that the rates between the SH2 domain of the adapter and a phosphorylated tyrosine kinase fragment is identical to the adapter–VEGFR rates. (4) We assume a 1 pL cell volume, to convert rates from M to molecules/cell.

Adapter phosphorylation

(1) All adapter phosphorylation rates (kp) are 0.01/s, so adapter phosphorylation is only dependent on VEGFR interaction kinetics. (2) Adapters do not undergo auto-dephosphorylation and are only dephosphorylated by phosphatases. (3) A generalized phosphatase (PTPN) binds and dephosphorylates all adapters, with the same interaction kinetics and dephosphorylation rate.

Predicting cell response from adapter phosphorylation

(1) The degradation cell response is identical to c-Cbl phosphorylation; only c-Cbl contributes to a degradation cell response. (2) Proliferation and migration cell responses are determined by a weighed sum of adapter phosphorylation. (3) Weights are calculated by the contribution each adapter provides toward the specific cell response, as determined experimentally (SI Table S3).

Tyrosine site specificity

(1) Multiple adapters can bind a single receptor if the combined size of the adapters is smaller than the available space between tyrosine sites (SI Table S4-S5). (2) Adapters bind the receptor in one-dimension (the y-direction). (3) Total adapter sizes are determined by measuring the maximal space the adapter crystal structure occupies in the y-direction. (4) The center of an adapter binds a VEGFR tyrosine site; thus, the amount of space a receptor occupies between VEGFR tyrosine sites is half the total adapter size. (5) We measure the average distance between VEGFR amino acids and use that distance to determine the space between VEGFR tyrosine sites. For example, the distance between individual amino acids in VEGFR1 was measured as 0.171 Å/amino acid, so the distance between tyrosine sites Tyr1242 and Tyr1333 is 15.6 Å.

Experimental methods

Experiments were performed in RAW 264.7 macrophages due to their high VEGFR1 expression (SI Fig S3), making them an ideal cell line to study VEGFR1 signaling.

Reagents and cell culture

Murine RAW 264.7 macrophages were cultured in Dulbecco’s Modified Eagle’s Medium supplemented with 10% fetal bovine serum and 1% penicillin–streptomycin. Cells were maintained in a humidified incubator at 37 °C and 5% CO2. Murine VEGF-A164 was purchased from BioLegend, and all inhibitors (Wortmannin, U73122, and Imatinib Mesylate) were purchased from Selleckchem. Enzyme-linked immunosorbent assay (ELISA) kits were purchased from Assay Biotechnology. The 3-[4,5-dimethylthiazol-2-yl]-2,5 diphenyl tetrazolium bromide (MTT) Cell Proliferation Assay Kit was purchased from Thermo Fisher Scientific.

Quantifying protein phosphorylation

RAWs were seeded into a 96-well plate, stimulated with VEGF or any inhibitors for specified times, and the phosphorylated and total proteins of interest (PLCγ, PI3K, and Abl) were measured using ELISAs. See SI Materials and Methods for details.

Cell migration assays

RAWs were seeded into a 12-well plate, scratched with a pipette tip, treated with VEGF or any inhibitors, and imaged at 0 and 24 h to characterize migration. See SI Materials and Methods for details.

Cell proliferation assays

RAWS were seeded into a 96-well plate, stimulated with VEGF or any inhibitors, and cell proliferation was measured after 24 h using a MTT assay. See SI Materials and Methods for details.

Flow cytometry

RAWs were labeled with Phycoerythrin (PE)-conjugated monoclonal antibodies specific to VEGFR1 or VEGFR2. Fluorescence given off by PE was captured in flow cytometry and converted to VEGFR level per cell (SI Fig S3). See SI Materials and Methods for details.

Data availability

Supplementary information includes a more detailed description of computational and experimental methods. Further data are available upon request.

Electronic supplementary material

Acknowledgements

This work was supported by the American Cancer Society (grant 282802) and the American Heart Association (grant 16SDG26940002). We would like to acknowledge Mathew Yang, Caroline Blassick, and Sunil Kumar for helping conduct literature searches to obtain initial concentrations, rates, and adapter-VEGFR1-binding motifs. We would also like to acknowledge Spencer Mamer and Ali Ansari for many helpful discussions with model building and data representation. This work was supported in part by grants from the American Cancer Society (282802), National Science Foundation (1512598), and American Heart Association (16SDG26940002).

Author contributions

Author contributions

J.C.W. and P.I.I. designed the study; J.C.W. developed the computational model, conducted model simulations, and performed the validation studies; S.C. performed the flow cytometric experiments and quantified receptor concentrations; all the authors analyzed the results and prepared the manuscript; P.I.I. is the guarantor.

Notes

Competing interests

The authors declare no competing financial interests.

Footnotes

Electronic supplementary material

Supplementary information accompanies the paper on the npj Systems Biology and Applications website (10.1038/s41540-017-0037-9).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References

1. Carmeliet P, Jain RK. Angiogenesis in cancer and other diseases. Nature. 2000;407:249–257. doi: 10.1038/35025220. [PubMed] [Cross Ref]
2. Hughes GC, Annex BH. Angiogenic therapy for coronary artery and peripheral arterial disease. Expert Rev. Cardiovasc. Ther. 2005;3:521–535. doi: 10.1586/14779072.3.3.521. [PubMed] [Cross Ref]
3. NHLBI. What Is Peripheral Artery Disease? (2016). Available at: https://www.nhlbi.nih.gov/health/healthtopics/topics/pad.
4. Bergers G, Benjamin LE. Tumorigenesis and the angiogenic switch. Nat. Rev. Cancer. 2003;3:401–410. doi: 10.1038/nrc1093. [PubMed] [Cross Ref]
5. Claesson-Welsh L, Welsh M. VEGFA and tumour angiogenesis. J. Intern. Med. 2013;273:114–127. doi: 10.1111/joim.12019. [PubMed] [Cross Ref]
6. Meadows, K. L. & Hurwitz, H. I. Anti-VEGF therapies in the clinic. Cold Spring Harb. Perspect. Med. 2, a006577–a006577 (2012). [PMC free article] [PubMed]
7. Carmeliet P, Jain RK. Molecular mechanisms and clinical applications of angiogenesis. Nature. 2011;473:298–307. doi: 10.1038/nature10144. [PMC free article] [PubMed] [Cross Ref]
8. Weddell JC, Imoukhuede PI. Quantitative characterization of cellular membrane-receptor heterogeneity through statistical and computational modeling. PLoS ONE. 2014;9:e97271. doi: 10.1371/journal.pone.0097271. [PMC free article] [PubMed] [Cross Ref]
9. Duda DG, et al. Plasma soluble VEGFR-1 is a potential dual biomarker of response and toxicity for bevacizumab with chemoradiation in locally advanced rectal cancer. Oncologist. 2010;15:577–583. doi: 10.1634/theoncologist.2010-0029. [PMC free article] [PubMed] [Cross Ref]
10. Lambrechts D, Lenz HJ, de Haas S, Carmeliet P, Scherer SJ. Markers of response for the antiangiogenic agent bevacizumab. J. Clin. Oncol. 2013;31:1219–1230. doi: 10.1200/JCO.2012.46.2762. [PubMed] [Cross Ref]
11. Shibuya M. Vascular endothelial growth factor receptor-1 (VEGFR-1/Flt-1): a dual regulator for angiogenesis. Angiogenesis. 2006;9:225–230. doi: 10.1007/s10456-006-9055-8. [PubMed] [Cross Ref]
12. Imoukhuede PI, Dokun AO, Annex BH, Popel AS. Endothelial cell-by-cell profiling reveals the temporal dynamics of VEGFR1 and VEGFR2 membrane localization after murine hindlimb ischemia. Am. J. Physiol. Heart Circ. Physiol. 2013;304:H1085–H1093. doi: 10.1152/ajpheart.00514.2012. [PubMed] [Cross Ref]
13. Imoukhuede PI, Popel AS. Quantitative fluorescent profiling of VEGFRs reveals tumor cell and endothelial cell heterogeneity in breast cancer xenografts. Cancer Med. 2014;3:225–244. doi: 10.1002/cam4.188. [PMC free article] [PubMed] [Cross Ref]
14. Carmeliet P, et al. Synergism between vascular endothelial growth factor and placental growth factor contributes to angiogenesis and plasma extravasation in pathological conditions. Nat. Med. 2001;7:575–583. doi: 10.1038/87904. [PubMed] [Cross Ref]
15. Fischer C, Mazzone M, Jonckx B, Carmeliet P. FLT1 and its ligands VEGFB and PlGF: drug targets for anti-angiogenic therapy? Nat. Rev. Cancer. 2008;8:942–956. doi: 10.1038/nrc2524. [PubMed] [Cross Ref]
16. Adini A, Kornaga T, Firoozbakht F, Benjamin LE. Placental growth factor is a survival factor for tumor endothelial cells and macrophages. Cancer Res. 2002;62:2749 LP–2752. [PubMed]
17. Fischer C, et al. Anti-PlGF inhibits growth of VEGF(R)-inhibitor-resistant tumors without affecting healthy vessels. Cell. 2007;131:463–475. doi: 10.1016/j.cell.2007.08.038. [PubMed] [Cross Ref]
18. Clegg LW, Mac Gabhann F. Site-specific phosphorylation of VEGFR2 is mediated by receptor trafficking: insights from a computational model. PLoS Comput. Biol. 2015;11:e1004158. doi: 10.1371/journal.pcbi.1004158. [PMC free article] [PubMed] [Cross Ref]
19. Matsumoto T, et al. VEGF receptor-2 Y951 signaling and a role for the adapter molecule TSAd in tumor angiogenesis. EMBO J. 2005;24:2342–2353. doi: 10.1038/sj.emboj.7600709. [PubMed] [Cross Ref]
20. Tong J, Taylor P, Peterman SM, Prakash A, Moran MF. Epidermal growth factor receptor phosphorylation sites Ser(991) and Tyr(998) are implicated in the regulation of receptor endocytosis and phosphorylations at Ser(1039) and Thr(1041) Mol. Cell. Proteomics. 2009;8:2131–2144. doi: 10.1074/mcp.M900148-MCP200. [PMC free article] [PubMed] [Cross Ref]
21. Tchaikovski V, Fellbrich G, Waltenberger J. The molecular basis of VEGFR-1 signal transduction pathways in primary human monocytes. Arterioscler. Thromb. Vasc. Biol. 2008;28:322–328. doi: 10.1161/ATVBAHA.107.158022. [PubMed] [Cross Ref]
22. Flynn DC. Adaptor proteins. Oncogene. 2001;20:6270–6272. doi: 10.1038/sj.onc.1204769. [PubMed] [Cross Ref]
23. Hunter T. Tyrosine phosphorylation: thirty years and counting. Curr. Opin. Cell Biol. 2009;21:140–146. doi: 10.1016/j.ceb.2009.01.028. [PMC free article] [PubMed] [Cross Ref]
24. Dephoure N, Gould KL, Gygi SP, Kellogg DR. Mapping and analysis of phosphorylation sites: a quick guide for cell biologists. Mol. Biol. Cell. 2013;24:535–542. doi: 10.1091/mbc.E12-09-0677. [PMC free article] [PubMed] [Cross Ref]
25. Duan H, Qu L, Shou C. Activation of EGFR-PI3K-AKT signaling is required for Mycoplasma hyorhinis-promoted gastric cancer cell migration. Cancer Cell Int. 2014;14:135. doi: 10.1186/s12935-014-0135-3. [PMC free article] [PubMed] [Cross Ref]
26. Germain RN, Meier-Schellersheim M, Nita-Lazar A, Fraser IDC. Systems biology in immunology – a computational modeling perspective. Annu. Rev. Immunol. 2011;29:527–585. doi: 10.1146/annurev-immunol-030409-101317. [PMC free article] [PubMed] [Cross Ref]
27. Williams AR, Timmis J, Qwarnstrom EE. Computational models of the NF-KB signalling pathway. Computation. 2014;2:131–158. doi: 10.3390/computation2040131. [Cross Ref]
28. Hogg, R. V., Tanis, E. A. & Zimmerman, D. Probability and Statistical Inference. (Pearson, Upper Saddle River, NJ, 2015).
29. Manning G, Whyte DB, Martinez R, Hunter T, Sudarsanam S. The protein kinase complement of the human genome. Science. 2002;298:1912–1934. doi: 10.1126/science.1075762. [PubMed] [Cross Ref]
30. Kholodenko BN, Demin OV, Moehren G, Hoek JB. Quantification of short term signaling by the epidermal growth factor receptor. J. Biol. Chem. 1999;274:30169–30181. doi: 10.1074/jbc.274.42.30169. [PubMed] [Cross Ref]
31. Mathew S, Banerjee I. Quantitative analysis of robustness of dynamic response and signal transfer in insulin mediated PI3K/Akt pathway. Comput. Chem. Eng. 2014;71:715–727. doi: 10.1016/j.compchemeng.2014.07.018. [PMC free article] [PubMed] [Cross Ref]
32. Tan WH, Popel AS, Mac Gabhann F. Computational model of Gab1/2-dependent VEGFR2 pathway to Akt activation. PLoS ONE. 2013;8:e67438. doi: 10.1371/journal.pone.0067438. [PMC free article] [PubMed] [Cross Ref]
33. Hsieh M, Yang S, Raymond-Stinz M, Edwards JS, Wilson BS. Spatio-temporal modeling of signaling protein recruitment to EGFR. BMC Syst. Biol. 2010;4:1–19. doi: 10.1186/1752-0509-4-57. [PMC free article] [PubMed] [Cross Ref]
34. Peirce SM, Van Gieson EJ, Skalak EC. Multicellular simulation predicts microvascular patterning and in silico tissue assembly. FASEB J. 2004;18:731–733. [PubMed]
35. Vivas J, Garzón-Alvarado D, Cerrolaza M. Modeling cell adhesion and proliferation: a cellular-automata based approach. Adv. Model Simul. Eng. Sci. 2015;2:1–12. doi: 10.1186/s40323-015-0053-5. [PMC free article] [PubMed] [Cross Ref]
36. Tian D, Kreeger PK. Analysis of the quantitative balance between insulin-like growth factor (IGF)-1 ligand, receptor, and binding protein levels to predict cell sensitivity and therapeutic efficacy. BMC Syst. Biol. 2014;8:98. doi: 10.1186/s12918-014-0098-y. [PMC free article] [PubMed] [Cross Ref]
37. Finley SD, Angelikopoulos P, Koumoutsakos P, Popel AS. Pharmacokinetics of anti-VEGF agent aflibercept in cancer predicted by data-driven, molecular-detailed model. CPT Pharmacomet. Syst. Pharmacol. 2015;4:641–649. doi: 10.1002/psp4.12040. [PMC free article] [PubMed] [Cross Ref]
38. Finley SD, Chu L, Popel AS. Computational systems biology approaches to anti-angiogenic cancer therapeutics. Drug Discov. Today. 2014;20:187–197. doi: 10.1016/j.drudis.2014.09.026. [PMC free article] [PubMed] [Cross Ref]
39. Finley SD, Dhar M, Popel AS. Compartment model predicts VEGF secretion and investigates the effects of VEGF trap in tumor-bearing mice. Front. Oncol. 2013;3:196. doi: 10.3389/fonc.2013.00196. [PMC free article] [PubMed] [Cross Ref]
40. Walpole J, Papin JA, Peirce SM. Multiscale computational models of complex biological systems. Annu. Rev. Biomed. Eng. 2013;15:137–154. doi: 10.1146/annurev-bioeng-071811-150104. [PMC free article] [PubMed] [Cross Ref]
41. Yu JS, Bagheri N. Multi-class and multi-scale models of complex biological phenomena. Curr. Opin. Biotechnol. 2016;39:167–173. doi: 10.1016/j.copbio.2016.04.002. [PubMed] [Cross Ref]
42. Zhang Z, Neiva KG, Lingen MW, Ellis LM, Nor JE. VEGF-dependent tumor angiogenesis requires inverse and reciprocal regulation of VEGFR1 and VEGFR2. Cell Death. Differ. 2010;17:499–512. doi: 10.1038/cdd.2009.152. [PMC free article] [PubMed] [Cross Ref]
43. Jakobsson L, et al. Endothelial cells dynamically compete for the tip cell position during angiogenic sprouting. Nat. Cell Biol. 2010;12:943–953. doi: 10.1038/ncb2103. [PubMed] [Cross Ref]
44. Blanco R, Gerhardt H. VEGF and Notch in tip and stalk cell selection. Cold Spring Harb. Perspect. Med. 2013;3:1–20. doi: 10.1101/cshperspect.a006569. [PMC free article] [PubMed] [Cross Ref]
45. Ou JM, et al. Knockdown of VEGFR2 inhibits proliferation and induces apoptosis in hemangioma-derived endothelial cells. Eur. J. Histochem. 2014;58:2263. doi: 10.4081/ejh.2014.2263. [PMC free article] [PubMed] [Cross Ref]
46. Zhang E, et al. Roles of PI3K/Akt and c-Jun signaling pathways in human papillomavirus type 16 oncoprotein-induced HIF-1α, VEGF, and IL-8 expression and in vitro angiogenesis in non-small cell lung cancer cells. PLoS ONE. 2014;9:e103440. doi: 10.1371/journal.pone.0103440. [PMC free article] [PubMed] [Cross Ref]
47. Zhao WF, Wang J, Tony To SS. The phosphatidylinositol 3-kinase/Akt and c-Jun N-terminal kinase signaling in cancer: alliance or contradiction? (Review) Int. J. Oncol. 2015;47:429–436. [PubMed]
48. Gharibi B, Ghuman MS, Hughes FJ. Akt- and Erk-mediated regulation of proliferation and differentiation during PDGFRβ-induced MSC self-renewal. J. Cell. Mol. Med. 2012;16:2789–2801. doi: 10.1111/j.1582-4934.2012.01602.x. [PMC free article] [PubMed] [Cross Ref]
49. Eichmann A, Simons M. VEGF signaling inside vascular endothelial cells and beyond. Curr. Opin. Cell Biol. 2012;24:188–193. doi: 10.1016/j.ceb.2012.02.002. [PMC free article] [PubMed] [Cross Ref]
50. Abhinand CS, Raju R, Soumya SJ, Arya PS, Sudhakaran PR. VEGF-A/VEGFR2 signaling network in endothelial cells relevant to angiogenesis. J. Cell Commun. Signal. 2016;10:347–354. doi: 10.1007/s12079-016-0352-8. [PMC free article] [PubMed] [Cross Ref]
51. Koch S, Claesson-Welsh L. Signal transduction by vascular endothelial growth factor receptors. Cold Spring Harb. Perspect. Med. 2012;2:a006502. doi: 10.1101/cshperspect.a006502. [PMC free article] [PubMed] [Cross Ref]
52. Le Novere N. Quantitative and logic modelling of molecular and gene networks. Nat. Rev. Genet. 2015;16:146–158. doi: 10.1038/nrg3885. [PMC free article] [PubMed] [Cross Ref]
53. Brodland GW. How computational models can help unlock biological systems. Semin. Cell Dev. Biol. 2015;47–48:62–73. doi: 10.1016/j.semcdb.2015.07.001. [PubMed] [Cross Ref]
54. Janes KA, Lauffenburger DA. Models of signalling networks – what cell biologists can gain from them and give to them. J. Cell Sci. 2013;126:1913 LP–1921. doi: 10.1242/jcs.112045. [PubMed] [Cross Ref]
55. Pannu KK, et al. Performance evaluation of quantiBRITE phycoerythrin beads. Cytometry. 2001;45:250–258. doi: 10.1002/1097-0320(20011201)45:4<250::AID-CYTO10021>3.0.CO;2-T. [PubMed] [Cross Ref]
56. Wang L, Abbasi F, Gaigalas AK, Vogt RF, Marti GE. Comparison of fluorescein and phycoerythrin conjugates for quantifying CD20 expression on normal and leukemic B-cells. Cytom. B Clin. Cytom. 2006;70:410–415. doi: 10.1002/cyto.b.20140. [PubMed] [Cross Ref]
57. Yen P, Finley SD, Engel-Stefanini MO, Popel AS. A two-compartment model of VEGF distribution in the mouse. PLoS ONE. 2011;6:e27514. doi: 10.1371/journal.pone.0027514. [PMC free article] [PubMed] [Cross Ref]
58. Finley SD, Engel-Stefanini MO, Imoukhuede P, Popel AS. Pharmacokinetics and pharmacodynamics of VEGF-neutralizing antibodies. BMC Syst. Biol. 2011;5:193. doi: 10.1186/1752-0509-5-193. [PMC free article] [PubMed] [Cross Ref]
59. Finley SD, Popel AS. Predicting the effects of anti-angiogenic agents targeting specific VEGF isoforms. AAPS J. 2012;14:500–509. doi: 10.1208/s12248-012-9363-4. [PMC free article] [PubMed] [Cross Ref]
60. Tan WH, Popel AS, Mac Gabhann F. Computational model of VEGFR2 pathway to ERK activation and modulation through receptor trafficking. Cell Signal. 2013;25:2496–2510. doi: 10.1016/j.cellsig.2013.08.015. [PMC free article] [PubMed] [Cross Ref]
61. Ko CL, Voit E, Wang FS. Estimating parameters for generalized mass action models with connectivity information. BMC Bioinformatics. 2009;10:140. doi: 10.1186/1471-2105-10-140. [PMC free article] [PubMed] [Cross Ref]
62. Bose AK, Janes KA. A high-throughput assay for phosphoprotein-specific phosphatase activity in cellular extracts. Mol. Cell. Proteomics. 2013;12:797–806. doi: 10.1074/mcp.O112.024059. [PMC free article] [PubMed] [Cross Ref]
63. McLaughlin AP, De Vries GW. Role of PLCγ and Ca2+ in VEGF- and FGF-induced choroidal endothelial cell proliferation. Am. J. Physiol. Cell. Physiol. 2001;281:C1448–C1456. [PubMed]
64. Putney JW. PLC-gamma: an old player has a new role. Nat. Cell Biol. 2002;4:E280–E281. doi: 10.1038/ncb1202-e280. [PubMed] [Cross Ref]
65. Purvis JE, Chatterjee MS, Brass LF, Diamond SL. A molecular signaling model of platelet phosphoinositide and calcium regulation during homeostasis and P2Y1 activation. Blood. 2008;112:4069–4079. doi: 10.1182/blood-2008-05-157883. [PubMed] [Cross Ref]
66. Jovic A, et al. Phase-locked signals elucidate circuit architecture of an oscillatory pathway. PLoS Comput. Biol. 2010;6:e1001040. doi: 10.1371/journal.pcbi.1001040. [PMC free article] [PubMed] [Cross Ref]
67. Jovic A, Wade SM, Neubig RR, Linderman JJ, Takayama S. Microfluidic interrogation and mathematical modeling of multi-regime calcium signaling dynamics. Integr. Biol. 2013;5:932–939. doi: 10.1039/c3ib40032h. [PMC free article] [PubMed] [Cross Ref]
68. Dolan AT, Diamond SL. Systems modeling of Ca2+ homeostasis and mobilization in platelets mediated by IP3 and store-operated Ca2+ entry. Biophys. J. 2015;106:2049–2060. doi: 10.1016/j.bpj.2014.03.028. [PubMed] [Cross Ref]
69. Marks PW, Maxfield FR. Transient increases in cytosolic free calcium appear to be required for the migration of adherent human neutrophils. J. Cell Biol. 1990;110:43–52. doi: 10.1083/jcb.110.1.43. [PMC free article] [PubMed] [Cross Ref]
70. Wei C, et al. Calcium flickers steer cell migration. Nature. 2009;457:901–905. doi: 10.1038/nature07577. [PMC free article] [PubMed] [Cross Ref]
71. Minton K. Cell migration: coordinating calcium signalling. Nat. Rev. Mol. Cell Biol. 2014;15:152. doi: 10.1038/nrm3761. [PubMed] [Cross Ref]
72. Cronshaw DG, et al. Evidence that phospholipase C-dependent, calcium-independent mechanisms are required for directional migration of T lymphocytes in response to the CCR4 ligands CCL17 and CCL22. J. Leukoc. Biol. 2006;79:1369–1380. doi: 10.1189/jlb.0106035. [PubMed] [Cross Ref]
73. Noren DP, et al. Endothelial cells decode VEGF-mediated Ca2+ signaling patterns to produce distinct functional responses. Sci. Signal. 2016;9:ra20–ra20. doi: 10.1126/scisignal.aad3188. [PMC free article] [PubMed] [Cross Ref]
74. Tsai F, Kuo G, Chang S, Tsai P. Ca2+ signaling in cytoskeletal reorganization, cell migration, and cancer metastasis. Biomed. Res. Int. 2015;2015:13. [PMC free article] [PubMed]
75. Mikoshiba K. IP3 receptor/Ca2+ channel: from discovery to new signaling concepts. J. Neurochem. 2007;102:1426–1446. doi: 10.1111/j.1471-4159.2007.04825.x. [PubMed] [Cross Ref]
76. Stefonek-Puccinelli TJ, Masters KS. Co-immobilization of gradient-patterned growth factors for directed cell migration. Ann. Biomed. Eng. 2008;36:2121–2133. doi: 10.1007/s10439-008-9581-1. [PMC free article] [PubMed] [Cross Ref]
77. Vanhaesebroeck B, Stephens L, Hawkins P. PI3K signalling: the path to discovery and understanding. Nat. Rev. Mol. Cell Biol. 2012;13:195–203. doi: 10.1038/nrm3290. [PubMed] [Cross Ref]
78. Hemmings BA, Restuccia DF. PI3K-PKB/Akt pathway. Cold Spring Harb. Perspect. Biol. 2012;4:a011189. doi: 10.1101/cshperspect.a011189. [PMC free article] [PubMed] [Cross Ref]
79. Viard P, et al. PI3K promotes voltage-dependent calcium channel trafficking to the plasma membrane. Nat. Neurosci. 2004;7:939–946. doi: 10.1038/nn1300. [PubMed] [Cross Ref]
80. Dawson NS, Zawieja DC, Wu MH, Granger HJ. Signaling pathways mediating VEGF165-induced calcium transients and membrane depolarization in human endothelial cells. FASEB J. 2006;20:991–993. doi: 10.1096/fj.05-3923fje. [PubMed] [Cross Ref]
81. Tauzin S, et al. The naturally processed CD95L elicits a c-Yes/calcium/PI3K-driven cell migration pathway. PLoS Biol. 2011;9:e1001090. doi: 10.1371/journal.pbio.1001090. [PMC free article] [PubMed] [Cross Ref]
82. Rameh LE, et al. Phosphoinositide 3-kinase regulates phospholipase Cgamma-mediated calcium signaling. J. Biol. Chem. 1998;273:23750–23757. doi: 10.1074/jbc.273.37.23750. [PubMed] [Cross Ref]
83. Cattaneo MG, Lucci G, Vicentini LM. Oxytocin stimulates in vitro angiogenesis via a Pyk-2/Src-dependent mechanism. Exp. Cell Res. 2009;315:3210–3219. doi: 10.1016/j.yexcr.2009.06.022. [PubMed] [Cross Ref]
84. He H, et al. Vascular endothelial growth factor signals endothelial cell production of nitric oxide and prostacyclin through Flk-1/KDR activation of c-Src. J. Biol. Chem. 1999;274:25130–25135. doi: 10.1074/jbc.274.35.25130. [PubMed] [Cross Ref]
85. Sun M, et al. Activation of phosphatidylinositol 3-kinase/Akt pathway by androgen through interaction of p85α, androgen receptor, and Src. J. Biol. Chem. 2003;278:42992–43000. doi: 10.1074/jbc.M306295200. [PubMed] [Cross Ref]
86. Liu G, et al. ICAM-1–activated Src and eNOS signaling increase endothelial cell surface PECAM-1 adhesivity and neutrophil transmigration. Blood. 2012;120:1942–1952. doi: 10.1182/blood-2011-12-397430. [PubMed] [Cross Ref]
87. Maniatis NA, et al. Novel mechanism of endothelial nitric oxide synthase activation mediated by caveolae internalization in endothelial cells. Circ. Res. 2006;99:870–877. doi: 10.1161/01.RES.0000245187.08026.47. [PubMed] [Cross Ref]
88. Irmisch A, Huelsken J. Metastasis: new insights into organ-specific extravasation and metastatic niches. Exp. Cell Res. 2013;319:1604–1610. doi: 10.1016/j.yexcr.2013.02.012. [PubMed] [Cross Ref]
89. Quail DF, Joyce JA. Microenvironmental regulation of tumor progression and metastasis. Nat. Med. 2013;19:1423–1437. doi: 10.1038/nm.3394. [PMC free article] [PubMed] [Cross Ref]
90. Kaplan RN, et al. VEGFR1-positive haematopoietic bone marrow progenitors initiate the pre-metastatic niche. Nature. 2005;438:820–827. doi: 10.1038/nature04186. [PMC free article] [PubMed] [Cross Ref]
91. Park KY, Li G, Platt MO. Monocyte-derived macrophage assisted breast cancer cell invasion as a personalized, predictive metric to score metastatic risk. Sci. Rep. 2015;5:13855. doi: 10.1038/srep13855. [PMC free article] [PubMed] [Cross Ref]
92. Mayer B. SH3 domains: complexity in moderation. J. Cell Sci. 2001;114:1253–1263. [PubMed]
93. Hu Q, Milfay D, Williams LT. Binding of NCK to SOS and activation of ras-dependent gene expression. Mol. Cell. Biol. 1995;15:1169–1174. doi: 10.1128/MCB.15.3.1169. [PMC free article] [PubMed] [Cross Ref]
94. Lemmon MA, Schlessinger J. Cell signaling by receptor tyrosine kinases. Cell. 2010;141:1117–1134. doi: 10.1016/j.cell.2010.06.011. [PMC free article] [PubMed] [Cross Ref]
95. Naegle KM, White FM, Lauffenburger DA, Yaffe MB. Robust co-regulation of tyrosine phosphorylation sites on proteins reveals novel protein interactions. Mol. Biosyst. 2012;8:2771–2782. doi: 10.1039/c2mb25200g. [PMC free article] [PubMed] [Cross Ref]
96. Böhmer F, Szedlacsek S, Tabernero L, Östman A, den Hertog J. Protein tyrosine phosphatase structure–function relationships in regulation and pathogenesis. FEBS J. 2013;280:413–431. doi: 10.1111/j.1742-4658.2012.08655.x. [PubMed] [Cross Ref]

Articles from NPJ Systems Biology and Applications are provided here courtesy of Nature Publishing Group