Increase of plasma VEGF after intravenous administration of bevacizumab is predicted by a pharmacokinetic model

Translation of Monoclonal Antibodies Pharmacokinetics from Animal to Human Using Physiologically Based Modeling in Open Systems Pharmacology (OSP) Suite: A Retrospective Analysis of Bevacizumab

Interspecies translation of monoclonal antibodies (mAbs) pharmacokinetics (PK) in presence of target-mediated drug disposition (TMDD) is particularly challenging. Incorporation of TMDD in physiologically based PK (PBPK) modeling is recent and needs to be consolidated and generalized to provide better prediction of TMDD regarding inter-species translation during preclinical and clinical development steps of mAbs. The objective of this study was to develop a generic PBPK translational approach for mAbs using the open-source software (PK-Sim® and Mobi®). The translation of bevacizumab based on data in non-human primates (NHP), healthy volunteers (HV), and cancer patients was used as a case example for model demonstration purpose. A PBPK model for bevacizumab concentration-time data was developed using data from literature and the Open Systems Pharmacology (OSP) Suite version 10. PK-sim® was used to build the linear part of bevacizumab PK (mainly FcRn-mediated), whereas MoBi® was used to develop the target-mediated part. The model was first developed for NHP and used for a priori PK prediction in HV. Then, the refined model obtained in HV was used for a priori prediction in cancer patients. A priori predictions were within 2-fold prediction error (predicted/observed) for both area under the concentration-time curve (AUC) and maximum concentration (Cmax) and all the predicted concentrations were within 2-fold average fold error (AFE) and average absolute fold error (AAFE). Sensitivity analysis showed that FcRn-mediated distribution and elimination processes must be accounted for at all mAb concentration levels, whereas the lower the mAb concentration, the more significant the target-mediated elimination. This project is the first step to generalize the full PBPK translational approach in Model-Informed Drug Development (MIDD) of mAbs using OSP Suite.

Basal VEGF-A and ACE Plasma Levels of Metastatic Colorectal Cancer Patients Have Prognostic Value for First-Line Treatment with Chemotherapy Plus Bevacizumab

Simple Summary

Molecular biology knowledge has enabled the incorporation of targeted therapies, such as the anti-angiogenic drug bevacizumab, into combined chemotherapy regimens for the treatment of metastatic colorectal cancer. However, to date, there are no reliable useful biomarkers to predict the efficacy of this anti-angiogenic therapy. The objective of this prospective study was to evaluate potential circulating plasma biomarkers in mCRC patients prior to the start of first-line treatment with chemotherapy plus bevacizumab. We found that high VEGF-A and low ACE plasma levels were associated with poor OS after treatment. Moreover, a simple scoring system combining both biomarkers efficiently stratified patients into high- or low-risk groups, which allows the selection of patients for anti-angiogenic therapy.

Abstract

The identification of factors that respond to anti-angiogenic therapy would represent a significant advance in the therapeutic management of metastatic-colorectal-cancer (mCRC) patients. We previously reported the relevance of VEGF-A and some components of the renin–angiotensin-aldosterone system (RAAS) in the response to anti-angiogenic therapy in cancer patients. Therefore, this prospective study aims to evaluate the prognostic value of basal plasma levels of VEGF-A and angiotensin-converting enzyme (ACE) in 73 mCRC patients who were to receive bevacizumab-based therapies as a first-line treatment. We found that high basal VEGF-A plasma levels were significantly associated with worse overall survival (OS) and progression-free survival (FPS). On the other hand, low ACE levels were significantly associated with poor OS. Importantly, a simple scoring system combining the basal plasma levels of VEGF-A and ACE efficiently stratified mCRC patients, according to OS, into high-risk or low-risk groups, prior to their treatment with bevacizumab. In conclusion, our study supports that VEGF-A and ACE may be potential biomarkers for selecting those mCRC patients who will most benefit from receiving chemotherapy plus bevacizumab treatment in first-line therapy. Additionally, our data reinforce the notion of a close association between the RAAS and the anti-angiogenic response in cancer.

Mechanistic characterization of endothelial sprouting mediated by pro-angiogenic signaling

Objective

We aim to quantitatively characterize the crosstalk between VEGF‐ and FGF‐mediated angiogenic signaling and endothelial sprouting, to gain mechanistic insights and identify novel therapeutic strategies.

Methods

We constructed an experimentally validated hybrid agent‐based mathematical model that characterizes endothelial sprouting driven by FGF‐ and VEGF‐mediated signaling. We predicted the total sprout length, number of sprouts, and average length by the mono‐ and co‐stimulation of FGF and VEGF.

Results

The experimentally fitted and validated model predicts that FGF induces stronger angiogenic responses in the long‐term compared with VEGF stimulation. Also, FGF plays a dominant role in the combination effects in endothelial sprouting. Moreover, the model suggests that ERK and Akt pathways and cellular responses contribute differently to the sprouting process. Last, the model predicts that the strategies to modulate endothelial sprouting are context‐dependent, and our model can identify potential effective pro‐ and anti‐angiogenic targets under different conditions and study their efficacy.

Conclusions

The model provides detailed mechanistic insight into VEGF and FGF interactions in sprouting angiogenesis. More broadly, this model can be utilized to identify targets that influence angiogenic signaling leading to endothelial sprouting and to study the effects of pro‐ and anti‐angiogenic therapies.

Cetuximab-Mediated Protection from Hypoxia- Induced Cell Death: Implications for Therapy Sequence in Colorectal Cancer

Simple Summary

Therapeutic antibodies are an integral part of treatment regimens for metastasized colorectal cancer. In KRAS wildtype tumors both bevacizumab and cetuximab are active. While bevacizumab has previously been shown to induce tumor hypoxia, we here report that EGFR inhibition by cetuximab protects colon cancer cells from hypoxia-induced cell death. This effect appears to be responsible for the inferior efficacy of a treatment sequence of bevacizumab followed by cetuximab versus an inverse sequence that we observed in a colorectal cancer mouse model. It also offers a mechanistic explanation for effects observed in clinical trials such as underadditive or even detrimental effects when combining bevacizumab and cetuximab (CAIRO2 trial) and the superior efficacy of first line cetuximab (FIRE-3 trial) under chemotherapy backbones in colorectal cancer.

Abstract

