Search

SignificanceMathematical modeling yields general principles for optimization of TRT in mouse models of multiple myeloma that can be extrapolated to other cancer models and clinical settings.

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.

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.

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 .

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.

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.

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.

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.

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.

ObjectiveAmyotrophic lateral sclerosis (ALS) is a heterogeneous disease with a complex etiology that lacks biomarkers predicting disease progression. The objective of this study was to use longitudinal cerebrospinal fluid (CSF) samples to identify biomarkers that distinguish fast progression (FP) from slow progression (SP) and assess their temporal response.MethodsWe utilized mass spectrometry (MS)-based proteomics to identify candidate biomarkers using longitudinal CSF from a discovery cohort of SP and FP ALS patients. Immunoassays were used to quantify and validate levels of the top biomarkers. A state-transition mathematical model was created using the longitudinal MS data that also predicted FP versus SP.ResultsWe identified a total of 1148 proteins in the CSF of all ALS patients. Pathway analysis determined enrichment of pathways related to complement and coagulation cascades in FPs and synaptogenesis and glucose metabolism in SPs. Longitudinal analysis revealed a panel of 59 candidate markers that could segregate FP and SP ALS. Based on multivariate analysis, we identified three biomarkers (F12, RBP4, and SERPINA4) as top candidates that segregate ALS based on rate of disease progression. These proteins were validated in the discovery and a separate validation cohort. Our state-transition model determined that the overall variance of the proteome over time was predictive of the disease progression rate.InterpretationWe identified pathways and protein biomarkers that distinguish rate of ALS disease progression. A mathematical model of the CSF proteome determined that the change in entropy of the proteome over time was predictive of FP versus SP.

Chimeric antigen receptor (CAR) T-cell based immunotherapy has shown its potential in treating blood cancers, and its application to solid tumors is currently being extensively investigated. For glioma brain tumors, various CAR T-cell targets include IL13Rα2, EGFRvIII, HER2, EphA2, GD2, B7-H3, and chlorotoxin. In this work, we are interested in developing a mathematical model of IL13Rα2 targeting CAR T-cells for treating glioma. We focus on extending the work of Kuznetsov et al. (1994) by considering binding of multiple CAR T-cells to a single glioma cell, and the dynamics of these multi-cellular conjugates. Our model more accurately describes experimentally observed CAR T-cell killing assay data than the models which do not consider multi-cellular conjugates. Moreover, we derive conditions in the CAR T-cell expansion rate that determines treatment success or failure. Finally, we show that our model captures distinct CAR T-cell killing dynamics from low to high antigen receptor densities in patient-derived brain tumor cells.

Efficacious stem cell-based therapies for traumatic brain injury (TBI) depend on successful delivery, migration, and engraftment of stem cells to induce neuroprotection. L-myc expressing human neural stem cells (LMNSC008) demonstrate an inherent tropism to injury sites after intranasal (IN) administration. We hypothesize that IN delivered LMNSC008 cells migrate to primary and secondary injury sites and modulate biomarkers associated with neuroprotection and tissue regeneration. To test this hypothesis, immunocompetent adult female rats received either controlled cortical impact injury or sham surgery. LMNSC008 cells or a vehicle were administered IN on postoperative days 7, 9, 11, 13, 15, and 17. The distribution and migration of eGFP-expressing LMNSC008 cells were quantified over 1 mm-thick optically cleared (CLARITY) coronal brain sections from TBI and SHAM controls. NSC migration was observed along white matter tracts projecting toward the hippocampus and regions of TBI. ELISA and Nanostring assays revealed a shift in tissue gene expression in LMNSC008 treated rats relative to controls. LMNSC008 treatment reduced expression of genes and pathways involved in inflammatory response, microglial function, and various cytokines and receptors. Our proof-of-concept studies, although preliminary, support the rationale of using intranasal delivery of LMNSC008 cells for functional studies in preclinical models of TBI and provide support for potential translatability in TBI patients.

