Search

A method is presented to select the optimal time points at which to measure DCE-MRI signal intensities, leaving time in the MR exam for high-spatial resolution image acquisition.Simplicial complexes are generated from the Kety-Tofts model pharmacokinetic parameters KThe DCE-MRI data acquired in women with invasive breast cancer (N = 27) were used to retrospectively compare parameter maps fit to full and subsampled time courses. Simplicial complexes were generated for a fixed range of Kety-Tofts model parameters and for the parameter ranges weighted by estimates from the fully sampled data. The largest-area manifolds determined the optimal three time points for each case. Simulations were performed along with retrospectively subsampled data fits. The agreement was computed between the model parameters fit to three points and those fit to all points.The optimal three-point sample times were from the data-informed simplicial complex analysis and determined to be 65, 204, and 393 s after arrival of the contrast agent to breast tissue. In the patient data, tumor-median parameter values fit using all points and the three selected time points agreed with concordance correlation coefficients of 0.97 for KIt is possible to accurately estimate pharmacokinetic parameters from three properly selected time points inserted into a clinical DCE-MRI breast exam. This technique can provide guidance on when to capture images for quantitative data between high-spatial-resolution DCE-MRI images.

Triple-negative breast cancer (TNBC) is persistently refractory to therapy, and methods to improve targeting and evaluation of responses to therapy in this disease are needed. Here, we integrate quantitative MRI data with biologically based mathematical modeling to accurately predict the response of TNBC to neoadjuvant systemic therapy (NAST) on an individual basis. Specifically, 56 patients with TNBC enrolled in the ARTEMIS trial (NCT02276443) underwent standard-of-care doxorubicin/cyclophosphamide (A/C) and then paclitaxel for NAST, where dynamic contrast-enhanced MRI and diffusion-weighted MRI were acquired before treatment and after two and four cycles of A/C. A biologically based model was established to characterize tumor cell movement, proliferation, and treatment-induced cell death. Two evaluation frameworks were investigated using: (i) images acquired before and after two cycles of A/C for calibration and predicting tumor status after A/C, and (ii) images acquired before, after two cycles, and after four cycles of A/C for calibration and predicting response following NAST. For Framework 1, the concordance correlation coefficients between the predicted and measured patient-specific, post-A/C changes in tumor cellularity and volume were 0.95 and 0.94, respectively. For Framework 2, the biologically based model achieved an area under the receiver operator characteristic curve of 0.89 (sensitivity/specificity = 0.72/0.95) for differentiating pathological complete response (pCR) from non-pCR, which is statistically superior (P < 0.05) to the value of 0.78 (sensitivity/specificity = 0.72/0.79) achieved by tumor volume measured after four cycles of A/C. Overall, this model successfully captured patient-specific, spatiotemporal dynamics of TNBC response to NAST, providing highly accurate predictions of NAST response. Significance: Integrating MRI data with biologically based mathematical modeling successfully predicts breast cancer response to chemotherapy, suggesting digital twins could facilitate a paradigm shift from simply assessing response to predicting and optimizing therapeutic efficacy.

Triple-negative breast cancer (TNBC) is persistently refractory to therapy, and methods to improve targeting and evaluation of responses to therapy in this disease are needed. Here, we integrate quantitative MRI data with biologically based mathematical modeling to accurately predict the response of TNBC to neoadjuvant systemic therapy (NAST) on an individual basis. Specifically, 56 patients with TNBC enrolled in the ARTEMIS trial (NCT02276443) underwent standard-of-care doxorubicin/cyclophosphamide (A/C) and then paclitaxel for NAST, where dynamic contrast-enhanced MRI and diffusion-weighted MRI were acquired before treatment and after two and four cycles of A/C. A biologically based model was established to characterize tumor cell movement, proliferation, and treatment-induced cell death. Two evaluation frameworks were investigated using: (i) images acquired before and after two cycles of A/C for calibration and predicting tumor status after A/C, and (ii) images acquired before, after two cycles, and after four cycles of A/C for calibration and predicting response following NAST. For Framework 1, the concordance correlation coefficients between the predicted and measured patient-specific, post-A/C changes in tumor cellularity and volume were 0.95 and 0.94, respectively. For Framework 2, the biologically based model achieved an area under the receiver operator characteristic curve of 0.89 (sensitivity/specificity = 0.72/0.95) for differentiating pathological complete response (pCR) from non-pCR, which is statistically superior (P < 0.05) to the value of 0.78 (sensitivity/specificity = 0.72/0.79) achieved by tumor volume measured after four cycles of A/C. Overall, this model successfully captured patient-specific, spatiotemporal dynamics of TNBC response to NAST, providing highly accurate predictions of NAST response.Significance:Integrating MRI data with biologically based mathematical modeling successfully predicts breast cancer response to chemotherapy, suggesting digital twins could facilitate a paradigm shift from simply assessing response to predicting and optimizing therapeutic efficacy.