Monoclonal antibodies like cetuximab, targeting the epidermal growth factor receptor (EGFR), and bevacizumab, targeting the vascular endothelial growth factor (VEGF), are an integral part of treatment regimens for metastasized colorectal cancer. However, inhibition of the EGFR has been shown to protect human glioma cells from cell death under hypoxic conditions. In colon carcinoma cells, the consequences of EGFR blockade in hypoxia (e.g., induced by bevacizumab) have not been evaluated yet. LIM1215 and SW948 colon carcinoma and LNT-229 glioblastoma cells were treated with cetuximab, PD153035, and erlotinib and analyzed for cell density and viability. The sequential administration of either cetuximab followed by bevacizumab (CET->BEV) or bevacizumab followed by cetuximab (BEV->CET) was investigated in a LIM1215 (KRAS wildtype) and SW948 (KRAS mutant) xenograft mouse model. In vitro, cetuximab protected from hypoxia. In the LIM1215 model, a survival benefit with cetuximab and bevacizumab monotherapy was observed, but only the sequence CET->BEV showed an additional benefit. This effect was confirmed in the SW948 model. Our observations support the hypothesis that bevacizumab modulates the tumor microenvironment (e.g., by inducing hypoxia) where cetuximab could trigger protective effects when administered later on. The sequence CET->BEV therefore seems to be superior as possible mutual adverse effects are bypassed.

Early Tumor Shrinkage and Depth of Response in the Second-Line Treatment for KRAS exon2 Wild-Type Metastatic Colorectal Cancer: An Exploratory Analysis of the Randomized Phase 2 Trial Comparing Panitumumab and Bevacizumab in Combination with FOLFIRI (WJOG6210G)

BACKGROUND

Predictive markers for the clinical outcomes of second-line treatment in patients with metastatic colorectal cancer (mCRC) remain unclear.

OBJECTIVE

This retrospective biomarker study was conducted to explore predictive markers for patients with KRAS exon 2 wild-type mCRC who were treated with FOLFIRI plus panitumumab (Pani) or bevacizumab (Bev) in the WJOG6210G trial.

PATIENTS AND METHODS

The associations of early tumor shrinkage (ETS), tumor location, and VEGF-D with progression-free survival (PFS) and overall survival (OS) were analyzed using a Cox proportional hazards model. Spearman's correlation coefficient was used to analyze the association of depth of response (DpR) with PFS and OS. Serum VEGF-D levels were measured in samples collected before treatment using magnetic bead panel Milliplex xMAP kits.

RESULTS

In total, 101 patients (Pani, n = 49; Bev, n = 52) were enrolled in this study. Patients with ETS had longer PFS (Pani: hazard ratio (HR) 0.40, P = 0.009; Bev: HR 0.078, P = 0.0002) and OS (Pani: HR 0.49, P = 0.044; Bev: HR 0.35, P = 0.048) than patients without ETS. The DpR was moderately correlated with PFS and OS in Pani (rs = 0.75, P < 0.001; rs = 0.60, P < 0.001) and Bev groups (rs = 0.68, P < 0.001; rs = 0.44, P = 0.002). No significant differences were observed in PFS and OS between the two treatment groups even if in left-sided tumors. No significant interaction between VEGF-D levels and treatment was observed in PFS and OS.

CONCLUSIONS

ETS and DpR serve as surrogate markers of PFS and OS in the second-line treatment with FOLFIRI plus targeted agent for mCRC.

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

Background

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

Methods

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

Results

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

Conclusions

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

The association between intravitreal ranibizumab therapy and serum cytokine concentrations in patients with diabetic macular edema

Aim

To investigate the association between intravitreal ranibizumab therapy and serum cytokine concentrations in patients with diabetic macular edema (DME).

Methods

Twenty-five patients with center-involved DME were recruited prospectively. Serum samples were collected from the patients before and 4 weeks after two ranibizumab injections. The levels of 32 cytokines at these two time points were assessed using a multiplex array assay.

Results

Following two ranibizumab injections, there was a statistically significant decrease in the median [interquartile range] levels of Interleukin 1-1beta (IL-1β) from 5.56 [3.6, 8.75] to 2.33 [1.51, 2.89], Interleukin 13 (IL-13) from 4.30 [1.84, 18.55] to 0.38 [0.38, 0.78], granulocyte-colony stimulating factor (G-CSF) from 64.65 [42.9, 108] to 37.8 [27.3, 46.37], Interferon gamma (IFN-γ) from 241 [103.33, 753.4] to 94.4626 [42.04, 118.58], Interferon gamma-induced protein 10 (IP-10) from 234.68 [144.16, 285.98] to 158.73 [94.71, 198.64], Macrophage Inflammatory Protein-1 alpha (MIP-1α) from 3.65 [2.62, 11.02] to 1.41 [0.94, 1.88], and Tumor necrosis factor- alpha (TNF-α) from 131.09 [100.68,28 240.27] to 45.19 [24.04, 68.55]. There was a statistically significant increase in the levels of Interleukin 9 (IL-9) from 0.76 [0.76, 7.03] to 19.67 [5.36 27.76], Macrophage Inflammatory Protein-1 beta (MIP-1β) from 0.28 [0.28, 30 0.28] to 6.79 [I3.74, 14.16], Vascular endothelial growth factor (VEGF) from 2.55 [2.55, 2.55] to 25.24 [14.51, 41.73], and soluble vascular endothelial growth factor -1 (sVEGFR-1) from 333.92 [204.99, 440.43] to 500.12 [38.7, 786.91]. A Bonferroni-corrected p value of 0.00156 was considered statistically significant.

Conclusions

In patients with DME, intravitreal ranibizumab therapy appears to influence the serum levels of a range of cytokines. After two injections, intravitreal ranibizumab therapy appears to be associated with a significant decrease in inflammatory mediators and a rise in VEGF and sVEGFR1.

Aluminum Metasurface with Hybrid Multipolar Plasmons for 1000-Fold Broadband Visible Fluorescence Enhancement and Multiplexed Biosensing

Aluminum (Al)-based nanoantennae traditionally suffer from weak plasmonic performance in the visible range, necessitating the application of more expensive noble metal substrates for rapidly expanding biosensing opportunities. We introduce a metasurface comprising Al nanoantennae of nanodisks-in-cavities that generate hybrid multipolar lossless plasmonic modes to strongly enhance local electromagnetic fields and increase the coupled emitter's local density of states throughout the visible regime. This results in highly efficient electromagnetic field confinement in visible wavelengths by these nanoantennae, favoring real-world plasmonic applications of Al over other noble metals. Additionally, we demonstrate spontaneous localization and concentration of target molecules at metasurface hotspots, leading to further improved on-chip detection sensitivity and a broadband fluorescence-enhancement factor above 1000 for visible wavelengths with respect to glass chips commonly used in bioassays. Using the metasurface and a multiplexing technique involving three visible wavelengths, we successfully detected three biomarkers, insulin, vascular endothelial growth factor, and thrombin relevant to diabetes, ocular and cardiovascular diseases, respectively, in a single 10 μL droplet containing only 1 fmol of each biomarker.