Non-small cell lung cancer (NSCLC) has a high rate of brain metastasis. The purpose of this study was to assess the differential distribution of brain metastases from primary NSCLC based on mutation status. Brain MRI scans of patients with brain metastases from primary NSCLC were retrospectively analyzed. Brain metastatic tumors were grouped according to mutation status of their primary NSCLC and the neuroimaging features of these brain metastases were analyzed. A total of 110 patients with 1386 brain metastases from primary NSCLC were included in this study. Gray matter density at the tumor center peaked at ~0.6 for all mutations. The median depths of tumors were 7.9 mm, 8.7 mm and 9.1 mm for EGFR, ALK and KRAS mutation groups, respectively (p = 0.044). Brain metastases for the EGFR mutation-positive group were more frequently located in the left cerebellum, left cuneus, left precuneus and right precentral gyrus. In the ALK mutation-positive group, brain metastases were more frequently located in the right middle occipital gyrus, right posterior cingulate, right precuneus, right precentral gyrus and right parietal lobe. In the KRAS mutation-positive patient group, brain metastases were more frequently located in the posterior left cerebellum. Our study showed differential spatial distribution of brain metastases in patients with NSCLC according to their mutation status. Information regarding distribution of brain metastases is clinically relevant as it could be helpful to guide treatment planning for targeted therapy, and for predicting prognosis.

In the development of cell-based cancer therapies, quantitative mathematical models of cellular interactions are instrumental in understanding treatment efficacy. Efforts to validate and interpret mathematical models of cancer cell growth and death hinge first on proposing a precise mathematical model, then analyzing experimental data in the context of the chosen model. In this work, we present the first application of the sparse identification of non-linear dynamics (SINDy) algorithm to a real biological system in order discover cell-cell interaction dynamics in in vitro experimental data, using chimeric antigen receptor (CAR) T-cells and patient-derived glioblastoma cells. By combining the techniques of latent variable analysis and SINDy, we infer key aspects of the interaction dynamics of CAR T-cell populations and cancer. Importantly, we show how the model terms can be interpreted biologically in relation to different CAR T-cell functional responses, single or double CAR T-cell-cancer cell binding models, and density-dependent growth dynamics in either of the CAR T-cell or cancer cell populations. We show how this data-driven model-discovery based approach provides unique insight into CAR T-cell dynamics when compared to an established model-first approach. These results demonstrate the potential for SINDy to improve the implementation and efficacy of CAR T-cell therapy in the clinic through an improved understanding of CAR T-cell dynamics.

High-grade (WHO grades III-IV) glioma remains one of the most lethal human cancers. Adoptive transfer of tumor-targeting chimeric antigen receptor (CAR)-redirected T cells for high-grade glioma has revealed promising indications of anti-tumor activity, but objective clinical responses remain elusive for most patients. A significant challenge to effective immunotherapy is the highly heterogeneous structure of these tumors, including large variations in the magnitudes and distributions of target antigen expression, observed both within individual tumors and between patients. To obtain a more detailed understanding of immunotherapy target antigens within patient tumors, we immunochemically mapped at single cell resolution three clinically-relevant targets, IL13Rα2, HER2 and EGFR, on tumor samples drawn from a 43-patient cohort. We observed that within individual tumor samples, expression of these antigens was neither random nor uniform, but rather that they mapped into local neighborhoods - phenotypically similar cells within regions of cellular tumor - reflecting not well understood properties of tumor cells and their milieu. Notably, tumor cell neighborhoods of high antigen expression were not arranged independently within regions. For example, in cellular tumor regions, neighborhoods of high IL13Rα2 and HER2 expression appeared to be reciprocal to those of EGFR, while in areas of pseudopalisading necrosis, expression of IL13Rα2 and HER2, but not EGFR, appeared to reflect the radial organization of tumor cells around hypoxic cores. Other structural features affecting expression of immunotherapy target antigens remain to be elucidated. This structured but heterogeneous organization of antigen expression in high grade glioma is highly permissive for antigen escape, and combinatorial antigen targeting is a commonly suggested potential mitigating strategy. Deeper understanding of antigen expression within and between patient tumors will enhance optimization of combination immunotherapies, the most immediate clinical application of the observations presented here being the importance of including (wild-type) EGFR as a target antigen.

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?

MicroRNAs (miRNAs) have been shown to hold prognostic value in acute myeloid leukemia (AML); however, the temporal dynamics of miRNA expression in AML are poorly understood. Using serial samples from a mouse model of AML to generate time-series miRNA sequencing data, we are the first to show that the miRNA transcriptome undergoes state-transition during AML initiation and progression. We modeled AML state-transition as a particle undergoing Brownian motion in a quasi-potential and validated the AML state-space and state-transition model to accurately predict time to AML in an independent cohort of mice. The critical points of the model provided a framework to align samples from mice that developed AML at different rates. Our mathematical approach allowed discovery of dynamic processes involved during AML development and, if translated to humans, has the potential to predict an individual's disease trajectory.

