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

Mechanistic computational modeling of sFLT1 secretion dynamics

Constitutively secreted by endothelial cells, soluble FLT1 (sFLT1 or sVEGFR1) binds and sequesters extracellular vascular endothelial growth factors (VEGF), thereby reducing VEGF binding to VEGF receptor tyrosine kinases and their downstream signaling. In doing so, sFLT1 plays an important role in vascular development and in the patterning of new blood vessels in angiogenesis. Here, we develop multiple mechanistic models of sFLT1 secretion and identify a minimal mechanistic model that recapitulates key qualitative and quantitative features of temporal experimental datasets of sFLT1 secretion from multiple studies. We show that the experimental data on sFLT1 secretion is best represented by a delay differential equation (DDE) system including a maturation term, reflecting the time required between synthesis and secretion. Using optimization to identify appropriate values for the key mechanistic parameters in the model, we show that two model parameters (extracellular degradation rate constant and maturation time) are very strongly constrained by the experimental data, and that the remaining parameters are related by two strongly constrained constants. Thus, only one degree of freedom remains, and measurements of the intracellular levels of sFLT1 would fix the remaining parameters. Comparison between simulation predictions and additional experimental data of the outcomes of chemical inhibitors and genetic perturbations suggest that intermediate values of the secretion rate constant best match the simulation with experiments, which would completely constrain the model. However, some of the inhibitors tested produce results that cannot be reproduced by the model simulations, suggesting that additional mechanisms not included here are required to explain those inhibitors. Overall, the model reproduces most available experimental data and suggests targets for further quantitative investigation of the sFLT1 system.

Author Summary

Proteins that are typically found outside cells are initially made inside cells, and later secreted into extracellular space. Many of these secreted proteins have important functions outside the cell that are well-studied; however, usually much less is known about the pre-secretion life of these molecules. Many computational models only represent the extracellular versions of secreted proteins, reducing all production and secretion steps into a single modeled process. Here, we develop a mechanistic model of the production and secretion of a specific secreted protein, sFLT1, which inhibits blood vessel growth by acting as an extracellular sponge for another set of secreted proteins, the vascular endothelial growth factors. We compare several models to existing experimentally-measured sFLT1 data, and we show that the data are most simply explained by including a delay between intracellular sFLT1 production and sFLT1 transport or degradation. This is consistent with the biology of the cell's secretory pathway, where immature proteins are gradually processed into mature forms over minutes to hours. Our approach could be incorporated into improved models for any pathway involving secreted proteins, including sFLT1-regulated models of blood vessel biology.

VEGF/VEGFR axis and its signaling in melanoma: Current knowledge toward therapeutic targeting agents and future perspectives

Melanoma is responsible for most skin cancer-associated deaths globally. The progression of melanoma is influenced by a number of pathogenic processes. Understanding the VEGF/VEGFR axis, which includes VEGF-A, PlGF, VEGF-B, VEGF-C, and VEGF-D and their receptors, VEGFR-1, VEGFR-2, and VEGFR-3, is of great importance in melanoma due to its crucial role in angiogenesis. This axis generates multifactorial and complex cellular signaling, engaging the MAPK/ERK, PI3K/AKT, PKC, PLC-γ, and FAK signaling pathways. Melanoma cell growth and proliferation, migration and metastasis, survival, and acquired resistance to therapy are influenced by this axis. The VEGF/VEGFR axis was extensively examined for their potential as diagnostic/prognostic biomarkers in melanoma patients and results showed that VEGF overexpression can be associated with unfavorable prognosis, higher level of tumor invasion and poor response to therapy. MicroRNAs linking to the VEGF/VEGFR axis were identified and, in this review, divided into two categories according to their functions, some of them promote melanoma angiogenesis (promotive group) and some restrict melanoma angiogenesis (protective group). In addition, the approach of treating melanoma by targeting the VEGF/VEGFR axis has garnered significant interest among researchers. These agents can be divided into two main groups: anti-VEGF and VEGFR inhibitors. These therapeutic options may be a prominent step along with the modern targeting and immune therapies for better coverage of pathological processes leading to melanoma progression and therapy resistance.

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.