Supplementary Data from MRI-Based Digital Models Forecast Patient-Specific Treatment Responses to Neoadjuvant Chemotherapy in Triple-Negative Breast Cancer

The early assessment of neoadjuvant therapy (NAT) response in triple-negative breast cancer (TNBC) would enable a treating oncologist to adjust a therapeutic plan of a non-responding patient, and thereby enhance outcomes while preventing unnecessary toxicities. To address this challenge, we propose leveraging personalized, in silico forecasts of tumor response to therapeutic regimens via a mechanistic mathematical model calibrated with patient-specific longitudinal magnetic resonance imaging (MRI) data acquired early during NAT. Here, we focus on identifying the driving mechanisms involved in the model formulation through a global sensitivity analysis. Our model describes the dynamics of TNBC cell density as a combination of mobility, which is formulated as a diffusion process constrained by local tumor-induced mechanical stress, and net proliferation, which is represented with a logistic term. To model TNBC response to NAT drug combinations, we adjust the tumor cell proliferation rate with a recent pharmacodynamic model, MuSyC, which also accounts for synergy of potency and efficacy. Tumor cell density is estimated from diffusion-weighted MRI data, while tissue mechanical properties are defined from segmented contrast-enhanced T1-weighted MRI data. To model the heterogeneous intratumoral delivery of drugs, we use perfusion maps estimated from dynamic contrast-enhanced MRI data. NAT drug pharmacokinetics are approximated with a linear model, which reasonably represents their temporal decay during NAT. Sobol’s method is used for the global sensitivity analysis of the NAT response of two different tumors (one well-perfused and one poorly-perfused) on a 3D, tissue-scale domain. This allows us to assess the total effect (ST) of each model parameter on tumor volume and global cellularity. Here, we focus on two standard NAT regimens: doxorubicin plus cyclophosphamide, and paclitaxel plus carboplatin. The parameter space is constructed by integrating three approaches. First, we use prior patient-specific in silico estimates of tumor cell mobility and proliferation. Second, we experimentally constrain the parameters accounting for potential synergistic drug activity by using time-resolved, high-throughput, automated microscopy assays that capture drug-induced changes in proliferation rates of various TNBC lines (HCC1143, SUM149, MDAMB231, and MDAMB468; perfosfamide was used in lieu of the pro-drug cyclophosphamide). Third, we scale the resulting in vitro parameter ranges to be clinically-relevant through an in silico study with our mechanistic model. Our results show that out of the 15 model parameters considered in the sensitivity analysis only a minority exhibited a driving role (ST > 0.1) in representing the dynamics of TNBC response to NAT, namely: the baseline tumor cell proliferation rate along with the effect and ratio of peak concentration to EC50 of doxorubicin, paclitaxel, and carboplatin. The other parameters have a marginal effect and can thus be fixed to any value within the parameter space. We select these constant parameters such that they contribute to simplifying the model formulation, leading to a surrogate reduced model. We further show that the reduced and the original model produce distributions of tumor volume and global cellularity that are in remarkable agreement both qualitatively (Dice similarity coefficient > 0.90 and > 0.94, respectively) and parameter-wise (concordance correlation coefficient > 0.85 and > 0.89, respectively). Thus, we conclude that our reduced model constitutes a feasible surrogate for future clinical calibration-forecasting studies, thereby facilitating personalized NAT response forecasting with patient-specific imaging datasets acquired in vivo. Citation Format: Guillermo Lorenzo, Angela M. Jarrett, Christian T. Meyer, Darren R. Tyson, Vito Quaranta, Thomas E. Yankeelov. In silico analysis of a novel mathematical model integrating in vitro and in vivo imaging data reveals driving mechanisms of breast cancer response to NAT for personalized tumor forecasting [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr P1-08-20.

Introduction:. Patients with locally advanced triple-negative breast cancer (TNBC) typically receive neoadjuvant therapy (NAT) to downstage the tumor and to improve the outcome of the subsequent breast conservation surgery. A critical unmet need is the lack of a method to accurately predict how a patient with TNBC will respond to NAT before surgery. In this work, we applied a clinical-computational framework to predict response of TNBC early in the course of NAT, by integrating quantitative MRI with mechanism-based mathematical modeling. Methods:. Patients and Data. Multiparametric quantitative MRI was acquired in patients (n = 46) before, and after 2 and 4 cycles of Adriamycin/Cyclophosphamide (A/C) regimen as part of the MD Anderson Cancer Center TNBC Moonshot Program. Within each imaging session, dynamic contrast-enhanced (DCE-), diffusion-weighted imaging (DWI), and a pre-contrast T1-map were acquired. Image processing. The processing pipeline consisted of three components. First, the images within each visit were registered to account for patient motion, and the parametric maps from the DCE and DWI images were computed. Second, inter-visit image registration was achieved by a non-rigid registration applied on breast, with a rigid penalty applied on the tumor region to preserve its size and shape. Third, post-processing was performed for preparation of modeling, including segmentation of the breast contour and tissues, and calculation of voxel-wise cellularity within tumors. Mathematical modeling. A predictive model was developed based on a reaction-diffusion equation (Eq. 1). The mobility of tumor cells is represented by diffusion coupled to mechanical properties of the tissue (Eq. 2), and the proliferation of the tumor is described with logistic growth. The injection and decay of administered therapies, inducing tumor cell death, is also represented in the model (Eq. 3). The variables and parameters used are listed in Table 1. Eq. 1: ∂N(x,t)/∂t = ∇⋅(D(x,t) ∇N(x,t)) + k(x) (1 - N(x,t)/θ)N(x,t) - (λ1(x,t) + λ2(x,t))N(x,t). Eq. 2: D(x,t) = D0 e-γσ(x,t). Eq. 3: λn(x,t) = αne-βn t C(x,t), n = 1, 2. For each patient, the domain and initial condition were generated from the pre-treatment images, and the images acquired during NAT were used for patient-specific calibration of parameters. The calibrated model was then used to predict the response to be observed at the end of NAT. We evaluated the model by comparing its predictions of tumor volume, longest axis, voxel-wise cellularity, and total tumor cellularity to the imaging measurements at the end of A/C. Results:. Our model predicted the tumor volume, total cellularity, and longest axis with a Pearson correlation coefficient (PCC) of 0.85, 0.80, and 0.60, respectively. The accuracy of voxel-wise cellularity achieved a PCC with the median (range) of 0.89 (0.77 - 0.93) between the prediction and the actual measurement. Moreover, we set criteria of 70% shrinkage of tumor volume to define response versus non-response cases, with which our model achieved a differentiation sensitivity/specificity of 0.90/0.73. Discussion:. Preliminary results of our study demonstrate the potential of the clinical-computational framework as a powerful tool for predicting response to NAT. Once validated, the method could also assist in optimizing treatment plans on a patient specific basis, or guiding patient selection in trials for novel NAT regimens. Table 1. Summary of the variables and parameters in the modelQuantitiesDefinition AssignmentDomainsΩbreast tissue domainGenerated from pre-treatment MRITEnd time point of NAT procedureDetermined from NAT schedulexCoordinate in breast tissueAssociated with spatial domain, ΩttimeAssociated with temporal domain, [0, T]VariablesN(x,t)Tumor cell numberInitialized from pre-treatment ADC, computed via Eq. 1D(x,t)Diffusive mobility of tumor cellsComputed via Eq. 2λn(x,t)Death rate induced by nth type of drugComputed via Eq. 3, n = 1 and 2 for A/Cσ(x,t)Von Mises stressComputed from gradient of N(x,t), based on Hormuth et al., 2018C(x,t)Spatiotemporal distribution of drugsAssigned based on NAT schedule and DCE imagesParametersk(x)Proliferation rate of tumor cellsLocally calibratedθTumor cells carry capacityGlobally calibratedαnEfficacy rate of nth type of drugGlobally calibratedβnDecay rate of of nth type of drugGlobally calibratedD0Diffusion coefficient of tumor cells in the absence of mechanical restrictionsGlobally calibratedγStress-tumor cell diffusion coupling constantAssigned based on Hormuth et al., 2018 Citation Format: Chengyue Wu, Angela M. Jarrett, Zijian Zhou, Nabil Elshafeey, Beatriz E. Adrada, Rosalind P. Candelaria, Rania M. Mohamed, Medine Boge, Lei Huo, Jason White, Debu Tripathy, Vicente Valero, Jennifer Litton, Stacy Moulder, Clinton Yam, Jong Bum Son, Jingfei Ma, Gaiane M. Rauch, Thomas E. Yankeelov. Forecasting treatment response to neoadjuvant systemic therapy in triple negative breast cancer viamathematical modeling and quantitative MRI [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr P1-08-08.

This protocol describes a complete data acquisition, analysis and computational forecasting pipeline for employing quantitative MRI data to predict the response of locally advanced breast cancer to neoadjuvant therapy in a community-based care setting. The methodology has previously been successfully applied to a heterogeneous patient population. The protocol details how to acquire the necessary images followed by registration, segmentation, quantitative perfusion and diffusion analysis, model calibration, and prediction. The data collection portion of the protocol requires ~25 min of scanning, postprocessing requires 2–3 h, and the model calibration and prediction components require ~10 h per patient depending on tumor size. The response of individual breast cancer patients to neoadjuvant therapy is forecast by application of a biophysical, reaction–diffusion mathematical model to these data. Successful application of the protocol results in coregistered MRI data from at least two scan visits that quantifies an individual tumor’s size, cellularity and vascular properties. This enables a spatially resolved prediction of how a particular patient’s tumor will respond to therapy. Expertise in image acquisition and analysis, as well as the numerical solution of partial differential equations, is required to carry out this protocol.

Background The purpose of this study was to determine whether advanced quantitative magnetic resonance imaging (MRI) can be deployed outside of large, research-oriented academic hospitals and into community care settings to predict eventual pathological complete response (pCR) to neoadjuvant therapy (NAT) in patients with locally advanced breast cancer. Methods Patients with stage II/III breast cancer (N = 28) were enrolled in a multicenter study performed in community radiology settings. Dynamic contrast-enhanced (DCE) and diffusion-weighted (DW)-MRI data were acquired at four time points during the course of NAT. Estimates of the vascular perfusion and permeability, as assessed by the volume transfer rate (Ktrans) using the Patlak model, were generated from the DCE-MRI data while estimates of cell density, as assessed by the apparent diffusion coefficient (ADC), were calculated from DW-MRI data. Tumor volume was calculated using semi-automatic segmentation and combined with Ktrans and ADC to yield bulk tumor blood flow and cellularity, respectively. The percent change in quantitative parameters at each MRI scan was calculated and compared to pathological response at the time of surgery. The predictive accuracy of each MRI parameter at different time points was quantified using receiver operating characteristic curves. Results Tumor size and quantitative MRI parameters were similar at baseline between groups that achieved pCR (n = 8) and those that did not (n = 20). Patients achieving a pCR had a larger decline in volume and cellularity than those who did not achieve pCR after one cycle of NAT (p Ktrans, ADC, cellularity, and bulk tumor flow from baseline (pre-treatment) were all significantly greater (p Conclusions Quantitative analysis of DCE-MRI and DW-MRI can be implemented in the community care setting to accurately predict the response of breast cancer to NAT. Dissemination of quantitative MRI into the community setting allows for the incorporation of these parameters into the standard of care and increases the number of clinical community sites able to participate in novel drug trials that require quantitative MRI.

The goal of this study is to calibrate a multiscale model of tumor angiogenesis with time-resolved data to allow for systematic testing of mathematical predictions of vascular sprouting. The multi-scale model consists of an agent-based description of tumor and endothelial cell dynamics coupled to a continuum model of vascular endothelial growth factor concentration. First, we calibrate ordinary differential equation models to time-resolved protein expression data to estimate the rates of secretion and consumption of vascular endothelial growth factor by endothelial and tumor cells, respectively. These parameters are then input into the multiscale tumor angiogenesis model, and the remaining model parameters are then calibrated to time resolved confocal microscopy images obtained within a 3D vascularized microfluidic platform. The microfluidic platform mimics a functional blood vessel with a surrounding collagen matrix seeded with inflammatory breast cancer cells, which induce tumor angiogenesis. Once the multi-scale model is fully parameterized, we forecast the spatiotemporal distribution of vascular sprouts at future time points and directly compare the predictions to experimentally measured data. We assess the ability of our model to globally recapitulate angiogenic vasculature density, resulting in an average relative calibration error of 17.7% ± 6.3% and an average prediction error of 20.2% ± 4% and 21.7% ± 3.6% using one and four calibrated parameters, respectively. We then assess the model’s ability to predict local vessel morphology (individualized vessel structure as opposed to global vascular density), initialized with the first time point and calibrated with two intermediate time points. To the best of our knowledge, this represents the first study to integrate well-controlled, experimental data into a mechanism-based, multiscale, mathematical model of angiogenic sprouting to make specific, testable predictions.

Introduction: Patients with locally advanced, triple-negative breast cancer (TNBC) typically receive neoadjuvant therapy (NAT) to downstage the tumor and to improve the outcome of subsequent breast conservation surgery. There are currently no methods to accurately predict how a TNBC patient will respond to NAT before surgery. In this work, we applied a digital twin framework to address this unmet clinical need, by integrating quantitative magnetic resonance imaging (MRI) data with mechanism-based mathematical modeling. Methods: Multiparametric MRI was acquired in patients (N = 50) before, after 2 and 4 cycles of Adriamycin/Cyclophosphamide (A/C), and again after 12 cycles of Paclitaxel as part of the ARTEMIS (NCT02276433) trial. Within each imaging session, dynamic contrast-enhanced (DCE) MRI, diffusion-weighted imaging (DWI), and a pre-contrast T1-map were acquired. The images were processed by a pipeline consisting of motion correction, multiparametric image alignment, inter-visit image registration to align the tumor and surrounding breast tissue, tissue segmentation, and estimation of tumor cellularity from DWI. A mechanism-based mathematical model, a reaction-diffusion equation, is used to characterize the mobility of tumor cells via diffusion damped by mechanical tissue properties, tumor proliferation via logistic growth, and treatment-induced cell death via the delivery and decay of therapies. For each patient, pre-treatment images were used for model initialization. The model calibration and prediction were implemented with two strategies: 1) using images acquired during the A/C for calibration and predicting up to the end of A/C, and 2) using images acquired during and after the A/C for calibration and predicting up to the end of NAT. For strategy 1), we evaluated the model by comparing its predicted tumor volume and total tumor cellularity to the imaging measurements at the end of A/C. For strategy 2), we evaluated the model by comparing its predicted final response to the post-surgical pathological findings. Results: For strategy 1), our framework predicted the change of tumor volume and total tumor cellularity with Pearson correlation coefficients of 0.91 and 0.89, respectively. Regarding strategy 2), our framework achieved an area under the receiver operator characteristic curve of 0.88 for distinguishing pCR from non-pCR. As a comparison, imaging measurement of tumor volume at the end of A/C achieved an AUC of 0.79. Conclusion: Our approach successfully captures the patient-specific dynamics of TNBC response to NAT and provides an improved prediction of final response, which demonstrates the potential of a digital twin framework to be a powerful tool for predicting response to NAT. Once validated, the method will provide a unique opportunity for optimizing treatment plans on a patient-specific basis. Citation Format: Chengyue Wu, Angela M. Jarrett, Zijian Zhou, Nabil Elshafeey, Beatriz E. Adrada, Rosalind P. Candelaria, Rania Mohamed, Medine Boge, Lei Huo, Jason White, Debu Tripathy, Vicente Valero, Jennifer Litton, Clinton Yam, Jong Bum Son, Jingfei Ma, Gaiane M. Rauch, Thomas E. Yankeelov. Forecasting treatment response to neoadjuvant therapy in triple-negative breast cancer via an image-guided digital twin [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 2736.

Background: Early and accurate prediction of response to neoadjuvant therapy (NAT) would empower personalization of breast cancer treatment regimens based on expected response. Noninvasive, quantitative dynamic contrast-enhanced (DCE-MRI) and diffusion-weighted magnetic resonance imaging (DW-MRI), when performed during the course of NAT, can accurately predict the ultimate pathological response. While these techniques have been incorporated into clinical trials in the academic setting, they have not yet been translated to the community setting, where most cancer patients receive their care. Implementation of quantitative MRI in the community setting widely expands the potential impact it can provide by 1) allowing community settings to participate in clinical trials that require quantitative MRI, and 2) advancing quantitative MRI towards standard-of-care for prediction of response in breast cancer. Methods: Women with locally advanced breast cancer (N = 28) were imaged four times during the course of NAT: 1) prior to the start of NAT, 2) after 1 cycle of NAT, 3) after 2-4 cycles of NAT, and 4) 1 cycle after MRI #3. Imaging data was acquired on 3T Siemens Skyra scanners equipped with breast coils and sited in a community hospital and radiology clinic, respectively. DW-MRI and DCE-MRI were acquired over 10 slices of 5 mm thickness. DW-MRI was acquired with diagonal monopolar diffusion-encoding gradients with b-values of 0, 200, and 800 s/mm2 in a total scan time of 1 minute 39 seconds. Voxel wise tumor cellularity was quantified using the apparent diffusion coefficient (ADC). For DCE-MRI, a gadolinium-based contrast agent was administered intravenously at 2 mL/sec after the acquisition of baseline scans. DCE-MRI data was acquired dynamically with a temporal resolution of 7.27 seconds for a total acquisition time of 8 minutes. The volume transfer constant Ktrans was calculated using Patlak analysis of DCE-MRI data to characterize the tumor vasculature. The tumor was semi-automatically segmented using a manually drawn region of interest followed by fuzzy c-means clustering of DCE-MRI data to identify a functional tumor volume. Measurements of tumor volume were combined with both ADC and Ktrans to yield metrics of tumor cellularity and bulk tumor flow, respectively. Results: Women who achieved pathological complete response at the time of surgery (pCR; n=8) displayed significantly different treatment-induced changes in MRI-derived tumor parameters versus women who did not achieve pCR (non-pCR, n=20). After 1 cycle of NAT, women who achieved pCR had smaller functional tumor volume and lower cellularity (p < 0.05) than non-pCR study participants. At the third and fourth MRI, tumor volume, ADC, Ktrans, cellularity, and bulk tumor flow were all significantly different between the pCR and non-pCR cohorts (p < 0.05). Of note, longest tumor diameter was not predictive of pCR at any time point in this study. Conclusions: This study demonstrates that quantitative DCE- and DW-MRI can be implemented successfully in community care facilities within standard-of-care settings for imaging locally advanced breast cancer. Metrics extracted from the change in DW-MRI (ADC) and DCE-MRI (Ktrans) can accurately predict pathological complete response to neoadjuvant therapy and may be more sensitive to tumor response than the RECIST criteria. Furthermore, incorporating quantitative metrics with tumor volume further increases the ability to predict pathological response to NAT in locally advanced breast cancer. While there are still challenges to address to effectively implement these quantitative metrics into the clinical workflow, this study is first in its kind to transition a decade’s worth of quantitative MRI advancements from academic settings into standard-of-care. Citation Format: John Virostko, Anna G Sorace, Kalina P Slavkova, Anum S Kazerouni, Angela M Jarrett, Julie C DiCarlo, Stefanie Woodard, Sarah Avery, Boone W Goodgame, Debra Patt, Thomas E Yankeelov. Quantitative multiparametric MRI predicts response to neoadjuvant therapy in the community setting [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr P3-03-03.

