In silico mouse study identifies tumour growth kinetics as biomarkers for the outcome of anti-angiogenic treatment

Combining Multikinase Tyrosine Kinase Inhibitors Targeting the Vascular Endothelial Growth Factor and Cluster of Differentiation 47 Signaling Pathways Is Predicted to Increase the Efficacy of Antiangiogenic Combination Therapies

Angiogenesis is a critical step in tumor growth, development, and invasion. Nascent tumor cells secrete vascular endothelial growth factor (VEGF) that significantly remodels the tumor microenvironment through interaction with multiple receptors on vascular endothelial cells, including type 2 VEGF receptor (VEGFR2). The complex pathways initiated by VEGF binding to VEGFR2 lead to enhanced proliferation, survival, and motility of vascular endothelial cells and formation of a new vascular network, enabling tumor growth. Antiangiogenic therapies that inhibit VEGF signaling pathways were among the first drugs that targeted stroma rather than tumor cells. Despite improvements in progression-free survival and higher response rates relative to chemotherapy in some types of solid tumors, the impact on overall survival (OS) has been limited, with the majority of tumors eventually relapsing due to resistance or activation of alternate angiogenic pathways. Here, we developed a molecularly detailed computational model of endothelial cell signaling and angiogenesis-driven tumor growth to investigate combination therapies targeting different nodes of the endothelial VEGF/VEGFR2 signaling pathway. Simulations predicted a strong threshold-like behavior in extracellular signal-regulated kinases 1/2 (ERK1/2) activation relative to phosphorylated VEGFR2 levels, as continuous inhibition of at least 95% of receptors was necessary to abrogate phosphorylated ERK1/2 (pERK1/2). Combinations with mitogen-activated protein kinase/ERK kinase (MEK) and spingosine-1-phosphate inhibitors were found to be effective in overcoming the ERK1/2 activation threshold and abolishing activation of the pathway. Modeling results also identified a mechanism of resistance whereby tumor cells could reduce pERK1/2 sensitivity to inhibitors of VEGFR2 by upregulation of Raf, MEK, and sphingosine kinase 1 (SphK1), thus highlighting the need for deeper investigation of the dynamics of the crosstalk between VEGFR2 and SphK1 pathways. Inhibition of VEGFR2 phosphorylation was found to be more effective at blocking protein kinase B, also known as AKT, activation; however, to effectively abolish AKT activation, simulations identified Axl autophosphorylation or the Src kinase domain as potent targets. Simulations also supported activating cluster of differentiation 47 (CD47) on endothelial cells as an effective combination partner with tyrosine kinase inhibitors to inhibit angiogenesis signaling and tumor growth. Virtual patient simulations supported the effectiveness of CD47 agonism in combination with inhibitors of VEGFR2 and SphK1 pathways. Overall, the rule-based system model developed here provides new insights, generates novel hypothesis, and makes predictions regarding combinations that may enhance the OS with currently approved antiangiogenic therapies.

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

Targeting a cell surface vitamin D receptor on tumor-associated macrophages in triple-negative breast cancer

Triple-negative breast cancer (TNBC) is an aggressive tumor with limited treatment options and poor prognosis. We applied the in vivo phage display technology to isolate peptides homing to the immunosuppressive cellular microenvironment of TNBC as a strategy for non-malignant target discovery. We identified a cyclic peptide (CSSTRESAC) that specifically binds to a vitamin D receptor, protein disulfide-isomerase A3 (PDIA3) expressed on the cell surface of tumor-associated macrophages (TAM), and targets breast cancer in syngeneic TNBC, non-TNBC xenograft, and transgenic mouse models. Systemic administration of CSSTRESAC to TNBC-bearing mice shifted the cytokine profile toward an antitumor immune response and delayed tumor growth. Moreover, CSSTRESAC enabled ligand-directed theranostic delivery to tumors and a mathematical model confirmed our experimental findings. Finally, in silico analysis showed PDIA3-expressing TAM in TNBC patients. This work uncovers a functional interplay between a cell surface vitamin D receptor in TAM and antitumor immune response that could be therapeutically exploited.