Neural stem cells (NSCs) offer a potential solution to treating brain tumors. This is because NSCs can circumvent the blood-brain barrier and migrate to areas of damage in the central nervous system, including tumors, stroke, and wound injuries. However, for successful clinical application of NSC treatment, a sufficient number of viable cells must reach the diseased or damaged area(s) in the brain, and evidence suggests that it may be affected by the paths the NSCs take through the brain, as well as the locations of tumors. To study the NSC migration in brain, we develop a mathematical model of therapeutic NSC migration towards brain tumor, that provides a low cost platform to investigate NSC treatment efficacy. Our model is an extension of the model developed in Rockne et al. (PLoS ONE 13, e0199967, 2018) that considers NSC migration in non-tumor bearing naive mouse brain. Here we modify the model in Rockne et al. in three ways: (i) we consider three-dimensional mouse brain geometry, (ii) we add chemotaxis to model the tumor-tropic nature of NSCs into tumor sites, and (iii) we model stochasticity of migration speed and chemosensitivity. The proposed model is used to study migration patterns of NSCs to sites of tumors for different injection strategies, in particular, intranasal and intracerebral delivery. We observe that intracerebral injection results in more NSCs arriving at the tumor site(s), but the relative fraction of NSCs depends on the location of injection relative to the target site(s). On the other hand, intranasal injection results in fewer NSCs at the tumor site, but yields a more even distribution of NSCs within and around the target tumor site(s).

Single-cell sequencing technologies have revolutionized molecular and cellular biology and stimulated the development of computational tools to analyze the data generated from these technology platforms. However, despite the recent explosion of computational analysis tools, relatively few mathematical models have been developed to utilize these data. Here we compare and contrast two cell state geometries for building mathematical models of cell state-transitions with single-cell RNA-sequencing data with hematopoeisis as a model system; (i) by using partial differential equations on a graph representing intermediate cell states between known cell types, and (ii) by using the equations on a multi-dimensional continuous cell state-space. As an application of our approach, we demonstrate how the calibrated models may be used to mathematically perturb normal hematopoeisis to simulate, predict, and study the emergence of novel cell states during the pathogenesis of acute myeloid leukemia. We particularly focus on comparing the strength and weakness of the graph model and multi-dimensional model.

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.

Cell fate commitment is driven by dynamic changes in chromatin architecture and activity of lineage-specific transcription factors (TFs). The chromatin assembly factor-1 (CAF-1) is a histone chaperone that regulates chromatin architecture by facilitating nucleosome assembly during DNA replication. Accumulating evidence supports a substantial role of CAF-1 in cell fate maintenance, but the mechanisms by which CAF-1 restricts lineage choice remain poorly understood. Here, we investigate how CAF-1 influences chromatin dynamics and TF activity during lineage differentiation. We show that CAF-1 suppression triggers rapid differentiation of myeloid stem and progenitor cells into a mixed lineage state. We find that CAF-1 sustains lineage fidelity by controlling chromatin accessibility at specific loci, and limiting the binding of ELF1 TF at newly-accessible diverging regulatory elements. Together, our findings decipher key traits of chromatin accessibility that sustain lineage integrity and point to a powerful strategy for dissecting transcriptional circuits central to cell fate commitment.

BackgroundDuring acute myeloid leukemia (AML) growth, the bone marrow (BM) niche acquires significant vascular changes that can be offset by therapeutic blast cytoreduction. The molecular mechanisms of this vascular plasticity remain to be fully elucidated. Herein, we report on the changes that occur in the vascular compartment of the FLT3-ITD+ AML BM niche pre and post treatment and their impact on leukemic stem cells (LSCs).MethodsBM vasculature was evaluated in FLT3-ITD+ AML models (MllPTD/WT/Flt3ITD/ITD mouse and patient-derived xenograft) by 3D confocal imaging of long bones, calvarium vascular permeability assays, and flow cytometry analysis. Cytokine levels were measured by Luminex assay and miR-126 levels evaluated by Q-RT-PCR and miRNA staining. Wild-type (wt) and MllPTD/WT/Flt3ITD/ITD mice with endothelial cell (EC) miR-126 knockout or overexpression served as controls. The impact of treatment-induced BM vascular changes on LSC activity was evaluated by secondary transplantation of BM cells after administration of tyrosine kinase inhibitors (TKIs) to MllPTD/WT/Flt3ITD/ITD mice with/without either EC miR-126 KO or co-treatment with tumor necrosis factor alpha (TNFα) or anti-miR-126 miRisten.ResultsIn the normal BM niche, CD31+Sca-1high ECs lining arterioles have miR-126 levels higher than CD31+Sca-1low ECs lining sinusoids. We noted that during FLT3-ITD+ AML growth, the BM niche lost arterioles and gained sinusoids. These changes were mediated by TNFα, a cytokine produced by AML blasts, which induced EC miR-126 downregulation and caused depletion of CD31+Sca-1high ECs and gain in CD31+Sca-1low ECs. Loss of miR-126high ECs led to a decreased EC miR-126 supply to LSCs, which then entered the cell cycle and promoted leukemia growth. Accordingly, antileukemic treatment with TKI decreased the BM blast-produced TNFα and increased miR-126high ECs and the EC miR-126 supply to LSCs. High miR-126 levels safeguarded LSCs, as shown by more severe disease in secondary transplanted mice. Conversely, EC miR-126 deprivation via genetic or pharmacological EC miR-126 knock-down prevented treatment-induced BM miR-126high EC expansion and in turn LSC protection.ConclusionsTreatment-induced CD31+Sca-1high EC re-vascularization of the leukemic BM niche may represent a LSC extrinsic mechanism of treatment resistance that can be overcome with therapeutic EC miR-126 deprivation.

