Towards a Framework for Predictive Mathematical Modeling of the Biomechanical Forces Causing Brain Tumor Mass-Effect

GBMs present with different growth phenotypes, ranging from invasive lesions without notable mass-effect to strongly displacing lesions that induce mechanical stresses and result in healthy-tissue deformation, midline shift or herniation. Biomechanical forces, such as those resulting from displacive tumor growth, are recognized to shape the tumor environment and to contribute to tumor progression. We therefore expect that biomechanical forces exerted by lesions on the brain parenchyma have implications on the biophysical level, and that they may affect treatment response and outcome. To better understand the role of biomechanics in the formation of different GBM phenotypes we started developing a framework for the predictive mathematical modeling of mechanical tumor-healthy tissue interaction on the macroscopic level. The tumor’s mass-effect is represented by a solid-mechanics model of brain tissue that computes tumor-induced strain based on local tumor cell concentration. The framework allows to seed tumors at multiple locations in a human brain atlas. It simulates tumor evolution over time and across different brain regions using literature-based parameter estimates for tumor cell proliferation, as well as isotropic motility, and mechanical tissue properties. Despite its simplicity, the mathematical model yielded realistic estimates of the mechanical impact of a growing tumor on intra-cranial pressure. However, comparison to publicly available GBM imaging data showed that asymmetric shapes could not be reproduced by isotropic growth assumptions. Here we present and evaluate an extended version of this mechanically-coupled reaction-diffusion model that takes into account tissue anisotropies based on MRI diffusion tensor imaging (MR-DTI). Structural anisotropies in brain tissue have been found to affect the directionality of tumor cell migration and are critical to mechanical behavior. This makes them likely to play a role also in the development of GBM phenotypes

Mathematical deconvolution of CAR T-cell proliferation and exhaustion from real-time killing assay data.

Chimeric antigen receptor (CAR) T-cell therapy has shown promise in the treatment of haematological cancers and is currently being investigated for solid tumours, including high-grade glioma brain tumours. There is a desperate need to quantitatively study the factors that contribute to the efficacy of CAR T-cell therapy in solid tumours. In this work, we use a mathematical model of predator-prey dynamics to explore the kinetics of CAR T-cell killing in glioma: the Chimeric Antigen Receptor T-cell treatment Response in GliOma (CARRGO) model. The model includes rates of cancer cell proliferation, CAR T-cell killing, proliferation, exhaustion, and persistence. We use patient-derived and engineered cancer cell lines with an in vitro real-time cell analyser to parametrize the CARRGO model. We observe that CAR T-cell dose correlates inversely with the killing rate and correlates directly with the net rate of proliferation and exhaustion. This suggests that at a lower dose of CAR T-cells, individual T-cells kill more cancer cells but become more exhausted when compared with higher doses. Furthermore, the exhaustion rate was observed to increase significantly with tumour growth rate and was dependent on level of antigen expression. The CARRGO model highlights nonlinear dynamics involved in CAR T-cell therapy and provides novel insights into the kinetics of CAR T-cell killing. The model suggests that CAR T-cell treatment may be tailored to individual tumour characteristics including tumour growth rate and antigen level to maximize therapeutic benefit

Early Changes in Tumor Perfusion from T1-Weighted Dynamic Contrast-Enhanced MRI following Neural Stem Cell-Mediated Therapy of Recurrent High-Grade Glioma Correlate with Overall Survival

Background. The aim of this study was to correlate T1-weighted dynamic contrast-enhanced MRI- (DCE-MRI-) derived perfusion parameters with overall survival of recurrent high-grade glioma patients who received neural stem cell- (NSC-) mediated enzyme/prodrug gene therapy. Methods. A total of 12 patients were included in this retrospective study. All patients were enrolled in a first-in-human study (NCT01172964) of NSC-mediated therapy for recurrent high-grade glioma. DCE-MRI data from all patients were collected and analyzed at three time points: MRI#1—day 1 postsurgery/treatment, MRI#2— day 7 ± 3 posttreatment, and MRI#3—one-month follow-up. Plasma volume (Vp), permeability (Ktr), and leakage (λtr) perfusion parameters were calculated by fitting a pharmacokinetic model to the DCE-MRI data. The contrast-enhancing (CE) volume was measured from the last dynamic phase acquired in the DCE sequence. Perfusion parameters and CE at each MRI time point were recorded along with their relative change between MRI#2 and MRI#3 (Δ32). Cox regression was used to analyze patient survival. Results. At MRI#1 and at MRI#3, none of the parameters showed a significant correlation with overall survival (OS). However, at MRI#2, CE and λtr were significantly associated with OS (p<0.05). The relative λtr and Vp from timepoint 2 to timepoint 3 (Δ32λtr and Δ32Vp) were each associated with a higher hazard ratio (p<0.05). All parameters were highly correlated, resulting in a multivariate model for OS including only CE at MRI#2 and Δ32Vp, with an R2 of 0.89. Conclusion. The change in perfusion parameter values from 1 week to 1 month following NSC-mediated therapy combined with contrast-enhancing volume may be a useful biomarker to predict overall survival in patients with recurrent high-grade glioma

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

IntroductionCancer 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.MethodsUsing 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.ResultsTherapy 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.DiscussionMathematical 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&nbsp;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

Differentiating Peripherally-Located Small Cell Lung Cancer From Non-small Cell Lung Cancer Using a CT Radiomic Approach

