Structural and practical identifiability of contrast transport models for DCE-MRI.

Contrast transport models are widely used to quantify blood flow and transport in dynamic contrast-enhanced magnetic resonance imaging. These models analyze the time course of the contrast agent concentration, providing diagnostic and prognostic value for many biological systems. Thus, ensuring accuracy and repeatability of the model parameter estimation is a fundamental concern. In this work, we analyze the structural and practical identifiability of a class of nested compartment models pervasively used in analysis of MRI data. We combine artificial and real data to study the role of noise in model parameter estimation. We observe that although all the models are structurally identifiable, practical identifiability strongly depends on the data characteristics. We analyze the impact of increasing data noise on parameter identifiability and show how the latter can be recovered with increased data quality. To complete the analysis, we show that the results do not depend on specific tissue characteristics or the type of enhancement patterns of contrast agent signal.

Dose-dependent thresholds of dexamethasone destabilize CAR T-cell treatment efficacy.

Chimeric antigen receptor (CAR) T-cell therapy is potentially an effective targeted immunotherapy for glioblastoma, yet there is presently little known about the efficacy of CAR T-cell treatment when combined with the widely used anti-inflammatory and immunosuppressant glucocorticoid, dexamethasone. Here we present a mathematical model-based analysis of three patient-derived glioblastoma cell lines treated in vitro with CAR T-cells and dexamethasone. Advanced in vitro experimental cell killing assay technologies allow for highly resolved temporal dynamics of tumor cells treated with CAR T-cells and dexamethasone, making this a valuable model system for studying the rich dynamics of nonlinear biological processes with translational applications. We model the system as a nonautonomous, two-species predator-prey interaction of tumor cells and CAR T-cells, with explicit time-dependence in the clearance rate of dexamethasone. Using time as a bifurcation parameter, we show that (1) dexamethasone destabilizes coexistence equilibria between CAR T-cells and tumor cells in a dose-dependent manner and (2) as dexamethasone is cleared from the system, a stable coexistence equilibrium returns in the form of a Hopf bifurcation. With the model fit to experimental data, we demonstrate that high concentrations of dexamethasone antagonizes CAR T-cell efficacy by exhausting, or reducing the activity of CAR T-cells, and by promoting tumor cell growth. Finally, we identify a critical threshold in the ratio of CAR T-cell death to CAR T-cell proliferation rates that predicts eventual treatment success or failure that may be used to guide the dose and timing of CAR T-cell therapy in the presence of dexamethasone in patients.

Roadmap on plasticity and epigenetics in cancer.

The role of plasticity and epigenetics in shaping cancer evolution and response to therapy has taken center stage with recent technological advances including single cell sequencing. This roadmap article is focused on state-of-the-art mathematical and experimental approaches to interrogate plasticity in cancer, and addresses the following themes and questions: is there a formal overarching framework that encompasses both non-genetic plasticity and mutation-driven somatic evolution? How do we measure and model the role of the microenvironment in influencing/controlling non-genetic plasticity? How can we experimentally study non-genetic plasticity? Which mathematical techniques are required or best suited? What are the clinical and practical applications and implications of these concepts?

From cells to tissue: How cell scale heterogeneity impacts glioblastoma growth and treatment response.

Glioblastomas are aggressive primary brain tumors known for their inter- and intratumor heterogeneity. This disease is uniformly fatal, with intratumor heterogeneity the major reason for treatment failure and recurrence. Just like the nature vs nurture debate, heterogeneity can arise from intrinsic or environmental influences. Whilst it is impossible to clinically separate observed behavior of cells from their environmental context, using a mathematical framework combined with multiscale data gives us insight into the relative roles of variation from different sources. To better understand the implications of intratumor heterogeneity on therapeutic outcomes, we created a hybrid agent-based mathematical model that captures both the overall tumor kinetics and the individual cellular behavior. We track single cells as agents, cell density on a coarser scale, and growth factor diffusion and dynamics on a finer scale over time and space. Our model parameters were fit utilizing serial MRI imaging and cell tracking data from ex vivo tissue slices acquired from a growth-factor driven glioblastoma murine model. When fitting our model to serial imaging only, there was a spectrum of equally-good parameter fits corresponding to a wide range of phenotypic behaviors. When fitting our model using imaging and cell scale data, we determined that environmental heterogeneity alone is insufficient to match the single cell data, and intrinsic heterogeneity is required to fully capture the migration behavior. The wide spectrum of in silico tumors also had a wide variety of responses to an application of an anti-proliferative treatment. Recurrent tumors were generally less proliferative than pre-treatment tumors as measured via the model simulations and validated from human GBM patient histology. Further, we found that all tumors continued to grow with an anti-migratory treatment alone, but the anti-proliferative/anti-migratory combination generally showed improvement over an anti-proliferative treatment alone. Together our results emphasize the need to better understand the underlying phenotypes and tumor heterogeneity present in a tumor when designing therapeutic regimens.

