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.

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.

PURPOSE: It has been reported that diffusion-weighted imaging (DWI) with ultrahigh b-value increases the diagnostic power of prostate cancer. DWI with higher b-value increases the diagnostic power of prostate cancer. DWI with higher b-value increases the diagnostic power of prostate cancer. DWI with higher b-value increases the diagnostic power of prostate cancer. DWI with higher Materials and Methods. Fifteen patients (7 malignant and 8 benign) were included in this study retrospectively with the institutional ethical committee approval. All images were acquired at a 3T MR scanner. The ADC values were calculated using a monoexponential model. Synthetic ADC (sADC) for higher b-value increases the diagnostic power of prostate cancer. DWI with higher. RESULTS: No significant difference was observed between actual ADC and sADC for b-value increases the diagnostic power of prostate cancer. DWI with higher p=0.002, paired t-test) in sDWI as compared to DWI. Malignant lesions showed significantly lower sADC as compared to benign lesions (p=0.002, paired t-test) in sDWI as compared to DWI. Malignant lesions showed significantly lower sADC as compared to benign lesions (Discussion/. CONCLUSION: Our initial investigation suggests that the ADC values corresponding to higher b-value can be computed using log-linear relationship derived from lower b-values (b ≤ 1000). Our method might help clinicians to decide the optimal b-value for prostate lesion identification.b-value increases the diagnostic power of prostate cancer. DWI with higher b-value increases the diagnostic power of prostate cancer. DWI with higher b-value increases the diagnostic power of prostate cancer. DWI with higher b-value increases the diagnostic power of prostate cancer. DWI with higher.

Background

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

Methods

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

Results

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

Conclusions

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

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

Cell-based therapies to the brain are promising for the treatment of multiple brain disorders including neurodegeneration and cancers. In order to access the brain parenchyma, there are multiple physiological barriers that must be overcome depending on the route of delivery. Specifically, the blood-brain barrier has been a major difficulty in drug delivery for decades, and it still presents a challenge for the delivery of therapeutic cells. Other barriers, including the blood-cerebrospinal fluid barrier and lymphatic-brain barrier, are less explored, but may offer specific challenges or opportunities for therapeutic delivery. Here we discuss the barriers to the brain and the strategies currently in place to deliver cell-based therapies, including engineered T cells, dendritic cells, and stem cells, to treat diseases. With a particular focus on cancers, we also highlight the current ongoing clinical trials that use cell-based therapies to treat disease, many of which show promise at treating some of the deadliest illnesses.

Targeted therapies for multiple myeloma (MM) include the anti-CD38 antibody daratumumab, which, in addition to its inherent cytotoxicity, can be radiolabeled with tracers for imaging and with β- and α-emitter radionuclides for radioimmunotherapy. Methods: We have compared the potential therapeutic efficacy of β- versus α-emitter radioimmunotherapy using radiolabeled DOTA-daratumumab in a preclinical model of disseminated multiple myeloma. Multiple dose levels were investigated to find the dose with the highest efficacy and lowest toxicity. Results: In a dose–response study with the β-emitter 177Lu-DOTA-daratumumab, the lowest tested dose, 1.85 MBq, extended survival from 37 to 47 d but did not delay tumor growth. Doses of 3.7 and 7.4 MBq extended survival to 55 and 58 d, respectively, causing a small equivalent delay in tumor growth, followed by regrowth. The higher dose, 11.1 MBq, eradicated the tumor but had no effect on survival compared with untreated controls, because of whole-body toxicity. In contrast, the α-emitter 225Ac-DOTA-daratumumab had a dose-dependent effect, in which 0.925, 1.85, and 3.7 kBq increased survival, compared with untreated controls (35 d), to 47, 52, and 73 d, respectively, with a significant delay in tumor growth for all 3 doses. Higher doses of 11.1 and 22.2 kBq resulted in equivalent survival to 82 d but with significant whole-body toxicity. Parallel studies with untargeted 225Ac-DOTA-trastuzumab conferred no improvement over untreated controls and resulted in whole-body toxicity. Conclusion: We conclude, and mathematic modeling confirms, that maximal biologic doses were achieved by targeted α-therapy and demonstrated 225Ac to be superior to 177Lu in delaying tumor growth and decreasing whole-body toxicity.

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.

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.

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.