We report here on a novel pro-leukemogenic role of FMS-like tyrosine kinase 3-internal tandem duplication (FLT3-ITD) that interferes with microRNAs (miRNAs) biogenesis in acute myeloid leukemia (AML) blasts. We showed that FLT3-ITD interferes with the canonical biogenesis of intron-hosted miRNAs such as miR-126, by phosphorylating SPRED1 protein and inhibiting the "gatekeeper" Exportin 5 (XPO5)/RAN-GTP complex that regulates the nucleus-to-cytoplasm transport of pre-miRNAs for completion of maturation into mature miRNAs. Of note, despite the blockage of "canonical" miRNA biogenesis, miR-155 remains upregulated in FLT3-ITD+ AML blasts, suggesting activation of alternative mechanisms of miRNA biogenesis that circumvent the XPO5/RAN-GTP blockage. MiR-155, a BIC-155 long noncoding (lnc) RNA-hosted oncogenic miRNA, has previously been implicated in FLT3-ITD+ AML blast hyperproliferation. We showed that FLT3-ITD upregulates miR-155 by inhibiting DDX3X, a protein implicated in the splicing of lncRNAs, via p-AKT. Inhibition of DDX3X increases unspliced BIC-155 that is then shuttled by NXF1 from the nucleus to the cytoplasm, where it is processed into mature miR-155 by cytoplasmic DROSHA, thereby bypassing the XPO5/RAN-GTP blockage via "non-canonical" mechanisms of miRNA biogenesis.

As the success of stem cell-based therapies is contingent on efficient cell delivery to damaged areas, neural stem cells (NSCs) have promising therapeutic potential because they inherently migrate to sites of central nervous system (CNS) damage. To explore the possibility of NSC-based therapy after traumatic brain injury (TBI), isoflurane-anesthetized adult male rats received a controlled cortical impact (CCI) of moderate severity (2.8 mm deformation at 4 m/s) or sham injury (i.e., no cortical impact). Beginning 1-week post-injury, the rats were immunosuppressed and 1 × 106 human NSCs (LM-NS008.GFP.fLuc) or vehicle (VEH) (2% human serum albumen) were administered intranasally (IN) on post-operative days 7, 9, 11, 13, 15, and 17. To evaluate the spatial distributions of the LM-NSC008 cells, half of the rats were euthanized on day 25, one day after completion of the cognitive task, and the other half were euthanized on day 46. 1 mm thick brain sections were optically cleared (CLARITY), and volumes were imaged by confocal microscopy. In addition, LM-NSC008 cell migration to the TBI site by immunohistochemistry for human-specific Nestin was observed at day 39. Acquisition of spatial learning was assessed in a well-established Morris water maze task on six successive days beginning on post-injury day 18. IN administration of LM-NSC008 cells after TBI (TBI + NSC) significantly facilitated spatial learning relative to TBI + VEH rats (p < 0.05) and had no effect on sham + NSC rats. Overall, these data indicate that IN-administered LM-NSC008 cells migrate to sites of TBI damage and that their presence correlates with cognitive improvement. Future studies will expand on these preliminary findings by evaluating other LM-NSC008 cell dosing paradigms and evaluating mechanisms by which LM-NSC008 cells contribute to cognitive recovery.