Dissecting Response to Cancer Immunotherapy by Applying Bayesian Network Analysis to Flow Cytometry Data.

Cancer immunotherapy, specifically immune checkpoint blockade, has been found to be effective in the treatment of metastatic cancers. However, only a subset of patients achieve clinical responses. Elucidating pretreatment biomarkers predictive of sustained clinical response is a major research priority. Another research priority is evaluating changes in the immune system before and after treatment in responders vs. nonresponders. Our group has been studying immune networks as an accurate reflection of the global immune state. Flow cytometry (FACS, fluorescence-activated cell sorting) data characterizing immune cell panels in peripheral blood mononuclear cells (PBMC) from gastroesophageal adenocarcinoma (GEA) patients were used to analyze changes in immune networks in this setting. Here, we describe a novel computational pipeline to perform secondary analyses of FACS data using systems biology/machine learning techniques and concepts. The pipeline is centered around comparative Bayesian network analyses of immune networks and is capable of detecting strong signals that conventional methods (such as FlowJo manual gating) might miss. Future studies are planned to validate and follow up the immune biomarkers (and combinations/interactions thereof) associated with clinical responses identified with this computational pipeline.

Intrinsic brain activity changes associated with adjuvant chemotherapy in older women with breast cancer: a pilot longitudinal study.

PURPOSE: Older cancer patients are at increased risk of cancer-related cognitive impairment. The purpose of this study was to assess the alterations in intrinsic brain activity associated with adjuvant chemotherapy in older women with breast cancer. METHODS: Chemotherapy treatment (CT) group included sixteen women aged ≥ 60 years (range 60-82 years) with stage I-III breast cancers, who underwent both resting-state functional magnetic resonance imaging (rs-fMRI) and neuropsychological testing with NIH Toolbox for Cognition before adjuvant chemotherapy, at time point 1 (TP1), and again within 1 month after completing chemotherapy, at time point 2 (TP2). Fourteen age- and sex-matched healthy controls (HC) underwent the same assessments at matched intervals. Three voxel-wise rs-fMRI parameters: amplitude of low-frequency fluctuation (ALFF), fractional ALFF (fALFF), and regional homogeneity, were computed at each time point. The changes in rs-fMRI parameters from TP1 to TP2 for each group, the group differences in changes (the CT group vs. the HC group), and the group difference in the baseline rs-fMRI parameters were assessed. In addition, correlative analysis between the rs-fMRI parameters and neuropsychological testing scores was also performed. RESULTS: In the CT group, one brain region, which included parts of the bilateral subcallosal gyri and right anterior cingulate gyrus, displayed increased ALFF from TP1 to TP2 (cluster p-corrected = 0.024); another brain region in the left precuneus displayed decreased fALFF from TP1 to TP2 (cluster level p-corrected = 0.025). No significant changes in the rs-fMRI parameters from TP1 to TP2 were observed in the HC group. Although ALFF and fALFF alterations were observed only in the CT group, none of the between-group differences in rs-fMRI parameter changes reached statistical significance. CONCLUSIONS: Our study results of ALFF and fALFF alterations in the chemotherapy-treated women suggest that adjuvant chemotherapy may affect intrinsic brain activity in older women with breast cancer.

The 2019 mathematical oncology roadmap.