Neuroprotective Effects of VEGF in the Enteric Nervous System

Although the enteric nervous system (ENS) functions largely autonomously as part of the peripheral nervous system (PNS), it is connected to the central nervous system (CNS) via the gut–brain axis. In many neurodegenerative diseases, pathological changes occur in addition to gastrointestinal symptoms, such as alpha-synuclein aggregates in Parkinson’s disease, which are found early in the ENS. In both the CNS and PNS, vascular endothelial growth factor (VEGF) mediates neuroprotective and neuroregenerative effects. Since the ENS with its close connection to the microbiome and the immune system is discussed as the origin of neurodegenerative diseases, it is necessary to investigate the possibly positive effects of VEGF on enteric neurons. Using laser microdissection and subsequent quantitative RT-PCR as well as immunohistochemistry, for the first time we were able to detect and localize VEGF receptor expression in rat myenteric neurons of different ages. Furthermore, we demonstrate direct neuroprotective effects of VEGF in the ENS in cell cultures. Thus, our results suggest a promising approach regarding neuroprotection, as the use of VEGF (may) prevent neuronal damage in the ENS.

Systems biology of angiogenesis signaling: Computational models and omics

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

This article is categorized under:

Cardiovascular Diseases > Computational Models

Cancer > Computational Models

VEGFR-1/Flt-1 inhibition increases angiogenesis and improves muscle function in a mouse model of Duchenne muscular dystrophy

Duchenne muscular dystrophy is characterized by structural degeneration of muscle, which is exacerbated by localized functional ischemia due to loss of nitric oxide synthase-induced vasodilation. Treatment strategies aimed at increasing vascular perfusion have been proposed. Toward this end, we have developed monoclonal antibodies (mAbs) that bind to the vascular endothelial growth factor (VEGF) receptor VEGFR-1 (Flt-1) and its soluble splice variant isoform (sFlt-1) leading to increased levels of free VEGF and proangiogenic signaling. The lead chimeric mAb, 21B3, had high affinity and specificity for both human and mouse sFlt-1 and inhibited VEGF binding to sFlt-1 in a competitive manner. Proof-of-concept studies in the mdx mouse model of Duchenne muscular dystrophy showed that intravenous administration of 21B3 led to elevated VEGF levels, increased vascularization and blood flow to muscles, and decreased fibrosis after 6-12 weeks of treatment. Greater muscle strength was also observed after 4 weeks of treatment. A humanized form of the mAb, 27H6, was engineered and demonstrated a comparable pharmacologic effect. Overall, administration of anti-Flt-1 mAbs in mdx mice inhibited the VEGF:Flt-1 interaction, promoted angiogenesis, and improved muscle function. These studies suggest a potential therapeutic benefit of Flt-1 inhibition for patients with Duchenne muscular dystrophy.

Disabling VEGF-Response of Purkinje Cells by Downregulation of KDR via miRNA-204-5p

The vascular endothelial growth factor (VEGF) is well known for its wide-ranging functions, not only in the vascular system, but also in the central (CNS) and peripheral nervous system (PNS). To study the role of VEGF in neuronal protection, growth and maturation processes have recently attracted much interest. These effects are mainly mediated by VEGF receptor 2 (VEGFR-2). Current studies have shown the age-dependent expression of VEGFR-2 in Purkinje cells (PC), promoting dendritogenesis in neonatal, but not in mature stages. We hypothesize that microRNAs (miRNA/miR) might be involved in the regulation of VEGFR-2 expression during the development of PC. In preliminary studies, we performed a miRNA profiling and identified miR204-5p as a potential regulator of VEGFR-2 expression. In the recent study, organotypic slice cultures of rat cerebella (postnatal day (p) 1 and 9) were cultivated and VEGFR-2 expression in PC was verified via immunohistochemistry. Additionally, PC at age p9 and p30 were isolated from cryosections by laser microdissection (LMD) to analyse VEGFR-2 expression by quantitative RT-PCR. To investigate the influence of miR204-5p on VEGFR-2 levels in PC, synthetic constructs including short hairpin (sh)-miR204-5p cassettes (miRNA-mimics), were microinjected into PC. The effects were analysed by confocal laser scanning microscopy (CLSM) and morphometric analysis. For the first time, we could show that miR204-5p has a negative effect on VEGF sensitivity in juvenile PC, resulting in a significant decrease of dendritic growth compared to untreated juvenile PC. In mature PC, the overexpression of miR204-5p leads to a shrinkage of dendrites despite VEGF treatment. The results of this study illustrate, for the first time, which miR204-5p expression has the potential to play a key role in cerebellar development by inhibiting VEGFR-2 expression in PC.