Resolvin D1 reduces cancer growth stimulating a protective neutrophil-dependent recruitment of anti-tumor monocytes

Background

Innovative therapies to target tumor-associated neutrophils (PMN) are of clinical interest, since these cells are centrally involved in cancer inflammation and tumor progression. Resolvin D1 (RvD1) is a lipid autacoid that promotes resolution of inflammation by regulating the activity of distinct immune and non-immune cells. Here, using human papilloma virus (HPV) tumorigenesis as a model, we investigated whether RvD1 modulates PMN to reduce tumor progression.

Methods

Growth-curve assays with multiple cell lines and in vivo grafting of two distinct HPV-positive cells in syngeneic mice were used to determine if RvD1 reduced cancer growth. To investigate if and how RvD1 modulates PMN activities, RNA sequencing and multiplex cytokine ELISA of human PMN in co-culture with HPV-positive cells, coupled with pharmacological depletion of PMN in vivo, were performed. The mouse intratumoral immune cell composition was evaluated through FACS analysis. Growth-curve assays and in vivo pharmacological depletion were used to evaluate anti-tumor activities of human and mouse monocytes, respectively. Bioinformatic analysis of The Cancer Genome Atlas (TCGA) database was exploited to validate experimental findings in patients.

Results

RvD1 decreased in vitro and in vivo proliferation of human and mouse HPV-positive cancer cells through stimulation of PMN anti-tumor activities. In addition, RvD1 stimulated a PMN-dependent recruitment of classical monocytes as key determinant to reduce tumor growth in vivo. In human in vitro systems, exposure of PMN to RvD1 increased the production of the monocyte chemoattractant protein-1 (MCP-1), and enhanced transmigration of classical monocytes, with potent anti-tumor actions, toward HPV-positive cancer cells. Consistently, mining of immune cells infiltration levels in cervical cancer patients from the TCGA database evidenced an enhanced immune reaction and better clinical outcomes in patients with higher intratumoral monocytes as compared to patients with higher PMN infiltration.

Conclusions

RvD1 reduces cancer growth by activating PMN anti-cancer activities and encouraging a protective PMN-dependent recruitment of anti-tumor monocytes. These findings demonstrate efficacy of RvD1 as an innovative therapeutic able to stimulate PMN reprogramming to an anti-cancer phenotype that restrains tumor growth.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-021-01937-3.

Quantitative Systems Pharmacology Modeling of PBMC-Humanized Mouse to Facilitate Preclinical Immuno-oncology Drug Development

Progress in immunotherapy has resulted in explosively increased new therapeutic interventions and they have shown promising results in the treatment of cancer. Animal testing is performed to provide preliminary efficacy and safety data for drugs under development prior to clinical trials. However, translational challenges remain for preclinical studies such as study design and the relevance of animal models to humans. Hence, only a small fraction of cancer patients showed response. The explosion of drug candidates and therapies makes preclinical assessment of every plausible option impossible, but it can be easily tested using Quantitative System Pharmacology (QSP) models. Here, we developed a QSP model for humanized mice. Tumor growth dynamics, T cell dynamics, cytokine release, immune checkpoint expression, and drug administration were modeled and calibrated using experimental data. Tumor growth inhibition data were used for model validation. Pharmacokinetics of T cell engager (TCE), tumor growth profile, T cell expansion in the blood and infiltration into tumor, T cell dissemination from primary tumor, cytokine release profile, and expression of additional PD-L1 induced by IFN-γ were modeled and calibrated using a variety of experimental data and showed good consistency. Mouse-specific response to T cell engager monotherapy also showed the key features of in vivo efficacy of TCE. This novel QSP model, designed for human peripheral blood mononuclear cells (PBMC) engrafted xenograft mice, incorporating the most critical components of the mouse model with key cancer and immune cells, can become an integral part of preclinical drug development.

Mathematical Model Predicts Effective Strategies to Inhibit VEGF-eNOS Signaling

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