1
|
Miniere HJM, Lima EABF, Lorenzo G, Hormuth II DA, Ty S, Brock A, Yankeelov TE. A mathematical model for predicting the spatiotemporal response of breast cancer cells treated with doxorubicin. Cancer Biol Ther 2024; 25:2321769. [PMID: 38411436 PMCID: PMC11057790 DOI: 10.1080/15384047.2024.2321769] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2023] [Accepted: 02/18/2024] [Indexed: 02/28/2024] Open
Abstract
Tumor heterogeneity contributes significantly to chemoresistance, a leading cause of treatment failure. To better personalize therapies, it is essential to develop tools capable of identifying and predicting intra- and inter-tumor heterogeneities. Biology-inspired mathematical models are capable of attacking this problem, but tumor heterogeneity is often overlooked in in-vivo modeling studies, while phenotypic considerations capturing spatial dynamics are not typically included in in-vitro modeling studies. We present a data assimilation-prediction pipeline with a two-phenotype model that includes a spatiotemporal component to characterize and predict the evolution of in-vitro breast cancer cells and their heterogeneous response to chemotherapy. Our model assumes that the cells can be divided into two subpopulations: surviving cells unaffected by the treatment, and irreversibly damaged cells undergoing treatment-induced death. MCF7 breast cancer cells were previously cultivated in wells for up to 1000 hours, treated with various concentrations of doxorubicin and imaged with time-resolved microscopy to record spatiotemporally-resolved cell count data. Images were used to generate cell density maps. Treatment response predictions were initialized by a training set and updated by weekly measurements. Our mathematical model successfully calibrated the spatiotemporal cell growth dynamics, achieving median [range] concordance correlation coefficients of > .99 [.88, >.99] and .73 [.58, .85] across the whole well and individual pixels, respectively. Our proposed data assimilation-prediction approach achieved values of .97 [.44, >.99] and .69 [.35, .79] for the whole well and individual pixels, respectively. Thus, our model can capture and predict the spatiotemporal dynamics of MCF7 cells treated with doxorubicin in an in-vitro setting.
Collapse
Affiliation(s)
- Hugo J. M. Miniere
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, USA
| | - Ernesto A. B. F. Lima
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, USA
| | - Guillermo Lorenzo
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, USA
- Department of Civil Engineering and Architecture, University of Pavia, Lombardy, Italy
| | - David A. Hormuth II
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, USA
- Livestrong Cancer Institutes, The University of Texas at Austin, Austin, USA
| | - Sophia Ty
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, USA
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, USA
| | - Amy Brock
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, USA
- Livestrong Cancer Institutes, The University of Texas at Austin, Austin, USA
| | - Thomas E. Yankeelov
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, USA
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, USA
- Livestrong Cancer Institutes, The University of Texas at Austin, Austin, USA
- Department of Diagnostic Medicine, The University of Texas at Austin, Austin, USA
- Department of Oncology, The University of Texas at Austin, Austin, USA
- Division of Diagnostic Imaging, The University of Texas M.D. Anderson Cancer Center, Houston, USA
| |
Collapse
|
2
|
Lorenzo G, Ahmed SR, Hormuth DA, Vaughn B, Kalpathy-Cramer J, Solorio L, Yankeelov TE, Gomez H. Patient-Specific, Mechanistic Models of Tumor Growth Incorporating Artificial Intelligence and Big Data. Annu Rev Biomed Eng 2024; 26:529-560. [PMID: 38594947 DOI: 10.1146/annurev-bioeng-081623-025834] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/11/2024]
Abstract
Despite the remarkable advances in cancer diagnosis, treatment, and management over the past decade, malignant tumors remain a major public health problem. Further progress in combating cancer may be enabled by personalizing the delivery of therapies according to the predicted response for each individual patient. The design of personalized therapies requires the integration of patient-specific information with an appropriate mathematical model of tumor response. A fundamental barrier to realizing this paradigm is the current lack of a rigorous yet practical mathematical theory of tumor initiation, development, invasion, and response to therapy. We begin this review with an overview of different approaches to modeling tumor growth and treatment, including mechanistic as well as data-driven models based on big data and artificial intelligence. We then present illustrative examples of mathematical models manifesting their utility and discuss the limitations of stand-alone mechanistic and data-driven models. We then discuss the potential of mechanistic models for not only predicting but also optimizing response to therapy on a patient-specific basis. We describe current efforts and future possibilities to integrate mechanistic and data-driven models. We conclude by proposing five fundamental challenges that must be addressed to fully realize personalized care for cancer patients driven by computational models.
Collapse
Affiliation(s)
- Guillermo Lorenzo
- Oden Institute for Computational Engineering and Sciences, University of Texas, Austin, Texas, USA
- Department of Civil Engineering and Architecture, University of Pavia, Pavia, Italy
| | - Syed Rakin Ahmed
- Geisel School of Medicine, Dartmouth College, Hanover, New Hampshire, USA
- Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
- Harvard Graduate Program in Biophysics, Harvard Medical School, Harvard University, Cambridge, Massachusetts, USA
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Boston, Massachusetts, USA
| | - David A Hormuth
- Livestrong Cancer Institutes, University of Texas, Austin, Texas, USA
- Oden Institute for Computational Engineering and Sciences, University of Texas, Austin, Texas, USA
| | - Brenna Vaughn
- Weldon School of Biomedical Engineering, Purdue University, West Lafayette, Indiana, USA;
| | | | - Luis Solorio
- Weldon School of Biomedical Engineering, Purdue University, West Lafayette, Indiana, USA;
| | - Thomas E Yankeelov
- Department of Imaging Physics, MD Anderson Cancer Center, Houston, Texas, USA
- Department of Biomedical Engineering, Department of Oncology, and Department of Diagnostic Medicine, University of Texas, Austin, Texas, USA
- Livestrong Cancer Institutes, University of Texas, Austin, Texas, USA
- Oden Institute for Computational Engineering and Sciences, University of Texas, Austin, Texas, USA
| | - Hector Gomez
- School of Mechanical Engineering and Purdue Center for Cancer Research, Purdue University, West Lafayette, Indiana, USA
- Weldon School of Biomedical Engineering, Purdue University, West Lafayette, Indiana, USA;
| |
Collapse
|
3
|
Bhandari A, Gu B, Kashkooli FM, Zhan W. Image-based predictive modelling frameworks for personalised drug delivery in cancer therapy. J Control Release 2024; 370:721-746. [PMID: 38718876 DOI: 10.1016/j.jconrel.2024.05.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2024] [Revised: 04/11/2024] [Accepted: 05/02/2024] [Indexed: 05/19/2024]
Abstract
Personalised drug delivery enables a tailored treatment plan for each patient compared to conventional drug delivery, where a generic strategy is commonly employed. It can not only achieve precise treatment to improve effectiveness but also reduce the risk of adverse effects to improve patients' quality of life. Drug delivery involves multiple interconnected physiological and physicochemical processes, which span a wide range of time and length scales. How to consider the impact of individual differences on these processes becomes critical. Multiphysics models are an open system that allows well-controlled studies on the individual and combined effects of influencing factors on drug delivery outcomes while accommodating the patient-specific in vivo environment, which is not economically feasible through experimental means. Extensive modelling frameworks have been developed to reveal the underlying mechanisms of drug delivery and optimise effective delivery plans. This review provides an overview of currently available models, their integration with advanced medical imaging modalities, and code packages for personalised drug delivery. The potential to incorporate new technologies (i.e., machine learning) in this field is also addressed for development.
Collapse
Affiliation(s)
- Ajay Bhandari
- Biofluids Research Lab, Department of Mechanical Engineering, Indian Institute of Technology (Indian School of Mines), Dhanbad, India
| | - Boram Gu
- School of Chemical Engineering, Chonnam National University, Gwangju, Republic of Korea
| | | | - Wenbo Zhan
- School of Engineering, University of Aberdeen, Aberdeen, UK.
| |
Collapse
|
4
|
Lorenzo G, Heiselman JS, Liss MA, Miga MI, Gomez H, Yankeelov TE, Reali A, Hughes TJ. A Pilot Study on Patient-specific Computational Forecasting of Prostate Cancer Growth during Active Surveillance Using an Imaging-informed Biomechanistic Model. CANCER RESEARCH COMMUNICATIONS 2024; 4:617-633. [PMID: 38426815 PMCID: PMC10906139 DOI: 10.1158/2767-9764.crc-23-0449] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/10/2023] [Revised: 12/15/2023] [Accepted: 02/09/2024] [Indexed: 03/02/2024]
Abstract
Active surveillance (AS) is a suitable management option for newly diagnosed prostate cancer, which usually presents low to intermediate clinical risk. Patients enrolled in AS have their tumor monitored via longitudinal multiparametric MRI (mpMRI), PSA tests, and biopsies. Hence, treatment is prescribed when these tests identify progression to higher-risk prostate cancer. However, current AS protocols rely on detecting tumor progression through direct observation according to population-based monitoring strategies. This approach limits the design of patient-specific AS plans and may delay the detection of tumor progression. Here, we present a pilot study to address these issues by leveraging personalized computational predictions of prostate cancer growth. Our forecasts are obtained with a spatiotemporal biomechanistic model informed by patient-specific longitudinal mpMRI data (T2-weighted MRI and apparent diffusion coefficient maps from diffusion-weighted MRI). Our results show that our technology can represent and forecast the global tumor burden for individual patients, achieving concordance correlation coefficients from 0.93 to 0.99 across our cohort (n = 7). In addition, we identify a model-based biomarker of higher-risk prostate cancer: the mean proliferation activity of the tumor (P = 0.041). Using logistic regression, we construct a prostate cancer risk classifier based on this biomarker that achieves an area under the ROC curve of 0.83. We further show that coupling our tumor forecasts with this prostate cancer risk classifier enables the early identification of prostate cancer progression to higher-risk disease by more than 1 year. Thus, we posit that our predictive technology constitutes a promising clinical decision-making tool to design personalized AS plans for patients with prostate cancer. SIGNIFICANCE Personalization of a biomechanistic model of prostate cancer with mpMRI data enables the prediction of tumor progression, thereby showing promise to guide clinical decision-making during AS for each individual patient.
Collapse
Affiliation(s)
- Guillermo Lorenzo
- Department of Civil Engineering and Architecture, University of Pavia, Pavia, Italy
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas
| | - Jon S. Heiselman
- Department of Biomedical Engineering, Vanderbilt University, Nashville, Tennessee
- Department of Surgery, Memorial Sloan-Kettering Cancer Center, New York, New York
| | - Michael A. Liss
- Department of Urology, The University of Texas Health Science Center at San Antonio, San Antonio, Texas
| | - Michael I. Miga
- Department of Biomedical Engineering, Vanderbilt University, Nashville, Tennessee
- Vanderbilt Institute for Surgery and Engineering, Vanderbilt University, Nashville, Tennessee
- Department of Neurological Surgery, Radiology, and Otolaryngology-Head and Neck Surgery, Vanderbilt University Medical Center, Nashville, Tennessee
| | - Hector Gomez
- School of Mechanical Engineering, Weldon School of Biomedical Engineering, and Purdue Institute for Cancer Research, Purdue University, West Lafayette, Indiana
| | - Thomas E. Yankeelov
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas
- Livestrong Cancer Institutes and Departments of Biomedical Engineering, Diagnostic Medicine, and Oncology, The University of Texas at Austin, Austin, Texas
- Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Alessandro Reali
- Department of Civil Engineering and Architecture, University of Pavia, Pavia, Italy
| | - Thomas J.R. Hughes
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas
| |
Collapse
|
5
|
Ren Z, Pineda FD, Howard FM, Fan X, Nanda R, Abe H, Kulkarni K, Karczmar GS. Bilateral asymmetry of quantitative parenchymal kinetics at ultrafast DCE-MRI predict response to neoadjuvant chemotherapy in patients with HER2+ breast cancer. Magn Reson Imaging 2023; 104:9-15. [PMID: 37611646 PMCID: PMC10879456 DOI: 10.1016/j.mri.2023.08.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2023] [Revised: 08/08/2023] [Accepted: 08/20/2023] [Indexed: 08/25/2023]
Abstract
PURPOSE To assess whether measurement of the bilateral asymmetry of semiquantitative and quantitative perfusion parameters from ultrafast dynamic contrast-enhanced MRI (DCE-MRI), allows early prediction of pathologic response after neoadjuvant chemotherapy (NAC) in patients with HER2+ breast cancer. MATERIALS AND METHODS Twenty-eight female patients with HER2+ breast cancer treated with NAC who underwent pre-NAC ultrafast DCE-MRI (3-9 s/phase) were enrolled for this study. Four semiquantitative and two quantitative parenchymal parameters were calculated for each patient. Ipsilateral/contralateral (I/C) ratio (for four parameters) and the difference between (for two parameters) ipsi- and contra-lateral parenchymal kinetic parameters (kBPE) were compared for patients with pathologic complete response (pCR) and those having residual disease. Lasso regression with leave-one-out cross validation was used to determine the optimal combination of parameters for a regression model and multivariable logistic regression was used to identify independent predictors for pCR. Chi-squared test, two-sided t-test and Kruskal-Wallis test were used. RESULTS The Ktrans I/C ratio cutoff value of 1.11 had a sensitivity of 83.3% and specificity of 75% for pCR. The ve I/C ratio cutoff value of 1.1 had a sensitivity of 75% and specificity of 81.3% for pCR. The area under the receiver operating characteristic curve of the three-kBPE parameter model, including initial area under the enhancement curve (AUC30) I/C ratio, KtransI/C ratio and ve I/C ratio, was 0.89 with sensitivity of 91.7% at specificity of 81.3%. CONCLUSION Quantitative assessment of bilateral asymmetry kBPE from pre-NAC ultrafast DCE-MRI can predict pCR in patients with HER2+ breast cancer.
Collapse
Affiliation(s)
- Zhen Ren
- Department of Radiology, The University of Chicago, Chicago, IL 60637, United States of America.
| | - Federico D Pineda
- Department of Radiology, University of Pittsburgh, Pittsburgh, PA 15260, United States of America.
| | - Frederick M Howard
- Section of Hematology and Oncology, Department of Medicine, The University of Chicago, Chicago, IL 60637, United States of America.
| | - Xiaobing Fan
- Department of Radiology, The University of Chicago, Chicago, IL 60637, United States of America.
| | - Rita Nanda
- Section of Hematology and Oncology, Department of Medicine, The University of Chicago, Chicago, IL 60637, United States of America.
| | - Hiroyuki Abe
- Department of Radiology, The University of Chicago, Chicago, IL 60637, United States of America.
| | - Kirti Kulkarni
- Department of Radiology, The University of Chicago, Chicago, IL 60637, United States of America.
| | - Gregory S Karczmar
- Department of Radiology, The University of Chicago, Chicago, IL 60637, United States of America.
| |
Collapse
|
6
|
Wagner A, Schlicke P, Fritz M, Kuttler C, Oden JT, Schumann C, Wohlmuth B. A phase-field model for non-small cell lung cancer under the effects of immunotherapy. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2023; 20:18670-18694. [PMID: 38052574 DOI: 10.3934/mbe.2023828] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/07/2023]
Abstract
Formulating mathematical models that estimate tumor growth under therapy is vital for improving patient-specific treatment plans. In this context, we present our recent work on simulating non-small-scale cell lung cancer (NSCLC) in a simple, deterministic setting for two different patients receiving an immunotherapeutic treatment. At its core, our model consists of a Cahn-Hilliard-based phase-field model describing the evolution of proliferative and necrotic tumor cells. These are coupled to a simplified nutrient model that drives the growth of the proliferative cells and their decay into necrotic cells. The applied immunotherapy decreases the proliferative cell concentration. Here, we model the immunotherapeutic agent concentration in the entire lung over time by an ordinary differential equation (ODE). Finally, reaction terms provide a coupling between all these equations. By assuming spherical, symmetric tumor growth and constant nutrient inflow, we simplify this full 3D cancer simulation model to a reduced 1D model. We can then resort to patient data gathered from computed tomography (CT) scans over several years to calibrate our model. Our model covers the case in which the immunotherapy is successful and limits the tumor size, as well as the case predicting a sudden relapse, leading to exponential tumor growth. Finally, we move from the reduced model back to the full 3D cancer simulation in the lung tissue. Thereby, we demonstrate the predictive benefits that a more detailed patient-specific simulation including spatial information as a possible generalization within our framework could yield in the future.
Collapse
Affiliation(s)
- Andreas Wagner
- School of Computation, Information and Technology, Technical University of Munich, Munich, Bavaria, Germany
| | - Pirmin Schlicke
- School of Computation, Information and Technology, Technical University of Munich, Munich, Bavaria, Germany
| | - Marvin Fritz
- Computational Methods for PDEs, Johann Radon Institute for Computational and Applied Mathematics, Linz, Upper Austria, Austria
| | - Christina Kuttler
- School of Computation, Information and Technology, Technical University of Munich, Munich, Bavaria, Germany
| | - J Tinsley Oden
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas, United States of America
| | - Christian Schumann
- Clinic of Pneumology, Thoracic Oncology, Sleep and Respiratory Critical Care, Klinikverbund Allgäu, Kempten, Bavaria, Germany
| | - Barbara Wohlmuth
- School of Computation, Information and Technology, Technical University of Munich, Munich, Bavaria, Germany
| |
Collapse
|
7
|
Lima EA, Wyde RA, Sorace AG, Yankeelov TE. Optimizing combination therapy in a murine model of HER2+ breast cancer. COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING 2022; 402:115484. [PMID: 37800167 PMCID: PMC10552906] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Subscribe] [Scholar Register] [Indexed: 10/07/2023]
Abstract
Human epidermal growth factor receptor 2 positive (HER2+) breast cancer is frequently treated with drugs that target the HER2 receptor, such as trastuzumab, in combination with chemotherapy, such as doxorubicin. However, an open problem in treatment design is to determine the therapeutic regimen that optimally combines these two treatments to yield optimal tumor control. Working with data quantifying temporal changes in tumor volume due to different trastuzumab and doxorubicin treatment protocols in a murine model of human HER2+ breast cancer, we propose a complete framework for model development, calibration, selection, and treatment optimization to find the optimal treatment protocol. Through different assumptions for the drug-tumor interactions, we propose ten different models to characterize the dynamic relationship between tumor volume and drug availability, as well as the drug-drug interaction. Using a Bayesian framework, each of these models are calibrated to the dataset and the model with the highest Bayesian information criterion weight is selected to represent the biological system. The selected model captures the inhibition of trastuzumab due to pre-treatment with doxorubicin, as well as the increase in doxorubicin efficacy due to pre-treatment with trastuzumab. We then apply optimal control theory (OCT) to this model to identify two optimal treatment protocols. In the first optimized protocol, we fix the maximum dosage for doxorubicin and trastuzumab to be the same as the maximum dose delivered experimentally, while trying to minimize tumor burden. Within this constraint, optimal control theory indicates the optimal regimen is to first deliver two doses of trastuzumab on days 35 and 36, followed by two doses of doxorubicin on days 37 and 38. This protocol predicts an additional 45% reduction in tumor burden compared to that achieved with the experimentally delivered regimen. In the second optimized protocol we fix the tumor control to be the same as that obtained experimentally, and attempt to reduce the doxorubicin dose. Within this constraint, the optimal regimen is the same as the first optimized protocol but uses only 43% of the doxorubicin dose used experimentally. This protocol predicts tumor control equivalent to that achieved experimentally. These results strongly suggest the utility of mathematical modeling and optimal control theory for identifying therapeutic regimens maximizing efficacy and minimizing toxicity.
Collapse
Affiliation(s)
- Ernesto A.B.F. Lima
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, United States of America
- Texas Advanced Computing Center, The University of Texas at Austin, United States of America
| | - Reid A.F. Wyde
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, United States of America
| | - Anna G. Sorace
- Department of Radiology, The University of Alabama at Birmingham, United States of America
- Department of Biomedical Engineering, The University of Alabama at Birmingham, United States of America
- O’Neal Comprehensive Cancer Center, The University of Alabama at Birmingham, United States of America
| | - Thomas E. Yankeelov
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, United States of America
- Department of Biomedical Engineering, The University of Texas at Austin, United States of America
- Department of Diagnostic Medicine, The University of Texas at Austin, United States of America
- Department of Oncology, The University of Texas at Austin, United States of America
- Livestrong Cancer Institutes, Dell Medical School, The University of Texas at Austin, United States of America
- Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, United States of America
| |
Collapse
|
8
|
Wu C, Jarrett AM, Zhou Z, Elshafeey N, Adrada BE, Candelaria RP, Mohamed RMM, Boge M, Huo L, White JB, Tripathy D, Valero V, Litton JK, Yam C, Son JB, Ma J, Rauch GM, Yankeelov TE. MRI-Based Digital Models Forecast Patient-Specific Treatment Responses to Neoadjuvant Chemotherapy in Triple-Negative Breast Cancer. Cancer Res 2022; 82:3394-3404. [PMID: 35914239 PMCID: PMC9481712 DOI: 10.1158/0008-5472.can-22-1329] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2022] [Revised: 06/14/2022] [Accepted: 07/26/2022] [Indexed: 02/07/2023]
Abstract
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.
Collapse
Affiliation(s)
- Chengyue Wu
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas
| | - Angela M Jarrett
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas.,Livestrong Cancer Institutes, The University of Texas at Austin, Austin, Texas
| | - Zijian Zhou
- Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Nabil Elshafeey
- Department of Abdominal Imaging, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Beatriz E Adrada
- Department of Breast Imaging, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Rosalind P Candelaria
- Department of Breast Imaging, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Rania M M Mohamed
- Department of Breast Imaging, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Medine Boge
- Department of Breast Imaging, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Lei Huo
- Department of Pathology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Jason B White
- Department of Breast Medical Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Debu Tripathy
- Department of Breast Medical Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Vicente Valero
- Department of Breast Medical Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Jennifer K Litton
- Department of Breast Medical Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Clinton Yam
- Department of Breast Medical Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Jong Bum Son
- Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Jingfei Ma
- Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Gaiane M Rauch
- Department of Abdominal Imaging, The University of Texas MD Anderson Cancer Center, Houston, Texas.,Department of Breast Imaging, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Thomas E Yankeelov
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas.,Livestrong Cancer Institutes, The University of Texas at Austin, Austin, Texas.,Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, Houston, Texas.,Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas.,Department of Diagnostic Medicine, The University of Texas at Austin, Austin, Texas.,Department of Oncology, The University of Texas at Austin, Austin, Texas
| |
Collapse
|
9
|
Wu C, Lorenzo G, Hormuth DA, Lima EABF, Slavkova KP, DiCarlo JC, Virostko J, Phillips CM, Patt D, Chung C, Yankeelov TE. Integrating mechanism-based modeling with biomedical imaging to build practical digital twins for clinical oncology. BIOPHYSICS REVIEWS 2022; 3:021304. [PMID: 35602761 PMCID: PMC9119003 DOI: 10.1063/5.0086789] [Citation(s) in RCA: 22] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/28/2022] [Accepted: 04/29/2022] [Indexed: 12/11/2022]
Abstract
Digital twins employ mathematical and computational models to virtually represent a physical object (e.g., planes and human organs), predict the behavior of the object, and enable decision-making to optimize the future behavior of the object. While digital twins have been widely used in engineering for decades, their applications to oncology are only just emerging. Due to advances in experimental techniques quantitatively characterizing cancer, as well as advances in the mathematical and computational sciences, the notion of building and applying digital twins to understand tumor dynamics and personalize the care of cancer patients has been increasingly appreciated. In this review, we present the opportunities and challenges of applying digital twins in clinical oncology, with a particular focus on integrating medical imaging with mechanism-based, tissue-scale mathematical modeling. Specifically, we first introduce the general digital twin framework and then illustrate existing applications of image-guided digital twins in healthcare. Next, we detail both the imaging and modeling techniques that provide practical opportunities to build patient-specific digital twins for oncology. We then describe the current challenges and limitations in developing image-guided, mechanism-based digital twins for oncology along with potential solutions. We conclude by outlining five fundamental questions that can serve as a roadmap when designing and building a practical digital twin for oncology and attempt to provide answers for a specific application to brain cancer. We hope that this contribution provides motivation for the imaging science, oncology, and computational communities to develop practical digital twin technologies to improve the care of patients battling cancer.
Collapse
Affiliation(s)
- Chengyue Wu
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas 78712, USA
| | | | | | | | - Kalina P. Slavkova
- Department of Physics, The University of Texas at Austin, Austin, Texas 78712, USA
| | | | | | - Caleb M. Phillips
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas 78712, USA
| | - Debra Patt
- Texas Oncology, Austin, Texas 78731, USA
| | - Caroline Chung
- Department of Radiation Oncology, MD Anderson Cancer Center, University of Texas, Houston, Texas 77030, USA
| | | |
Collapse
|