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.

A real-time IR navigation system for pleural photodynamic therapy with a 3D surface acquisition system.

Photodynamic therapy (PDT) has been used intraoperatively to treat patients with malignant pleural mesothelioma. For the efficiency of PDT, it is crucial to deliver light doses uniformly. The current procedure utilizes eight light detectors placed inside the pleural cavity to monitor the light. An updated navigation system, combined with a novel scanning system, is developed to provide real-time guidance for physicians during pleural PDT to improve light delivery. The scanning system consists of two handheld three-dimensional (3D) scanners to capture the pleural cavity's surface topographies quickly and precisely before PDT so that the target surface can be identified for real-time light fluence distribution calculation during PDT. An algorithm is developed to further process the scanned volume to denoise for accurate light fluence calculation and rotate the local coordinate system into any desired direction for a clear visualization during the real-time guidance. The navigation coordinate system is registered to the patient coordinate system utilizing at least three markers to track the light source point position within the pleural cavity throughout the treatment. During PDT, the light source position, the scanned pleural cavity, and the light fluence distribution for the cavity's surface will be displayed in 3D and 2D, respectively. For validation, this novel system is tested using phantom studies with a large chest phantom and 3D-printed lung phantoms of different volumes based on a personal CT scan, immersed in a liquid tissue-simulating phantom with different optical properties, and treated with eight isotropic detectors and the navigation system.

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/cm2 when the irradiance is ≥8.6 mW/cm2 (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.

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

We compare previously reported benzoporphyrin derivative (BPD)-mediated photodynamic therapy (PDT) results for reactive singlet oxygen concentration (also called singlet oxygen dose) on mice with simulations using a computational device, Dosie™, that calculates light transport and photokinetics for PDT in near real-time. The two sets of results are consistent and validate the use of the device in PDT treatment planning to predict BPD-mediated PDT outcomes in mice animal studies based on singlet oxygen dose, which showed a much better correlation with the cure index than the conventional light dose.

Light Fluence Dosimetry in Lung-simulating Cavities.

Accurate light dosimery is critical to ensure consistent outcome for pleural photodynamic therapy (pPDT). Ellipsoid shaped cavities with different sizes surrounded by turbid medium are used to simulate the intracavity lung geometry. An isotropic light source is introduced and surrounded by turbid media. Direct measurements of light fluence rate were compared to Monte Carlo simulated values on the surface of the cavities for various optical properties. The primary component of the light was determined by measurements performed in air in the same geometry. The scattered component was found by submerging the air-filled cavity in scattering media (Intralipid) and absorbent media (ink). The light source was located centrally with the azimuthal angle, but placed in two locations (vertically centered and 2 cm below the center) for measurements. Light fluence rate was measured using isotropic detectors placed at various angles on the ellipsoid surface. The measurements and simulations show that the scattered dose is uniform along the surface of the intracavity ellipsoid geometries in turbid media. One can express the light fluence rate empirically as ϕ =4S/As *Rd/(1 - Rd), where Rd is the diffuse reflectance, As is the surface area, and S is the source power. The measurements agree with this empirical formula to within an uncertainty of 10% for the range of optical properties studied. GPU voxel-based Monte-Carlo simulation is performed to compare with measured results. This empirical formula can be applied to arbitrary geometries, such as the pleural or intraperitoneal cavity.

Interstitial Photodynamic Therapy-A Focused Review.

Multiple clinical studies have shown that interstitial photodynamic therapy (I-PDT) is a promising modality in the treatment of locally-advanced cancerous tumors. However, the utilization of I-PDT has been limited to several centers. The objective of this focused review is to highlight the different approaches employed to administer I-PDT with photosensitizers that are either approved or in clinical studies for the treatment of prostate cancer, pancreatic cancer, head and neck cancer, and brain cancer. Our review suggests that I-PDT is a promising treatment in patients with large-volume or thick tumors. Image-based treatment planning and real-time dosimetry are required to optimize and further advance the utilization of I-PDT. In addition, pre- and post-imaging using computed tomography (CT) with contrast may be utilized to assess the response.

Effects of excited-state absorption on two-photon absorbing materials.

Many chromophores with a large two-photon absorptive cross section are hybrid materials where the two-photon absorption (TPA) is coupled to an excited-state absorption (ESA). We develop a numerical technique to investigate hybrid two-photon processes in nonlinear absorbers. Our numerical method compares well with published results. In addition to customary calculation of the transmission curve, we demonstrate the importance of the ESA following the TPA, which may cause significant temporal and radial distortion. We also show that improvements in the transmission can result in significant radial and temporal pulse distortion, which may actually reduce the material effect.