Background: Standard of care (SOC) breast MRI exams typically acquire 4-7 frames of dynamic contrast-enhanced MRI (DCE-MRI) for cancer screening and staging. Post-contrast images depict lesion spiculations and boundaries to identify and characterize tumors. Pharmacokinetic (PK) analysis of DCE-MRI involves modeling blood flow to the lesion and surrounding tissue and has shown promise in diagnosis and prediction of therapeutic response. Currently, SOC DCE-MRI requires ~60-90 seconds per volume for images with sufficient quality and spatial resolution. However, PK analysis of DCE-MRI requires faster time course sampling. For this reason, PK modeling is limited to research scans with lower spatial resolution and higher temporal resolution. PK modeling would improve feedback of treatment response, and implementation in the SOC exam would increase imaging trial participation. In this study, we tested the estimation of Ktrans, a mixed perfusion and permeability PK parameter, from three images at optimal time points after contrast agent (CA) injection, and compared it to the Ktrans estimation from analysis of the full-length time course.. Methods: Women (N=23) with newly diagnosed invasive breast cancers who were eligible for neoadjuvant therapy (NAT) were scanned with a research MRI protocol as part of a treatment-monitoring study. Images acquired prior to the start of NAT were used. MRI was performed on 3.0T Siemens Skyra scanners at two sites with bilateral breast coils. The research protocol included ten sagittal slices centered about the primary tumor. The DCE-MRI images came from a fast sequence with 1.3 × 1.3 × 5.0 mm resolution acquired at 7.3 seconds per frame (66 frames total,) with a gadolinium-based CA injected one minute into the scan. A population arterial input function was used to implement a mathematical graph-based search of possible tissue CA concentration curves from the expected range of PK parameters. The search results gave a set of three optimized sub-sampled timepoints, Topt, from the full set of sample times, Tfull, at which to best sample the CA concentration curves to optimally estimate PK values. The imaging data was analyzed to find one parameter map from image times Tfull, and another from the subset of images at times Topt. The difference in Ktrans was computed at each parameter map voxel, and the concordance correlation coefficient (CCC) was computed per patient to determine agreement. The median Ktrans values were also compared for each patient. Results: The graph-based search of CA concentration curves found optimal times Topt of 37, 66, and 153 seconds after injection. The averaged values over all patients for median and maximum Ktrans from the original Tfull image set were 0.07 and 0.5 (min)-1. The average difference in Ktrans values between the Topt and Tfull sets was 0.02 (min)-1. When the median Ktrans values for each patient were compared, the average difference in median Ktrans values was 15% +/- 9%. The concordance correlation coefficients comparing the Topt and Tfull -sampled parameter maps for each patient were 0.89 +/- 0.12, showing high agreement. Discussion: This retrospective analysis suggests that it is possible to estimate PK parameters from a few properly selected post-contrast images inserted into a SOC DCE-MRI exam. The combination of optimal timing with fast acquisition techniques for high-resolution imaging could be used to provide quantitative data while preserving post-contrast images with the necessary spatial resolution for clinical reading. Importantly, the test images were acquired in the community setting with widely available MRI hardware, further indicating the potential for integration with SOC exams. Funding: NCI U24 CA226110, NCI U01 CA174706, NCI U01 CA142565, CPRIT RR160005 Citation Format: Julie C DiCarlo, Angela M Jarrett, Anum S Kazerouni, John Virostko, Anna G Sorace, Kalina P Slavkova, Debra Patt, Boone W Goodgame, Sarah Avery, Thomas E Yankeelov. Three timepoint pharmacokinetic modeling to incorporate within standard of care MRI breast exams [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr P3-02-09.