Lung cancer can be classified into two main categories: small cell lung cancer (SCLC) and non-small cell lung cancer (NSCLC), which are different in treatment strategy and survival probability. The lung CT images of SCLC and NSCLC are similar such that their subtle differences are hardly visually discernible by the human eye through conventional imaging evaluation. We hypothesize that SCLC/NSCLC differentiation could be achieved via computerized image feature analysis and classification in feature space, as termed a radiomic model. The purpose of this study was to use CT radiomics to differentiate SCLC from NSCLC adenocarcinoma. Patients with primary lung cancer, either SCLC or NSCLC adenocarcinoma, were retrospectively identified. The post-diagnosis pre-treatment lung CT images were used to segment the lung cancers. Radiomic features were extracted from histogram-based statistics, textural analysis of tumor images and their wavelet transforms. A minimal-redundancy-maximal-relevance method was used for feature selection. The predictive model was constructed with a multilayer artificial neural network. The performance of the SCLC/NSCLC adenocarcinoma classifier was evaluated by the area under the receiver operating characteristic curve (AUC). Our study cohort consisted of 69 primary lung cancer patients with SCLC (n = 35; age mean ± SD = 66.91± 9.75 years), and NSCLC adenocarcinoma (n = 34; age mean ± SD = 58.55 ± 11.94 years). The SCLC group had more male patients and smokers than the NSCLC group (P \u3c 0.05). Our SCLC/NSCLC classifier achieved an overall performance of AUC of 0.93 (95% confidence interval = [0.85, 0.97]), sensitivity = 0.85, and specificity = 0.85). Adding clinical data such as smoking history could improve the performance slightly. The top ranking radiomic features were mostly textural features. Our results showed that CT radiomics could quantitatively represent tumor heterogeneity and therefore could be used to differentiate primary lung cancer subtypes with satisfying results. CT image processing with the wavelet transformation technique enhanced the radiomic features for SCLC/NSCLC classification. Our pilot study should motivate further investigation of radiomics as a non-invasive approach for early diagnosis and treatment of lung cancer

Quantitative Evaluation of Intraventricular Delivery of Therapeutic Neural Stem Cells to Orthotopic Glioma

Neural stem cells (NSCs) are inherently tumor-tropic, which allows them to migrate through normal tissue and selectively localize to invasive tumor sites in the brain. We have engineered a clonal, immortalized allogeneic NSC line (HB1.F3.CD21; CD-NSCs) that maintains its stem-like properties, a normal karyotype and is HLA Class II negative. It is genetically and functionally stable over time and multiple passages, and has demonstrated safety in phase I glioma trials. These properties enable the production of an “off-the-shelf” therapy that can be readily available for patient treatment. There are multiple factors contributing to stem cell tumor-tropism, and much remains to be elucidated. The route of NSC delivery and the distribution of NSCs at tumor sites are key factors in the development of effective cell-based therapies. Stem cells can be engineered to deliver and/or produce many different therapeutic agents, including prodrug activating enzymes (which locally convert systemically administered prodrugs to active chemotherapeutic agents); oncolytic viruses; tumor-targeted antibodies; therapeutic nanoparticles; and extracellular vesicles that contain therapeutic oligonucleotides. By targeting these therapeutics selectively to tumor foci, we aim to minimize toxicity to normal tissues and maximize therapeutic benefits. In this manuscript, we demonstrate that NSCs administered via intracerebral/ventricular (IVEN) routes can migrate efficiently toward single or multiple tumor foci. IVEN delivery will enable repeat administrations for patients through an Ommaya reservoir, potentially resulting in improved therapeutic outcomes. In our preclinical studies using various glioma lines, we have quantified NSC migration and distribution in mouse brains and have found robust migration of our clinically relevant HB1.F3.CD21 NSC line toward invasive tumor foci, irrespective of their origin. These results establish proof-of-concept and demonstrate the potential of developing a multitude of therapeutic options using modified NSCs

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)

Patient-Specific Metrics of Invasiveness Reveal Significant Prognostic Benefit of Resection in a Predictable Subset of Gliomas

Object Malignant gliomas are incurable, primary brain neoplasms noted for their potential to extensively invade brain parenchyma. Current methods of clinical imaging do not elucidate the full extent of brain invasion, making it difficult to predict which, if any, patients are likely to benefit from gross total resection. Our goal was to apply a mathematical modeling approach to estimate the overall tumor invasiveness on a patient-by-patient basis and determine whether gross total resection would improve survival in patients with relatively less invasive gliomas. Methods In 243 patients presenting with contrast-enhancing gliomas, estimates of the relative invasiveness of each patient's tumor, in terms of the ratio of net proliferation rate of the glioma cells to their net dispersal rate, were derived by applying a patient-specific mathematical model to routine pretreatment MR imaging. The effect of varying degrees of extent of resection on overall survival was assessed for cohorts of patients grouped by tumor invasiveness. Results We demonstrate that patients with more diffuse tumors showed no survival benefit (P = 0.532) from gross total resection over subtotal/biopsy, while those with nodular (less diffuse) tumors showed a significant benefit (P = 0.00142) with a striking median survival benefit of over eight months compared to sub-totally resected tumors in the same cohort (an 80% improvement in survival time for GTR only seen for nodular tumors). Conclusions These results suggest that our patient-specific, model-based estimates of tumor invasiveness have clinical utility in surgical decision making. Quantification of relative invasiveness assessed from routinely obtained pre-operative imaging provides a practical predictor of the benefit of gross total resection

MultiCellDS: a standard and a community for sharing multicellular data

Cell biology is increasingly focused on cellular heterogeneity and multicellular systems. To make the fullest use of experimental, clinical, and computational efforts, we need standardized data formats, community-curated "public data libraries", and tools to combine and analyze shared data. To address these needs, our multidisciplinary community created MultiCellDS (MultiCellular Data Standard): an extensible standard, a library of digital cell lines and tissue snapshots, and support software. With the help of experimentalists, clinicians, modelers, and data and library scientists, we can grow this seed into a community-owned ecosystem of shared data and tools, to the benefit of basic science, engineering, and human health