Exploring the Extracellular Regulation of the Tumor Angiogenic Interaction Network Using a Systems Biology Model

Tumor angiogenesis is regulated by pro- and anti-angiogenic factors. Anti-angiogenic agents target the interconnected network of angiogenic factors to inhibit neovascularization, which subsequently impedes tumor growth. Due to the complexity of this network, optimizing anti-angiogenic cancer treatments requires detailed knowledge at a systems level. In this study, we constructed a tumor tissue-based model to better understand how the angiogenic network is regulated by opposing mediators at the extracellular level. We consider the network comprised of two pro-angiogenic factors: vascular endothelial growth factor (VEGF) and basic fibroblast growth factor (FGF2), and two anti-angiogenic factors: thrombospondin-1 (TSP1) and platelet factor 4 (PF4). The model's prediction of angiogenic factors' distribution in tumor tissue reveals the localization of different factors and indicates the angiogenic state of the tumor. We explored how the distributions are affected by the secretion of the pro- and anti-angiogenic factors, illustrating how the angiogenic network is regulated in the extracellular space. Interestingly, we identified a counterintuitive result that the secretion of the anti-angiogenic factor PF4 can enhance pro-angiogenic signaling by elevating the levels of the interstitial and surface-level pro-angiogenic species. This counterintuitive situation is pertinent to the clinical setting, such as the release of anti-angiogenic factors in platelet activation or the administration of exogenous PF4 for anti-angiogenic therapy. Our study provides mechanistic insights into this counterintuitive result and highlights the role of heparan sulfate proteoglycans in regulating the interactions between angiogenic factors. This work complements previous studies aimed at understanding the formation of angiogenic complexes in tumor tissue and helps in the development of anti-cancer strategies targeting angiogenesis.

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

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

Predictive Approach Identifies Molecular Targets and Interventions to Restore Angiogenesis in Wounds With Delayed Healing

Impaired angiogenesis is a hallmark of wounds with delayed healing, and currently used therapies to restore angiogenesis have limited efficacy. Here, we employ a computational simulation-based approach to identify influential molecular and cellular processes, as well as protein targets, whose modulation may stimulate angiogenesis in wounds. We developed a mathematical model that captures the time courses for platelets, 9 cell types, 29 proteins, and oxygen, which are involved in inflammation, proliferation, and angiogenesis during wound healing. We validated our model using previously published experimental data. By performing global sensitivity analysis on thousands of simulated wound-healing scenarios, we identified six processes (among the 133 modeled in total) whose modulation may improve angiogenesis in wounds. By simulating knockouts of 25 modeled proteins and by simulating different wound-oxygenation levels, we identified four proteins [namely, transforming growth factor (TGF)-β, vascular endothelial growth factor (VEGF), fibroblast growth factor-2 (FGF-2), and angiopoietin-2 (ANG-2)], as well as oxygen, as therapeutic targets for stimulating angiogenesis in wounds. Our modeling results indicated that simultaneous inhibition of TGF-β and supplementation of either FGF-2 or ANG-2 could be more effective in stimulating wound angiogenesis than the modulation of either protein alone. Our findings suggest experimentally testable intervention strategies to restore angiogenesis in wounds with delayed healing.

Sequential cetuximab/bevacizumab therapy is associated with improved outcomes in patients with wild-type KRAS exon 2 metastatic colorectal cancer

Purpose

Combination of biological therapy and chemotherapy improves the survival of patients with metastatic colorectal cancer (mCRC). However, the optimal biological therapy sequence remains unclear. In this retrospective study, we evaluated the clinical outcomes of patients with mCRC treated with different sequences of biological therapies as first‐ and third‐line therapy.

Methods

We only included patients with wild‐type KRAS exon 2 mCRC who had received cetuximab, bevacizumab, and standard chemotherapy. The patients were treated with cetuximab or bevacizumab as first‐ or third‐line therapy combined with a similar chemotherapy backbone.

Results

In total, 102 patients were included. Forty‐six patients received first‐line cetuximab therapy followed by third‐line bevacizumab therapy (cetuximab → bevacizumab group) and 56 patients received first‐line bevacizumab therapy followed by third‐line cetuximab therapy (bevacizumab → cetuximab group). The cetuximab → bevacizumab group was associated with increased survival (OS) compared with the bevacizumab → cetuximab group (median OS: 30.4 months vs 25.7 months, hazard ratio (HR): 0.55, 95% confidence interval (CI): 0.36‐0.86). When calculated from the start of second‐ and third‐line therapies, OS was also higher in the cetuximab → bevacizumab group (second‐line: 20.6 months vs 14.8 months, HR: 0.54, 95% CI: 0.34‐0.81; third‐line: 12.5 months vs 9.9 months, HR: 0.53, 95% CI: 0.35‐0.83). The cetuximab → bevacizumab group was also associated with better progression‐free survival than the bevacizumab → cetuximab group (8.8 vs 4.5 months, HR: 0.43, 95% CI: 0.25‐0.58) in the third‐line setting, but not in the first‐ or second‐line settings.

Conclusions

Our study demonstrated that first‐line cetuximab therapy followed by third‐line bevacizumab therapy was associated with favorable clinical outcomes as compared to the reverse sequence.

In silico simulation of a clinical trial with anti-CTLA-4 and anti-PD-L1 immunotherapies in metastatic breast cancer using a systems pharmacology model

The low response rate of immune checkpoint blockade in breast cancer has highlighted the need for predictive biomarkers to identify responders. While a number of clinical trials are ongoing, testing all possible combinations is not feasible. In this study, a quantitative systems pharmacology model is built to integrate immune-cancer cell interactions in patients with breast cancer, including central, peripheral, tumour-draining lymph node (TDLN) and tumour compartments. The model can describe the immune suppression and evasion in both TDLN and the tumour microenvironment due to checkpoint expression, and mimic the tumour response to checkpoint blockade therapy. We investigate the relationship between the tumour response to checkpoint blockade therapy and composite tumour burden, PD-L1 expression and antigen intensity, including their individual and combined effects on the immune system, using model-based simulations. The proposed model demonstrates the potential to make predictions of tumour response of individual patients given sufficient clinical measurements, and provides a platform that can be further adapted to other types of immunotherapy and their combination with molecular-targeted therapies. The patient predictions demonstrate how this systems pharmacology model can be used to individualize immunotherapy treatments. When appropriately validated, these approaches may contribute to optimization of breast cancer treatment.

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

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