Whether the nom de guerre is Mathematical Oncology, Computational or Systems Biology, Theoretical Biology, Evolutionary Oncology, Bioinformatics, or simply Basic Science, there is no denying that mathematics continues to play an increasingly prominent role in cancer research. Mathematical Oncology-defined here simply as the use of mathematics in cancer research-complements and overlaps with a number of other fields that rely on mathematics as a core methodology. As a result, Mathematical Oncology has a broad scope, ranging from theoretical studies to clinical trials designed with mathematical models. This Roadmap differentiates Mathematical Oncology from related fields and demonstrates specific areas of focus within this unique field of research. The dominant theme of this Roadmap is the personalization of medicine through mathematics, modelling, and simulation. This is achieved through the use of patient-specific clinical data to: develop individualized screening strategies to detect cancer earlier; make predictions of response to therapy; design adaptive, patient-specific treatment plans to overcome therapy resistance; and establish domain-specific standards to share model predictions and to make models and simulations reproducible. The cover art for this Roadmap was chosen as an apt metaphor for the beautiful, strange, and evolving relationship between mathematics and cancer.

Assessing brain volume changes in older women with breast cancer receiving adjuvant chemotherapy: a brain magnetic resonance imaging pilot study.

BACKGROUND: Cognitive decline is among the most feared treatment-related outcomes of older adults with cancer. The majority of older patients with breast cancer self-report cognitive problems during and after chemotherapy. Prior neuroimaging research has been performed mostly in younger patients with cancer. The purpose of this study was to evaluate longitudinal changes in brain volumes and cognition in older women with breast cancer receiving adjuvant chemotherapy. METHODS: Women aged ≥ 60 years with stage I-III breast cancer receiving adjuvant chemotherapy and age-matched and sex-matched healthy controls were enrolled. All participants underwent neuropsychological testing with the US National Institutes of Health (NIH) Toolbox for Cognition and brain magnetic resonance imaging (MRI) prior to chemotherapy, and again around one month after the last infusion of chemotherapy. Brain volumes were measured using Neuroreader™ software. Longitudinal changes in brain volumes and neuropsychological scores were analyzed utilizing linear mixed models. RESULTS: A total of 16 patients with breast cancer (mean age 67.0, SD 5.39 years) and 14 age-matched and sex-matched healthy controls (mean age 67.8, SD 5.24 years) were included: 7 patients received docetaxel and cyclophosphamide (TC) and 9 received chemotherapy regimens other than TC (non-TC). There were no significant differences in segmented brain volumes between the healthy control group and the chemotherapy group pre-chemotherapy (p > 0.05). Exploratory hypothesis generating analyses focusing on the effect of the chemotherapy regimen demonstrated that the TC group had greater volume reduction in the temporal lobe (change = - 0.26) compared to the non-TC group (change = 0.04, p for interaction = 0.02) and healthy controls (change = 0.08, p for interaction = 0.004). Similarly, the TC group had a decrease in oral reading recognition scores (change = - 6.94) compared to the non-TC group (change = - 1.21, p for interaction = 0.07) and healthy controls (change = 0.09, p for interaction = 0.02). CONCLUSIONS: There were no significant differences in segmented brain volumes between the healthy control group and the chemotherapy group; however, exploratory analyses demonstrated a reduction in both temporal lobe volume and oral reading recognition scores among patients on the TC regimen. These results suggest that different chemotherapy regimens may have differential effects on brain volume and cognition. Future, larger studies focusing on older adults with cancer on different treatment regimens are needed to confirm these findings. TRIAL REGISTRATION: ClinicalTrials.gov, NCT01992432 . Registered on 25 November 2013. Retrospectively registered.

Gray matter density reduction associated with adjuvant chemotherapy in older women with breast cancer.

