Validation of combined Monte Carlo and photokinetic simulations for the outcome correlation analysis of benzoporphyrin derivative-mediated photodynamic therapy on mice

Computational Optimization of Irradiance and Fluence for Interstitial Photodynamic Therapy Treatment of Patients with Malignant Central Airway Obstruction

There are no effective treatments for patients with extrinsic malignant central airway obstruction (MCAO). In a recent clinical study, we demonstrated that interstitial photodynamic therapy (I-PDT) is a safe and potentially effective treatment for patients with extrinsic MCAO. In previous preclinical studies, we reported that a minimum light irradiance and fluence should be maintained within a significant volume of the target tumor to obtain an effective PDT response. In this paper, we present a computational approach to personalized treatment planning of light delivery in I-PDT that simultaneously optimizes the delivered irradiance and fluence using finite element method (FEM) solvers of either Comsol Multiphysics^® or Dosie™ for light propagation. The FEM simulations were validated with light dosimetry measurements in a solid phantom with tissue-like optical properties. The agreement between the treatment plans generated by two FEMs was tested using typical imaging data from four patients with extrinsic MCAO treated with I-PDT. The concordance correlation coefficient (CCC) and its 95% confidence interval (95% CI) were used to test the agreement between the simulation results and measurements, and between the two FEMs treatment plans. Dosie with CCC = 0.994 (95% CI, 0.953-0.996) and Comsol with CCC = 0.999 (95% CI, 0.985-0.999) showed excellent agreement with light measurements in the phantom. The CCC analysis showed very good agreement between Comsol and Dosie treatment plans for irradiance (95% CI, CCC: 0.996-0.999) and fluence (95% CI, CCC: 0.916-0.987) in using patients' data. In previous preclinical work, we demonstrated that effective I-PDT is associated with a computed light dose of ≥45 J/cm² when the irradiance is ≥8.6 mW/cm² (i.e., the effective rate-based light dose). In this paper, we show how to use Comsol and Dosie packages to optimize rate-based light dose, and we present Dosie's newly developed domination sub-maps method to improve the planning of the delivery of the effective rate-based light dose. We conclude that image-based treatment planning using Comsol or Dosie FEM-solvers is a valid approach to guide the light dosimetry in I-PDT of patients with MCAO.

Integrating clinical access limitations into iPDT treatment planning with PDT-SPACE

PDT-SPACE is an open-source software tool that automates interstitial photodynamic therapy treatment planning by providing patient-specific placement of light sources to destroy a tumor while minimizing healthy tissue damage. This work extends PDT-SPACE in two ways. The first enhancement allows specification of clinical access constraints on light source insertion to avoid penetrating critical structures and to minimize surgical complexity. Constraining fiber access to a single burr hole of adequate size increases healthy tissue damage by 10%. The second enhancement generates an initial placement of light sources as a starting point for refinement, rather than requiring entry of a starting solution by the clinician. This feature improves productivity and also leads to solutions with 4.5% less healthy tissue damage. The two features are used in concert to perform simulations of various surgery options of virtual glioblastoma multiforme brain tumors.

A Monte Carlo simulation for Moving Light Source in Intracavity PDT

We developed a simulation method for modeling the light fluence delivery in intracavity Photodynamic Therapy (icav-PDT) for pleural lung cancer using a moving light source. Due to the large surface area of the pleural lung cavity, the light source needs to be moved to deliver a uniform dose around the entire cavity. While multiple fixed detectors are used for dosimetry at a few locations, an accurate simulation of light fluence and fluence rate is still needed for the rest of the cavity. We extended an existing Monte Carlo (MC) based light propagation solver to support moving light sources by densely sampling the continuous light source trajectory and assigning the proper number of photon packages launched along the way. The performance of Simphotek GPU CUDA-based implementation of the method - PEDSy-MC - has been demonstrated on a life-size lung-shaped phantom, custom printed for testing icav-PDT navigation system at the Perlman School of Medicine (PSM) - calculations completed under a minute (for some cases) and within minutes have been achieved. We demonstrate results within a 5% error of the analytic solution for multiple detectors in the phantom. PEDSy-MC is accompanied by a dose-cavity visualization tool that allows real-time inspection of dose values of the treated cavity in 2D and 3D, which will be expanded to ongoing clinical trials at PSM. PSM has developed a technology to measure 8-detectors in a pleural cavity phantom using Photofrin-mediated PDT that has been used during validation.

