Abstract
We recently developed a computational model of cisplatin pharmacodynamics in an endobronchial lung tumor following ultrasound-guided transbronchial needle injection (EBUS-TBNI). The model suggests that it is more efficacious to apportion the cisplatin dose between injections at different sites rather than giving it all in a single central injection, but the model was calibrated only on blood cisplatin data from a single patient. Accordingly, we applied a modified version of our original model in a set of 32 patients undergoing EBUS-TBNI for non-small cell lung cancer (NSCLC). We used the model to predict clinical responses and compared them retrospectively to actual patient outcomes. The model correctly predicted the clinical response in 72% of cases, with 80% accuracy for adenocarcinomas and 62.5% accuracy for squamous-cell lung cancer. We also found a power-law relationship between tumor volume and the minimal dose needed to induce a response, with the power-law exponent depending on the number of injections administered. Our results suggest that current injection strategies may be significantly over- or under-dosing the agent depending on tumor size, and that computational modeling can be a useful planning tool for EBUS-TBNI of cisplatin in lung cancer.
Similar content being viewed by others
Introduction
Lung cancer is the leading cause of cancer-related mortality in the world, with 1.76 million deaths in 2018. This is more than the next three deadliest types of cancer combined 1. Moreover, even though lung cancer accounts for 13% of new cancer cases, it is responsible for 22% of the deaths2. There is a significant interest in therapies that could potentially improve lung cancer response rates.
Endobronchial ultrasound-guided transbronchial needle injection (EBUS-TBNI) of cisplatin has recently emerged as a safe alternative to systemic delivery for treating recurrent centrally located non-small cell lung cancer (NSCLC)3,4. The potential advantage of intratumoral delivery of chemotherapy is that it can achieve high concentrations of cytotoxic agent within a tumor while reducing the off-target tissue burden that can lead to adverse side effects 5,6. However, there is still no consensus as to exactly how intratumoral cisplatin should be delivered to NSCLC tumors. In fact, only two strategies have been employed clinically to date: (1) administering the entire cisplatin dose as a single central injection7, and (2) dividing the dose among five distinct injection sites throughout the tumor8. In addition, the dose itself has been chosen empirically, most commonly being up to 40 mg9. There would thus seem to be significant opportunity to advance the efficacy of EBUS-TBNI of cisplatin by determining the dose that achieves the best tradeoff between cytotoxicity and systemic side-effects, and by employing an injection strategy that ensures all tumor cells receive a lethal concentration of agent. Indeed, our previous computational modeling work10 suggests that there are enormous gains to be made in distributing a given cisplatin dose at optimally chosen injection sites throughout the tumor as opposed to depositing it all at a single location.
The goal of the present work was therefore to investigate how various doses and injection strategies influence treatment response to EBUS-TBNI of cisplatin. As we are not yet at the point where this question can be studied prospectively, we examined the responses achieved in a case series of patients who had received EBUS-TBNI of cisplatin as salvage therapy for late-stage NSCLC. We adapted our previously developed computational model to predict the likelihood of tumor reduction in each patient based on the dose and number of injections they were given, and then compared these predictions to the clinically observed responses. We also used the model to estimate the minimal cisplatin dose required to induce a positive response in each tumor, as well as how this dose depends on the number of spatially distributed injections, with the goal of establishing an initial foundation upon which to set rationalized guidelines for cisplatin dose as a function of tumor volume and delivery strategy.
Materials and methods
Patient data
The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). This study was approved by the University of Vermont (UVM) Committee on Human Research in the Medical Sciences (CHRMS 17-075) and by The University of Florida (UFL) Institutional Review Board (201700864). Informed consent was not required by these oversight committees due to the retrospective nature of the study. A cohort of patients that underwent EBUS-TBNI from January 2009 to November 2018 was retrospectively assessed.
Computational modeling
We implemented our previously developed computational model of cisplatin pharmacodynamics following intratumoral injection10. Briefly, this model considers a tumor to have a boundary defined by segmentation of the tumor in a 3D CT scan11. The voxel size reflects the resolution of the 3-D CT scan so that tumor volume and shape can be accurately represented in the model. Volume and shape are the characteristics that distinguish the different tumors from the perspective of our model. Each voxel in the segmented tumor is assumed to be comprised of the superposition of two distinct compartments corresponding to extracellular and intracellular spaces. Cisplatin diffuses within the extracellular compartment \(\left({\varphi }_{e}\right)\) with a spatially uniform diffusion constant \(D\). It also moves irreversibly into the intracellular compartment \(\left({\varphi }_{i}\right)\) with rate-constant \({k}_{i}\), where it exerts its cytotoxic effects. Cisplatin is cleared with rate-constant \({k}_{f}\) by the tumor perfusion and deposited within a systemic fluid compartment \(\left({\varphi }_{f}\right)\) having volume \({V}_{f}\) from which it is eliminated via the kidneys with rate-constant \({k}_{r}\). The model is represented diagrammatically in Fig. 1 and mathematically in Eqs. (1–3). \({\varphi }_{e}\), \({\varphi }_{i}\) and \({\varphi }_{f}\) represent the cisplatin concentrations in the extracellular, intracellular and fluid compartments, respectively.
Based on empirical observations12, the cisplatin diffusion coefficient within each tissue space was estimated as a function of molecular weight to be \(2.47\times {10}^{-6}\) cm2/s. The intracellular uptake rate constant was taken from the literature13,14 to be \(1.05\times {10}^{-4}\)/s. The total volume of the fluid space \(\left({V}_{f}\right)\) and the cisplatin clearance rate from the tumor to the fluid space \(\left({k}_{f}\right)\) were set to be \(12.2\) L and \(1.46\times {10}^{-4}\)/s, respectively. These values are based on the fit of the model to the cisplatin concentrations in the blood from a single patient monitored for 2 h following intratumoral injection10.
The model is initialized by loading the initial dose of cisplatin into the extracellular compartment \(\left({m}_{j}\right)\), apportioned between 2 and 5 locations \(\left(N\right)\) spaced according to our previously described algorithm10. Briefly, this algorithm identifies those injection sites that will give rise to the most homogeneous cisplatin distribution pattern within the tumor. Because we assume cisplatin to bind irreversibly to DNA in the intracellular compartment, its concentration in that compartment increases monotonically with time toward an asymptote, \({\varphi }_{\infty }\), defined by the solution of Eqs. (1–3)10:
Even though several different features play a role in determining tumor response to treatment, for the sake of simplicity we assume that apoptosis of cancer cells occurs when the intracellular cisplatin concentration exceeds a threshold value, \({\varphi }_{t}\) that was empirically set to be \(0.5\times {10}^{-7}\) mg/ml in the original version of the model. In the updated version of the model used in the present study, \({\varphi }_{t}\) is set using the IC50 values estimated in vitro from experiments in immortalized human cancer cell lines (Genomics of Drug Sensitivity in Cancer)15. We assume that treatment response is related to a 30% reduction in tumor diameter, following RECIST v.1.1 criteria16. This is roughly equivalent to a 66% volume reduction for a spherical tumor. Therefore, treatment response will be achieved in the model if at least 66% of the tumor volume reaches an intracellular concentration above \({\varphi }_{t}\). We used the computational model to determine whether or not this was the case for each tumor using the injection strategy that was employed in the tumor in question.
Results
We screened 37 patients treated for lung cancer with EBUS-TBNI of cisplatin, the inclusion criteria being7 age 18–80 years, pathologically confirmed non-small cell lung cancer, histologic or cytologic recurrence of cancer following therapy at initial diagnosis, recurrence in an EBUS-accessible site, and computed tomography (CT) scans performed less than five weeks prior to treatment. The characteristics of the cohort are shown in Table 1.
Out of the 37 patients in the dataset, 33 met the inclusion criteria. We excluded one additional patient with squamous cell lung cancer who passed away prior to assessment of therapeutic response. Of the remaining 32 patients, 15 had confirmed adenocarcinoma, 16 had confirmed squamous-cell carcinoma, and 1 had large-cell carcinoma.
All patients received a total cisplatin dose, determined empirically by the treating physician, ranging from 10 to 40 mg per tumor. Treated lesion response to therapy was determined based on RECIST v1.1 critera16. Partial or complete responses by RECIST were considered positive, so the patients were divided into a group that responded to treatment (complete or partial remission) and a group that did not respond (stable or progressive disease). Of the 22 patients treated at the University of Florida, 1 received 2 injections and 21 received 4 injections according to the injected dose, which was 20 or 40 mg, respectively. Of the 10 patients treated at the University of Vermont, all were administered 5 injections regardless of the dose that ranged from 10 to 40 mg. Multiple injections were administered at sites distributed roughly uniformly throughout the tumor as determined empirically by the treating physician, and dose was selected without formal regard to tumor size.
All CT scans were evaluated by an observer who was blinded to the tumor response. Semi-automated segmentation of the treated volume was performed as previously described using the freely available Chest Imaging Platform (https://chestimagingplatform.org)17.
Treatment response was observed in 25 patients (78%), with an 80% response rate in patients with adenocarcinoma, a 68.7% response rate in those with squamous-cell carcinoma, and a 100% response rate in patients with large-cell carcinoma.
The dataset from the Genomics of Drug Sensitivity in Cancer platform15 showed that the 72-h IC50 for cisplatin in 75 immortalized human NSCLC cell lines (adenocarcinoma and squamous-cell) ranged from 0.66 to 155.1 μM (\(1.99\times {10}^{-4} \,\mathrm{to} \, 4.65\times {10}^{-2}\) mg/ml). The logarithm of the IC50 follows a normal distribution (p > 0.05, Shapiro–Wilk test). Figure 2 shows the distribution of the log-transformed IC50 values and the fitted Gaussian curve \(\left(\mu =-2.59, \sigma =0.50\right)\). In order to embrace the heterogeneity observed in the IC50 values for different cell lines, we drew the threshold intracellular concentration for cell apoptosis in each voxel \(\left({\varphi }_{t}\right)\) randomly from the Gaussian distribution shown in Fig. 2.
The model correctly predicted treatment response in 72% of the 32 patients in the dataset. Among the adenocarcinomas, the model performed correctly in 80% of the cases (Table 2). Among the squamous-cell carcinomas, the model performed correctly in 62.5% of the cases (Table 3). The model also correctly classified the outcome for the single case of large-cell carcinoma (Table 4). The confusion matrices for all patients, those with adenocarcinoma, and those with squamous-cell carcinoma are given in Fig. 3A–C, respectively.
Interestingly, the two tumors associated with the false negatives in the adenocarcinoma group (ADN-1 and ADN-15) were substantially larger (152.6 and 229.9 cm3, respectively) than the group mean of 14.1 cm3 (SD 19.2 cm3) (p = 0.019, Mann–Whitney U test). Among the six small squamous-cell tumors, which all had volumes lower than 4 cm3, there were three false positives and three true positives. In the group with tumor volumes > 4 cm3 there were two false negatives and one false positive.
We performed a model sensitivity analysis by reassessing its prediction accuracy following an order of magnitude variation in each of the four parameters in turn. These parameters are (1) the diffusion coefficient \(\left(D\right)\), (2) the flux rate to the intracellular space \(\left({k}_{i}\right)\), (3) the clearance rate to the fluid space \(\left({k}_{f}\right)\) and (4) the mean value of the IC50 distribution \(\left(\mu \right)\) that determines the intracellular threshold concentration for cell apoptosis \(\left({\varphi }_{t}\right)\).
The results are presented in Table 5 and show that the model predictions of therapeutic outcome were essentially unaffected by large variations in model parameter values, with the threshold concentration for apoptosis in the intracellular space being the most sensitive parameter. Indeed, the only meaningful alteration in model prediction occurred with a tenfold reduction in the mean value of the IC50 distribution, implying that particular attention needs to be given to the determination of this quantity10.
The minimum cisplatin dose estimated by the model to induce a treatment response is non-linearly dependent on tumor volume (Fig. 4A). For tumors < 40 cm3 in volume, this relationship is accurately approximated by a power-law of the form \(dose=\alpha \cdot {Volume}^{\beta }\), where \(\alpha \) and \(\beta \) are constants that depend on the number of injections. The estimated values of \(\alpha \) and \(\beta \) are presented in Table 6. Of note, the exponent of the power-law relationship, \(\beta \), decreased systematically as the number of injections increased from 1 to 5 (Fig. 4B–F).
Discussion
Our computational model of cisplatin pharmacodynamics following intratumoral injection in lung cancer is, like all models, based on a set of simplifying assumptions10. In the present study, we put these assumptions to the test by comparing the predictions of the model to clinical outcomes in a cohort of non-small cell lung cancer patients receiving EBUS-TBNI of cisplatin (Fig. 3, Tables 2, 3, 4). The predictive accuracy of the model was 62.5% for squamous-cell lung cancers and 80% for adenocarcinoma. Moreover, a sensitivity analysis showed that the model is robust to large changes in the values of its key parameters (Table 5). We believe these findings should be validated, of course, given the limited size of the patient cohort. This reflects the fact that EBUS-TBNI of cisplatin for lung cancer is still an innovative therapy, its sole current indication being as salvage therapy for patients with few other therapeutic options. In fact, the data we used were collected retrospectively from the only two centers in the US that have reported on the use of this therapy, so our cohort remains the largest described. Nevertheless, our findings suggest that the model, despite its many simplifying assumptions, manages to capture important features of cisplatin pharmacodynamics following injection and thus has the potential to serve as a useful adjunct for treatment planning.
Using a literature-based estimate of the IC50 range for cisplatin (Fig. 2), we also estimated the minimum cisplatin dose to achieve a therapeutic response for different delivery strategies as a function of tumor volume (Fig. 4A). Interestingly, we found a strong nonlinear relationship between cisplatin dose and tumor volume that can be well represented by a power-law (Fig. 4B–F) in the 0–40 cm3 range, which encompasses 87% of our dataset. This is a potentially important finding that challenges the previously established empirical guideline of 2 mg of cisplatin per ml of tumor18. Moreover, a 10–40 mg strategy regardless of volume and number of injections does not seem to capture the myriad of nuances that drive therapeutic response; it is clearly an overdose for small tumors to receive 40 mg in 5 injections, while delivering only 10 mg in a single injection might not be nearly enough for a large tumor.
Due to the nonlinear nature of the power-law relationships identified in Fig. 4, it is apparent that treatment efficacy depends strongly on the number of injections, especially for larger tumors, and that even small variations in the exponent of the power law relationship can result in substantial effects on treatment outcome. For example, for a 40 cm3 tumor, the model predicts that 5 injections can reduce the total dose required for efficacy by 3–5 orders of magnitude. On the other hand, dose reduction is only 1–2 orders of magnitude for tumors smaller than 4 cm3. This makes intuitive sense given the smaller diffusion distances involved in having cisplatin reach all corners of a smaller tumor, and is also fortuitous given the greater practical challenges of accurately spacing injection sites throughout a smaller tumor. Our model simulations also suggest a formula for calculating cisplatin dosing, up to 40 mg, for a given number of injections according to a power-law relationship with parameters listed in Table 6.
Evaluation of the cases misclassified by the model provides further insight into potential areas for improvement. For adenocarcinomas, two out of the three misclassifications were false negative in tumors significantly larger than the rest of the cohort. Large tumors can develop high interstitial fluid pressures, permanently shutting down small blood vessels that perfuse the tumor. This limits oxygenation in effected areas and causes local necrosis19,20. Although necrosis is less common in adenocarcinomas, it can occur in very large tumors21. The absence of blood vessels in parts of a tumor limits vascular clearance and so facilitates cisplatin dissemination around the extracellular space, potentially explaining a positive clinical outcome at odds with model predictions for large tumors. Large tumor size may thus suggest the presence of poorly oxygenated areas as a result of necrosis, but this could only be confirmed by histological assessment.
This study has a number of limitations. For example, we did not have the data necessary to calibrate model parameters to the individual patients in our cohort, so the flux of cisplatin to the extratumoral space, and the volume of this space, was based on previous data from a single patient in whom serial plasma concentrations of cisplatin following injection were measured10. We used a fixed cisplatin diffusion coefficient estimated from the molecular weight of the drug according to an empirical law for normal tissues at 37 °C12. The flux of cisplatin to the intracellular space was based on data from head, neck, and gastric carcinomas13,14. Threshold intracellular concentration for cell apoptosis was pulled from immortalized human cell lines (Genomics of Drug Sensitivity in Cancer platform15), and this was used to gauge treatment response according to the RECIST v1.1 criteria16. Finally, the model does not take into account regional heterogeneities of structure and/or biophysical properties within a tumor, such as may arise spontaneously or as a result of prior chemotherapy or radiation therapy. Nevertheless, the model remained consistent in predicting the therapeutic outcomes of the cohort when parameters were modified by an order of magnitude, demonstrating a degree of robustness in the face of parameter value uncertainties. Going forward, as we develop a greater understanding of how these various factors vary within a tumor and between patients, and develop methods for assessing them in individual tumors, we expect the accuracy of our model predictions to improve accordingly.
The computational model itself also has limitations because of the numerous simplifying assumptions that had to be made in its construction, mostly reflective of a lack of biologic data. Thus, even though we have been able to demonstrate promising performance, we recognize that the ability of the model to predict clinical outcomes could probably be improved substantially with better and more detailed information about the tumor microenvironment. As cytotoxic cisplatin levels within a tumor depend heavily on the competition between drug diffusion and clearance, spatial variations in the cisplatin diffusion coefficient and clearance rate likely impact outcomes significantly. The tumor microenvironment is highly heterogeneous, with densely fibrotic areas22 where cisplatin diffusion from sites of injection to distant malignant cells could be particularly limited, yet our model assumes fixed diffusivity throughout the tumor. Relating local diffusivity to Hounsfield unit (HU) values on CT might offer a means of refining the model in this regard, although it remains unclear as to the precise link between diffusivity and HU. Estimating local tissue clearance due to perfusion is even more problematic, as CT does not have the resolution to quantify the microvasculature. In principle, we can estimate global cisplatin clearance from a tumor from the kinetics of cisplatin levels in the blood, as we have shown previously10, but unfortunately these data were not available for the patients of the present study. Accordingly, we used an estimated clearance parameter obtained previously in a single patient, but it is likely that this parameter varies significantly between patients, and indeed can even vary over the course of treatment in single patient due to angiogenesis on the time-scale of weeks.
In conclusion, we have demonstrated that our previously developed computational model of cisplatin pharmacodynamics following intratumoral injection has promise as a treatment planning tool for EBUS-TBNI cisplatin of lung tumors, particularly adenocarcinomas. Even in its current form the model suggests that 1–2 injections might be the optimal choice for tumors smaller than 4 cm3, while for larger tumors we advocate for as many injections as possible subject to tumor accessibility and patient tolerance. Furthermore, we found a power law relationship between the optimal dose to elicit a clinical response and tumor volume, with the exponent of the power law depending on the number of injections.
References
Organization WH. Fact Sheet: Cancer (WHO, 2018).
Howlader, N., Noone, A. & Krapcho, M. SEER Cancer Statistics Review, 1975–2014, Based on November 2016 SEER Data Submission (National Cancer Institute, 2017).
Khan, F., Anker, C.J., Garrison, G., Kinsey, C.M. Endobronchial ultrasound–guided transbronchial needle injection for local control of recurrent non–small cell lung cancer. Ann. Am. Thorac. Soc. 12, 101–104 (2015).
Cappuzzo, F. et al. Erlotinib as maintenance treatment in advanced non-small-cell lung cancer: A multicentre, randomised, placebo-controlled phase 3 study. Lancet Oncol. 11, 521–529 (2010).
Goldberg, E.P., Hadba, A.R., Almond, B.A., Marotta, J.S. Intratumoral cancer chemotherapy and immunotherapy: Opportunities for nonsystemic preoperative drug delivery. J. Pharm. Pharmacol. 54, 159–180 (2002).
Walter, K. A., Tamargo, R. J., Olivi, A., Burger, P. C. & Brem, H. Intratumoral chemotherapy. Neurosurgery 37, 1129–1145 (1995).
Mehta, H. J. et al. Treatment of isolated mediastinal and hilar recurrence of lung cancer with bronchoscopic endobronchial ultrasound guided intratumoral injection of chemotherapy with cisplatin. Lung Cancer 90, 542–547 (2015).
Kinsey, C. M. et al. Tumor density is associated with response to endobronchial ultrasound-guided transbronchial needle injection of cisplatin. J. Thorac. Dis. 12, 4825 (2020).
Celikoglu, F., Celikoglu, S. I. & Goldberg, E. P. Bronchoscopic intratumoral chemotherapy of lung cancer. Lung Cancer 61, 1–12 (2008).
Mori, V., Roy, G. S., Bates, J. H. & Kinsey, C. M. Cisplatin pharmacodynamics following endobronchial ultrasound-guided transbronchial needle injection into lung tumors. Sci. Rep. 9, 1–10 (2019).
Estepar, R., Jose, S., Ross, J., et al. Chest imaging platform: An open-source library and workstation for quantitative chest imaging. Am. J. Respirat. Crit. Care Med. 191, A4975 (2015).
Swabb, E.A., Wei, J., Gullino, P.M. Diffusion and convection in normal and neoplastic tissues. Cancer Res. 34, 2814–2822 (1974).
Troger, V., Fischel, J., Formento, P., Gioanni, J. & Milano, G. Effects of prolonged exposure to cisplatin on cytotoxicity and intracellular drug concentration. Eur. J. Cancer 28, 82–86 (1992).
Kurihara, N., Kubota, T., Hoshiya, Y., et al. Antitumor activity of cis-diamminedichloroplatinum II against human tumor xenografts depends on its area under the curve in nude mice. J. Surg. Oncol. 61, 138–142 (1996).
Yang, W. et al. Genomics of Drug Sensitivity in Cancer (GDSC): A resource for therapeutic biomarker discovery in cancer cells. Nucleic Acids Res. 41, D955–D961 (2012).
Eisenhauer, E. A. et al. New response evaluation criteria in solid tumours: Revised RECIST guideline (version 1.1). Eur. J. Cancer. 45, 228–247 (2009).
Kinsey, C. M., San José Estépar, R., Wei, Y., Washko, G. R. & Christiani, D. C. Regional emphysema of a non-small cell tumor is associated with larger tumors and decreased survival rates. Ann. Am. Thorac. Soc. 12, 1197–1205 (2015).
Celikoglu, S. I., Celikoglu, F. & Goldberg, E. P. Intratumoral cancer chemotherapy through a flexible bronchoscope. J. Bronchol. Intervent. Pulmonol. 11, 260–265 (2004).
Milross, C. G. et al. The effect of tumor size on necrosis and polarographically measured pO2. Acta Oncol. 36, 183–189 (1997).
Forster, J. C., Harriss-Phillips, W. M., Douglass, M. J. & Bezak, E. A review of the development of tumor vasculature and its effects on the tumor microenvironment. Hypoxia 5, 21 (2017).
Caruso, R. et al. Histologic coagulative tumour necrosis as a prognostic indicator of aggressiveness in renal, lung, thyroid and colorectal carcinomas: A brief review. Oncol. Lett. 3, 16–18 (2012).
Salmon, H. et al. Matrix architecture defines the preferential localization and migration of T cells into the stroma of human lung tumors. J. Clin. Investig. 122, 899–910 (2012).
Acknowledgements
Vitor Mori is a Damon Runyon Fellow supported by the Damon Runyon Cancer Research Foundation (DRQ-0620). C. Matthew Kinsey is funded by the University of Vermont Department of Medicine Grant.
Author information
Authors and Affiliations
Contributions
Conception and design: V.M., J.H.T.B., C.M.K. Administrative support: J.H.T.B., C.M.K. Provision of study materials or patients: M.J., H.M., C.M.K. Collection and assembly of data: M.J., H.M., C.M.K. Data analysis and interpretation: V.M., J.H.T.B., C.M.K. Manuscript writing: V.M., J.H.T.B., C.M.K. Final approval of manuscript: V.M., J.H.T.B., M.J., H.M., C.M.K.
Corresponding author
Ethics declarations
Competing interests
JHTB is a consultant for Johnson & Johnson on matters related to this study, and is co-applicant on the patent application “Methods for computational modeling to guide intratumoral therapy.” U.S. Patent Application No. 62/542,623. Filed: August 8, 2017. CMK is a consultant for Olympus America, Johnson and Johnson, and consultant and equity holder for Quantitative Imaging Solutions. He serves on the steering committee for Nuvaira and the Scientific Advisory Board for Galvanize Therapeutics. He reports grants from the NIH, clinical trial funding through Johnson and Johnson, and is a co-applicant on the patent application "Methods for computational modeling to guide intratumoral therapy". U.S. Patent Application No. 62/542,623. Filed: August 8, 2017.” VM, MJ and HM declare that there are no conflict of interest.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Mori, V., Bates, J.H.T., Jantz, M. et al. A computational modeling approach for dosing endoscopic intratumoral chemotherapy for advanced non-small cell lung cancer. Sci Rep 12, 44 (2022). https://doi.org/10.1038/s41598-021-03849-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-021-03849-w
This article is cited by
-
Functional biomaterials for biomimetic 3D in vitro tumor microenvironment modeling
In vitro models (2023)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.