Disease Control With FOLFIRI Plus Ziv-aflibercept (zFOLFIRI) Beyond FOLFIRI Plus Bevacizumab: Case Series in Metastatic Colorectal Cancer (mCRC)

Background: The prognosis of patients with metastatic colorectal cancer (mCRC) is poor, especially after failure of initial systemic therapy. The VELOUR study showed modestly prolonged overall survival (OS) with ziv-aflibercept plus 5-fluorouracil, leucovorin, and irinotecan (zFOLFIRI) vs. placebo+FOLFIRI after progression on 5-fluoruracil, leucovorin, and oxaliplatin (FOLFOX) ± bevacizumab. The utility of zFOLFIRI after bevacizumab+FOLFIRI is unknown and not recommended in NCCN guidelines. We explored whether zFOLFIRI may be active beyond progression on bevacizumab+FOLFIRI. Methods: We undertook a retrospective analysis of patients treated in routine clinical practice. A chart review was conducted for a cohort (N = 19) of advanced cancer patients (18 mCRC) who received zFOLFIRI from 2014 to 2018 at Fox Chase Cancer Center (FCCC). Analysis included time on zFOLFIRI, PFS, OS, CEA trends and adverse events. A second mCRC cohort (N = 26) from the Flatiron Health EHR-derived database treated with zFOLFIRI after prior bevacizumab+FOLFOX and bevacizumab+FOLFIRI was analyzed for time-on-treatment and overall survival. Results: Median age of mCRC cohort at zFOLFIRI treatment was 54 (FCCC; N = 18) and 62 (Flatiron Health-cohort; N = 26). Of 18 FCCC mCRC patients, 1 patient had prior bevacizumab+FOLFOX and ramucirumab+irinotecan prior to zFOLFIRI for 8.5 months. Of 17 FCCC mCRC patients with prior bevacizumab+FOLFIRI who received zFOLFIRI, 13 had mutant-KRAS, 3 WT-KRAS, and one BRAF-V600E. The patient with BRAF-V600E mutation achieved stable disease on zFOLFIRI after multiple BRAF-targeted therapies. One patient (WT-KRAS mCRC) remained on zFOLFIRI for 14 months. Of 14 patients with mutated-KRAS, 8 remained on zFOLFIRI for >5 months including 3 for >15 months. The rate-of-change in CEA measures on zFOLFIRI was significantly different (p = 0.004) between rapid progressors and those with PFS>4 months. For mCRC patients treated with zFOLFIRI in the 3rd line or greater (N = 18), median PFS was 7.1 months (214 days) and median OS was 13.8 months (416 days). Median time-on-treatment with zFOLFIRI in the Flatiron Health cohort was 4.4 months, median OS was 7.8 months, and longest time-on-treatment with zFOLFIRI was 266 days. Conclusions: In these small real-world cohorts, clinical meaningful stable disease and overall survival on zFOLFIRI beyond progression on bevacizumab+FOLFIRI was observed in patients with mCRC. Further exploration of this approach is warranted.

Treatment sequencing in metastatic colorectal cancer

Metastatic colorectal cancer (mCRC) remains incurable in most cases, but survival has improved with advances in cytotoxic chemotherapy and targeted agents. However, the optimal use and sequencing of these agents across multiple lines of treatment is unclear. Here, we review current treatment approaches and optimal treatment sequencing across the first-, second- and third-line settings in mCRC, including biological aspects affecting sequencing and rechallenge. Effective first-line therapy is a key determinant of treatment outcomes and should be selected after considering both clinical factors and biological markers, notably RAS and BRAF. The second-line regimen choice depends on the systemic therapies given in first-line. Anti-angiogenic agents (e.g. bevacizumab, ramucirumab and aflibercept) are indicated for most patients, whereas epidermal growth factor receptor (EGFR) inhibitors do not improve survival in the second-line setting. Molecular profiling is important in third-line treatment, with options in RAS wild-type patients including EGFR inhibitors (cetuximab or panitumumab), regorafenib and trifluridine/tipiracil. Immunotherapy with pembrolizumab or nivolumab ± ipilimumab may be considered for patients with high microsatellite instability disease. Targeting HER2/neu amplification shows promise for the subset of CRC tumours displaying this abnormality. Sequencing decisions are complicated by the potential for any treatment break or de-escalation to evoke a distinct clinical progression type. Ongoing trials are investigating the optimal sequencing and timing of therapies for mCRC. Molecular profiling has established new targets, and increasing knowledge of tumour evolution under drug pressure will possibly impact on sequencing.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

The VEGF rise in blood of bevacizumab patients is not based on tumor escape but a host-blockade of VEGF clearance

Vascular endothelial growth factor (VEGF) has become a major target in cancer treatment as it promotes tumor angiogenesis. Therapy with anti-VEGF antibody bevacizumab reportedly induces high levels of circulating VEGF which may potentially contribute to resistance. Based on animal or computational models, mechanisms of VEGF induction by bevacizumab have been proposed but not verified in the clinical setting. Hence, we evaluated sixty patients with colorectal cancer metastases for changes in plasma VEGF during neoadjuvant/conversion and adjuvant chemotherapy with or without bevacizumab. VEGF expression was assessed in tissue sections of liver metastases. The VEGF source was investigated with in vitro cultures of tumor, endothelial cells, fibroblasts and platelets, and potential protein stabilization due to anti-VEGF therapy was addressed. A VEGF rise was observed in blood of bevacizumab patients but not in chemotherapy controls, and VEGF was found to be largely complexed by the antibody. A comparable VEGF increase occurred in the presence (neoadjuvant) and absence of the tumor (adjuvant). Accordingly, VEGF expression in tumor tissue was not determined by bevacizumab treatment. Investigations with isolated cell types did not reveal VEGF production in response to bevacizumab. However, antibody addition to endothelial cultures led to a dose-dependent blockade of VEGF internalization and hence stabilized VEGF in the supernatant. In conclusion, the VEGF rise in cancer patients treated with bevacizumab is not originating from the tumor. The accumulation of primarily host-derived VEGF in circulation can be explained by antibody interference with receptor-mediated endocytosis and protein degradation. Thus, the VEGF increase in response to bevacizumab therapy should not be regarded as a tumor escape mechanism.

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.