Photodynamic therapy in a pleural cavity using monte carlo simulations with 2D/3D Graphical Visualization

Cancer therapy using Photodynamic Therapy (PDT) has been investigated for some time [1,2] and now it is a growing area of interest in clinical trials [3]. Monte Carlo (MC) simulations were used for early laboratory studies [4,5] for analysis in PDT. Various improvements in the MC method have advanced the field in recent years.

In Silico, Combined Plasmonic Photothermal and Photodynamic Therapy in Mice

Plasmonic photothermal and photodynamic therapy (PPTT and PDT, respectively) are two cancer treatments that have the potential to be combined in a synergistic scheme. The aim of this study is to optimize the PPTT treatment part, in order to account for the PDT lack of coverage in the hypoxic tumor volume and in cancer areas laying in deep sites. For the needs of this study, a mouse was modeled, subjected to PDT and its necrotic area was estimated by using the MATLAB software. The same procedure was repeated for PPTT, using COMSOL Multiphysics. PPTT treatment parameters, namely laser power and irradiation time, were optimized in order to achieve the optimum therapeutic effect of the combined scheme. The PDT alone resulted in 54.8% tumor necrosis, covering the upper cancer layers. When the PPTT was also applied, the total necrosis percentage raised up to 99.3%, while all of the surrounding studied organs (skin, heart, lungs and trachea, ribs, liver and spleen) were spared. The optimized values of the PPTT parameters were 550 mW of laser power and 70 s of irradiation time. Hence, the PPTT–PDT combination shows great potential in achieving high levels of tumor necrosis while sparing the healthy tissues.

Machine learning for real-time optical property recovery in interstitial photodynamic therapy: a stimulation-based study

With the continued development of non-toxic photosensitizer drugs, interstitial photodynamic therapy (iPDT) is showing more favorable outcomes in recent clinical trials. IPDT planning is crucial to further increase the treatment efficacy. However, it remains a major challenge to generate a high-quality, patient-specific plan due to uncertainty in tissue optical properties (OPs), µ _a and µ _s . These parameters govern how light propagates inside tissues, and any deviation from the planning-assumed values during treatment could significantly affect the treatment outcome. In this work, we increase the robustness of iPDT against OP variations by using machine learning models to recover the patient-specific OPs from light dosimetry measurements and then re-optimizing the diffusers' optical powers to adapt to these OPs in real time. Simulations on virtual brain tumor models show that reoptimizing the power allocation with the recovered OPs significantly reduces uncertainty in the predicted light dosimetry for all tissues involved.

FullMonteCUDA: a fast, flexible, and accurate GPU-accelerated Monte Carlo simulator for light propagation in turbid media

Optimizing light delivery for photodynamic therapy, quantifying tissue optical properties or reconstructing 3D distributions of sources in bioluminescence imaging and absorbers in diffuse optical imaging all involve solving an inverse problem. This can require thousands of forward light propagation simulations to determine the parameters to optimize treatment, image tissue or quantify tissue optical properties, which is time-consuming and computationally expensive. Addressing this problem requires a light propagation simulator that produces results quickly given modelling parameters. In previous work, we developed FullMonteSW: currently the fastest, tetrahedral-mesh, Monte Carlo light propagation simulator written in software. Additional software optimizations showed diminishing performance improvements, so we investigated hardware acceleration methods. This work focuses on FullMonteCUDA: a GPU-accelerated version of FullMonteSW which targets NVIDIA GPUs. FullMonteCUDA has been validated across several benchmark models and, through various GPU-specific optimizations, achieves a 288-936x speedup over the single-threaded, non-vectorized version of FullMonteSW and a 4-13x speedup over the highly optimized, hand-vectorized and multi-threaded version. The increase in performance allows inverse problems to be solved more efficiently and effectively.