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.

MicroRNA networks in FLT3-ITD acute myeloid leukemia.

MiR-126 and miR-155 are key microRNAs (miRNAs) that regulate, respectively, hematopoietic cell quiescence and proliferation. Herein we showed that in acute myeloid leukemia (AML), the biogenesis of these two miRNAs is interconnected through a network of regulatory loops driven by the FMS-like tyrosine kinase 3-internal tandem duplication (FLT3-ITD). In fact, FLT3-ITD induces the expression of miR-155 through a noncanonical mechanism of miRNA biogenesis that implicates cytoplasmic Drosha ribonuclease III (DROSHA). In turn, miR-155 down-regulates SH2-containing inositol phosphatase 1 (SHIP1), thereby increasing phosphor-protein kinase B (AKT) that in turn serine-phosphorylates, stabilizes, and activates Sprouty related EVH1 domain containing 1 (SPRED1). Activated SPRED1 inhibits the RAN/XPO5 complex and blocks the nucleus-to-cytoplasm transport of pre-miR-126, which cannot then complete the last steps of biogenesis. The net result is aberrantly low levels of mature miR-126 that allow quiescent leukemia blasts to be recruited into the cell cycle and proliferate. Thus, miR-126 down-regulation in proliferating AML blasts is downstream of FLT3-ITD­dependent miR-155 expression that initiates a complex circuit of concatenated regulatory feedback (i.e., miR-126/SPRED1, miR-155/human dead-box protein 3 [DDX3X]) and feed-forward (i.e., miR-155/SHIP1/AKT/miR-126) regulatory loops that eventually converge into an output signal for leukemic growth.

Sequential CAR T cell and targeted alpha immunotherapy in disseminated multiple myeloma.

Multiple myeloma (MM) is still an incurable disorder despite improved antibody and cellular therapies against different MM antigens. Single targeted antigens have so far been ineffective against MM with most patients relapsing after initial response. Hence, sequential immunotherapies directed at different targets are expected to perform better than monotherapy alone. Here, we optimized and established in preclinical studies the therapeutic rationale of using targeted alpha therapy (TAT) directed against CD38 antigen (225Ac-DOTA-daratumumab) with CAR T cell therapy directed at CS1 antigen in a systemic MM model. The sequential therapies compared CAR T therapy followed by TAT to TAT followed by CAR T therapy. CAR T cell monotherapy increased median survival from 49 days (d) in untreated controls to 71d with a modest improvement to 89d for 3.7 kBq of TAT given 14d later. When CAR T was followed by 7.4 kBq of TAT 29d later, sequential therapy increased median survival from 47d in untreated controls to 106d, compared to 68d for CAR T monotherapy. When CAR T therapy was followed by untargeted alpha immunotherapy using 7.4 kBq of 225Ac-DOTA-trastuzumab (anti-HER2) antibody 29d later, there was only a slight improvement in response over CAR T monotherapy demonstrating the role of tumor targeting. TAT (7.4 kBq) followed by CAR T therapy was also effective when CAR T therapy was delayed for 21d vs 14d or 28d post TAT, highlighting the importance of timing sequential therapies. Sequential targeted therapies using CS1 CAR T or 225Ac-DOTA-CD38 TAT in either order shows promise over monotherapies alone.

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.

Transcriptome free energy can serve as a dynamic patient-specific biomarker in acute myeloid leukemia.

Acute myeloid leukemia (AML) is prevalent in both adult and pediatric patients. Despite advances in patient categorization, the heterogeneity of AML remains a challenge. Recent studies have explored the use of gene expression data to enhance AML diagnosis and prognosis, however, alternative approaches rooted in physics and chemistry may provide another level of insight into AML transformation. Utilizing publicly available databases, we analyze 884 human and mouse blood and bone marrow samples. We employ a personalized medicine strategy, combining state-transition theory and surprisal analysis, to assess the RNA transcriptome of individual patients. The transcriptome is transformed into physical parameters that represent each sample's steady state and the free energy change (FEC) from that steady state, which is the state with the lowest free energy.We found the transcriptome steady state was invariant across normal and AML samples. FEC, representing active molecular processes, varied significantly between samples and was used to create patient-specific barcodes to characterize the biology of the disease. We discovered that AML samples that were in a transition state had the highest FEC. This disease state may be characterized as the most unstable and hence the most therapeutically targetable since a change in free energy is a thermodynamic requirement for disease progression. We also found that distinct sets of ongoing processes may be at the root of otherwise similar clinical phenotypes, implying that our integrated analysis of transcriptome profiles may facilitate a personalized medicine approach to cure AML and restore a steady state in each patient.

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.

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.

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.

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.