Generative mathematical modelling to demonstrate virtual simulations of neovascular age related macular degeneration

Purpose

To develop a generative mathematical model of wet age-related macular degeneration (AMD) and model the impact of injections of anti-vascular endothelial growth factor to virtual patients with the condition.

Methods

We isolated key pathophysiological components of macular degeneration in terms of macular edema development and response to anti-vascular endothelial growth factor (VEGF) agents. We developed mathematical models for each of these components using constants determined from published biological experimentation. Consequently, we combined the mathematical models of the separate components to arrive at an end-to-end model of the evolution of macular edema size and its response to treatment.

Results

We present a series of simulations based upon our idealised model. Initially, we demonstrate the theoretical change in macular edema height in wet macular degeneration over time without and with anti-VEGF interventions. In our final simulation, we demonstrate the powerful possibilities of virtual clinical trials by simulating a virtual model of a landmark study using our existing mathematical AMD model.

Conclusions

Using our mathematical modelling based upon known pathological and pharmacological processes we have been able to model the effect of intravitreal injection of an anti-VEGF agent on macular edema from age related macular degeneration. We were subsequently able to mathematically simulate a major clinical trial with results that mirror many key features of the clinical established study. We anticipate that the generative model presented here can evolve to be a useful supportive tool in the challenge to deliver optimal therapy for patients with wet macular degeneration.

Discovery of High-Affinity PDGF-VEGFR Interactions: Redefining RTK Dynamics

Nearly all studies of angiogenesis have focused on uni-family ligand-receptor binding, e.g., VEGFs bind to VEGF receptors, PDGFs bind to PDGF receptors, etc. The discovery of VEGF-PDGFRs binding challenges this paradigm and calls for investigation of other ligand-receptor binding possibilities. We utilized surface plasmon resonance to identify and measure PDGF-to-VEGFR binding rates, establishing cut-offs for binding and non-binding interactions. We quantified the kinetics of the recent VEGF-A:PDGFRβ interaction for the first time with K_D = 340 pM. We discovered new PDGF:VEGFR2 interactions with PDGF-AA:R2 K_D = 530 nM, PDGF-AB:R2 K_D = 110 pM, PDGF-BB:R2 K_D = 40 nM, and PDGF-CC:R2 K_D = 70 pM. We computationally predict that cross-family PDGF binding could contribute up to 96% of VEGFR2 ligation in healthy conditions and in cancer. Together the identification, quantification, and simulation of these novel cross-family interactions posits new mechanisms for understanding anti-angiogenic drug resistance and presents an expanded role of growth factor signaling with significance in health and disease.

A phase I trial of intraperitoneal GEN-1, an IL-12 plasmid formulated with PEG-PEI-cholesterol lipopolymer, administered with pegylated liposomal doxorubicin in patients with recurrent or persistent epithelial ovarian, fallopian tube or primary peritoneal cancers: An NRG Oncology/Gynecologic Oncology Group study

OBJECTIVE

The study's purpose was to assess safety and efficacy of escalating doses of weekly GEN-1 with pegylated liposomal doxorubicin (PLD) in patients with recurrent or persistent epithelial ovarian, fallopian tube or primary peritoneal cancers (EOC).

METHODS

Patients had persistent or recurrent platinum-resistant EOC. The trial was a standard 3+3 phase I dose escalation design with patients receiving intravenous PLD 40mg/m2 (dose level 1 and 2) or 50mg/m2 (dose level 3) every 28days and intraperitoneal GEN-1 at 24mg/m2 (dose level 1) or 36mg/m2 (dose level 2 and 3) on days 1, 8, 15, and 22 of a 28day cycle. Cycles were repeated every 28days until disease progression. Patients were monitored for toxicity, clinical efficacy, and evidence of systemic and intraperitoneal immunologic effect.

RESULTS

Sixteen evaluable patients received a median of 4cycles (range 1-8). No dose limiting toxicities were found. The adverse side effects were 4 grade 3 anemia, 2 grade 3 abdominal pain, 7 grade 3 neutropenia, and 2 grade 4 neutropenia. A clinical benefit of 57.1% (PR=21.4%; SD=35.7%) was found in the 14 patients with measurable disease. The highest number of partial responses (28.6%) and stable disease (57.1%) were found at dose level 3. The maximum tolerated dose was not reached. Increases in IL-12, IFN-γ, and TNF-α levels were found in peritoneal fluid following GEN-1 treatment.

CONCLUSIONS

GEN-1 in combination with PLD has encouraging clinical benefit and biological activity in recurrent or persistent EOC and warrants further investigation with escalating doses of GEN-1.

Multiscale modeling methods in biomechanics

More and more frequently, computational biomechanics deals with problems where the portion of physical reality to be modeled spans over such a large range of spatial and temporal dimensions, that it is impossible to represent it as a single space-time continuum. We are forced to consider multiple space-time continua, each representing the phenomenon of interest at a characteristic space-time scale. Multiscale models describe a complex process across multiple scales, and account for how quantities transform as we move from one scale to another. This review offers a set of definitions for this emerging field, and provides a brief summary of the most recent developments on multiscale modeling in biomechanics. Of all possible perspectives, we chose that of the modeling intent, which vastly affect the nature and the structure of each research activity. To the purpose we organized all papers reviewed in three categories: 'causal confirmation,' where multiscale models are used as materializations of the causation theories; 'predictive accuracy,' where multiscale modeling is aimed to improve the predictive accuracy; and 'determination of effect,' where multiscale modeling is used to model how a change at one scale manifests in an effect at another radically different space-time scale. Consistent with how the volume of computational biomechanics research is distributed across application targets, we extensively reviewed papers targeting the musculoskeletal and the cardiovascular systems, and covered only a few exemplary papers targeting other organ systems. The review shows a research subdomain still in its infancy, where causal confirmation papers remain the most common. WIREs Syst Biol Med 2017, 9:e1375. doi: 10.1002/wsbm.1375 For further resources related to this article, please visit the WIREs website.

Immune modulation associated with vascular endothelial growth factor (VEGF) blockade in patients with glioblastoma

BACKGROUND