Recent successes of immune-modulating therapies for cancer have stimulated research on information flow within the immune system and, in turn, clinical applications of concepts from information theory. Through information theory, one can describe and formalize, in a mathematically rigorous fashion, the function of interconnected components of the immune system in health and disease. Specifically, using concepts including entropy, mutual information, and channel capacity, one can quantify the storage, transmission, encoding, and flow of information within and between cellular components of the immune system on multiple temporal and spatial scales. To understand, at the quantitative level, immune signaling function and dysfunction in cancer, we present a methodology-oriented review of information-theoretic treatment of biochemical signal transduction and transmission coupled with mathematical modeling.

Over 90% of breast cancer is cured; yet there remain highly aggressive breast cancers that develop rapidly and are extremely difficult to treat, much less prevent. Breast cancers that rapidly develop between breast image screening are called "interval cancers." The efforts of our team focus on identifying multiscale integrated strategies to identify biologically aggressive precancerous breast lesions. Our goal is to identify spatiotemporal changes that occur prior to development of interval breast cancers. To accomplish this requires integration of new technology. Our team has the ability to perform single cell in situ transcriptional profiling, noncontrast biological imaging, mathematical analysis, and nanoscale evaluation of receptor organization and signaling. These technological innovations allow us to start to identify multidimensional spatial and temporal relationships that drive the transition from biologically aggressive precancer to biologically aggressive interval breast cancer. This article is categorized under: Cancer > Computational Models Cancer > Molecular and Cellular Physiology Cancer > Genetics/Genomics/Epigenetics.

Acute myeloid leukemia (AML) harboring inv(16)(p13q22) expresses high levels of miR-126. Here we show that the CBFB-MYH11 (CM) fusion gene upregulates miR-126 expression through aberrant miR-126 transcription and perturbed miR-126 biogenesis via the HDAC8/RAN-XPO5-RCC1 axis. Aberrant miR-126 upregulation promotes survival of leukemia-initiating progenitors and is critical for initiating and maintaining CM-driven AML. We show that miR-126 enhances MYC activity through the SPRED1/PLK2-ERK-MYC axis. Notably, genetic deletion of miR-126 significantly reduces AML rate and extends survival in CM knock-in mice. Therapeutic depletion of miR-126 with an anti-miR-126 (miRisten) inhibits AML cell survival, reduces leukemia burden and leukemia stem cell (LSC) activity in inv(16) AML murine and xenograft models. The combination of miRisten with chemotherapy further enhances the anti-leukemia and anti-LSC activity. Overall, this study provides molecular insights for the mechanism and impact of miR-126 dysregulation in leukemogenesis and highlights the potential of miR-126 depletion as a therapeutic approach for inv(16) AML.

BackgroundDynamic contrast-enhanced MRI (DCE-MRI) parameters have been shown to be biomarkers for treatment response in glioblastoma (GBM). However, variations in analysis and measurement methodology complicate determination of biological changes measured via DCE. The aim of this study is to quantify DCE-MRI variations attributable to analysis methodology and image quality in GBM patients.MethodsThe Extended Tofts model (eTM) and Leaky Tracer Kinetic Model (LTKM), with manually and automatically segmented vascular input functions (VIFs), were used to calculate perfusion kinetic parameters from 29 GBM patients with double-baseline DCE-MRI data. DCE-MRI images were acquired 2-5 days apart with no change in treatment. Repeatability of kinetic parameters was quantified with Bland-Altman and percent repeatability coefficient (%RC) analysis.ResultsThe perfusion parameter with the least RC was the plasma volume fraction (v p ), with a %RC of 53%. The extra-cellular extra-vascular volume fraction (v e ) %RC was 82% and 81%, for extended Tofts-Kety Model (eTM) and LTKM respectively. The %RC of the volume transfer rate constant (K trans ) was 72% for the eTM, and 82% for the LTKM, respectively. Using an automatic VIF resulted in smaller %RCs for all model parameters, as compared to manual VIF.ConclusionsAs much as 72% change in K trans (eTM, autoVIF) can be attributable to non-biological changes in the 2-5 days between double-baseline imaging. Poor K trans repeatability may result from inferior temporal resolution and short image acquisition time. This variation suggests DCE-MRI repeatability studies should be performed institutionally, using an automatic VIF method and following quantitative imaging biomarkers alliance guidelines.

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.