Injured Axons Instruct Schwann Cells to Build Constricting Actin Spheres to Accelerate Axonal Disintegration

After a peripheral nerve lesion, distal ends of injured axons disintegrate into small fragments that are subsequently cleared by Schwann cells and later by macrophages. Axonal debris clearing is an early step of the repair process that facilitates regeneration. We show here that Schwann cells promote distal cut axon disintegration for timely clearing. By combining cell-based and in vivo models of nerve lesion with mouse genetics, we show that this mechanism is induced by distal cut axons, which signal to Schwann cells through PlGF mediating the activation and upregulation of VEGFR1 in Schwann cells. In turn, VEGFR1 activates Pak1, leading to the formation of constricting actomyosin spheres along unfragmented distal cut axons to mediate their disintegration. Interestingly, oligodendrocytes can acquire a similar behavior as Schwann cells by enforced expression of VEGFR1. These results thus identify controllable molecular cues of a neuron-glia crosstalk essential for timely clearing of damaged axons.

Systems Biology Will Direct Vascular-Targeted Therapy for Obesity

Healthy adipose tissue expansion and metabolism during weight gain require coordinated angiogenesis and lymphangiogenesis. These vascular growth processes rely on the vascular endothelial growth factor (VEGF) family of ligands and receptors (VEGFRs). Several studies have shown that controlling vascular growth by regulating VEGF:VEGFR signaling can be beneficial for treating obesity; however, dysregulated angiogenesis and lymphangiogenesis are associated with several chronic tissue inflammation symptoms, including hypoxia, immune cell accumulation, and fibrosis, leading to obesity-related metabolic disorders. An ideal obesity treatment should minimize adipose tissue expansion and the advent of adverse metabolic consequences, which could be achieved by normalizing VEGF:VEGFR signaling. Toward this goal, a systematic investigation of the interdependency of vascular and metabolic systems in obesity and tools to predict personalized treatment ranges are necessary to improve patient outcomes through vascular-targeted therapies. Systems biology can identify the critical VEGF:VEGFR signaling mechanisms that can be targeted to regress adipose tissue expansion and can predict the metabolic consequences of different vascular-targeted approaches. Establishing a predictive, biologically faithful platform requires appropriate computational models and quantitative tissue-specific data. Here, we discuss the involvement of VEGF:VEGFR signaling in angiogenesis, lymphangiogenesis, adipogenesis, and macrophage specification – key mechanisms that regulate adipose tissue expansion and metabolism. We then provide useful computational approaches for simulating these mechanisms, and detail quantitative techniques for acquiring tissue-specific parameters. Systems biology, through computational models and quantitative data, will enable an accurate representation of obese adipose tissue that can be used to direct the development of vascular-targeted therapies for obesity and associated metabolic disorders.

Multifaceted Role of the Placental Growth Factor (PlGF) in the Antitumor Immune Response and Cancer Progression

