Predicting breast cancer response to neoadjuvant chemotherapy based on tumor vascular features in needle biopsies

An ultrasound-based nomogram for predicting axillary node pathologic complete response after neoadjuvant chemotherapy in breast cancer: Modeling and external validation

INTRODUCTION

The staging and treatment of axillary nodes in breast cancer have become a focus of research. For breast cancer patients with fine-needle aspiration-or core needle biopsy-confirmed positive nodes, axillary lymph node dissection (ALND) after neoadjuvant chemotherapy (NAC) is still a standard treatment. However, some patients achieve an axillary pathologic complete response (pCR) after NAC. In this study, the authors sought to construct a model to predict axillary pCR in patients with positive axillary lymph nodes (cN+) breast cancer.

METHODS

Data from patients with pathologically proven cN+ breast cancer treated with NAC followed by ALND between January 2010 and April 2019 at the Peking University Cancer Hospital were reviewed. Axillary lymph node status was assessed using ultrasonography before and after NAC. The patient cohort was assigned to the construction and internal validation cohorts according to admission time. A nomogram was constructed based on the significant factors associated with axillary pCR. The predictive performance of the model was externally validated using data from Peking University First Hospital.

RESULTS

This study included 953 and 267 patients from Peking University Cancer Hospital and Peking University First Hospital, respectively. In the construction cohort, 39.7% (238 of 600) of patients achieved axillary pCR after NAC. The result of multivariate logistic regression analysis showed that tumor grade, clinical nodal response, NAC regimen, tumor pCR, lymphovascular invasion, and tumor biologic subtype were significant independent predictors of ypN0 (p < 0.05). The areas under the receiver operating characteristic curves for the construction, validation, and independent testing cohorts were 0.87 (95% confidence interval [CI], 0.84-0.90), 0.83 (95% CI, 0.79-0.87), and 0.84 (0.79-0.89), respectively.

CONCLUSIONS

A nomogram was constructed to predict the pCR of axillary lymph nodes after NAC for breast cancer. Validation of both the internal and external cohorts achieved good predictive performance, indicating that the model has preliminary clinical application prospects.

Prediction of complete response after neoadjuvant chemotherapy for invasive breast cancers: The utility of shear wave elastography and superb microvascular imaging in pretreatment breast ultrasound

PURPOSE

To investigate whether multiparametric parameters of pretreatment breast ultrasound (US) and clinicopathologic factors are associated with pathologic complete response (pCR) after neoadjuvant chemotherapy (NAC) for breast cancer.

METHODS

Between November 2018 and September 2022, 88 patients who underwent NAC and subsequent surgery were included in this study (median age, 55 years; interquartile range [IQR], 45, 59.3). Multiparametric breast US including grayscale, shear wave elastography (SWE) and superb microvascular imaging (SMI) of pathologically proven invasive breast cancers were retrospectively reviewed. Clinicopathological and multiparametric parameters of breast US, including size, SWEmax, SWEratio and vascular index on SMI (SMIVI) were compared between the groups. Univariate and multivariate logistic regression analyses were performed to determine factors predicting pCR after NAC. AUROC curve analysis was performed to determine the predictors' optimal cut-off values and diagnostic performance.

RESULTS

The pCR group (n = 24) showed a significantly smaller tumor size, lower SWEmax, higher Ki-67 index, higher hormone receptor negativity and negative axillary lymph node metastasis compared to the non-pCR group (n = 64). Multivariate regression analysis showed that SWEmax (adjusted odds ratio[aOR] = 0.956, 95 % confidence interval [CI] = 0.919-0.994, P = 0.025) and Ki-67 index (aOR = 1.083, 95 % CI = 1.012-1.159, P = 0.021) were independently associated with pathologically complete response. The optimal cut-off values for predicting pCR were 27.5 % for Ki-67 with an AUC of 0.743 and 134.8 kPa for SWEmax with an AUC of 0.779. A combination model including clinical factors and SWEmax showed the best diagnostic performance with an AUC of 0.876.

CONCLUSION

A higher Ki-67 index and lower SWEmax measured on pretreatment breast US were independently associated with pCR in invasive breast cancer after NAC.

Protocol for mathematical prediction of patient response and survival to immune checkpoint inhibitor immunotherapy

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)¹ and Butner et al. (2021).².

Dedifferentiation-mediated stem cell niche maintenance in early-stage ductal carcinoma in situ progression: insights from a multiscale modeling study

We present a multiscale agent-based model of ductal carcinoma in situ (DCIS) to study how key phenotypic and signaling pathways are involved in the early stages of disease progression. The model includes a phenotypic hierarchy, and key endocrine and paracrine signaling pathways, and simulates cancer ductal growth in a 3D lattice-free domain. In particular, by considering stochastic cell dedifferentiation plasticity, the model allows for study of how dedifferentiation to a more stem-like phenotype plays key roles in the maintenance of cancer stem cell populations and disease progression. Through extensive parameter perturbation studies, we have quantified and ranked how DCIS is sensitive to perturbations in several key mechanisms that are instrumental to early disease development. Our studies reveal that long-term maintenance of multipotent stem-like cell niches within the tumor are dependent on cell dedifferentiation plasticity, and that disease progression will become arrested due to dilution of the multipotent stem-like population in the absence of dedifferentiation. We have identified dedifferentiation rates necessary to maintain biologically relevant multipotent cell populations, and also explored quantitative relationships between dedifferentiation rates and disease progression rates, which may potentially help to optimize the efficacy of emerging anti-cancer stem cell therapeutics.

Early prediction of clinical response to checkpoint inhibitor therapy in human solid tumors through mathematical modeling

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.