Vascular endothelial growth factor (VEGF), in addition to being pro-angiogenic, is an immunomodulatory cytokine systemically and in the tumor microenvironment. We previously reported the immunomodulatory effects of radiation and temozolomide (TMZ) in newly diagnosed glioblastoma. This study aimed to assess changes in peripheral blood mononuclear cell (PBMC) populations, plasma cytokines, and growth factor concentrations following treatment with radiation, TMZ, and bevacizumab (BEV).

METHODS

Eleven patients with newly diagnosed glioblastoma were treated with radiation, TMZ, and BEV, following surgery. We measured immune-related PBMC subsets using multi-parameter flow cytometry and plasma cytokine and growth factor concentrations using electrochemiluminescence-based multiplex analysis at baseline and after 6 weeks of treatment.

RESULTS

The absolute number of peripheral blood regulatory T cells (Tregs) decreased significantly following treatment. The lower number of peripheral Tregs was associated with a CD4+ lymphopenia, and thus, the ratio of Tregs to PBMCs was unchanged. The addition of bevacizumab to standard radiation and temozolomide led to the decrease in the number of circulating Tregs when compared with our prior study. There was a significant decrease in CD8+ cytotoxic and CD4+ recent thymic emigrant T cells, but no change in the number of myeloid-derived suppressor cells. Significant increases in plasma VEGF and placental growth factor (PlGF) concentrations were observed.

CONCLUSIONS

Treatment with radiation, TMZ, and BEV decreased the number but not the proportion of peripheral Tregs and increased the concentration of circulating VEGF. This shift in the peripheral immune cell profile may modulate the tumor environment and have implications for combining immunotherapy with anti-angiogenic therapy.

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

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

Does bevacizumab impact anti-EGFR therapy efficacy in metastatic colorectal cancer?

Anti-EGFR therapy and antiangiogenic therapies are used alone or in combination with chemotherapies to improve survival in metastatic colorectal cancer. However, it is unknown whether pretreatment with antiangiogenic therapy could impact on the efficacy of anti-EGFR therapy.

We selected one hundred and twenty eight patients diagnosed with advanced colorectal cancer with a KRAS and NRAS unmutated tumor. These patients were treated with cetuximab or panitumumab alone or with chemotherapy as second or third-line. Univariate and multivariate Cox model analysis were performed to estimate the effect of a previous bevacizumab regimen on progression free survival and on overall survival during anti-EGFR therapy. In vitro studies using wild type KRAS and NRAS colon cancer cells were performed to evaluate the impact of VEGF-A on cetuximab-induced cell death.

The median progression free survival (PFS) during anti-EGFR treatment was significantly different between the bevacizumab group and the non-bevacizumab group (2.8 and 4 months respectively; p = 0.003). The median overall survival from the beginning of the metastatic disease was similar in the two groups (41.3 and 42 months respectively; p = 0.7). In vitro, VEGF-A induced a resistance toward cetuximab cytotoxicity on three KRAS and NRAS wild type colon cancer cell lines in a VEGFR2 and Stat-3-dependent manner.

All in all, our clinical data, supported by in vitro procedures, suggest that a previous anti-VEGF therapy decreases anti-EGFR efficacy. Although these results are observed in a limited cohort, they could be taken into consideration for a better strategy of care for patient suffering from metastatic colorectal cancer.

A functional bioassay to determine the activity of anti-VEGF antibody therapy in blood of patients with cancer

Background:

Only a small proportion of patients respond to anti-VEGF therapy, pressing the need for a reliable biomarker that can identify patients who will benefit. We studied the biological activity of anti-VEGF antibodies in patients' blood during anti-VEGF therapy by using the Ba/F3-VEGFR2 cell line, which is dependent on VEGF for its growth.

Methods:

Serum samples from 22 patients with cancer before and during treatment with bevacizumab were tested for their effect on proliferation of Ba/F3-VEGFR2 cells. Vascular endothelial growth factor as well as bevacizumab concentrations in serum samples from these patients were determined by enzyme linked immunosorbent assay (ELISA).

Results:

The hVEGF-driven cell proliferation was effectively blocked by bevacizumab (IC₅₀ 3.7 μg ml⁻¹; 95% CI 1.7–8.3 μg ml⁻¹). Cell proliferation was significantly reduced when patients' serum during treatment with bevacizumab was added (22–103% inhibition compared with pre-treatment). Although bevacizumab levels were not related, on-treatment serum VEGF levels were correlated with Ba/F3-VEGFR2 cell proliferation.

Conclusions:

We found that the neutralising effect of anti-VEGF antibody therapy on the biological activity of circulating VEGF can be accurately determined with a Ba/F3-VEGFR2 bioassay. The value of this bioassay to predict clinical benefit of anti-VEGF antibody therapy needs further clinical evaluation in a larger randomised cohort.

Blood-based biomarkers for monitoring antiangiogenic therapy in non-small cell lung cancer

Tumor angiogenesis pathways have been identified as important therapeutic targets in non-small cell lung cancer. However, no biomarkers have been described as predictors of response to antiangiogenic therapy in these patients. In this study, plasma levels of VEGF, bFGF, E-selectin, and S-ICAM and gene expression profiles of peripheral blood mononuclear cells from non-small cell lung cancer patients treated with chemotherapy plus bevacizumab were analyzed before and after treatment. Values were correlated with clinicopathological characteristics and treatment response. Plasma factor levels were measured using commercially available ELISA kits. The TaqMan(®) human angiogenesis array was used to investigate the effect of treatment on gene expression profiles. Kyoto Encyclopedia of Genes and Genomes and Gene Ontology enrichment analysis was performed for differentially expressed genes using WEB-based GEne SeT AnaLysis Toolkit. Our results suggest a benefit for patients with increased plasma levels of VEGF, E-selectin, and S-ICAM in the course of bevacizumab treatment. Also, we identified differentially expressed genes between paired blood samples from patients before and after treatment, and significantly perturbed pathways were predicted. These changes in gene expression and levels of plasma factors could be used to assess the effectiveness of antiangiogenic therapy, in addition to standard clinical and radiological evaluations.

A bi-functional antibody-receptor domain fusion protein simultaneously targeting IGF-IR and VEGF for degradation