The sharing of molecules function that affects both tumor growth and neoangiogenesis with cells of the immune system creates a mutual interplay that impairs the host’s immune response against tumor progression. Increasing evidence shows that tumors are able to create an immunosuppressive microenvironment by recruiting specific immune cells. Moreover, molecules produced by tumor and inflammatory cells in the tumor microenvironment create an immunosuppressive milieu able to inhibit the development of an efficient immune response against cancer cells and thus fostering tumor growth and progression. In addition, the immunoediting could select cancer cells that are less immunogenic or more resistant to lysis. In this review, we summarize recent findings regarding the immunomodulatory effects and cancer progression of the angiogenic growth factor namely placental growth factor (PlGF) and address the biological complex effects of this cytokine. Different pathways of the innate and adaptive immune response in which, directly or indirectly, PlGF is involved in promoting tumor immune escape and metastasis will be described. PlGF is important for building up vascular structures and functions. Although PlGF effects on vascular and tumor growth have been widely summarized, its functions in modulating the immune intra-tumoral microenvironment have been less highlighted. In agreement with PlGF functions, different antitumor strategies can be envisioned.

Context-dependent regulation of receptor tyrosine kinases: Insights from systems biology approaches

Receptor tyrosine kinases (RTKs) are cell membrane proteins that provide cells with the ability to sense proteins in their environments. Many RTKs are essential to development and organ growth. Derangement of RTKs-by mutation or by overexpression-is central to several developmental and adult disorders including cancer, short stature, and vascular pathologies. The mechanism of action of RTKs is complex and is regulated by contextual components, including the existence of multiple competing ligands and receptors in many families, the intracellular location of the RTK, the dynamic and cell-specific coexpression of other RTKs, and the commonality of downstream signaling pathways. This means that both the state of the cell and the microenvironment outside the cell play a role, which makes sense given the pivotal location of RTKs as the nexus linking the extracellular milieu to intracellular signaling and modification of cell behavior. In this review, we describe these different contextual components through the lens of systems biology, in which both computational modeling and experimental "omics" approaches have been used to better understand RTK networks. The complexity of these networks is such that using these systems biology approaches is necessary to get a handle on the mechanisms of pathology and the design of therapeutics targeting RTKs. In particular, we describe in detail three concrete examples (involving ErbB3, VEGFR2, and AXL) that illustrate how systems approaches can reveal key mechanistic and therapeutic insights. This article is categorized under: Biological Mechanisms > Cell Signaling Models of Systems Properties and Processes > Mechanistic Models Translational, Genomic, and Systems Medicine > Therapeutic Methods.

Cell migration and growth induced by photo-immobilised vascular endothelial growth factor (VEGF) isoforms

VEGF isoforms immobilised by photo-reactive gelatin (AzPhe-gelatin) enhance cell migration and proliferation.

Apatinib for salvage treatment of advanced malignant pleural mesothelioma: A case report

RATIONALE

Malignant Pleural Mesothelioma (MPM) is rare cancer and has a poor prognosis with resistance to chemotherapy or radiotherapy. Until now there is no standard third-line treatment for patients who have failed second-line therapy.

PATIENT CONCERNS

A 58-year-old non-smoking female peasant of ethnic Han was admitted to the oncology department of the 363 Hospital with a primary complaint of chest tightness and breathlessness from 3 months ago.

DIAGNOSES

Positron emission tomography-computed tomography (PET/CT) examination showed "dirty" pleural and parietal pleural involvement as well as mediastinal and pulmonary hilar lymph node enlargement. Finally, cancer cells were seen after repeated pleural effusion cell examination. Immunohistochemistry confirmed epithelioid of pleural mesothelioma.

INTERVENTIONS

Apatinib as a third-line treatment after failure from pemetrexed/cisplatin (PC) as the first-line chemotherapy and gemcitabine/cisplatin (GP) as the second-line chemotherapy. At first, 250 mg/day was given and 1 week later, the dose was increased to 500 mg/day.

OUTCOMES

A 5-month progression-free survival was achieved and toxicity included severe hand-foot syndrome, mild proteinuria, and hypertension.

LESSONS

Apatinib may be a potential therapeutic drug for MPM, particularly as a third-line treatment in cases resistant to chemotherapeutic options.

The role of placental growth factor (PlGF) and its receptor system in retinal vascular diseases