PURPOSE: The purpose of this study was to evaluate longitudinal changes in brain gray matter density (GMD) before and after adjuvant chemotherapy in older women with breast cancer. METHODS: We recruited 16 women aged ≥ 60 years with stage I-III breast cancers receiving adjuvant chemotherapy (CT) and 15 age- and sex-matched healthy controls (HC). The CT group underwent brain MRI and the NIH Toolbox for Cognition testing prior to adjuvant chemotherapy (time point 1, TP1) and within 1 month after chemotherapy (time point 2, TP2). The HC group underwent the same assessments at matched intervals. GMD was evaluated with the voxel-based morphometry. RESULTS: The mean age was 67 years in the CT group and 68.5 years in the HC group. There was significant GMD reduction within the chemotherapy group from TP1 to TP2. Compared to the HC group, the CT group displayed statistically significantly greater GMD reductions from TP1 to TP2 in the brain regions involving the left anterior cingulate gyrus, right insula, and left middle temporal gyrus (pFWE(family-wise error)-corrected < 0.05). The baseline GMD in left insula was positively correlated with the baseline list-sorting working memory score in the HC group (pFWE-corrected < 0.05). No correlation was observed for the changes in GMD with the changes in cognitive testing scores from TP1 to TP2 (pFWE-corrected < 0.05). CONCLUSIONS: Our findings indicate that GMD reductions were associated with adjuvant chemotherapy in older women with breast cancer. Future studies are needed to understand the clinical significance of the neuroimaging findings. This study is registered on ClinicalTrials.gov (NCT01992432).

Simulating PDGF-Driven Glioma Growth and Invasion in an Anatomically Accurate Brain Domain.

Gliomas are the most common of all primary brain tumors. They are characterized by their diffuse infiltration of the brain tissue and are uniformly fatal, with glioblastoma being the most aggressive form of the disease. In recent years, the over-expression of platelet-derived growth factor (PDGF) has been shown to produce tumors in experimental rodent models that closely resemble this human disease, specifically the proneural subtype of glioblastoma. We have previously modeled this system, focusing on the key attribute of these experimental tumors-the "recruitment" of oligodendroglial progenitor cells (OPCs) to participate in tumor formation by PDGF-expressing retrovirally transduced cells-in one dimension, with spherical symmetry. However, it has been observed that these recruitable progenitor cells are not uniformly distributed throughout the brain and that tumor cells migrate at different rates depending on the material properties in different regions of the brain. Here we model the differential diffusion of PDGF-expressing and recruited cell populations via a system of partial differential equations with spatially variable diffusion coefficients and solve the equations in two spatial dimensions on a mouse brain atlas using a flux-differencing numerical approach. Simulations of our in silico model demonstrate qualitative agreement with the observed tumor distribution in the experimental animal system. Additionally, we show that while there are higher concentrations of OPCs in white matter, the level of recruitment of these plays little role in the appearance of "white matter disease," where the tumor shows a preponderance for white matter. Instead, simulations show that this is largely driven by the ratio of the diffusion rate in white matter as compared to gray. However, this ratio has less effect on the speed of tumor growth than does the degree of OPC recruitment in the tumor. It was observed that tumor simulations with greater degrees of recruitment grow faster and develop more nodular tumors than if there is no recruitment at all, similar to our prior results from implementing our model in one dimension. Combined, these results show that recruitment remains an important consideration in understanding and slowing glioma growth.

Patient-specific mathematical neuro-oncology: using a simple proliferation and invasion tumor model to inform clinical practice.

Glioblastoma multiforme (GBM) is the most common malignant primary brain tumor associated with a poor median survival of 15-18 months, yet there is wide heterogeneity across and within patients. This heterogeneity has been the source of significant clinical challenges facing patients with GBM and has hampered the drive toward more precision or personalized medicine approaches to treating these challenging tumors. Over the last two decades, the field of Mathematical Neuro-oncology has grown out of desire to use (often patient-specific) mathematical modeling to better treat GBMs. Here, we will focus on a series of clinically relevant results using patient-specific mathematical modeling. The core model at the center of these results incorporates two hallmark features of GBM, proliferation [Formula: see text] and invasion (D), as key parameters. Based on routinely obtained magnetic resonance images, each patient's tumor can be characterized using these two parameters. The Proliferation-Invasion (PI) model uses [Formula: see text] and D to create patient-specific growth predictions. The PI model, its predictions, and parameters have been used in a number of ways to derive biological insight. Beyond predicting growth, the PI model has been utilized to identify patients who benefit from different surgery strategies, to prognosticate response to radiation therapy, to develop a treatment response metric, and to connect clinical imaging features and genetic information. Demonstration of the PI model's clinical relevance supports the growing role for it and other mathematical models in routine clinical practice.

Mathematical Modeling Unveils Optimization Strategies for Targeted Radionuclide Therapy of Blood Cancers.