Targeting a cell surface vitamin D receptor on tumor-associated macrophages in triple-negative breast cancer

Triple-negative breast cancer (TNBC) is an aggressive tumor with limited treatment options and poor prognosis. We applied the in vivo phage display technology to isolate peptides homing to the immunosuppressive cellular microenvironment of TNBC as a strategy for non-malignant target discovery. We identified a cyclic peptide (CSSTRESAC) that specifically binds to a vitamin D receptor, protein disulfide-isomerase A3 (PDIA3) expressed on the cell surface of tumor-associated macrophages (TAM), and targets breast cancer in syngeneic TNBC, non-TNBC xenograft, and transgenic mouse models. Systemic administration of CSSTRESAC to TNBC-bearing mice shifted the cytokine profile toward an antitumor immune response and delayed tumor growth. Moreover, CSSTRESAC enabled ligand-directed theranostic delivery to tumors and a mathematical model confirmed our experimental findings. Finally, in silico analysis showed PDIA3-expressing TAM in TNBC patients. This work uncovers a functional interplay between a cell surface vitamin D receptor in TAM and antitumor immune response that could be therapeutically exploited.

A Mathematical Model to Estimate Chemotherapy Concentration at the Tumor-Site and Predict Therapy Response in Colorectal Cancer Patients with Liver Metastases

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.

A mathematical model for the quantification of a patient's sensitivity to checkpoint inhibitors and long-term tumour burden

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.

Investigating the Effect of Aging on the Pharmacokinetics and Tumor Delivery of Nanomaterials using Mathematical Modeling

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.

Therapeutic potential of FLANC, a novel primate-specific long non-coding RNA in colorectal cancer

OBJECTIVE

To investigate the function of a novel primate-specific long non-coding RNA (lncRNA), named FLANC, based on its genomic location (co-localised with a pyknon motif), and to characterise its potential as a biomarker and therapeutic target.

DESIGN

FLANC expression was analysed in 349 tumours from four cohorts and correlated to clinical data. In a series of multiple in vitro and in vivo models and molecular analyses, we characterised the fundamental biological roles of this lncRNA. We further explored the therapeutic potential of targeting FLANC in a mouse model of colorectal cancer (CRC) metastases.

RESULTS

FLANC, a primate-specific lncRNA feebly expressed in normal colon cells, was significantly upregulated in cancer cells compared with normal colon samples in two independent cohorts. High levels of FLANC were associated with poor survival in two additional independent CRC patient cohorts. Both in vitro and in vivo experiments demonstrated that the modulation of FLANC expression influenced cellular growth, apoptosis, migration, angiogenesis and metastases formation ability of CRC cells. In vivo pharmacological targeting of FLANC by administration of 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine nanoparticles loaded with a specific small interfering RNA, induced significant decrease in metastases, without evident tissue toxicity or pro-inflammatory effects. Mechanistically, FLANC upregulated and prolonged the half-life of phosphorylated STAT3, inducing the overexpression of VEGFA, a key regulator of angiogenesis.

CONCLUSIONS

Based on our findings, we discovered, FLANC as a novel primate-specific lncRNA that is highly upregulated in CRC cells and regulates metastases formation. Targeting primate-specific transcripts such as FLANC may represent a novel and low toxic therapeutic strategy for the treatment of patients.

Mathematical prediction of clinical outcomes in advanced cancer patients treated with checkpoint inhibitor immunotherapy

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.

Image-guided mathematical modeling for pharmacological evaluation of nanomaterials and monoclonal antibodies

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.

Crosstalk between HER2 and PD-1/PD-L1 in Breast Cancer: From Clinical Applications to Mathematical Models

Breast cancer is one of the major causes of mortality in women worldwide. The most aggressive breast cancer subtypes are human epidermal growth factor receptor-positive (HER2+) and triple-negative breast cancers. Therapies targeting HER2 receptors have significantly improved HER2+ breast cancer patient outcomes. However, several recent studies have pointed out the deficiency of existing treatment protocols in combatting disease relapse and improving response rates to treatment. Overriding the inherent actions of the immune system to detect and annihilate cancer via the immune checkpoint pathways is one of the important hallmarks of cancer. Thus, restoration of these pathways by various means of immunomodulation has shown beneficial effects in the management of various types of cancers, including breast. We herein review the recent progress in the management of HER2+ breast cancer via HER2-targeted therapies, and its association with the programmed death receptor-1 (PD-1)/programmed death ligand-1 (PD-L1) axis. In order to link research in the areas of medicine and mathematics and point out specific opportunities for providing efficient theoretical analysis related to HER2+ breast cancer management, we also review mathematical models pertaining to the dynamics of HER2+ breast cancer and immune checkpoint inhibitors.

A mathematical model to predict nanomedicine pharmacokinetics and tumor delivery

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.

Development of a Physiologically-Based Mathematical Model for Quantifying Nanoparticle Distribution in Tumors

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.

Mathematical Modeling to Address Challenges in Pancreatic Cancer

Pancreatic Ductal Adenocarcinoma (PDAC) is regarded as one of the most lethal cancer types for its challenges associated with early diagnosis and resistance to standard chemotherapeutic agents, thereby leading to a poor five-year survival rate. The complexity of the disease calls for a multidisciplinary approach to better manage the disease and improve the status quo in PDAC diagnosis, prognosis, and treatment. To this end, the application of quantitative tools can help improve the understanding of disease mechanisms, develop biomarkers for early diagnosis, and design patient-specific treatment strategies to improve therapeutic outcomes. However, such approaches have only been minimally applied towards the investigation of PDAC, and we review the current status of mathematical modeling works in this field.