Placental growth factor (PlGF) is a member of the vascular endothelial growth factor (VEGF) family. Upon binding to VEGF- and neuropilin-receptor sub-types, PlGF modulates a range of neural, glial and vascular cell responses that are distinct from VEGF-A. As PlGF expression is selectively associated with pathological angiogenesis and inflammation, its blockade does not affect the healthy vasculature. PlGF actions have been extensively described in tumor biology but more recently there has been accumulating preclinical evidence that indicates that this growth factor could have an important role in retinal diseases. High levels of PlGF have been found in aqueous humor, vitreous and/or retina of patients exhibiting retinopathies, especially those with diabetic retinopathy (DR) and neovascular age-related macular degeneration (nvAMD). Expression of this growth factor seems to correlate closely with many of the key pathogenic features of early and late retinopathy in preclinical models. For example, studies using genetic modification and/or pharmacological treatment to block PlGF in the laser-induced choroidal neovascularization (CNV) model, oxygen-induced retinopathy model, as well as various murine diabetic models, have shown that PlGF deletion or inhibition can reduce neovascularization, retinal leakage, inflammation and gliosis, without affecting vascular development or inducing neuronal degeneration. Moreover, an inhibitory effect of PlGF blockade on retinal scarring in the mouse CNV model has also been recently demonstrated and was found to be unique for PlGF inhibition, as compared to various VEGF inhibition strategies. Together, these preclinical results suggest that anti-PlGF therapy might have advantages over anti-VEGF treatment, and that it may have clinical applications as a standalone treatment or in combination with anti-VEGF. Additional clinical studies are clearly needed to further elucidate the role of PlGF and its potential as a therapeutic target in ocular diseases.

Adrenergic receptor blockade attenuates placental ischemia-induced hypertension

Preeclampsia (PE), a disorder of new-onset maternal hypertension and vascular dysfunction during pregnancy, is thought to be linked to placental ischemia-induced release of prohypertensive factors and reductions of vasoprotective factors in the maternal circulation. Although markers of sympathetic nervous activity are elevated in experimental models of placental ischemia-induced hypertension and women with PE compared with their normal pregnant counterparts, the importance of adrenergic receptor signaling in the development of hypertension in PE is unknown. Therefore, we tested the hypothesis that adrenergic receptor blockade attenuates the development of placental ischemia-induced hypertension in rats. Wistar Hannover rats underwent reduced uterine perfusion pressure (RUPP) or Sham surgeries on gestational day 14. By day 19, mean arterial blood pressure (MAP) was increased in RUPP over Sham rats. Groups of RUPP and Sham pregnant rats received terazosin and propranolol (3 mg/kg per day of each via subcutaneous osmotic minipump) to block α1- and β-adrenergic receptors, respectively, beginning on gestational day 14. Adrenergic blockade significantly attenuated the development of hypertension in the RUPP rats with a slight blood pressure-lowering response in the Sham, normal pregnant rats by day 19. In conclusion, these data implicate that placental ischemia-induced hypertension involves adrenergic receptor signaling to promote increases in blood pressure during PE.

Real-Time Ligand Binding of Fluorescent VEGF-A Isoforms that Discriminate between VEGFR2 and NRP1 in Living Cells

Fluorescent VEGF-A isoforms have been evaluated for their ability to discriminate between VEGFR2 and NRP1 in real-time ligand binding studies in live cells using BRET. To enable this, we synthesized single-site (N-terminal cysteine) labeled versions of VEGF₁₆₅a, VEGF₁₆₅b, and VEGF₁₂₁a. These were used in combination with N-terminal NanoLuc-tagged VEGFR2 or NRP1 to evaluate the selectivity of VEGF isoforms for these two membrane proteins. All fluorescent VEGF-A isoforms displayed high affinity for VEGFR2. Only VEGF₁₆₅a-TMR bound to NanoLuc-NRP1 with a similar high affinity (4.4 nM). Competition NRP1 binding experiments yielded a rank order of potency of VEGF₁₆₅a > VEGF₁₈₉a > VEGF₁₄₅a. VEGF₁₆₅b, VEGF-Ax, VEGF₁₂₁a, and VEGF₁₁₁a were unable to bind to NRP1. There were marked differences in the kinetic binding profiles of VEGF₁₆₅a-TMR for NRP1 and VEGFR2. These data emphasize the importance of the kinetic aspects of ligand binding to VEGFR2 and its co-receptors in the dynamics of VEGF signaling.