Targeted radionuclide therapy is based on injections of cancer-specific molecules conjugated with radioactive nuclides. Despite the specificity of this treatment, it is not devoid of side-effects limiting its use and is especially harmful for rapidly proliferating organs well perfused by blood, like bone marrow. Optimization of radioconjugates administration accounting for toxicity constraints can increase treatment efficacy. Based on our experiments on disseminated multiple myeloma mouse model treated by 225Ac-DOTA-daratumumab, we developed a mathematical model which investigation highlighted the following principles for optimization of targeted radionuclide therapy. 1) Nuclide to antibody ratio importance. The density of radioconjugates on cancer cells determines the density of radiation energy deposited in them. Low labeling ratio as well as accumulation of unlabeled antibodies and antibodies attached to decay products in the bloodstream can mitigate cancer radiation damage due to excessive occupation of specific receptors by antibodies devoid of radioactive nuclides. 2) Cancer binding capacity-based dosing. The total number of specific receptors on cancer cells is a critical factor for treatment optimization, which estimation may allow increasing treatment efficacy close to its theoretical limit. Injection of doses significantly exceeding cancer binding capacity should be avoided since radioconjugates remaining in the bloodstream have negligible efficacy to toxicity ratio. 3) Particle range-guided multi-dosing. The use of short-range particle emitters and high-affinity antibodies can allow for robust treatment optimization via initial saturation of cancer binding capacity, enabling redistribution of further injected radioconjugates and deposited dose towards still viable cells that continue expressing specific receptors.

Mathematical Modeling Unveils Optimization Strategies for Targeted Radionuclide Therapy of Blood Cancers.

Targeted radionuclide therapy is based on injections of cancer-specific molecules conjugated with radioactive nuclides. Despite the specificity of this treatment, it is not devoid of side-effects limiting its use and is especially harmful for rapidly proliferating organs well perfused by blood, like bone marrow. Optimization of radioconjugates administration accounting for toxicity constraints can increase treatment efficacy. Based on our experiments on disseminated multiple myeloma mouse model treated by 225Ac-DOTA-daratumumab, we developed a mathematical model which investigation highlighted the following principles for optimization of targeted radionuclide therapy. 1) Nuclide to antibody ratio importance. The density of radioconjugates on cancer cells determines the density of radiation energy deposited in them. Low labeling ratio as well as accumulation of unlabeled antibodies and antibodies attached to decay products in the bloodstream can mitigate cancer radiation damage due to excessive occupation of specific receptors by antibodies devoid of radioactive nuclides. 2) Cancer binding capacity-based dosing. The rate of binding of drug to cancer cells depends on the total number of their specific receptors, which therefore can be estimated from the pharmacokinetic curve of diagnostic radioconjugates. Injection of doses significantly exceeding cancer binding capacity should be avoided since radioconjugates remaining in the bloodstream have negligible efficacy to toxicity ratio. 3) Particle range-guided multi-dosing. The use of short-range particle emitters and high-affinity antibodies allows for robust treatment optimization via initial saturation of cancer binding capacity, enabling redistribution of further injected radioconjugates and deposited dose towards still viable cells that continue expressing specific receptors.

Model discovery approach enables noninvasive measurement of intra-tumoral fluid transport in dynamic MRI.

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) is a routine method to noninvasively quantify perfusion dynamics in tissues. The standard practice for analyzing DCE-MRI data is to fit an ordinary differential equation to each voxel. Recent advances in data science provide an opportunity to move beyond existing methods to obtain more accurate measurements of fluid properties. Here, we developed a localized convolutional function regression that enables simultaneous measurement of interstitial fluid velocity, diffusion, and perfusion in 3D. We validated the method computationally and experimentally, demonstrating accurate measurement of fluid dynamics in situ and in vivo. Applying the method to human MRIs, we observed tissue-specific differences in fluid dynamics, with an increased fluid velocity in breast cancer as compared to brain cancer. Overall, our method represents an improved strategy for studying interstitial flows and interstitial transport in tumors and patients. We expect that our method will contribute to the better understanding of cancer progression and therapeutic response.

State-transition modeling of blood transcriptome predicts disease evolution and treatment response in chronic myeloid leukemia.