Bi-specific antibodies (BsAbs), which can simultaneously block 2 tumor targets, have emerged as promising therapeutic alternatives to combinations of individual monoclonal antibodies. Here, we describe the engineering and development of a novel, human bi-functional antibody-receptor domain fusion molecule with ligand capture (bi-AbCap) through the fusion of the domain 2 of human vascular endothelial growth factor receptor 1 (VEGFR1) to an antibody directed against insulin-like growth factor - type I receptor (IGF-IR). The bi-AbCap possesses excellent stability and developability, and is the result of minimal engineering. Beyond potent neutralizing activities against IGF-IR and VEGF, the bi-AbCap is capable of cross-linking VEGF to IGF-IR, leading to co-internalization and degradation of both targets by tumor cells. In multiple mouse xenograft tumor models, the bi-AbCap improves anti-tumor activity over individual monotherapies. More importantly, it exhibits superior inhibition of tumor growth, compared with the combination of anti-IGF-IR and anti-VEGF therapies, via powerful blockade of both direct tumor cell growth and tumor angiogenesis. The unique "capture-for-degradation" mechanism of the bi-AbCap is informative for the design of next-generation bi-functional anti-cancer therapies directed against independent signaling pathways. The bi-AbCap design represents an alternative approach to the creation of dual-targeting antibody fusion molecules by taking advantage of natural receptor-ligand interactions.

The Best. First. Anti-EGFR before anti-VEGF, in the first-line treatment of RAS wild-type metastatic colorectal cancer: from bench to bedside

BACKGROUND

Since 2013, informative trials exploring the optimal use of available biologic agents in the first-line setting of metastatic colorectal cancer (mCRC) have been presented. These trials have opened a stimulating debate on the biological effect that first-line therapies may have on subsequent lines of treatment even long after the first-line progression.

MATERIALS AND METHODS

We reviewed available preclinical and clinical data on the effect of different sequences of the biological drugs approved for use in mCRC patients. The importance of molecular selection of patients based on RAS mutational status and toxicity and quality-of-life issues were also analyzed.

RESULTS

Convincing evidence exists on the optimal therapeutic effect obtained by using anti-EGFR agents in first-line treatment before anti-VEGF agents. On the contrary, up-front anti-VEGF agents' use seems to determine biological changes that increase the risk of acquired resistance to subsequent EGFR inhibitors. This hypothesis is confirmed by the scarce evidence of EGFR inhibitor activity in second-line treatment. Such a therapeutic optimum is subject to a fine molecular selection based on RAS mutational status.

CONCLUSION

There is accumulating evidence suggesting that, after precise and well-established molecular selection, anti-EGFR agents deliver their maximum efficacy in mCRC patients when given early in the treatment strategy. Their toxicity profile seems manageable under the supervision of experienced physicians. Large randomized trials prospectively confirming the impact of different sequencing strategies are eagerly awaited.

Dose - response relationship of bevacizumab in hereditary hemorrhagic telangiectasia

Hereditary hemorrhagic telangiectasia (HHT), a genetic vascular disorder associated with epistaxis and hepatic shunts, is responsible for high-output cardiac failure in rare cases. Bevacizumab, which targets vascular endothelial growth factor, was shown to decrease both cardiac index (CI) and epistaxis duration in HHT patients with severe liver involvement. The relationship between its serum concentration and change in both CI and epistaxis duration was investigated to design the bevacizumab maintenance dosing regimen of future therapeutic studies. Twenty-five HHT patients with dyspnea and high CI were included in a prospective non-comparative study. They received bevacizumab at a dose of 5 mg/kg per infusion every 14 days for a total of 6 injections. The relationships between bevacizumab serum concentration and both CI and epistaxis duration were described using transit compartments and direct inhibition pharmacokinetic-pharmacodynamic models. The performances of different maintenance regimens were evaluated using simulation. Infusions every 3, 2 and one months were predicted to maintain 41%, 45% and 50% of patients with CI <4 L/min/m(2) at 24 months, respectively. The fraction of patients with <20 min epistaxis per month was predicted to be 34%, 43% and 60%, with infusion every 3, 2 or one months, respectively. Simulations of the effects of different maintenance dosing regimens predict that monthly 5 mg/kg infusions of bevacizumab should allow sustained control of both cardiac index and epistaxis.

Intravitreally Injected Anti-VEGF Antibody Reduces Brown Fat in Neonatal Mice

Anti-vascular endothelial growth factor (VEGF) agents are the mainstay treatment for various angiogenesis-related retinal diseases. Currently, bevacizumab, a recombinant humanized anti-VEGF antibody, is trailed in retinopathy of prematurity, a vasoproliferative retinal disorder in premature infants. However, the risks of systemic complications after intravitreal injection of anti-VEGF antibody in infants are not well understood. In this study, we show that intravitreally injected anti-VEGF antibody is transported into the systemic circulation into the periphery where it reduces brown fat in neonatal C57BL/6 mice. A considerable amount of anti-VEGF antibody was detected in serum after intravitreal injection. Furthermore, in interscapular brown adipose tissue, we found lipid droplet accumulation, decreased VEGF levels, loss of vascular network, and decreased expression of mitochondria-related genes, Ppargc1a and Ucp1, all of which are characteristics of "whitening" of brown fat. With increasing age and body weight, brown fat restored its morphology and vascularity. Our results show that there is a transient, but significant impact of intravitreally administered anti-VEGF antibody on brown adipose tissue in neonatal mice. We suggest that more attention should be focused on the metabolic and developmental significance of brown adipose tissue in bevacizumab treated retinopathy of prematurity infants.

Molecular mechanism matters: Benefits of mechanistic computational models for drug development

Making drug development a more efficient and cost-effective process will have a transformative effect on human health. A key, yet underutilized, tool to aid in this transformation is mechanistic computational modeling. By incorporating decades of hard-won prior knowledge of molecular interactions, cellular signaling, and cellular behavior, mechanistic models can achieve a level of predictiveness that is not feasible using solely empirical characterization of drug pharmacodynamics. These models can integrate diverse types of data from cell culture and animal experiments, including high-throughput systems biology experiments, and translate the results into the context of human disease. This provides a framework for identification of new drug targets, measurable biomarkers for drug action in target tissues, and patient populations for which a drug is likely to be effective or ineffective. Additionally, mechanistic models are valuable in virtual screening of new therapeutic strategies, such as gene or cell therapy and tissue regeneration, identifying the key requirements for these approaches to succeed in a heterogeneous patient population. These capabilities, which are distinct from and complementary to those of existing drug development strategies, demonstrate the opportunity to improve success rates in the drug development pipeline through the use of mechanistic computational models.

Systems biology of the microvasculature

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

A pharmacokinetic binding model for bevacizumab and VEGF165 in colorectal cancer patients