•

VEGF₁₆₅a, VEGF₁₂₁a, and VEGF₁₆₅b were single-site labeled with tetramethylrhodamine

•

NanoBRET quantified that VEGF-A isoforms have similar binding properties at VEGFR2

•

NRP1 expressed in live cells does not bind VEGF₁₆₅b, VEGF₁₂₁a, VEGF-Ax, or VEGF₁₁₁a

•

VEGFR2 and NRP1 have markedly distinct kinetic profiles binding VEGF₁₆₅a-TMR

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.

Molecular Pharmacology of VEGF-A Isoforms: Binding and Signalling at VEGFR2

Vascular endothelial growth factor-A (VEGF-A) is a key mediator of angiogenesis, signalling via the class IV tyrosine kinase receptor family of VEGF Receptors (VEGFRs). Although VEGF-A ligands bind to both VEGFR1 and VEGFR2, they primarily signal via VEGFR2 leading to endothelial cell proliferation, survival, migration and vascular permeability. Distinct VEGF-A isoforms result from alternative splicing of the Vegfa gene at exon 8, resulting in VEGFxxxa or VEGFxxxb isoforms. Alternative splicing events at exons 5⁻7, in addition to recently identified posttranslational read-through events, produce VEGF-A isoforms that differ in their bioavailability and interaction with the co-receptor Neuropilin-1. This review explores the molecular pharmacology of VEGF-A isoforms at VEGFR2 in respect to ligand binding and downstream signalling. To understand how VEGF-A isoforms have distinct signalling despite similar affinities for VEGFR2, this review re-evaluates the typical classification of these isoforms relative to the prototypical, “pro-angiogenic” VEGF165a. We also examine the molecular mechanisms underpinning the regulation of VEGF-A isoform signalling and the importance of interactions with other membrane and extracellular matrix proteins. As approved therapeutics targeting the VEGF-A/VEGFR signalling axis largely lack long-term efficacy, understanding these isoform-specific mechanisms could aid future drug discovery efforts targeting VEGF receptor pharmacology.

A computational analysis of pro-angiogenic therapies for peripheral artery disease

Inducing therapeutic angiogenesis to effectively form hierarchical, non-leaky networks of perfused vessels in tissue engineering applications and ischemic disease remains an unmet challenge, despite extensive research and multiple clinical trials. Here, we use a previously-developed, multi-scale, computational systems pharmacology model of human peripheral artery disease to screen a diverse array of promising pro-angiogenic strategies, including gene therapy, biomaterials, and antibodies. Our previously-validated model explicitly accounts for VEGF immobilization, Neuropilin-1 binding, and weak activation of VEGF receptor 2 (VEGFR2) by the "VEGFxxxb" isoforms. First, we examine biomaterial-based delivery of VEGF engineered for increased affinity to the extracellular matrix. We show that these constructs maintain VEGF close to physiological levels and extend the duration of VEGFR2 activation. We demonstrate the importance of sub-saturating VEGF dosing to prevent angioma formation. Second, we examine the potential of ligand- or receptor-based gene therapy to normalize VEGF receptor signaling. Third, we explore the potential for antibody-based pro-angiogenic therapy. Our model supports recent observations that improvement in perfusion following treatment with anti-VEGF165b in mice is mediated by VEGF-receptor 1, not VEGFR2. Surprisingly, the model predicts that the approved anti-VEGF cancer drug, bevacizumab, may actually improve signaling of both VEGFR1 and VEGFR2 via a novel 'antibody swapping' effect that we demonstrate here. Altogether, this model provides insight into the mechanisms of action of several classes of pro-angiogenic strategies within the context of the complex molecular and physiological processes occurring in vivo. We identify molecular signaling similarities between promising approaches and key differences between promising and ineffective strategies.