BackgroundGlioblastoma (GBM) is the deadliest and most common brain tumor in adults, with poor survival and response to aggressive therapy. Limited access of drugs to tumor cells is one reason for such grim clinical outcomes. A driving force for therapeutic delivery is interstitial fluid flow (IFF), both within the tumor and in the surrounding brain parenchyma. However, convective and diffusive transport mechanisms are understudied. In this study, we examined the application of a novel image analysis method to measure fluid flow and diffusion in GBM patients.MethodsHere, we applied an imaging methodology that had been previously tested and validated in vitro, in silico, and in preclinical models of disease to archival patient data from the Ivy Glioblastoma Atlas Project (GAP) dataset. The analysis required the use of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI), which is readily available in the database. The analysis results, which consisted of IFF flow velocity and diffusion coefficients, were then compared to patient outcomes such as survival.ResultsWe characterized IFF and diffusion patterns in patients. We found strong correlations between flow rates measured within tumors and in the surrounding parenchymal space, where we hypothesized that velocities would be higher. Analyzing overall magnitudes indicated a significant correlation with both age and survival in this patient cohort. Additionally, we found that neither tumor size nor resection significantly altered the velocity magnitude. Lastly, we mapped the flow pathways in patient tumors and found a variability in the degree of directionality that we hypothesize may lead to information concerning treatment, invasive spread, and progression in future studies.ConclusionsAn analysis of standard DCE-MRI in patients with GBM offers more information regarding IFF and transport within and around the tumor, shows that IFF is still detected post-resection, and indicates that velocity magnitudes correlate with patient prognosis.

Radioimmunotherapy, an approach using radiolabeled antibodies, has had minimal success in the clinic with several β-emitting radionuclides for the treatment of ovarian cancer. Alternatively, radioimmunotherapy with α-emitters offers the advantage of depositing much higher energy over shorter distances but was thought to be inappropriate for the treatment of solid tumors, for which antibody penetration is limited to a few cell diameters around the vascular system. However, the deposition of high-energy α-emitters to tumor markers adjacent to a typical leaky tumor vascular system may have large antitumor effects at the tumor vascular level, and their reduced penetration in normal tissue would be expected to lower off-target toxicity. Methods: To evaluate this concept, DOTAylated-huCC49 was labeled with the α-emitter 225Ac to target tumor-associated glycoprotein 72-positive xenografts in a murine model of ovarian cancer. Results:225Ac-labeled DOTAylated-huCC49 radioimmunotherapy significantly reduced tumor growth in a dose-dependent manner (1.85, 3.7, and 7.4 kBq), with the 7.4-kBq dose extending survival by more than 3-fold compared with the untreated control. Additionally, a multitreatment regime (1.85 kBq followed by 5 weekly doses of 0.70 kBq for a total of 5.4 kBq) extended survival almost 3-fold compared with the untreated control group, without significant off-target toxicity. Conclusion: These results establish the potential for antibody-targeted α-radionuclide therapy for ovarian cancer, which may be generalized to α-radioimmunotherapy in other solid tumors.

The homeostasis of vascular endothelium is crucial for cardiovascular health and endothelial cell (EC) aging and dysfunction could negatively impact vascular function. Leveraging transcriptome profiles from ECs subjected to various stimuli, including time-series data obtained from ECs under physiological pulsatile flow vs. pathophysiological oscillatory flow, we performed principal component analysis (PCA) to identify key genes contributing to divergent transcriptional states of ECs. Through bioinformatics analysis, we identified that a long non-coding RNA (lncRNA) RAMP2-AS1 encoded on the antisense of RAMP2, a determinant of endothelial homeostasis and vascular integrity, is a novel regulator essential for EC homeostasis and function. Knockdown of RAMP2-AS1 suppressed RAMP2 expression and caused EC functional changes promoting aging, including impaired angiogenesis and increased senescence. Our study demonstrates an integrative approach to quantifying EC aging based on transcriptome changes, which also identified a number of novel regulators, including protein-coding genes and many lncRNAs involved EC functional modulation, exemplified by RAMP2-AS1.