PURPOSE

To characterize the population pharmacokinetics of bevacizumab, its binding properties to VEGF165 and the effect of demographic data and VEGF-A polymorphisms on the interplay between bevacizumab serum pharmacokinetics and VEGF165 serum concentrations in patients with colorectal cancer stage IV.

METHODS

Bevacizumab and VEGF165 data were collected from 19 adult patients with metastatic colorectal cancer enrolled in an observational clinical study. Bevacizumab was administered with one of the following combinations: 5-FU/Leucovorin/Irinotecan, 5-FU/Leucovorin/Oxaliplatin, Capecitabine/Irinotecan at doses ranging from 5 to 10 mg/kg every 2 or 3 weeks. Data analysis was performed using nonlinear mixed-effects modeling implemented in NONMEM 7.3.

RESULTS

A target-mediated drug disposition model adequately described bevacizumab concentration changes over time and its binding characteristics to VEGF165. The estimated clearance of bevacizumab was 0.18 L/day, the free VEGF165 levels at baseline were 212 ng/L, and the elimination rate constant of free VEGF165 was 0.401 day(-1). Body weight was allometrically included in all PK parameters.

CONCLUSION

The final model adequately described the pre- and post-dose concentrations of total bevacizumab and free VEGF165 in patients with colorectal cancer. Model parameters were consistent with those previously reported for patients with solid tumors. Correlations between the binding affinity of bevacizumab and the VEGF-2578C/A and VEGF-634G/C polymorphisms were noticed.

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

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

Multiphasic changes in systemic VEGF following intravitreal injections of ranibizumab in a child

PURPOSE

To investigate whether intravitreal ranibizumab injections administered to a child alter systemic plasma levels of total and free VEGF 165.

METHODS

A 9-year-old child sustained a choroidal rupture from blunt trauma. He subsequently developed a secondary choroidal neovascular membrane, which was treated with five ranibizumab injections over a period of 8 months. Peripheral venous blood samples were taken at each visit over a period of 12 months and plasma was extracted. Plasma VEGF 165 levels were determined using enzyme-linked immunosorbent assay and were assayed both pre- and post-immunodepletion to remove complexed VEGF.

RESULTS

Plasma VEGF 165 levels proved labile following intravitreal injection of ranibizumab. Levels increased by 30% above baseline following the first intravitreal ranibizumab injection, but then returned to baseline despite two subsequent injections. There was then a rebound increase of 67% in total plasma VEGF levels following a further injection, which remained above baseline for 12 weeks despite two further intravitreal ranibizumab injections. Baseline levels were re-attained 26 weeks after the final injection.

CONCLUSIONS

These results suggest intravitreal ranibizumab injections can cause significant, multiphasic changes in systemic VEGF levels. This may be of particular clinical significance in children as VEGF is known to be vital in the development of major organs, in addition to its role in the maintenance of normal organ function in adults.

Decreased peritherapeutic VEGF expression could be a predictor of responsiveness to first-line FOLFIRI plus bevacizumab in mCRC patients

OBJECTIVE

Bevacizumab is the only anti-angiogenic agent approved in first-line therapy for metastatic colorectal cancer (mCRC). Although chemotherapy plus bevacizumab has led to improve outcomes for mCRC patients and is a common choice for first-line treatment of mCRC, previous research has established no prominent biomarker that can help to select patients who may benefit from bevacizumab in order to improve cost-effectiveness and therapeutic outcomes. The aim of this study was to compare pre- and post-therapeutic VEGF immunohistochemical (IHC) expression in mCRC patients treated with FOLFIRI plus bevacizumab to identify its potential role as a predictive biomarker.

METHODS

A total of 57 mCRC patients who underwent FOLFIRI combined with bevacizumab chemotherapy as a first-line neoadjuvant regimen were enrolled and clinical outcome data analyzed.

RESULTS

Low post-therapeutic VEGF expression (P < 0.001) and decreased peri-therapeutic VEGF expression (P < 0.001) were significantly predictive factors of responders. Furthermore, the 6-month progression-free survival (PFS) rate in mCRC patients with decreased peri-therapeutic VEGF expression was significantly better than the rate for those patients with no peri-therapeutic VEGF expression alterations (P = 0.033).

CONCLUSIONS

Decreased peri-therapeutic VEGF expression in mCRC patients could probably be used to predict responsiveness to bevacizumab and subsequent PFS in clinical practice.

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

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

Computational systems biology approaches to anti-angiogenic cancer therapeutics

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

Do endothelial cells dream of eclectic shape?

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

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

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

Quantitative fluorescent profiling of VEGFRs reveals tumor cell and endothelial cell heterogeneity in breast cancer xenografts

Plasma membrane-localized vascular endothelial growth factor receptors (VEGFR) play a critical role in transducing VEGF signaling toward pro and antiangiogenic outcomes and quantitative characterization of these receptors is critical toward identifying biomarkers for antiangiogenic therapies, understanding mechanisms of action of antiangiogenic drugs, and advancing predictive computational models. While in vitro analysis of cell surface-VEGFRs has been performed, little is known about the levels of cell surface-VEGFR on tumor cells. Therefore, we inoculate nude mice with the human triple-negative breast cancer, MDA-MB-231, cell line; isolate human tumor cells and mouse tumor endothelial cells from xenografts; and quantitatively characterize the VEGFR localization on these cells. We observe 15,000 surface-VEGFR1/tumor endothelial cell versus 8200 surface-VEGFR1/tumor endothelial cell at 3 and 6 weeks of tumor growth, respectively; and we quantify 1200-1700 surface-VEGFR2/tumor endothelial cell. The tumor cell levels of VEGFR1 and VEGFR2 are relatively constant between 3 and 6 weeks: 2000-2200 surface-VEGFR1/tumor cell and ~1000 surface-VEGFR2/tumor cell. Cell-by-cell analysis provides additional insight into tumor heterogeneity by identifying four cellular subpopulations based on size and levels of cell membrane-localized VEGFR. Furthermore, when these ex vivo data are compared to in vitro data, we observe little to no VEGFRs on MDA-MB-231 cells, and the MDA-MB-231 VEGFR surface levels are not regulated by a saturating dose of VEGF. Overall, the quantification of these dissimilarities for the first time in tumor provides insight into the balance of modulatory (VEGFR1) and proangiogenic (VEGFR2) receptors.

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

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

A systems biology view of blood vessel growth and remodelling

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