Integrating Transcriptomic Data with Mechanistic Systems Pharmacology Models for Virtual Drug Combination Trials

Monotherapy clinical trials with mutation-targeted kinase inhibitors, despite some success in other cancers, have yet to impact glioblastoma (GBM). Besides insufficient blood-brain barrier penetration, combinations are key to overcoming obstacles such as intratumoral heterogeneity, adaptive resistance, and the epistatic nature of tumor genomics that cause mutation-targeted therapies to fail. With now hundreds of potential drugs, exploring the combination space clinically and preclinically is daunting. We are building a simulation-based approach that integrates patient-specific data with a mechanistic computational model of pan-cancer driver pathways (receptor tyrosine kinases, RAS/RAF/ERK, PI3K/AKT/mTOR, cell cycle, apoptosis, and DNA damage) to prioritize drug combinations by their simulated effects on tumor cell proliferation and death. Here we illustrate a first step, tailoring the model to 14 GBM patients from The Cancer Genome Atlas defined by an mRNA-seq transcriptome, and then simulating responses to three promiscuous FDA-approved kinase inhibitors (bosutinib, ibrutinib, and cabozantinib) with evidence for blood-brain barrier penetration. The model captures binding of the drug to primary targets and off-targets based on published affinity data and simulates responses of 100 heterogeneous tumor cells within a patient. Single drugs are marginally effective or even counterproductive. Common copy number alterations (PTEN loss, EGFR amplification, and NF1 loss) have a negligible correlation with single-drug or combination efficacy, reinforcing the importance of postgenetic approaches that account for kinase inhibitor promiscuity to match drugs to patients. Drug combinations tend to be either cytostatic or cytotoxic, but seldom both, highlighting the need for considering targeted and nontargeted therapy. Although we focus on GBM, the approach is generally applicable.

Multi-physics interactions drive VEGFR2 relocation on endothelial cells

Vascular Endothelial Growth Factor Receptor-2 (VEGFR2) is a pro-angiogenic receptor, expressed on endothelial cells (ECs). Although biochemical pathways that follow the VEGFR2 activation are well established, knowledge about the dynamics of receptors on the plasma membrane remains limited. Ligand stimulation induces the polarization of ECs and the relocation of VEGFR2, either in cell protrusions or in the basal aspect in cells plated on ligand-enriched extracellular matrix (ECM). We develop a mathematical model in order to simulate the relocation of VEGFR2 on the cell membrane during the mechanical adhesion of cells onto a ligand-enriched substrate. Co-designing the in vitro experiments with the simulations allows identifying three phases of the receptor dynamics, which are controlled respectively by the high chemical reaction rate, by the mechanical deformation rate, and by the diffusion of free receptors on the membrane. The identification of the laws that regulate receptor polarization opens new perspectives toward developing innovative anti-angiogenic strategies through the modulation of EC activation.

Systems Pharmacology of VEGF165b in Peripheral Artery Disease

We built a whole‐body computational model to study the role of the poorly understood vascular endothelial growth factor (VEGF)_165b splice isoform in peripheral artery disease (PAD). This model was built and validated using published and new experimental data from cells, mice, and humans, and explicitly accounts for known properties of VEGF_165b: lack of extracellular matrix (ECM)‐binding and weak phosphorylation of vascular endothelial growth factor receptor‐2 (VEGFR2) in vitro. The resulting model captures all known information about VEGF_165b distribution and signaling in human PAD, and provides novel, nonintuitive insight into VEGF_165b mechanism of action in vivo. Although VEGF_165a and VEGF_165b compete for VEGFR2 in vitro, simulations show that these isoforms do not compete for VEGFR2 at much lower physiological concentrations. Instead, reduced VEGF_165a may drive impaired VEGFR2 signaling. The model predicts that VEGF_165b does compete for binding to VEGFR1, supporting a VEGFR1‐mediated response to anti‐VEGF_165b. The model predicts a key role for VEGF_165b in PAD, but in a different way than previously hypothesized.