Chronic myeloid leukemia (CML) is initiated and maintained by BCR::ABL which is clinically targeted using tyrosine kinase inhibitors (TKIs). TKIs can induce long-term remission but are also not curative. Thus, CML is an ideal system to test our hypothesis that transcriptome-based state-transition models accurately predict cancer evolution and treatment response. We collected time-sequential blood samples from tetracycline-off (Tet-Off) BCR::ABL-inducible transgenic mice and wild-type controls. From the transcriptome, we constructed a CML state-space and a three-well leukemogenic potential landscape. The potential's stable critical points defined observable disease states. Early states were characterized by anti-CML genes opposing leukemia; late states were characterized by pro-CML genes. Genes with expression patterns shaped similarly to the potential landscape were identified as drivers of disease transition. Re-introduction of tetracycline to silence the BCR::ABL gene returned diseased mice transcriptomes to a near healthy state, without reaching it, suggesting parts of the transition are irreversible. TKI only reverted the transcriptome to an intermediate disease state, without approaching a state of health; disease relapse occurred soon after treatment. Using only the earliest time-point as initial conditions, our state-transition models accurately predicted both disease progression and treatment response, supporting this as a potentially valuable approach to time clinical intervention, before phenotypic changes become detectable.

Designing combination therapies for cancer treatment: application of a mathematical framework combining CAR T-cell immunotherapy and targeted radionuclide therapy.

Introduction: Cancer combination treatments involving immunotherapies with targeted radiation therapy are at the forefront of treating cancers. However, dosing and scheduling of these therapies pose a challenge. Mathematical models provide a unique way of optimizing these therapies. Methods: Using a preclinical model of multiple myeloma as an example, we demonstrate the capability of a mathematical model to combine these therapies to achieve maximum response, defined as delay in tumor growth. Data from mice studies with targeted radionuclide therapy (TRT) and chimeric antigen receptor (CAR)-T cell monotherapies and combinations with different intervals between them was used to calibrate mathematical model parameters. The dependence of progression-free survival (PFS), overall survival (OS), and the time to minimum tumor burden on dosing and scheduling was evaluated. Different dosing and scheduling schemes were evaluated to maximize the PFS and optimize timings of TRT and CAR-T cell therapies. Results: Therapy intervals that were too close or too far apart are shown to be detrimental to the therapeutic efficacy, as TRT too close to CAR-T cell therapy results in radiation related CAR-T cell killing while the therapies being too far apart result in tumor regrowth, negatively impacting tumor control and survival. We show that splitting a dose of TRT or CAR-T cells when administered in combination is advantageous only if the first therapy delivered can produce a significant benefit as a monotherapy. Discussion: Mathematical models are crucial tools for optimizing the delivery of cancer combination therapy regimens with application along the lines of achieving cure, maximizing survival or minimizing toxicity.

Challenges with sirolimus experimental data to inform QSP model of post-transplantation cyclophosphamide regimens.

Dose optimization of sirolimus may further improve outcomes in allogeneic hematopoietic cell transplant (HCT) patients receiving post-transplantation cyclophosphamide (PTCy) to prevent graft-versus-host disease (GVHD). Sirolimus exposure-response association studies in HCT patients (i.e., the association of trough concentration with clinical outcomes) have been conflicting. Sirolimus has important effects on T-cells, including conventional (Tcons) and regulatory T-cells (Tregs), both of which have been implicated in the mechanisms by which PTCy prevents GVHD, but there is an absence of validated biomarkers of sirolimus effects on these cell subsets. Considering the paucity of existing biomarkers and the complexities of the immune system, we conducted a literature review to inform a quantitative systems pharmacology (QSP) model of GVHD. The published literature presented multiple challenges. The sirolimus pharmacokinetic models insufficiently describe sirolimus distribution to relevant physiological compartments. Despite multiple publications describing sirolimus effects on Tcons and Tregs in preclinical and human ex vivo models, consistent parameters relating sirolimus concentrations to circulating Tcons and Tregs could not be found. Each aspect presents a challenge in building a QSP model of sirolimus and its temporal effects on T-cell subsets and GVHD prevention. To optimize GVHD prevention regimens, phase I studies and systematic studies of immunosuppressant concentration-effect association are needed for QSP modeling.

Systems profiling reveals recurrently dysregulated cytokine signaling responses in ER+ breast cancer patients' blood.

