1
|
Butner JD, Dogra P, Chung C, Koay EJ, Welsh JW, Hong DS, Cristini V, Wang Z. Hybridizing mechanistic mathematical modeling with deep learning methods to predict individual cancer patient survival after immune checkpoint inhibitor therapy. Res Sq 2024:rs.3.rs-4151883. [PMID: 38586046 PMCID: PMC10996814 DOI: 10.21203/rs.3.rs-4151883/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/09/2024]
Abstract
We present a study where predictive mechanistic modeling is used in combination with deep learning methods to predict individual patient survival probabilities under immune checkpoint inhibitor (ICI) therapy. This hybrid approach enables prediction based on both measures that are calculable from mechanistic models (but may not be directly measurable in the clinic) and easily measurable quantities or characteristics (that are not always readily incorporated into predictive mechanistic models). The mechanistic model we have applied here can predict tumor response from CT or MRI imaging based on key mechanisms underlying checkpoint inhibitor therapy, and in the present work, its parameters were combined with readily-available clinical measures from 93 patients into a hybrid training set for a deep learning time-to-event predictive model. Analysis revealed that training an artificial neural network with both mechanistic modeling-derived and clinical measures achieved higher per-patient predictive accuracy based on event-time concordance, Brier score, and negative binomial log-likelihood-based criteria than when only mechanistic model-derived values or only clinical data were used. Feature importance analysis revealed that both clinical and model-derived parameters play prominent roles in neural network decision making, and in increasing prediction accuracy, further supporting the advantage of our hybrid approach. We anticipate that many existing mechanistic models may be hybridized with deep learning methods in a similar manner to improve predictive accuracy through addition of additional data that may not be readily implemented in mechanistic descriptions.
Collapse
Affiliation(s)
- Joseph D Butner
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA
- Institute for Data Science in Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA
- Master in Clinical Translation Management Program, The Cameron School of Business, University of St. Thomas, Houston, TX 77006, USA
| | - Prashant Dogra
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA
- Department of Physiology and Biophysics, Weill Cornell Medicine, New York, NY 10065, USA
| | - Caroline Chung
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA
- Institute for Data Science in Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA
| | - Eugene J Koay
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA
| | - James W Welsh
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA
| | - David S Hong
- Department of Investigational Cancer Therapeutics, University of Texas MD Anderson Cancer Center, Houston, Texas 77230, USA
| | - Vittorio Cristini
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA
- Neal Cancer Center, Houston Methodist Research Institute, Houston, TX 77030, USA
- Physiology, Biophysics, and Systems Biology Program, Graduate School of Medical Sciences, Weill Cornell Medicine, New York, NY 10065, USA
- Department of Imaging Physics, University of Texas MD Anderson Cancer Center, Houston, TX 77230, USA
| | - Zhihui Wang
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA
- Department of Physiology and Biophysics, Weill Cornell Medicine, New York, NY 10065, USA
- Neal Cancer Center, Houston Methodist Research Institute, Houston, TX 77030, USA
- Department of Medical Education, Texas A&M University School of Medicine, Bryan, TX 77807, USA
| |
Collapse
|
2
|
Dogra P, Shinglot V, Ruiz-Ramírez J, Cave J, Butner JD, Schiavone C, Duda DG, Kaseb AO, Chung C, Koay EJ, Cristini V, Ozpolat B, Calin GA, Wang Z. Translational modeling-based evidence for enhanced efficacy of standard-of-care drugs in combination with anti-microRNA-155 in non-small-cell lung cancer. medRxiv 2024:2024.03.14.24304306. [PMID: 38559070 PMCID: PMC10980136 DOI: 10.1101/2024.03.14.24304306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
Elevated microRNA-155 (miR-155) expression in non-small-cell lung cancer (NSCLC) promotes cisplatin resistance and negatively impacts treatment outcomes. However, miR-155 can also boost anti-tumor immunity by suppressing PD-L1 expression. We developed a multiscale mechanistic model, calibrated with in vivo data and then extrapolated to humans, to investigate the therapeutic effects of nanoparticle-delivered anti-miR-155 in NSCLC, alone or in combination with standard-of-care drugs. Model simulations and analyses of the clinical scenario revealed that monotherapy with anti-miR-155 at a dose of 2.5 mg/kg administered once every three weeks has substantial anti-cancer activity. It led to a median progression-free survival (PFS) of 6.7 months, which compared favorably to cisplatin and immune checkpoint inhibitors. Further, we explored the combinations of anti-miR-155 with standard-of-care drugs, and found strongly synergistic two- and three-drug combinations. A three-drug combination of anti-miR-155, cisplatin, and pembrolizumab resulted in a median PFS of 13.1 months, while a two-drug combination of anti-miR-155 and cisplatin resulted in a median PFS of 11.3 months, which emerged as a more practical option due to its simple design and cost-effectiveness. Our analyses also provided valuable insights into unfavorable dose ratios for drug combinations, highlighting the need for optimizing dose regimen to prevent antagonistic effects. Thus, this work bridges the gap between preclinical development and clinical translation of anti-miR-155 and unravels the potential of anti-miR-155 combination therapies in NSCLC.
Collapse
Affiliation(s)
- Prashant Dogra
- Mathematics in Medicine Program, Department of Medicine, Houston Methodist Research Institute, Houston, TX, USA
- Department of Physiology and Biophysics, Weill Cornell Medical College, New York, NY, USA
| | - Vrushaly Shinglot
- Mathematics in Medicine Program, Department of Medicine, Houston Methodist Research Institute, Houston, TX, USA
| | | | - Joseph Cave
- Mathematics in Medicine Program, Department of Medicine, Houston Methodist Research Institute, Houston, TX, USA
- Physiology, Biophysics, and Systems Biology Program, Graduate School of Medical Sciences, Weill Cornell Medicine, New York, NY, USA
| | - Joseph D. Butner
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - Carmine Schiavone
- Mathematics in Medicine Program, Department of Medicine, Houston Methodist Research Institute, Houston, TX, USA
- Department of Chemical, Materials and Industrial Production Engineering, University of Naples Federico II, Naples, Italy
| | - Dan G. Duda
- Edwin. L. Steele Laboratories for Tumor Biology, Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA
| | - Ahmed O. Kaseb
- Department of Gastrointestinal Medical Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - Caroline Chung
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - Eugene J. Koay
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - Vittorio Cristini
- Mathematics in Medicine Program, Department of Medicine, Houston Methodist Research Institute, Houston, TX, USA
- Physiology, Biophysics, and Systems Biology Program, Graduate School of Medical Sciences, Weill Cornell Medicine, New York, NY, USA
- Neal Cancer Center, Houston Methodist Research Institute, Houston, TX, USA
- Department of Imaging Physics, University of Texas M.D. Anderson Cancer Center, Houston, TX, USA
| | - Bulent Ozpolat
- Department of Nanomedicine, Houston Methodist Research Institute, Houston, TX, USA
| | - George A. Calin
- Department of Translational Molecular Pathology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - Zhihui Wang
- Mathematics in Medicine Program, Department of Medicine, Houston Methodist Research Institute, Houston, TX, USA
- Department of Physiology and Biophysics, Weill Cornell Medical College, New York, NY, USA
- Neal Cancer Center, Houston Methodist Research Institute, Houston, TX, USA
| |
Collapse
|
3
|
Cave J, Shinglot V, Butner JD, Cristini V, Ozpolat B, Calin GA, Dogra P, Wang Z. Mechanistic modeling of anti-microRNA-155 therapy combinations in lung cancer. Annu Int Conf IEEE Eng Med Biol Soc 2023; 2023:1-4. [PMID: 38083518 DOI: 10.1109/embc40787.2023.10341114] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/18/2023]
Abstract
To improve treatment outcomes in non-small cell lung cancer (NSCLC), it is crucial to identify treatment strategies with the potential to exhibit drug synergism. This can lower the required effective dose, reducing exposure to drugs and associated toxicities, while improving treatment efficacy. In previous studies, drugs targeting the microRNA-155 or PD-L1 have been promising in restraining NSCLC tumor growth. We have developed a mathematical model that simulates the in vivo pharmacokinetics and pharmacodynamics of the novel nanoparticle-delivered anti-microRNA-155 for potential use with standard-of-care drug atezolizumab for NSCLC. Through modeling and simulation, we identified possible drug synergism between the two drugs that holds promise to improve tumor response at reduced drug exposure.Clinical Relevance-Identifying the possibility of drug synergism for an anti-microRNA-155 based nanotherapeutic with standard-of-care immunotherapy to improve lung cancer treatment outcomes.
Collapse
|
4
|
Syed M, Cagely M, Dogra P, Hollmer L, Butner JD, Cristini V, Koay EJ. Immune-checkpoint inhibitor therapy response evaluation using oncophysics-based mathematical models. Wiley Interdiscip Rev Nanomed Nanobiotechnol 2023; 15:e1855. [PMID: 36148978 DOI: 10.1002/wnan.1855] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/07/2021] [Revised: 06/10/2022] [Accepted: 08/23/2022] [Indexed: 11/08/2022]
Abstract
The field of oncology has transformed with the advent of immunotherapies. The standard of care for multiple cancers now includes novel drugs that target key checkpoints that function to modulate immune responses, enabling the patient's immune system to elicit an effective anti-tumor response. While these immune-based approaches can have dramatic effects in terms of significantly reducing tumor burden and prolonging survival for patients, the therapeutic approach remains active only in a minority of patients and is often not durable. Multiple biological investigations have identified key markers that predict response to the most common form of immunotherapy-immune checkpoint inhibitors (ICI). These biomarkers help enrich patients for ICI but are not 100% predictive. Understanding the complex interactions of these biomarkers with other pathways and factors that lead to ICI resistance remains a major goal. Principles of oncophysics-the idea that cancer can be described as a multiscale physical aberration-have shown promise in recent years in terms of capturing the essence of the complexities of ICI interactions. Here, we review the biological knowledge of mechanisms of ICI action and how these are incorporated into modern oncophysics-based mathematical models. Building on the success of oncophysics-based mathematical models may help to discover new, rational methods to engineer immunotherapy for patients in the future. This article is categorized under: Therapeutic Approaches and Drug Discovery > Nanomedicine for Oncologic Disease.
Collapse
Affiliation(s)
- Mustafa Syed
- Department of Gastrointestinal Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas, USA
| | - Matthew Cagely
- Department of Gastrointestinal Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas, USA
| | - Prashant Dogra
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, Texas, USA.,Department of Physiology and Biophysics, Weill Cornell Medical College, New York, New York, USA
| | - Lauren Hollmer
- Department of Gastrointestinal Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas, USA
| | - Joseph D Butner
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, Texas, USA
| | - Vittorio Cristini
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, Texas, USA.,Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, Houston, Texas, USA.,Physiology, Biophysics, and Systems Biology Program, Graduate School of Medical Sciences, Weill Cornell Medicine, New York, New York, USA
| | - Eugene J Koay
- Department of Gastrointestinal Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas, USA
| |
Collapse
|
5
|
Butner JD, Dogra P, Chung C, Pasqualini R, Arap W, Lowengrub J, Cristini V, Wang Z. Mathematical modeling of cancer immunotherapy for personalized clinical translation. Nat Comput Sci 2022; 2:785-796. [PMID: 38126024 PMCID: PMC10732566 DOI: 10.1038/s43588-022-00377-z] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/19/2022] [Accepted: 11/14/2022] [Indexed: 12/23/2023]
Abstract
Encouraging advances are being made in cancer immunotherapy modeling, especially in the key areas of developing personalized treatment strategies based on individual patient parameters, predicting treatment outcomes and optimizing immunotherapy synergy when used in combination with other treatment approaches. Here we present a focused review of the most recent mathematical modeling work on cancer immunotherapy with a focus on clinical translatability. It can be seen that this field is transitioning from pure basic science to applications that can make impactful differences in patients' lives. We discuss how researchers are integrating experimental and clinical data to fully inform models so that they can be applied for clinical predictions, and present the challenges that remain to be overcome if widespread clinical adaptation is to be realized. Lastly, we discuss the most promising future applications and areas that are expected to be the focus of extensive upcoming modeling studies.
Collapse
Affiliation(s)
- Joseph D. Butner
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX, USA
| | - Prashant Dogra
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX, USA
| | - Caroline Chung
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - Renata Pasqualini
- Rutgers Cancer Institute of New Jersey, Newark, NJ, USA
- Department of Radiation Oncology, Division of Cancer Biology, Rutgers New Jersey Medical School, Newark, NJ, USA
| | - Wadih Arap
- Rutgers Cancer Institute of New Jersey, Newark, NJ, USA
- Department of Medicine, Division of Hematology/Oncology, Rutgers New Jersey Medical School, Newark, NJ, USA
| | - John Lowengrub
- Department of Mathematics, University of California at Irvine, Irvine, CA, USA
| | - Vittorio Cristini
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX, USA
- Neal Cancer Center, Houston Methodist Research Institute, Houston, TX, USA
- Department of Imaging Physics, University of Texas MD Anderson Cancer Center, Houston, TX, USA
- Physiology, Biophysics, and Systems Biology Program, Graduate School of Medical Sciences, Weill Cornell Medicine, New York, NY, USA
| | - Zhihui Wang
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX, USA
- Neal Cancer Center, Houston Methodist Research Institute, Houston, TX, USA
- Department of Physiology and Biophysics, Weill Cornell Medicine, New York, NY, USA
| |
Collapse
|
6
|
Butner JD, Farhat M, Cristini V, Chung C, Wang Z. Protocol for mathematical prediction of patient response and survival to immune checkpoint inhibitor immunotherapy. STAR Protoc 2022; 3:101886. [PMID: 36595890 PMCID: PMC9719106 DOI: 10.1016/j.xpro.2022.101886] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2022] [Revised: 10/03/2022] [Accepted: 11/03/2022] [Indexed: 12/03/2022] Open
Abstract
This protocol describes the application of a mechanistic mathematical model of immune checkpoint inhibitor (ICI) immunotherapy to patient tumor imaging data for predicting solid tumor response and patient survival under ICI intervention. We describe steps for data collection and processing, data pipelines, and approaches to increase precision. The protocol is highly predictive as early as the first restaging after treatment start and can be used with standard-of-care imaging measures. For complete details on the use and execution of this protocol, please refer to Butner et al. (2020)1 and Butner et al. (2021).2.
Collapse
Affiliation(s)
- Joseph D. Butner
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA,Corresponding author
| | - Maguy Farhat
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA
| | - Vittorio Cristini
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA,Neal Cancer Center, Houston Methodist Research Institute, Houston, TX 77030, USA,Department of Imaging Physics, University of Texas MD Anderson Cancer Center, Houston, TX 77230, USA,Physiology, Biophysics, and Systems Biology Program, Graduate School of Medical Sciences, Weill Cornell Medicine, New York, NY 10065, USA
| | - Caroline Chung
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA,Corresponding author
| | - Zhihui Wang
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA,Neal Cancer Center, Houston Methodist Research Institute, Houston, TX 77030, USA,Department of Physiology and Biophysics, Weill Cornell Medicine, New York, NY 10065, USA,Department of Medical Education, Texas A&M University School of Medicine, Bryan, TX 77807, USA,Corresponding author
| |
Collapse
|
7
|
Sun K, Xu Y, Zhang L, Niravath P, Darcourt J, Patel T, Teh BS, Farach AM, Guerrero C, Mathur S, Sultenfuss MA, Gupta N, Schwartz MR, Haley SL, Nair S, Li X, Nguyen TTA, Butner JD, Ensor J, Mejia JA, Mei Z, Butler EB, Chen SH, Bernicker EH, Chang JC. A Phase 2 Trial of Enhancing Immune Checkpoint Blockade by Stereotactic Radiation and In Situ Virus Gene Therapy in Metastatic Triple-Negative Breast Cancer. Clin Cancer Res 2022; 28:4392-4401. [PMID: 35877117 PMCID: PMC9561553 DOI: 10.1158/1078-0432.ccr-22-0622] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2022] [Revised: 05/06/2022] [Accepted: 07/21/2022] [Indexed: 01/07/2023]
Abstract
PURPOSE A Phase 2 trial of stereotactic radiotherapy and in situ cytotoxic virus therapy in patients with metastatic triple-negative breast cancer (mTNBC) followed by pembrolizumab (STOMP) was designed to evaluate dual approach of enhancing single-agent immune checkpoint blockade with adenovirus-mediated expression of herpes-simplex-virus thymidine-kinase (ADV/HSV-tk) plus valacyclovir gene therapy and stereotactic body radiotherapy (SBRT) in patients with mTNBC. PATIENTS AND METHODS In this single-arm, open-label Phase 2 trial, patients with mTNBC were treated with ADV/HSV-tk [5 × 1011 virus particles (vp)] intratumoral injection, followed by SBRT to the injected tumor site, then pembrolizumab (200 mg, every 3 weeks). The primary endpoint was clinical benefit rate [CBR; complete response (CR), partial response (PR), or stable disease (SD) ≥ 24 weeks per RECIST version1.1 at non-irradiated site]. Secondary endpoints included duration on treatment (DoT), overall survival (OS), and safety. Exploratory endpoints included immune response to treatment assessed by correlative tissue and blood-based biomarkers. RESULTS Twenty-eight patients were enrolled and treated. CBR was seen in 6 patients (21.4%), including 2 CR (7.1%), 1 PR (3.6%), and 3 SD (10.7%). Patients with clinical benefit had durable responses, with median DoT of 9.6 months and OS of 14.7 months. The median OS was 6.6 months in the total population. The combination was well tolerated. Correlative studies with Cytometry by Time of Flight (CyTOF) and imaging mass cytometry (IMC) revealed a significant increase of CD8 T cells in responders and of myeloid cells in non-responders. CONCLUSIONS The median OS increased by more than 2-fold in patients with clinical benefit. The therapy is a well-tolerated treatment in heavily pretreated patients with mTNBC. Early detection of increased effector and effector memory CD8 T cells and myeloids correlate with response and non-response, respectively.
Collapse
Affiliation(s)
- Kai Sun
- Houston Methodist Neal Cancer Center, Houston, Texas
| | - Yitian Xu
- Houston Methodist Research Institute, Center for Immunotherapy Research, Houston, Texas
| | - Licheng Zhang
- Houston Methodist Research Institute, Center for Immunotherapy Research, Houston, Texas
| | | | | | - Tejal Patel
- Houston Methodist Neal Cancer Center, Houston, Texas
| | - Bin S. Teh
- Department of Radiation Oncology, Houston Methodist Hospital, Houston, Texas
| | - Andrew M. Farach
- Department of Radiation Oncology, Houston Methodist Hospital, Houston, Texas
| | | | - Sunil Mathur
- Houston Methodist Neal Cancer Center, Houston, Texas
| | | | - Nakul Gupta
- Department of Radiology, Houston Methodist Hospital, Houston, Texas
| | - Mary R. Schwartz
- Department of Pathology and Genomic Medicine, Houston Methodist Hospital, Houston, Texas
| | - Susan L. Haley
- Department of Pathology and Genomic Medicine, Houston Methodist Hospital, Houston, Texas
| | - Sindhu Nair
- Houston Methodist Neal Cancer Center, Houston, Texas
| | - Xiaoxian Li
- Department of Pathology and Laboratory Medicine, Emory University, Atlanta, Georgia
| | - Thi Truc Anh Nguyen
- Department of Pathology and Laboratory Medicine, Emory University, Atlanta, Georgia
| | - Joseph D. Butner
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, Texas
| | - Joe Ensor
- Houston Methodist Neal Cancer Center, Houston, Texas
| | | | - Zhuyong Mei
- Center for Cell and Gene Therapy, Baylor College of Medicine, Houston, Texas
| | - E. Brian Butler
- Department of Radiation Oncology, Houston Methodist Hospital, Houston, Texas
| | - Shu-hsia Chen
- Houston Methodist Research Institute, Center for Immunotherapy Research, Houston, Texas
| | | | - Jenny C. Chang
- Houston Methodist Neal Cancer Center, Houston, Texas.,Corresponding Author: Jenny C. Chang, Houston Methodist Research Institute, 6445 Main Street, Floor 24, Houston, TX 77030. Phone: 713-441-9948; Fax: 713-441-8791; E-mail:
| |
Collapse
|
8
|
Wang CX, Elganainy D, Zaid MM, Butner JD, Agrawal A, Nizzero S, Minsky BD, Holliday EB, Taniguchi CM, Smith GL, Koong AC, Herman JM, Das P, Maitra A, Wang H, Wolff RA, Katz MHG, Crane CH, Cristini V, Koay EJ. Mass Transport Model of Radiation Response: Calibration and Application to Chemoradiation for Pancreatic Cancer. Int J Radiat Oncol Biol Phys 2022; 114:163-172. [PMID: 35643254 PMCID: PMC10042520 DOI: 10.1016/j.ijrobp.2022.04.044] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2021] [Revised: 03/22/2022] [Accepted: 04/28/2022] [Indexed: 11/25/2022]
Abstract
PURPOSE The benefit of radiation therapy for pancreatic ductal adenocarcinoma (PDAC) remains unclear. We hypothesized that a new mechanistic mathematical model of chemotherapy and radiation response could predict clinical outcomes a priori, using a previously described baseline measurement of perfusion from computed tomography scans, normalized area under the enhancement curve (nAUC). METHODS AND MATERIALS We simplified an existing mass transport model that predicted cancer cell death by replacing previously unknown variables with averaged direct measurements from randomly selected pathologic sections of untreated PDAC. This allowed using nAUC as the sole model input to approximate tumor perfusion. We then compared the predicted cancer cell death to the actual cell death measured from corresponding resected tumors treated with neoadjuvant chemoradiation in a calibration cohort (n = 80) and prospective cohort (n = 25). After calibration, we applied the model to 2 separate cohorts for pathologic and clinical associations: targeted therapy cohort (n = 101), cetuximab/bevacizumab + radiosensitizing chemotherapy, and standard chemoradiation cohort (n = 81), radiosensitizing chemotherapy to 50.4 Gy in 28 fractions. RESULTS We established the relationship between pretreatment computed v nAUC to pathologically verified blood volume fraction of the tumor (r = 0.65; P = .009) and fractional tumor cell death (r = 0.97-0.99; P < .0001) in the calibration and prospective cohorts. On multivariate analyses, accounting for traditional covariates, nAUC independently associated with overall survival in all cohorts (mean hazard ratios, 0.14-0.31). Receiver operator characteristic analyses revealed discrimination of good and bad prognostic groups in the cohorts with area under the curve values of 0.64 to 0.71. CONCLUSIONS This work presents a new mathematical modeling approach to predict clinical response from chemotherapy and radiation for PDAC. Our findings indicate that oxygen/drug diffusion strongly influences clinical responses and that nAUC is a potential tool to select patients with PDAC for radiation therapy.
Collapse
Affiliation(s)
- Charles X Wang
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, Texas; Department of Radiation Oncology, University of California Davis Medical Center, Sacramento, California
| | - Dalia Elganainy
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Mohamed M Zaid
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Joseph D Butner
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, Texas
| | - Anshuman Agrawal
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Sara Nizzero
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, Texas
| | - Bruce D Minsky
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Emma B Holliday
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Cullen M Taniguchi
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Grace L Smith
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Albert C Koong
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Joseph M Herman
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Prajnan Das
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, Texas
| | | | | | | | - Matthew H G Katz
- Surgical Oncology, University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Christopher H Crane
- Department of Radiation Oncology, Memorial Sloan Kettering Cancer Center, New York, New York
| | - Vittorio Cristini
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, Texas; Department of Imaging Physics, University of Texas MD Anderson Cancer Center, Houston, Texas; Physiology, Biophysics, and Systems Biology Program, Graduate School of Medical Sciences, Weill Cornell Medicine, New York, New York
| | - Eugene J Koay
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, Texas.
| |
Collapse
|
9
|
Dogra P, Ramírez JR, Butner JD, Peláez MJ, Chung C, Hooda-Nehra A, Pasqualini R, Arap W, Cristini V, Calin GA, Ozpolat B, Wang Z. Translational Modeling Identifies Synergy between Nanoparticle-Delivered miRNA-22 and Standard-of-Care Drugs in Triple-Negative Breast Cancer. Pharm Res 2022; 39:511-528. [PMID: 35294699 PMCID: PMC8986735 DOI: 10.1007/s11095-022-03176-3] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2021] [Accepted: 01/21/2022] [Indexed: 12/29/2022]
Abstract
Purpose Downregulation of miRNA-22 in triple-negative breast cancer (TNBC) is associated with upregulation of eukaryotic elongation 2 factor kinase (eEF2K) protein, which regulates tumor growth, chemoresistance, and tumor immunosurveillance. Moreover, exogenous administration of miRNA-22, loaded in nanoparticles to prevent degradation and improve tumor delivery (termed miRNA-22 nanotherapy), to suppress eEF2K production has shown potential as an investigational therapeutic agent in vivo. Methods To evaluate the translational potential of miRNA-22 nanotherapy, we developed a multiscale mechanistic model, calibrated to published in vivo data and extrapolated to the human scale, to describe and quantify the pharmacokinetics and pharmacodynamics of miRNA-22 in virtual patient populations. Results Our analysis revealed the dose-response relationship, suggested optimal treatment frequency for miRNA-22 nanotherapy, and highlighted key determinants of therapy response, from which combination with immune checkpoint inhibitors was identified as a candidate strategy for improving treatment outcomes. More importantly, drug synergy was identified between miRNA-22 and standard-of-care drugs against TNBC, providing a basis for rational therapeutic combinations for improved response Conclusions The present study highlights the translational potential of miRNA-22 nanotherapy for TNBC in combination with standard-of-care drugs. Supplementary Information The online version contains supplementary material available at 10.1007/s11095-022-03176-3.
Collapse
Affiliation(s)
- Prashant Dogra
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, Texas, 77030, USA
- Department of Physiology and Biophysics, Weill Cornell Medical College, New York, New York, 10065, USA
| | - Javier Ruiz Ramírez
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, Texas, 77030, USA
| | - Joseph D Butner
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, Texas, 77030, USA
| | - Maria J Peláez
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, Texas, 77030, USA
| | - Caroline Chung
- Department of Radiation Oncology, The University of Texas M.D. Anderson Cancer Center, Houston, Texas, 77030, USA
| | - Anupama Hooda-Nehra
- Rutgers Cancer Institute of New Jersey, Newark, New Jersey, 07101, USA
- Department of Medicine, Division of Hematology/Oncology, Rutgers New Jersey Medical School, Newark, New Jersey, 07103, USA
| | - Renata Pasqualini
- Rutgers Cancer Institute of New Jersey, Newark, New Jersey, 07101, USA
- Department of Radiation Oncology, Division of Cancer Biology, Rutgers New Jersey Medical School, Newark, New Jersey, 07103, USA
| | - Wadih Arap
- Rutgers Cancer Institute of New Jersey, Newark, New Jersey, 07101, USA
- Department of Medicine, Division of Hematology/Oncology, Rutgers New Jersey Medical School, Newark, New Jersey, 07103, USA
| | - Vittorio Cristini
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, Texas, 77030, USA
- Department of Imaging Physics, The University of Texas M.D. Anderson Cancer Center, Houston, Texas, 77230, USA
- Physiology, Biophysics, and Systems Biology Program, Graduate School of Medical Sciences, Weill Cornell Medicine, New York, New York, 10065, USA
| | - George A Calin
- Department of Translational Molecular Pathology, The University of Texas M.D. Anderson Cancer Center, Houston, Texas, 77030, USA
| | - Bulent Ozpolat
- Department of Experimental Therapeutics, The University of Texas M.D. Anderson Cancer Center, Houston, Texas, 77030, USA
| | - Zhihui Wang
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, Texas, 77030, USA.
- Department of Physiology and Biophysics, Weill Cornell Medical College, New York, New York, 10065, USA.
| |
Collapse
|
10
|
Butner JD, Martin GV, Wang Z, Corradetti B, Ferrari M, Esnaola N, Chung C, Hong DS, Welsh JW, Hasegawa N, Mittendorf EA, Curley SA, Chen SH, Pan PY, Libutti SK, Ganesan S, Sidman RL, Pasqualini R, Arap W, Koay EJ, Cristini V. Early prediction of clinical response to checkpoint inhibitor therapy in human solid tumors through mathematical modeling. eLife 2021; 10:70130. [PMID: 34749885 PMCID: PMC8629426 DOI: 10.7554/elife.70130] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2021] [Accepted: 10/25/2021] [Indexed: 11/13/2022] Open
Abstract
Background: Checkpoint inhibitor therapy of cancer has led to markedly improved survival of a subset of patients in multiple solid malignant tumor types, yet the factors driving these clinical responses or lack thereof are not known. We have developed a mechanistic mathematical model for better understanding these factors and their relations in order to predict treatment outcome and optimize personal treatment strategies. Methods: Here, we present a translational mathematical model dependent on three key parameters for describing efficacy of checkpoint inhibitors in human cancer: tumor growth rate (α), tumor-immune infiltration (Λ), and immunotherapy-mediated amplification of anti-tumor response (µ). The model was calibrated by fitting it to a compiled clinical tumor response dataset (n = 189 patients) obtained from published anti-PD-1 and anti-PD-L1 clinical trials, and then validated on an additional validation cohort (n = 64 patients) obtained from our in-house clinical trials. Results: The derived parameters Λ and µ were both significantly different between responding versus nonresponding patients. Of note, our model appropriately classified response in 81.4% of patients by using only tumor volume measurements and within 2 months of treatment initiation in a retrospective analysis. The model reliably predicted clinical response to the PD-1/PD-L1 class of checkpoint inhibitors across multiple solid malignant tumor types. Comparison of model parameters to immunohistochemical measurement of PD-L1 and CD8+ T cells confirmed robust relationships between model parameters and their underlying biology. Conclusions: These results have demonstrated reliable methods to inform model parameters directly from biopsy samples, which are conveniently obtainable as early as the start of treatment. Together, these suggest that the model parameters may serve as early and robust biomarkers of the efficacy of checkpoint inhibitor therapy on an individualized per-patient basis. Funding: We gratefully acknowledge support from the Andrew Sabin Family Fellowship, Center for Radiation Oncology Research, Sheikh Ahmed Center for Pancreatic Cancer Research, GE Healthcare, Philips Healthcare, and institutional funds from the University of Texas M.D. Anderson Cancer Center. We have also received Cancer Center Support Grants from the National Cancer Institute (P30CA016672 to the University of Texas M.D. Anderson Cancer Center and P30CA072720 the Rutgers Cancer Institute of New Jersey). This research has also been supported in part by grants from the National Science Foundation Grant DMS-1930583 (ZW, VC), the National Institutes of Health (NIH) 1R01CA253865 (ZW, VC), 1U01CA196403 (ZW, VC), 1U01CA213759 (ZW, VC), 1R01CA226537 (ZW, RP, WA, VC), 1R01CA222007 (ZW, VC), U54CA210181 (ZW, VC), and the University of Texas System STARS Award (VC). BC acknowledges support through the SER Cymru II Programme, funded by the European Commission through the Horizon 2020 Marie Skłodowska-Curie Actions (MSCA) COFUND scheme and the Welsh European Funding Office (WEFO) under the European Regional Development Fund (ERDF). EK has also received support from the Project Purple, NIH (U54CA210181, U01CA200468, and U01CA196403), and the Pancreatic Cancer Action Network (16-65-SING). MF was supported through NIH/NCI center grant U54CA210181, R01CA222959, DoD Breast Cancer Research Breakthrough Level IV Award W81XWH-17-1-0389, and the Ernest Cockrell Jr. Presidential Distinguished Chair at Houston Methodist Research Institute. RP and WA received serial research awards from AngelWorks, the Gillson-Longenbaugh Foundation, and the Marcus Foundation. This work was also supported in part by grants from the National Cancer Institute to SHC (R01CA109322, R01CA127483, R01CA208703, and U54CA210181 CITO pilot grant) and to PYP (R01CA140243, R01CA188610, and U54CA210181 CITO pilot grant). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Collapse
Affiliation(s)
- Joseph D Butner
- The Houston Methodist Research Institute, Houston, United States
| | - Geoffrey V Martin
- The University of Texas MD Anderson Cancer Center, Houston, United States
| | - Zhihui Wang
- The Houston Methodist Research Institute, Houston, United States
| | - Bruna Corradetti
- The Houston Methodist Research Institute, Houston, United States
| | - Mauro Ferrari
- The Houston Methodist Research Institute, Houston, United States
| | - Nestor Esnaola
- The Houston Methodist Research Institute, Houston, United States
| | - Caroline Chung
- The University of Texas MD Anderson Cancer Center, Houston, United States
| | - David S Hong
- The University of Texas MD Anderson Cancer Center, Houston, United States
| | - James W Welsh
- The Houston Methodist Research Institute, Houston, United States
| | - Naomi Hasegawa
- University of Texas Health Science Center, Houston, United States
| | | | | | - Shu-Hsia Chen
- The Houston Methodist Research Institute, Houston, United States
| | - Ping-Ying Pan
- The Houston Methodist Research Institute, Houston, United States
| | | | | | - Richard L Sidman
- Department of Neurology, Harvard Medical School, Boston, United States
| | | | - Wadih Arap
- Hematology and Oncology, Rutgers Cancer Institute of New Jersey, Newark, United States
| | - Eugene J Koay
- University of Texas MD Anderson Cancer Center, Houston, United States
| | | |
Collapse
|
11
|
Dogra P, Ramirez JR, Butner JD, Pelaez MJ, Cristini V, Wang Z. A Multiscale Model to Identify Limiting Factors in Nanoparticle-Based miRNA Delivery for Tumor Inhibition . Annu Int Conf IEEE Eng Med Biol Soc 2021; 2021:4230-4233. [PMID: 34892157 PMCID: PMC8712117 DOI: 10.1109/embc46164.2021.9630862] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
MicroRNA-based gene therapy for cancer treatment via nanoparticles (NPs) requires navigation of multiple physical and physiological barriers in order to efficiently deliver the miRNAs to the cancer cell cytoplasm. We here present a mathematical model to investigate the variability associated with tumor, NP, and miRNA characteristics, and identify the limiting factors in miRNA delivery to tumors. Through global parameter analysis, the miRNA release rate from NPs and NP degradability were found to have the most significant impact on cytosolic accumulation of miRNAs. These NP properties can be fine-tuned in order to optimize the delivery system for enhancing the efficacy of miRNA-based therapy.Clinical Relevance-Understanding the effect of nanoparticle, tumor, and miRNA characteristics in governing the efficacy of miRNA-based cancer therapy will support its clinical translation.
Collapse
|
12
|
Affiliation(s)
- Joseph D Butner
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA
| | - Zhihui Wang
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA.,Cancer Center, Houston Methodist Research Institute, Houston, TX 77030, USA.,Department of Physiology & Biophysics, Weill Cornell Medicine, New York, NY 10065, USA
| |
Collapse
|
13
|
Dogra P, Ruiz-Ramírez J, Sinha K, Butner JD, Peláez MJ, Rawat M, Yellepeddi VK, Pasqualini R, Arap W, Sostman HD, Cristini V, Wang Z. Innate Immunity Plays a Key Role in Controlling Viral Load in COVID-19: Mechanistic Insights from a Whole-Body Infection Dynamics Model. ACS Pharmacol Transl Sci 2021; 4:248-265. [PMID: 33615177 PMCID: PMC7805603 DOI: 10.1021/acsptsci.0c00183] [Citation(s) in RCA: 27] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2020] [Indexed: 12/18/2022]
Abstract
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a pathogen of immense public health concern. Efforts to control the disease have only proven mildly successful, and the disease will likely continue to cause excessive fatalities until effective preventative measures (such as a vaccine) are developed. To develop disease management strategies, a better understanding of SARS-CoV-2 pathogenesis and population susceptibility to infection are needed. To this end, mathematical modeling can provide a robust in silico tool to understand COVID-19 pathophysiology and the in vivo dynamics of SARS-CoV-2. Guided by ACE2-tropism (ACE2 receptor dependency for infection) of the virus and by incorporating cellular-scale viral dynamics and innate and adaptive immune responses, we have developed a multiscale mechanistic model for simulating the time-dependent evolution of viral load distribution in susceptible organs of the body (respiratory tract, gut, liver, spleen, heart, kidneys, and brain). Following parameter quantification with in vivo and clinical data, we used the model to simulate viral load progression in a virtual patient with varying degrees of compromised immune status. Further, we ranked model parameters through sensitivity analysis for their significance in governing clearance of viral load to understand the effects of physiological factors and underlying conditions on viral load dynamics. Antiviral drug therapy, interferon therapy, and their combination were simulated to study the effects on viral load kinetics of SARS-CoV-2. The model revealed the dominant role of innate immunity (specifically interferons and resident macrophages) in controlling viral load, and the importance of timing when initiating therapy after infection.
Collapse
Affiliation(s)
- Prashant Dogra
- Mathematics
in Medicine Program, Houston Methodist Research
Institute, Houston, Texas 77030, United States
| | - Javier Ruiz-Ramírez
- Mathematics
in Medicine Program, Houston Methodist Research
Institute, Houston, Texas 77030, United States
| | - Kavya Sinha
- DeBakey
Heart and Vascular Center, Houston Methodist
Hospital, Houston, Texas 77030, United States
| | - Joseph D. Butner
- Mathematics
in Medicine Program, Houston Methodist Research
Institute, Houston, Texas 77030, United States
| | - Maria J. Peláez
- Mathematics
in Medicine Program, Houston Methodist Research
Institute, Houston, Texas 77030, United States
| | - Manmeet Rawat
- Department
of Internal Medicine, University of New
Mexico School of Medicine, Albuquerque, New Mexico 87131, United States
| | - Venkata K. Yellepeddi
- Division
of Clinical Pharmacology, Department of Pediatrics, School of Medicine, University of Utah, Salt Lake City, Utah 84132, United States
- Department
of Pharmaceutics and Pharmaceutical Chemistry, College of Pharmacy, University of Utah, Salt Lake City, Utah 84112, United States
| | - Renata Pasqualini
- Rutgers
Cancer Institute of New Jersey, Newark, New Jersey 07101, United States
- Department
of Radiation Oncology, Division of Cancer Biology, Rutgers New Jersey Medical School, Newark, New Jersey 07103, United States
| | - Wadih Arap
- Rutgers
Cancer Institute of New Jersey, Newark, New Jersey 07101, United States
- Department
of Medicine, Division of Hematology/Oncology, Rutgers New Jersey Medical School, Newark, New Jersey 07103, United States
| | - H. Dirk Sostman
- Weill
Cornell Medicine, New York, New York 10065, United States
- Houston
Methodist Research Institute, Houston, Texas 77030, United States
- Houston
Methodist Academic Institute, Houston, Texas 77030, United States
| | - Vittorio Cristini
- Mathematics
in Medicine Program, Houston Methodist Research
Institute, Houston, Texas 77030, United States
- Weill
Cornell Medicine, New York, New York 10065, United States
| | - Zhihui Wang
- Mathematics
in Medicine Program, Houston Methodist Research
Institute, Houston, Texas 77030, United States
- Weill
Cornell Medicine, New York, New York 10065, United States
| |
Collapse
|
14
|
Anaya DA, Dogra P, Wang Z, Haider M, Ehab J, Jeong DK, Ghayouri M, Lauwers GY, Thomas K, Kim R, Butner JD, Nizzero S, Ramírez JR, Plodinec M, Sidman RL, Cavenee WK, Pasqualini R, Arap W, Fleming JB, Cristini V. A Mathematical Model to Estimate Chemotherapy Concentration at the Tumor-Site and Predict Therapy Response in Colorectal Cancer Patients with Liver Metastases. Cancers (Basel) 2021; 13:cancers13030444. [PMID: 33503971 PMCID: PMC7866038 DOI: 10.3390/cancers13030444] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2021] [Accepted: 01/21/2021] [Indexed: 12/22/2022] Open
Abstract
Simple Summary It is known that drug transport barriers in the tumor determine drug concentration at the tumor site, causing disparity from the systemic (plasma) drug concentration. However, current clinical standard of care still bases dosage and treatment optimization on the systemic concentration of drugs. Here, we present a proof of concept observational cohort study to accurately estimate drug concentration at the tumor site from mathematical modeling using biologic, clinical, and imaging/perfusion data, and correlate it with outcome in colorectal cancer liver metastases. We demonstrate that drug concentration at the tumor site, not in systemic circulation, can be used as a credible biomarker for predicting chemotherapy outcome, and thus our mathematical modeling approach can be applied prospectively in the clinic to personalize treatment design to optimize outcome. Abstract Chemotherapy remains a primary treatment for metastatic cancer, with tumor response being the benchmark outcome marker. However, therapeutic response in cancer is unpredictable due to heterogeneity in drug delivery from systemic circulation to solid tumors. In this proof-of-concept study, we evaluated chemotherapy concentration at the tumor-site and its association with therapy response by applying a mathematical model. By using pre-treatment imaging, clinical and biologic variables, and chemotherapy regimen to inform the model, we estimated tumor-site chemotherapy concentration in patients with colorectal cancer liver metastases, who received treatment prior to surgical hepatic resection with curative-intent. The differential response to therapy in resected specimens, measured with the gold-standard Tumor Regression Grade (TRG; from 1, complete response to 5, no response) was examined, relative to the model predicted systemic and tumor-site chemotherapy concentrations. We found that the average calculated plasma concentration of the cytotoxic drug was essentially equivalent across patients exhibiting different TRGs, while the estimated tumor-site chemotherapeutic concentration (eTSCC) showed a quadratic decline from TRG = 1 to TRG = 5 (p < 0.001). The eTSCC was significantly lower than the observed plasma concentration and dropped by a factor of ~5 between patients with complete response (TRG = 1) and those with no response (TRG = 5), while the plasma concentration remained stable across TRG groups. TRG variations were driven and predicted by differences in tumor perfusion and eTSCC. If confirmed in carefully planned prospective studies, these findings will form the basis of a paradigm shift in the care of patients with potentially curable colorectal cancer and liver metastases.
Collapse
Affiliation(s)
- Daniel A. Anaya
- Department of Gastrointestinal Oncology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL 33612, USA; (M.H.); (J.E.); (R.K.); (J.B.F.)
- Correspondence: (D.A.A.); (V.C.); Tel.: +1-813-745-1432 (D.A.A.); +1-505-934-1813 (V.C.); Fax: +1-813-745-7229 (D.A.A.)
| | - Prashant Dogra
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA; (P.D.); (Z.W.); (J.D.B.); (S.N.); (J.R.R.)
| | - Zhihui Wang
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA; (P.D.); (Z.W.); (J.D.B.); (S.N.); (J.R.R.)
| | - Mintallah Haider
- Department of Gastrointestinal Oncology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL 33612, USA; (M.H.); (J.E.); (R.K.); (J.B.F.)
| | - Jasmina Ehab
- Department of Gastrointestinal Oncology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL 33612, USA; (M.H.); (J.E.); (R.K.); (J.B.F.)
| | - Daniel K. Jeong
- Department of Diagnostic Imaging and Interventional Radiology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL 33612, USA; (D.K.J.); (M.G.); (G.Y.L.); (K.T.)
| | - Masoumeh Ghayouri
- Department of Diagnostic Imaging and Interventional Radiology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL 33612, USA; (D.K.J.); (M.G.); (G.Y.L.); (K.T.)
| | - Gregory Y. Lauwers
- Department of Diagnostic Imaging and Interventional Radiology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL 33612, USA; (D.K.J.); (M.G.); (G.Y.L.); (K.T.)
| | - Kerry Thomas
- Department of Diagnostic Imaging and Interventional Radiology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL 33612, USA; (D.K.J.); (M.G.); (G.Y.L.); (K.T.)
| | - Richard Kim
- Department of Gastrointestinal Oncology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL 33612, USA; (M.H.); (J.E.); (R.K.); (J.B.F.)
| | - Joseph D. Butner
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA; (P.D.); (Z.W.); (J.D.B.); (S.N.); (J.R.R.)
| | - Sara Nizzero
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA; (P.D.); (Z.W.); (J.D.B.); (S.N.); (J.R.R.)
| | - Javier Ruiz Ramírez
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA; (P.D.); (Z.W.); (J.D.B.); (S.N.); (J.R.R.)
| | - Marija Plodinec
- Biozentrum and the Swiss Nanoscience Institute & ARTIDIS AG, University of Basel, 4056 Basel, Switzerland;
| | - Richard L. Sidman
- Department of Neurology, Harvard Medical School, Boston, MA 02115, USA;
| | - Webster K. Cavenee
- Ludwig Institute for Cancer Research, University of California-San Diego, La Jolla, CA 92093, USA;
| | - Renata Pasqualini
- Rutgers Cancer Institute of New Jersey & Division of Cancer Biology, Department of Radiation Oncology, Rutgers New Jersey Medical School, Newark, NJ 07103, USA;
| | - Wadih Arap
- Rutgers Cancer Institute of New Jersey & Division of Hematology/Oncology, Department of Medicine Rutgers New Jersey Medical School, Newark, NJ 07103, USA;
| | - Jason B. Fleming
- Department of Gastrointestinal Oncology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL 33612, USA; (M.H.); (J.E.); (R.K.); (J.B.F.)
| | - Vittorio Cristini
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA; (P.D.); (Z.W.); (J.D.B.); (S.N.); (J.R.R.)
- Correspondence: (D.A.A.); (V.C.); Tel.: +1-813-745-1432 (D.A.A.); +1-505-934-1813 (V.C.); Fax: +1-813-745-7229 (D.A.A.)
| |
Collapse
|
15
|
Butner JD, Fuentes D, Ozpolat B, Calin GA, Zhou X, Lowengrub J, Cristini V, Wang Z. Corrections to "A multiscale Agent-Based Model of Ductal Carcinoma in Situ". IEEE Trans Biomed Eng 2021; 68:736. [PMID: 33476261 DOI: 10.1109/tbme.2020.3044508] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
16
|
Butner JD, Wang Z, Elganainy D, Al Feghali KA, Plodinec M, Calin GA, Dogra P, Nizzero S, Ruiz-Ramírez J, Martin GV, Tawbi HA, Chung C, Koay EJ, Welsh JW, Hong DS, Cristini V. A mathematical model for the quantification of a patient's sensitivity to checkpoint inhibitors and long-term tumour burden. Nat Biomed Eng 2021; 5:297-308. [PMID: 33398132 PMCID: PMC8669771 DOI: 10.1038/s41551-020-00662-0] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2019] [Accepted: 11/14/2020] [Indexed: 02/06/2023]
Abstract
A large proportion of patients with cancer are unresponsive to treatment with immune checkpoint blockade and other immunotherapies. Here, we report a mathematical model of the time-course of tumour responses to immune-checkpoint inhibitors. The model takes into account intrinsic tumour-growth rates, the rates of immune activation and of tumour–immune-cell interactions, and the efficacy of immune-mediated tumour killing. For 124 patients, four cancer types and two immunotherapy agents, the model reliably described the immune responses and final tumour burden across all different cancers and drug combinations examined. In validation cohorts from four clinical trials of checkpoint inhibitors (with a total of 177 patients), the model accurately stratified the patients according to reduced or increased long-term tumour burden. We also provide model-derived quantitative measures of treatment sensitivity for specific drug–cancer combinations. The model can be used to predict responses to therapy and to quantify specific drug–cancer sensitivities in individual patients.
Collapse
Affiliation(s)
- Joseph D Butner
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX, USA
| | - Zhihui Wang
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX, USA. .,Department of Imaging Physics, University of Texas MD Anderson Cancer Center, Houston, TX, USA.
| | - Dalia Elganainy
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - Karine A Al Feghali
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - Marija Plodinec
- Biozentrum and the Swiss Nanoscience Institute, University of Basel, Basel, Switzerland
| | - George A Calin
- Department of Translational Molecular Pathology, University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - Prashant Dogra
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX, USA
| | - Sara Nizzero
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX, USA
| | - Javier Ruiz-Ramírez
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX, USA
| | - Geoffrey V Martin
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - Hussein A Tawbi
- Department of Melanoma Medical Oncology, University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - Caroline Chung
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - Eugene J Koay
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - James W Welsh
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - David S Hong
- Department of Investigational Cancer Therapeutics, University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - Vittorio Cristini
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX, USA. .,Department of Imaging Physics, University of Texas MD Anderson Cancer Center, Houston, TX, USA. .,Physiology, Biophysics, and Systems Biology Program, Graduate School of Medical Sciences, Weill Cornell Medicine, New York, NY, USA.
| |
Collapse
|
17
|
Dogra P, Ruiz-Ramírez J, Sinha K, Butner JD, Peláez MJ, Rawat M, Yellepeddi VK, Pasqualini R, Arap W, Sostman HD, Cristini V, Wang Z. Innate immunity plays a key role in controlling viral load in COVID-19: mechanistic insights from a whole-body infection dynamics model. medRxiv 2020:2020.10.30.20215335. [PMID: 33173913 PMCID: PMC7654909 DOI: 10.1101/2020.10.30.20215335] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/29/2023]
Abstract
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a pathogen of immense public health concern. Efforts to control the disease have only proven mildly successful, and the disease will likely continue to cause excessive fatalities until effective preventative measures (such as a vaccine) are developed. To develop disease management strategies, a better understanding of SARS-CoV-2 pathogenesis and population susceptibility to infection are needed. To this end, physiologically-relevant mathematical modeling can provide a robust in silico tool to understand COVID-19 pathophysiology and the in vivo dynamics of SARS-CoV-2. Guided by ACE2-tropism (ACE2 receptor dependency for infection) of the virus, and by incorporating cellular-scale viral dynamics and innate and adaptive immune responses, we have developed a multiscale mechanistic model for simulating the time-dependent evolution of viral load distribution in susceptible organs of the body (respiratory tract, gut, liver, spleen, heart, kidneys, and brain). Following calibration with in vivo and clinical data, we used the model to simulate viral load progression in a virtual patient with varying degrees of compromised immune status. Further, we conducted global sensitivity analysis of model parameters and ranked them for their significance in governing clearance of viral load to understand the effects of physiological factors and underlying conditions on viral load dynamics. Antiviral drug therapy, interferon therapy, and their combination was simulated to study the effects on viral load kinetics of SARS-CoV-2. The model revealed the dominant role of innate immunity (specifically interferons and resident macrophages) in controlling viral load, and the importance of timing when initiating therapy following infection.
Collapse
Affiliation(s)
- Prashant Dogra
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA
| | - Javier Ruiz-Ramírez
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA
| | - Kavya Sinha
- DeBakey Heart and Vascular Center, Houston Methodist Hospital, Houston, TX 77030, USA
| | - Joseph D. Butner
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA
| | - Maria J Peláez
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA
| | - Manmeet Rawat
- Department of Internal Medicine, University of New Mexico School of Medicine, Albuquerque, NM 87131, USA
| | - Venkata K. Yellepeddi
- Division of Clinical Pharmacology, Department of Pediatrics, School of Medicine, University of Utah, Salt Lake City, UT 84132, USA
- Department of Pharmaceutics and Pharmaceutical Chemistry, College of Pharmacy, University of Utah, Salt Lake City, UT 84112, USA
| | - Renata Pasqualini
- Rutgers Cancer Institute of New Jersey, Newark, NJ, 07101, USA
- Department of Radiation Oncology, Division of Cancer Biology, Rutgers Cancer Institute of New Jersey, Rutgers New Jersey Medical School, Newark, NJ, 07103, USA
| | - Wadih Arap
- Rutgers Cancer Institute of New Jersey, Newark, NJ, 07101, USA
- Department of Medicine, Division of Hematology/Oncology, Rutgers New Jersey Medical School, Newark, NJ, 07103, USA
| | - H. Dirk Sostman
- Weill Cornell Medicine, New York, NY 10065, USA
- Houston Methodist Research Institute, Houston, TX 77030, USA
- Houston Methodist Academic Institute, Houston, TX 77030, USA
| | - Vittorio Cristini
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA
| | - Zhihui Wang
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA
| |
Collapse
|
18
|
Dogra P, Butner JD, Ramirez JR, Cristini V, Wang Z. Investigating the Effect of Aging on the Pharmacokinetics and Tumor Delivery of Nanomaterials using Mathematical Modeling. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2020; 2020:2447-2450. [PMID: 33018501 DOI: 10.1109/embc44109.2020.9175322] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
The application of nanomedicine for diagnosis and treatment of cancer has immense potential, but has witnessed only limited clinical success, in part due to insufficient understanding of the role of nanomaterial properties and physiological variables in governing nanoparticle (NP) pharmacology. Here, we present a multiscale mathematical model to examine the effects of physiological changes associated with patient age on the pharmacokinetics and tumor delivery efficiency of NPs. We show that physiological changes due to aging prolong the residence of NPs in the systemic circulation, thereby improving passive accumulation of NPs in tumors.Clinical Relevance - Understanding the effect of inter-individual variability on the pharmacological behavior of nanomaterials will improve their clinical translatability.
Collapse
|
19
|
Butner JD, Elganainy D, Wang CX, Wang Z, Chen SH, Esnaola NF, Pasqualini R, Arap W, Hong DS, Welsh J, Koay EJ, Cristini V. Mathematical prediction of clinical outcomes in advanced cancer patients treated with checkpoint inhibitor immunotherapy. Sci Adv 2020; 6:eaay6298. [PMID: 32426472 PMCID: PMC7190324 DOI: 10.1126/sciadv.aay6298] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/04/2019] [Accepted: 01/21/2020] [Indexed: 05/20/2023]
Abstract
We present a mechanistic mathematical model of immune checkpoint inhibitor therapy to address the oncological need for early, broadly applicable readouts (biomarkers) of patient response to immunotherapy. The model is built upon the complex biological and physical interactions between the immune system and cancer, and is informed using only standard-of-care CT. We have retrospectively applied the model to 245 patients from multiple clinical trials treated with anti-CTLA-4 or anti-PD-1/PD-L1 antibodies. We found that model parameters distinctly identified patients with common (n = 18) and rare (n = 10) malignancy types who benefited and did not benefit from these monotherapies with accuracy as high as 88% at first restaging (median 53 days). Further, the parameters successfully differentiated pseudo-progression from true progression, providing previously unidentified insights into the unique biophysical characteristics of pseudo-progression. Our mathematical model offers a clinically relevant tool for personalized oncology and for engineering immunotherapy regimens.
Collapse
Affiliation(s)
- Joseph D. Butner
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX, USA
| | - Dalia Elganainy
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - Charles X. Wang
- MD/PhD Program, Drexel University College of Medicine, Philadelphia, PA 19102, USA
| | - Zhihui Wang
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX, USA
- Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - Shu-Hsia Chen
- Immunotherapy Research Center, Houston Methodist Research Institute, Houston, TX, USA
- Cancer Center, Houston Methodist Research Institute, Houston, TX, USA
| | - Nestor F. Esnaola
- Cancer Center, Houston Methodist Research Institute, Houston, TX, USA
- Department of Surgical Oncology, Fox Chase Cancer Center—Temple Health, Philadelphia, PA, USA
| | - Renata Pasqualini
- Rutgers Cancer Institute of New Jersey at University Hospital, Newark, NJ, USA
- Division of Cancer Biology, Department of Radiation Oncology, Rutgers New Jersey Medical School, Newark, NJ, USA
| | - Wadih Arap
- Rutgers Cancer Institute of New Jersey at University Hospital, Newark, NJ, USA
- Division of Hematology/Oncology, Department of Medicine, Rutgers New Jersey Medical School, Newark, NJ, USA
| | - David S. Hong
- Department of Investigational Cancer Therapeutics, The University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - James Welsh
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - Eugene J. Koay
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA
- Corresponding author. (V.C.); (E.J.K.)
| | - Vittorio Cristini
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX, USA
- MD/PhD Program, Drexel University College of Medicine, Philadelphia, PA 19102, USA
- Corresponding author. (V.C.); (E.J.K.)
| |
Collapse
|
20
|
Abstract
OBJECTIVE we present a multiscale agent-based model of Ductal Carcinoma in Situ (DCIS) in order to gain a detailed understanding of the cell-scale population dynamics, phenotypic distributions, and the associated interplay of important molecular signaling pathways that are involved in DCIS ductal invasion into the duct cavity (a process we refer to as duct advance rate here). METHODS DCIS is modeled mathematically through a hybridized discrete cell-scale model and a continuum molecular scale model, which are explicitly linked through a bidirectional feedback mechanism. RESULTS we find that duct advance rates occur in two distinct phases, characterized by an early exponential population expansion, followed by a long-term steady linear phase of population expansion, a result that is consistent with other modeling work. We further found that the rates were influenced most strongly by endocrine and paracrine signaling intensity, as well as by the effects of cell density induced quiescence within the DCIS population. CONCLUSION our model analysis identified a complex interplay between phenotypic diversity that may provide a tumor adaptation mechanism to overcome proliferation limiting conditions, allowing for dynamic shifts in phenotypic populations in response to variation in molecular signaling intensity. Further, sensitivity analysis determined DCIS axial advance rates and calcification rates were most sensitive to cell cycle time variation. SIGNIFICANCE this model may serve as a useful tool to study the cell-scale dynamics involved in DCIS initiation and intraductal invasion, and may provide insights into promising areas of future experimental research.
Collapse
|
21
|
Dogra P, Butner JD, Nizzero S, Ruiz Ramírez J, Noureddine A, Peláez MJ, Elganainy D, Yang Z, Le AD, Goel S, Leong HS, Koay EJ, Brinker CJ, Cristini V, Wang Z. Image-guided mathematical modeling for pharmacological evaluation of nanomaterials and monoclonal antibodies. Wiley Interdiscip Rev Nanomed Nanobiotechnol 2020; 12:e1628. [PMID: 32314552 PMCID: PMC7507140 DOI: 10.1002/wnan.1628] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/15/2019] [Revised: 02/06/2020] [Accepted: 02/15/2020] [Indexed: 12/13/2022]
Abstract
While plasma concentration kinetics has traditionally been the predictor of drug pharmacological effects, it can occasionally fail to represent kinetics at the site of action, particularly for solid tumors. This is especially true in the case of delivery of therapeutic macromolecules (drug-loaded nanomaterials or monoclonal antibodies), which can experience challenges to effective delivery due to particle size-dependent diffusion barriers at the target site. As a result, disparity between therapeutic plasma kinetics and kinetics at the site of action may exist, highlighting the importance of target site concentration kinetics in determining the pharmacodynamic effects of macromolecular therapeutic agents. Assessment of concentration kinetics at the target site has been facilitated by non-invasive in vivo imaging modalities. This allows for visualization and quantification of the whole-body disposition behavior of therapeutics that is essential for a comprehensive understanding of their pharmacokinetics and pharmacodynamics. Quantitative non-invasive imaging can also help guide the development and parameterization of mathematical models for descriptive and predictive purposes. Here, we present a review of the application of state-of-the-art imaging modalities for quantitative pharmacological evaluation of therapeutic nanoparticles and monoclonal antibodies, with a focus on their integration with mathematical models, and identify challenges and opportunities. This article is categorized under: Therapeutic Approaches and Drug Discovery > Nanomedicine for Oncologic Disease Diagnostic Tools > in vivo Nanodiagnostics and Imaging Nanotechnology Approaches to Biology > Nanoscale Systems in Biology.
Collapse
Affiliation(s)
- Prashant Dogra
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, Texas, USA
| | - Joseph D Butner
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, Texas, USA
| | - Sara Nizzero
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, Texas, USA
| | - Javier Ruiz Ramírez
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, Texas, USA
| | - Achraf Noureddine
- Department of Chemical and Biological Engineering, University of New Mexico, Albuquerque, New Mexico, USA
| | - María J Peláez
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, Texas, USA.,Applied Physics Graduate Program, Rice University, Houston, Texas, USA
| | - Dalia Elganainy
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, Texas, USA
| | - Zhen Yang
- Center for Bioenergetics, Houston Methodist Research Institute, Houston, Texas, USA
| | - Anh-Dung Le
- Nanoscience and Microsystems Engineering, University of New Mexico, Albuquerque, New Mexico, USA
| | - Shreya Goel
- Cancer Systems Imaging, University of Texas MD Anderson Cancer Center, Houston, Texas, USA
| | - Hon S Leong
- Biological Sciences Platform, Sunnybrook Research Institute, Toronto, Ontario, Canada.,Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada
| | - Eugene J Koay
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, Texas, USA
| | - C Jeffrey Brinker
- Department of Chemical and Biological Engineering and UNM Comprehensive Cancer Center, University of New Mexico, Albuquerque, New Mexico, USA
| | - Vittorio Cristini
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, Texas, USA
| | - Zhihui Wang
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, Texas, USA
| |
Collapse
|
22
|
Dogra P, Butner JD, Ruiz Ramírez J, Chuang YL, Noureddine A, Jeffrey Brinker C, Cristini V, Wang Z. A mathematical model to predict nanomedicine pharmacokinetics and tumor delivery. Comput Struct Biotechnol J 2020; 18:518-531. [PMID: 32206211 PMCID: PMC7078505 DOI: 10.1016/j.csbj.2020.02.014] [Citation(s) in RCA: 47] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2019] [Revised: 02/14/2020] [Accepted: 02/22/2020] [Indexed: 02/07/2023] Open
Abstract
Towards clinical translation of cancer nanomedicine, it is important to systematically investigate the various parameters related to nanoparticle (NP) physicochemical properties, tumor characteristics, and inter-individual variability that affect the tumor delivery efficiency of therapeutic nanomaterials. Comprehensive investigation of these parameters using traditional experimental approaches is impractical due to the vast parameter space; mathematical models provide a more tractable approach to navigate through such a multidimensional space. To this end, we have developed a predictive mathematical model of whole-body NP pharmacokinetics and their tumor delivery in vivo, and have conducted local and global sensitivity analyses to identify the factors that result in low tumor delivery efficiency and high off-target accumulation of NPs. Our analyses reveal that NP degradation rate, tumor blood viscosity, NP size, tumor vascular fraction, and tumor vascular porosity are the key parameters in governing NP kinetics in the tumor interstitium. The impact of these parameters on tumor delivery efficiency of NPs is discussed, and optimal values for maximizing NP delivery are presented.
Collapse
Affiliation(s)
- Prashant Dogra
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA
| | - Joseph D. Butner
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA
| | - Javier Ruiz Ramírez
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA
| | - Yao-li Chuang
- Department of Mathematics, California State University, Northridge, CA 91330, USA
| | - Achraf Noureddine
- Department of Chemical and Biological Engineering, University of New Mexico, Albuquerque, NM 87106, USA
| | - C. Jeffrey Brinker
- Department of Chemical and Biological Engineering, University of New Mexico, Albuquerque, NM 87106, USA
- UNM Comprehensive Cancer Center, University of New Mexico, Albuquerque, NM 87102, USA
| | - Vittorio Cristini
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA
| | - Zhihui Wang
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX 77030, USA
- Corresponding author at: Mathematics in Medicine Program, The Houston Methodist Research Institute, HMRI R8-122, 6670 Bertner Ave, Houston, TX 77030, USA.
| |
Collapse
|
23
|
Dogra P, Chuang YL, Butner JD, Cristini V, Wang Z. Development of a Physiologically-Based Mathematical Model for Quantifying Nanoparticle Distribution in Tumors. Annu Int Conf IEEE Eng Med Biol Soc 2020; 2019:2852-2855. [PMID: 31946487 DOI: 10.1109/embc.2019.8856503] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Nanomedicine holds promise for the treatment of cancer, as it enables tumor-targeted drug delivery. However, reports on translation of most nanomedicine strategies to the clinic so far have been less than satisfactory, in part due to insufficient understanding of the effects of nanoparticle (NP) physiochemical properties and physiological variables on their pharmacological behavior. In this paper, we present a multiscale mathematical model to examine the efficacy of NP delivery to solid tumors; as a case example, we apply the model to a clinically detectable primary pancreatic ductal adenocarcinoma (PDAC) to assess tissue-scale spatiotemporal distribution profiles of NPs. We integrate NP systemic disposition kinetics with NP-cell interactions in PDAC abstractly described as a two-dimensional structure, which is then parameterized with human physiological data obtained from published literature. Through model analysis of delivery efficiency, we verify the multiscale approach by showing that NP concentration kinetics of interest in various compartments predicted by the whole-body scale model were in agreement with those obtained from the tissue-scale model. We also found that more NPs were trapped in the outer well-perfused tumor region than the inner semi-necrotic domain. Further development of the model may provide a useful tool for optimal NP design and physiological interventions.
Collapse
|
24
|
Butner JD, Cristini V, Wang Z. Multiscale Modeling of Ductal Carcinoma In Situ. Biophys J 2019. [DOI: 10.1016/j.bpj.2018.11.1750] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022] Open
|
25
|
Butner JD, Cristini V. Development of a three dimensional, multiscale agent-based model of ductal carcinoma in situ. Annu Int Conf IEEE Eng Med Biol Soc 2018; 2017:86-89. [PMID: 29059817 DOI: 10.1109/embc.2017.8036769] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
Ductal carcinoma in situ (DCIS) is the most commonly diagnosed form of non-invasive breast cancer, constituting 20% of all new breast cancer cases in the United States. Although non-invasive, DCIS is usually treated surgically through resection. Interestingly, long-term survival studies have shown that patient survival rates are not significantly impacted by the type or resection, indicating that increased breast conservation through minimized surgical resection may indeed be possible. This requires a greater understanding of disease development, so that clinicians may more accurately determine surgical margins which minimize patient impact while maintaining survival rates. To this end, we have developed a three-dimensional, lattice-free multiscale agent based model of DCIS designed to help quantify ductal invasion rates, in order to allow clinicians to better estimate disease age and extent of invasion, and to predict surgical margins based on parameters obtainable from non-invasive testing (i.e., mammography). Here, we present the model development to date, and discuss some preliminary results.
Collapse
|
26
|
Butner JD, Cristini V. Development of a three dimensional, lattice-free multiscale model of the mammary terminal end bud. Annu Int Conf IEEE Eng Med Biol Soc 2017; 2016:6134-6137. [PMID: 28269652 DOI: 10.1109/embc.2016.7592128] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
The terminal end bud (TEB) is a bulbous structure composed of highly proliferative cells that is responsible for mammary gland development during the pubertal stage. This is a highly organized process, involving cellular differentiation hierarchies regulated by endocrine and paracrine signaling. Here, we present development of a lattice-free, three dimensional multiscale agent based model of the TEB to study the effects of cellular phenotypic hierarchies, endocrine and paracrine signaling, and proliferation demographics on pubertal mammary gland development. Cells in the TEB experience complex physical interaction during the active growth involved in pubertal ductal elongation, which we represent mathematically based on the physical forces involved in cell-cell and cell-microenvironment interactions. We observe that maximum ductal elongation rates are achievable due to each progenitive phenotype cell only undergoing a couple proliferation cycles before losing the progenitive capability, and that molecular signaling is necessary to restrict ductal elongation to biologically relevant rates. Cellular proliferation and growth is sufficient to achieve these elongation rates in the absence of other cellular behaviors such as migration or conformational changes. This model serves as a valuable tool to gain insights into the cell population dynamics of mammary gland development, and can serve as a foundation to study the early stages of breast cancer development based on endocrine-mediated phenotypic population shifts.
Collapse
|
27
|
Wang Z, Kerketta R, Chuang YL, Dogra P, Butner JD, Brocato TA, Day A, Xu R, Shen H, Simbawa E, AL-Fhaid AS, Mahmoud SR, Curley SA, Ferrari M, Koay EJ, Cristini V. Theory and Experimental Validation of a Spatio-temporal Model of Chemotherapy Transport to Enhance Tumor Cell Kill. PLoS Comput Biol 2016; 12:e1004969. [PMID: 27286441 PMCID: PMC4902302 DOI: 10.1371/journal.pcbi.1004969] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2015] [Accepted: 05/09/2016] [Indexed: 12/14/2022] Open
Abstract
It has been hypothesized that continuously releasing drug molecules into the tumor over an extended period of time may significantly improve the chemotherapeutic efficacy by overcoming physical transport limitations of conventional bolus drug treatment. In this paper, we present a generalized space- and time-dependent mathematical model of drug transport and drug-cell interactions to quantitatively formulate this hypothesis. Model parameters describe: perfusion and tissue architecture (blood volume fraction and blood vessel radius); diffusion penetration distance of drug (i.e., a function of tissue compactness and drug uptake rates by tumor cells); and cell death rates (as function of history of drug uptake). We performed preliminary testing and validation of the mathematical model using in vivo experiments with different drug delivery methods on a breast cancer mouse model. Experimental data demonstrated a 3-fold increase in response using nano-vectored drug vs. free drug delivery, in excellent quantitative agreement with the model predictions. Our model results implicate that therapeutically targeting blood volume fraction, e.g., through vascular normalization, would achieve a better outcome due to enhanced drug delivery.
Collapse
Affiliation(s)
- Zhihui Wang
- Department of NanoMedicine and Biomedical Engineering, University of Texas Medical School at Houston, Houston, Texas, United States of America
- Brown Foundation Institute of Molecular Medicine, University of Texas Medical School at Houston, Houston, Texas, United States of America
- Department of Imaging Physics, University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
| | - Romica Kerketta
- Department of Pathology, University of New Mexico, Albuquerque, New Mexico, United States of America
| | - Yao-Li Chuang
- Department of Mathematics, California State University, Northridge, California, United States of America
| | - Prashant Dogra
- Department of Pathology, University of New Mexico, Albuquerque, New Mexico, United States of America
| | - Joseph D. Butner
- Department of Chemical and Biological Engineering and Center for Biomedical Engineering, University of New Mexico, Albuquerque, New Mexico, United States of America
| | - Terisse A. Brocato
- Department of Chemical and Biological Engineering and Center for Biomedical Engineering, University of New Mexico, Albuquerque, New Mexico, United States of America
| | - Armin Day
- Department of Pathology, University of New Mexico, Albuquerque, New Mexico, United States of America
| | - Rong Xu
- Department of Nanomedicine, Methodist Hospital Research Institute, Houston, Texas, United States of America
| | - Haifa Shen
- Department of Nanomedicine, Methodist Hospital Research Institute, Houston, Texas, United States of America
| | - Eman Simbawa
- Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah, Saudi Arabia
| | - A. S. AL-Fhaid
- Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah, Saudi Arabia
| | - S. R. Mahmoud
- Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Steven A. Curley
- Michael E. DeBakey Department of Surgery, Baylor College of Medicine, Houston, Texas, United States of America
| | - Mauro Ferrari
- Department of Nanomedicine, Methodist Hospital Research Institute, Houston, Texas, United States of America
| | - Eugene J. Koay
- Department of Radiation Oncology, University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
- * E-mail: (EJK); (VC)
| | - Vittorio Cristini
- Department of NanoMedicine and Biomedical Engineering, University of Texas Medical School at Houston, Houston, Texas, United States of America
- Brown Foundation Institute of Molecular Medicine, University of Texas Medical School at Houston, Houston, Texas, United States of America
- Department of Imaging Physics, University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
- Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah, Saudi Arabia
- * E-mail: (EJK); (VC)
| |
Collapse
|
28
|
Abstract
Terminal end buds (TEBs) are bulb-like structures at the growing tips of elongating mammary ducts, and the growth of a TEB is a complex, organized biological process. In this paper, we present a hybrid continuum-discrete agent-based model to provide quantitative insight into the properties of cell symmetric and asymmetric division on the spatial and developing cell rearrangement within the TEB during ductal elongation. An interplay of endocrine-paracrine signaling and cell lineage has been implemented in the model. Our results show that higher symmetric division rates resulted in more progenitor cells remaining in the TEB, while lower rates resulted in more differentiated cells in the TEB. Moreover, pure proliferation alone was enough to result in ductal elongation in the absence of any cellular migration, a result consistent with current experimental data. This model can also serve as a platform to study how mutation-induced phenotypic changes contribute to developmental defects in mammary gland development.
Collapse
|
29
|
Wang Z, Butner JD, Cristini V, Deisboeck TS. Integrated PK-PD and agent-based modeling in oncology. J Pharmacokinet Pharmacodyn 2015; 42:179-89. [PMID: 25588379 DOI: 10.1007/s10928-015-9403-7] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2014] [Accepted: 01/08/2015] [Indexed: 01/11/2023]
Abstract
Mathematical modeling has become a valuable tool that strives to complement conventional biomedical research modalities in order to predict experimental outcome, generate new medical hypotheses, and optimize clinical therapies. Two specific approaches, pharmacokinetic-pharmacodynamic (PK-PD) modeling, and agent-based modeling (ABM), have been widely applied in cancer research. While they have made important contributions on their own (e.g., PK-PD in examining chemotherapy drug efficacy and resistance, and ABM in describing and predicting tumor growth and metastasis), only a few groups have started to combine both approaches together in an effort to gain more insights into the details of drug dynamics and the resulting impact on tumor growth. In this review, we focus our discussion on some of the most recent modeling studies building on a combined PK-PD and ABM approach that have generated experimentally testable hypotheses. Some future directions are also discussed.
Collapse
Affiliation(s)
- Zhihui Wang
- Department of Pathology, University of New Mexico, Albuquerque, NM, 87131, USA
| | | | | | | |
Collapse
|
30
|
Wang Z, Butner JD, Kerketta R, Cristini V, Deisboeck TS. Simulating cancer growth with multiscale agent-based modeling. Semin Cancer Biol 2014; 30:70-8. [PMID: 24793698 DOI: 10.1016/j.semcancer.2014.04.001] [Citation(s) in RCA: 129] [Impact Index Per Article: 12.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2014] [Revised: 03/18/2014] [Accepted: 04/04/2014] [Indexed: 01/01/2023]
Abstract
There have been many techniques developed in recent years to in silico model a variety of cancer behaviors. Agent-based modeling is a specific discrete-based hybrid modeling approach that allows simulating the role of diversity in cell populations as well as within each individual cell; it has therefore become a powerful modeling method widely used by computational cancer researchers. Many aspects of tumor morphology including phenotype-changing mutations, the adaptation to microenvironment, the process of angiogenesis, the influence of extracellular matrix, reactions to chemotherapy or surgical intervention, the effects of oxygen and nutrient availability, and metastasis and invasion of healthy tissues have been incorporated and investigated in agent-based models. In this review, we introduce some of the most recent agent-based models that have provided insight into the understanding of cancer growth and invasion, spanning multiple biological scales in time and space, and we further describe several experimentally testable hypotheses generated by those models. We also discuss some of the current challenges of multiscale agent-based cancer models.
Collapse
Affiliation(s)
- Zhihui Wang
- Department of Pathology, University of New Mexico, Albuquerque, NM 87131, USA.
| | - Joseph D Butner
- Department of Chemical Engineering and Center for Biomedical Engineering, University of New Mexico, Albuquerque, NM 87131, USA
| | - Romica Kerketta
- Department of Pathology, University of New Mexico, Albuquerque, NM 87131, USA
| | - Vittorio Cristini
- Department of Pathology, University of New Mexico, Albuquerque, NM 87131, USA; Department of Chemical Engineering and Center for Biomedical Engineering, University of New Mexico, Albuquerque, NM 87131, USA; Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
| | | |
Collapse
|
31
|
Pascal J, Ashley CE, Wang Z, Brocato TA, Butner JD, Carnes EC, Koay EJ, Brinker CJ, Cristini V. Mechanistic modeling identifies drug-uptake history as predictor of tumor drug resistance and nano-carrier-mediated response. ACS Nano 2013; 7:11174-11182. [PMID: 24187963 PMCID: PMC3891887 DOI: 10.1021/nn4048974] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
A quantitative understanding of the advantages of nanoparticle-based drug delivery vis-à-vis conventional free drug chemotherapy has yet to be established for cancer or other diseases despite numerous investigations. Here, we employ first-principles cell biophysics, drug pharmaco-kinetics, and drug pharmaco-dynamics to model the delivery of doxorubicin (DOX) to hepatocellular carcinoma (HCC) tumor cells and predict the resultant experimental cytotoxicity data. The fundamental, mechanistic hypothesis of our mathematical model is that the integrated history of drug uptake by the cells over time of exposure, which sets the cell death rate parameter, and the uptake rate are the sole determinants of the dose response relationship. A universal solution of the model equations is capable of predicting the entire, nonlinear dose response of the cells to any drug concentration based on just two separate measurements of these cellular parameters. This analysis reveals that nanocarrier-mediated delivery overcomes resistance to the free drug because of improved cellular uptake rates, and that dose response curves to nanocarrier mediated drug delivery are equivalent to those for free-drug, but "shifted to the left;" that is, lower amounts of drug achieve the same cell kill. We then demonstrate the model's general applicability to different tumor and drug types, and cell-exposure time courses by investigating HCC cells exposed to cisplatin and 5-fluorouracil, breast cancer MCF-7 cells exposed to DOX, and pancreatic adenocarcinoma PANC-1 cells exposed to gemcitabine. The model will help in the optimal design of nanocarriers for clinical applications and improve the current, largely empirical understanding of in vivo drug transport and tumor response.
Collapse
Affiliation(s)
- Jennifer Pascal
- Department of Pathology, The University of New Mexico Health Sciences Center, Albuquerque, NM 87131, USA
| | - Carlee E. Ashley
- Biotechnology and Bioengineering Department, Sandia National Laboratories, Livermore, CA 94551-0969, USA
- Cancer Research and Treatment Center, The University of New Mexico Health Sciences Center, Albuquerque, NM 87131, USA
| | - Zhihui Wang
- Department of Pathology, The University of New Mexico Health Sciences Center, Albuquerque, NM 87131, USA
| | - Terisse A. Brocato
- Department of Chemical and Nuclear Engineering, and Center for Biomedical Engineering, The University of New Mexico, Albuquerque, NM 87131, USA
| | - Joseph D. Butner
- Department of Chemical and Nuclear Engineering, and Center for Biomedical Engineering, The University of New Mexico, Albuquerque, NM 87131, USA
| | - Eric C. Carnes
- Nanobiology Department, Sandia National Laboratories, Albuquerque, CA 87185-1349
- Department of Chemical and Nuclear Engineering, and Center for Biomedical Engineering, The University of New Mexico, Albuquerque, NM 87131, USA
- Cancer Research and Treatment Center, The University of New Mexico Health Sciences Center, Albuquerque, NM 87131, USA
| | - Eugene J. Koay
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX 77030
- Department of Nanomedicine, Houston Methodist Research Institute, Houston, TX 77030
| | - C. Jeffrey Brinker
- Department of Chemical and Nuclear Engineering, and Center for Biomedical Engineering, The University of New Mexico, Albuquerque, NM 87131, USA
- Center for Micro-Engineered Materials, The University of New Mexico, Albuquerque, NM 87131, USA
- Cancer Research and Treatment Center, The University of New Mexico Health Sciences Center, Albuquerque, NM 87131, USA
- Department of Molecular Genetics and Microbiology, the University of New Mexico Health Sciences Center, Albuquerque, NM 87131, USA
- Self-Assembled Materials Department, Sandia National Laboratories, Albuquerque, NM 87185-1349, USA
| | - Vittorio Cristini
- Department of Pathology, The University of New Mexico Health Sciences Center, Albuquerque, NM 87131, USA
- Department of Chemical and Nuclear Engineering, and Center for Biomedical Engineering, The University of New Mexico, Albuquerque, NM 87131, USA
- Cancer Research and Treatment Center, The University of New Mexico Health Sciences Center, Albuquerque, NM 87131, USA
| |
Collapse
|