Background: Brain metastases are associated with poor survival. Molecular genetic testing informs on targeted therapy and survival. The purpose of this study was to perform a MR imaging-based radiomic analysis of brain metastases from non-small cell lung cancer (NSCLC) to identify radiomic features that were important for predicting survival duration. Methods: We retrospectively identified our study cohort via an institutional database search for patients with brain metastases from EGFR, ALK, and/or KRAS mutation-positive NSCLC. We segmented the brain metastatic tumors on the brain MR images, extracted radiomic features, constructed radiomic scores from significant radiomic features based on multivariate Cox regression analysis (p < 0.05), and built predictive models for survival duration. Result: Of the 110 patients in the cohort (mean age 57.51 ± 12.32 years; range: 22-85 years, M:F = 37:73), 75, 26, and 15 had NSCLC with EGFR, ALK, and KRAS mutations, respectively. Predictive modeling of survival duration using both clinical and radiomic features yielded areas under the receiver operative characteristic curve of 0.977, 0.905, and 0.947 for the EGFR, ALK, and KRAS mutation-positive groups, respectively. Radiomic scores enabled the separation of each mutation-positive group into two subgroups with significantly different survival durations, i.e., shorter vs. longer duration when comparing to the median survival duration of the group. Conclusion: Our data supports the use of radiomic scores, based on MR imaging of brain metastases from NSCLC, as non-invasive biomarkers for survival duration. Future research with a larger sample size and external cohorts is needed to validate our results.

Temporal dynamics of gene expression inform cellular and molecular perturbations associated with disease development and evolution. Given the complexity of high-dimensional temporal genomic data, an analytic framework guided by a robust theory is needed to interpret time-sequential changes and to predict system dynamics. Here we model temporal dynamics of the transcriptome of peripheral blood mononuclear cells in a two-dimensional state-space representing states of health and leukemia using time-sequential bulk RNA-seq data from a murine model of acute myeloid leukemia (AML). The state-transition model identified critical points that accurately predict AML development and identifies stepwise transcriptomic perturbations that drive leukemia progression. The geometry of the transcriptome state-space provided a biological interpretation of gene dynamics, aligned gene signals that are not synchronized in time across mice, and allowed quantification of gene and pathway contributions to leukemia development. Our state-transition model synthesizes information from multiple cell types in the peripheral blood and identifies critical points in the transition from health to leukemia to guide interpretation of changes in the transcriptome as a whole to predict disease progression. SIGNIFICANCE: These findings apply the theory of state transitions to model the initiation and development of acute myeloid leukemia, identifying transcriptomic perturbations that accurately predict time to disease development.See related commentary by Kuijjer, p. 3072 GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/15/3157/F1.large.jpg.

Lung cancer metastases comprise most of all brain metastases in adults and most brain metastases are diagnosed by magnetic resonance (MR) scans. The purpose of this study was to conduct an MR imaging-based radiomic analysis of brain metastatic lesions from patients with primary lung cancer to classify mutational status of the metastatic disease. We retrospectively identified lung cancer patients with brain metastases treated at our institution between 2009 and 2017 who underwent genotype testing of their primary lung cancer. Brain MR Images were used for segmentation of enhancing tumors and peritumoral edema, and for radiomic feature extraction. The most relevant radiomic features were identified and used with clinical data to train random forest classifiers to classify the mutation status. Of 110 patients in the study cohort (mean age 57.51 ± 12.32 years; M: F = 37:73), 75 had an EGFR mutation, 21 had an ALK translocation, and 15 had a KRAS mutation. One patient had both ALK translocation and EGFR mutation. Majority of radiomic features most relevant for mutation classification were textural. Model building using both radiomic features and clinical data yielded more accurate classifications than using either alone. For classification of EGFR, ALK, and KRAS mutation status, the model built with both radiomic features and clinical data resulted in area-under-the-curve (AUC) values based on cross-validation of 0.912, 0.915, and 0.985, respectively. Our study demonstrated that MR imaging-based radiomic analysis of brain metastases in patients with primary lung cancer may be used to classify mutation status. This approach may be useful for devising treatment strategies and informing prognosis.

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.

While targeted therapies exist for human epidermal growth factor receptor 2 positive (HER2 +) breast cancer, HER2 + patients do not always respond to therapy. We present the results of utilizing a biophysical mathematical model to predict tumor response for two HER2 + breast cancer patients treated with the same therapeutic regimen but who achieved different treatment outcomes. Quantitative data from magnetic resonance imaging (MRI) and 64Cu-DOTA-trastuzumab positron emission tomography (PET) are used to estimate tumor density, perfusion, and distribution of HER2-targeted antibodies for each individual patient. MRI and PET data are collected prior to therapy, and follow-up MRI scans are acquired at a midpoint in therapy. Given these data types, we align the data sets to a common image space to enable model calibration. Once the model is parameterized with these data, we forecast treatment response with and without HER2-targeted therapy. By incorporating targeted therapy into the model, the resulting predictions are able to distinguish between the two different patient responses, increasing the difference in tumor volume change between the two patients by > 40%. This work provides a proof-of-concept strategy for processing and integrating PET and MRI modalities into a predictive, clinical-mathematical framework to provide patient-specific predictions of HER2 + treatment response.

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 < 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.