Cytokines operate in concert to maintain immune homeostasis and coordinate immune responses. In cases of ER+ breast cancer, peripheral immune cells exhibit altered responses to several cytokines, and these alterations are correlated strongly with patient outcomes. To develop a systems-level understanding of this dysregulation, we measured a panel of cytokine responses and receptor abundances in the peripheral blood of healthy controls and ER+ breast cancer patients across immune cell types. Using tensor factorization to model this multidimensional data, we found that breast cancer patients exhibited widespread alterations in response, including drastically reduced response to IL-10 and heightened basal levels of pSmad2/3 and pSTAT4. ER+ patients also featured upregulation of PD-L1, IL6Rα, and IL2Rα, among other receptors. Despite this, alterations in response to cytokines were not explained by changes in receptor abundances. Thus, tensor factorization helped to reveal a coordinated reprogramming of the immune system that was consistent across our cohort.

A novel class of inhibitors that disrupts the stability of integrin heterodimers identified by CRISPR-tiling-instructed genetic screens.

The plasma membrane is enriched for receptors and signaling proteins that are accessible from the extracellular space for pharmacological intervention. Here we conducted a series of CRISPR screens using human cell surface proteome and integrin family libraries in multiple cancer models. Our results identified ITGAV (integrin αV) and its heterodimer partner ITGB5 (integrin ß5) as the essential integrin α/ß pair for cancer cell expansion. High-density CRISPR gene tiling further pinpointed the integral pocket within the ß-propeller domain of ITGAV for integrin αVß5 dimerization. Combined with in silico compound docking, we developed a CRISPR-Tiling-Instructed Computer-Aided (CRISPR-TICA) pipeline for drug discovery and identified Cpd_AV2 as a lead inhibitor targeting the ß-propeller central pocket of ITGAV. Cpd_AV2 treatment led to rapid uncoupling of integrin αVß5 and cellular apoptosis, providing a unique class of therapeutic action that eliminates the integrin signaling via heterodimer dissociation. We also foresee the CRISPR-TICA approach to be an accessible method for future drug discovery studies.

Locoregional delivery of IL-13Rα2-targeting CAR-T cells in recurrent high-grade glioma: a phase 1 trial.

Chimeric antigen receptor T cell (CAR-T) therapy is an emerging strategy to improve treatment outcomes for recurrent high-grade glioma, a cancer that responds poorly to current therapies. Here we report a completed phase I trial evaluating IL-13Rα2-targeted CAR-T cells in 65 patients with recurrent high-grade glioma, the majority being recurrent glioblastoma (rGBM). Primary objectives were safety and feasibility, maximum tolerated dose/maximum feasible dose and a recommended phase 2 dose plan. Secondary objectives included overall survival, disease response, cytokine dynamics and tumor immune contexture biomarkers. This trial evolved to evaluate three routes of locoregional T cell administration (intratumoral (ICT), intraventricular (ICV) and dual ICT/ICV) and two manufacturing platforms, culminating in arm 5, which utilized dual ICT/ICV delivery and an optimized manufacturing process. Locoregional CAR-T cell administration was feasible and well tolerated, and as there were no dose-limiting toxicities across all arms, a maximum tolerated dose was not determined. Probable treatment-related grade 3+ toxicities were one grade 3 encephalopathy and one grade 3 ataxia. A clinical maximum feasible dose of 200 × 106 CAR-T cells per infusion cycle was achieved for arm 5; however, other arms either did not test or achieve this dose due to manufacturing feasibility. A recommended phase 2 dose will be refined in future studies based on data from this trial. Stable disease or better was achieved in 50% (29/58) of patients, with two partial responses, one complete response and a second complete response after additional CAR-T cycles off protocol. For rGBM, median overall survival for all patients was 7.7 months and for arm 5 was 10.2 months. Central nervous system increases in inflammatory cytokines, including IFNγ, CXCL9 and CXCL10, were associated with CAR-T cell administration and bioactivity. Pretreatment intratumoral CD3 T cell levels were positively associated with survival. These findings demonstrate that locoregional IL-13Rα2-targeted CAR-T therapy is safe with promising clinical activity in a subset of patients. ClinicalTrials.gov Identifier: NCT02208362 .