Abdelmoula, Walid M.

Abstract Motivation Mass spectrometry imaging (MSI) provides rich biochemical information in a label-free manner and therefore holds promise to substantially impact current practice in disease diagnosis. However, the complex nature of MSI data poses computational challenges in its analysis. The complexity of the data arises from its large size, high-dimensionality and spectral nonlinearity. Preprocessing, including peak picking, has been used to reduce raw data complexity; however, peak picking is sensitive to parameter selection that, perhaps prematurely, shapes the downstream analysis for tissue classification and ensuing biological interpretation. Results We propose a deep learning model, massNet, that provides the desired qualities of scalability, nonlinearity and speed in MSI data analysis. This deep learning model was used, without prior preprocessing and peak picking, to classify MSI data from a mouse brain harboring a patient-derived tumor. The massNet architecture established automatically learning of predictive features, and automated methods were incorporated to identify peaks with potential for tumor delineation. The model’s performance was assessed using cross-validation, and the results demonstrate higher accuracy and a substantial gain in speed compared to the established classical machine learning method, support vector machine. Availability and implementation https://github.com/wabdelmoula/massNet. The data underlying this article are available in the NIH Common Fund’s National Metabolomics Data Repository (NMDR) Metabolomics Workbench under project id (PR001292) with http://dx.doi.org/10.21228/M8Q70T. Supplementary information Supplementary data are available at Bioinformatics online.

Abstract INTRODUCTION Glioblastoma (GBM) is a diffusely invasive primary brain tumor with significant spread of tumor cells to the periphery of visible image abnormality. Enhancement of Gadolinium (Gd) contrast agent on magnetic resonance imaging (MRI) has historically been considered a confirmation of local breakdown of the blood brain barrier (BBB) and sufficient drug delivery to the bulk of tumors. In this work, we used GBM-derived xenograft (PDX) models to compare drug delivery in GBM brain for high and low BBB-permeable drugs. MATERIALS AND METHODS Five patient-derived orthotopic xenograft models from two GBM cell lines (GBM39 and GBM12) were co-dosed with erlotinib and osimertinib, two drugs with low and high BBB-permeability, respectively. T1Gd and T2-weighted MRIs were acquired from all animals prior to model sacrifice. Tumors were manually segmented on denoised and standardized MRIs and intensity patterns were captured using first and second order statistical features in the moving 3x3 kernel. We compared drug levels found in Matrix Assisted Laser Desorption Ionization (MALDI) in T1Gd enhancement, T2 enhancement, and normal brain. We also performed linear regression modeling to predict drug levels using MRI features. Model performance was measured using root mean squared error (RMSE). RESULTS Our analysis showed correlations between imaging features and MALDI drug levels. Osimertinib had a uniform distribution across the brain for all animals and all cell lines, consistent with our expectation for a high BBB-penetrant drug. Erlotinib showed the highest drug levels in T2 for GBM39 and in T1Gd for GBM12. Regression models showed promising results for predicting Erlotinib with a low RMSE of 0.037. CONCLUSION Our preliminary results suggest MRI can be predictive of drug levels for low-BBB penetrant drugs. Understanding the relationship between MRIs and drug distribution in diffuse tumors can be beneficial to developing effective treatment.

Abstract Background Response to targeted therapy varies between patients for largely unknown reasons. Here, we developed and applied an integrative platform using mass spectrometry imaging (MSI), phosphoproteomics, and multiplexed tissue imaging for mapping drug distribution, target engagement, and adaptive response to gain insights into heterogeneous response to therapy. Methods Patient-derived xenograft (PDX) lines of glioblastoma were treated with adavosertib, a Wee1 inhibitor, and tissue drug distribution was measured with MALDI-MSI. Phosphoproteomics was measured in the same tumors to identify biomarkers of drug target engagement and cellular adaptive response. Multiplexed tissue imaging was performed on sister sections to evaluate spatial co-localization of drug and cellular response. The integrated platform was then applied on clinical specimens from glioblastoma patients enrolled in the phase 1 clinical trial. Results PDX tumors exposed to different doses of adavosertib revealed intra- and inter-tumoral heterogeneity of drug distribution and integration of the heterogeneous drug distribution with phosphoproteomics and multiplexed tissue imaging revealed new markers of molecular response to adavosertib. Analysis of paired clinical specimens from patients enrolled in the phase 1 clinical trial informed the translational potential of the identified biomarkers in studying patient’s response to adavosertib. Conclusions The multimodal platform identified a signature of drug efficacy and patient-specific adaptive responses applicable to preclinical and clinical drug development. The information generated by the approach may inform mechanisms of success and failure in future early phase clinical trials, providing information for optimizing clinical trial design and guiding future application into clinical practice.

Abstract Glioblastoma (GBM) is increasingly recognized as a disease involving dysfunctional cellular metabolism. GBMs are known to be complex heterogeneous systems containing multiple distinct cell populations and are supported by an aberrant network of blood vessels. A better understanding of GBM metabolism, its variation with respect to the tumor microenvironment, and resulting regional changes in chemical composition is required. This may shed light on the observed heterogeneous drug distribution, which cannot be fully described by limited or uneven disruption of the blood–brain barrier. In this work, we used mass spectrometry imaging (MSI) to map metabolites and lipids in patient-derived xenograft models of GBM. A data analysis workflow revealed that distinctive spectral signatures were detected from different regions of the intracranial tumor model. A series of long-chain acylcarnitines were identified and detected with increased intensity at the tumor edge. A 3D MSI dataset demonstrated that these molecules were observed throughout the entire tumor/normal interface and were not confined to a single plane. mRNA sequencing demonstrated that hallmark genes related to fatty acid metabolism were highly expressed in samples with higher acylcarnitine content. These data suggest that cells in the core and the edge of the tumor undergo different fatty acid metabolism, resulting in different chemical environments within the tumor. This may influence drug distribution through changes in tissue drug affinity or transport and constitute an important consideration for therapeutic strategies in the treatment of GBM. Significance: GBM tumors exhibit a metabolic gradient that should be taken into consideration when designing therapeutic strategies for treatment. See related commentary by Tan and Weljie, p. 1231