Tumor recurrence in glioblastoma multiforme (GBM) is often attributed to acquired resistance to the standard chemotherapeutic agent, temozolomide (TMZ). Promoter methylation of the DNA repair gene MGMT (O6-methylguanine-DNA methyltransferase) has been associated with sensitivity to TMZ, whereas increased expression of MGMT has been associated with TMZ resistance. Clinical studies have observed a downward shift in MGMT methylation percentage from primary to recurrent stage tumors; however, the evolutionary processes that drive this shift and more generally the emergence and growth of TMZ-resistant tumor subpopulations are still poorly understood. Here, we develop a mathematical model, parameterized using clinical and experimental data, to investigate the role of MGMT methylation in TMZ resistance during the standard treatment regimen for GBM-surgery, chemotherapy, and radiation. We first found that the observed downward shift in MGMT promoter methylation status between detection and recurrence cannot be explained solely by evolutionary selection. Next, our model suggests that TMZ has an inhibitory effect on maintenance methylation of MGMT after cell division. Finally, incorporating this inhibitory effect, we study the optimal number of TMZ doses per adjuvant cycle for patients with GBM with high and low levels of MGMT methylation at diagnosis.

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.

Natural products have long been a viable source of therapeutic agents, providing unique structures and mechanisms that may be beneficial for cancer treatment. Herein we first report on the anticancer activity OST-01, a natural product from Baccharis Coridifolia, on breast cancer cells, including triple-negative breast cancer (TNBC). OST-01 significantly inhibited cell proliferation and oncogenic activities of TNBC cells in vitro. OST-01 also markedly inhibited TNBC tumor growth in vivo, with > 50% reduction in tumor size compared to vehicle control treatment in different in vivo models, i.e., cell line-derived (CDX), patient-derived (PDX), and mammary fat pad xenografts. Mechanistically, OST-01 induces ferroptosis by downregulating LRP8-regulated selenoproteins, i.e., GPX4. A shift from a basal-mesenchymal to a luminal-epithelial state of breast cancer stem cells (BCSCs) as supported by the downregulation of stemness (e.g., CD44) and mesenchymal (e.g., FN1 and vimentin) markers, along with the upregulation of differentiation markers (e.g., CD24) and luminal-epithelial markers (e.g., CK19), was also observed. [ABSTRACT FROM AUTHOR]

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.

Glioblastoma (GBM), a highly aggressive form of brain tumor, is a disease marked by extensive invasion into the surrounding brain. Interstitial fluid flow (IFF), or the movement of fluid within the spaces between cells, has been linked to increased invasion of GBM cells. Better characterization of IFF could elucidate underlying mechanisms driving this invasion in vivo. Here, we develop a technique to noninvasively measure interstitial flow velocities in the glioma microenvironment of mice using dynamic contrast-enhanced magnetic resonance imaging (MRI), a common clinical technique. Using our in vitro model as a phantom "tumor" system and in silico models of velocity vector fields, we show we can measure average velocities and accurately reconstruct velocity directions. With our combined MR and analysis method, we show that velocity magnitudes are similar across four human GBM cell line xenograft models and the direction of fluid flow is heterogeneous within and around the tumors, and not always in the outward direction. These values were not linked to the tumor size. Finally, we compare our flow velocity magnitudes and the direction of flow to a classical marker of vessel leakage and bulk fluid drainage, Evans blue. With these data, we validate its use as a marker of high and low IFF rates and IFF in the outward direction from the tumor border in implanted glioma models. These methods show, for the first time, the nature of interstitial fluid flow in models of glioma using a technique that is translatable to clinical and preclinical models currently using contrast-enhanced MRI.

BackgroundThe 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.MethodsA 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.ResultsAt 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.ConclusionThe 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.

Here we present a theoretical and mathematical perspective on the process of aging. We extend the concepts of physical space and time to an abstract, mathematically-defined space, which we associate with a concept of "biological space-time" in which biological dynamics may be represented. We hypothesize that biological dynamics, represented as trajectories in biological space-time, may be used to model and study different rates of biological aging. As a consequence of this hypothesis, we show how dilation or contraction of time analogous to relativistic corrections of physical time resulting from accelerated or decelerated biological dynamics may be used to study precipitous or protracted aging. We show specific examples of how these principles may be used to model different rates of aging, with an emphasis on cancer in aging. We discuss how this theory may be tested or falsified, as well as novel concepts and implications of this theory that may improve our interpretation of biological aging.