Clinical PDT dose dosimetry for pleural Photofrin-mediated photodynamic therapy.

Significance: Photodynamic therapy (PDT) is an established cancer treatment utilizing light-activated photosensitizers (PS). Effective treatment hinges on the PDT dose-dependent on PS concentration and light fluence-delivered over time. We introduce an innovative eight-channel PDT dose dosimetry system capable of concurrently measuring light fluence and PS concentration during treatment. Aim: We aim to develop and evaluate an eight-channel PDT dose dosimetry system for simultaneous measurement of light fluence and PS concentration. By addressing uncertainties due to tissue variations, the system enhances accurate PDT dosimetry for improved treatment outcomes. Approach: The study positions eight isotropic detectors strategically within the pleural cavity before PDT. These detectors are linked to bifurcated fibers, distributing signals to both a photodiode and a spectrometer. Calibration techniques are applied to counter tissue-related variations and improve measurement accuracy. The fluorescence signal is normalized using the measured light fluence, compensating for variations in tissue properties. Measurements were taken in 78 sites in the pleural cavities of 20 patients. Results: Observations reveal minimal Photofrin concentration variation during PDT at each site, juxtaposed with significant intra- and inter-patient heterogeneities. Across 78 treated sites in 20 patients, the average Photofrin concentration for all 78 sites is 4.98 µM, with a median concentration of 4.47 µM. The average PDT dose for all 78 sites is 493.17 µMJ/cm2, with a median dose of 442.79 µMJ/cm2. A significant variation in PDT doses is observed, with a maximum difference of 3.1 times among all sites within one patient and a maximum difference of 9.8 times across all patients. Conclusions: The introduced eight-channel PDT dose dosimetry system serves as a valuable real-time monitoring tool for light fluence and PS concentration during PDT. Its ability to mitigate uncertainties arising from tissue properties enhances dosimetry accuracy, thus optimizing treatment outcomes and bolstering the effectiveness of PDT in cancer therapy.

Fractionated Photofrin-Mediated Photodynamic Therapy Significantly Improves Long-Term Survival.

This study investigates the effect of fractionated (two-part) PDT on the long-term local control rate (LCR) using the concentration of reactive oxygen species ([ROS]rx) as a dosimetry quantity. Groups with different fractionation schemes are examined, including a 2 h interval between light delivery sessions to cumulative fluences of 135, 180, and 225 J/cm2. While the total treatment time remains constant within each group, the division of treatment time between the first and second fractionations are explored to assess the impact on long-term survival at 90 days. In all preclinical studies, Photofrin is intravenously administered to mice at a concentration of 5 mg/kg, with an incubation period between 18 and 24 h before the first light delivery session. Fluence rate is fixed at 75 mW/cm2. Treatment ensues via a collimated laser beam, 1 cm in diameter, emitting light at 630 nm. Dosimetric quantities are assessed for all groups along with long-term (90 days) treatment outcomes. This study demonstrated a significant improvement in long-term survival after fractionated treatment schemes compared to single-fraction treatment, with the optimal 90-day survival increasing to 63%, 86%, and 100% vs. 20%, 25%, and 50%, respectively, for the three cumulative fluences. The threshold [ROS]rx for the optimal scheme of fractionated Photofrin-mediated PDT, set at 0.78 mM, is significantly lower than that for the single-fraction PDT, at 1.08 mM.

Real-time PDT Dose Dosimetry for Pleural Photodynamic Therapy.

PDT dose is the product of the photosensitizer concentration and the light fluence in the target tissue. For improved dosimetry during plural photodynamic therapy (PDT), an eight-channel PDT dose dosimeter was developed to measure both the light fluence and the photosensitizer concentration simultaneously from eight different sites in the pleural cavity during PDT. An isotropic detector with bifurcated fibers was used for each channel to ensure detected light was split equally to the photodiode and spectrometer. The light fluence rate distribution is monitored using an IR navigation system. The navigation system allows 2D light fluence mapping throughout the whole pleural cavity rather than just the selected points. The fluorescence signal is normalized by the light fluence measured at treatment wavelength. We have shown that the absolute photosensitizer concentration can be obtained by applying optical properties correction and linear spectral fitting to the measured fluorescence data. The detection limit and the optical property correction factor of each channel were determined and validated using tissue-simulating phantoms with known varying concentration of Photofrin. Tissue optical properties are determined using an absorption spectroscopy probe immediately before PDT at the same sites. The combination of 8-channel PDT dosimeter system and IR navigation system, which can calculate light fluence rate in the pleural cavity in real-time, providing a mean to determine the distribution of PDT dose on the entire pleural cavity to investigate the heterogeneity of PDT dose on the pleural cavity.

Monte Carlo simulation of Cherenkov imaging for Total Skin Electron Treatment with CT DICOM realistic patient geometry.

This Monte Carlo (MC) simulation study provides an evaluation of dose uniformity in a patient and the difference between dose and Cherenkov distributions, which is invaluable in developing conversion factors to relate observed Cherenkov images to actual dose distributions for TSET patients. This MC simulations with TOPAS is performed using realistic patient geometries obtained with a 3D scanner during total skin electron treatments (TSET) at UPenn. For each treatment posture in the Stanford technique, the differences between Cherenkov photon distributions and dose distributions produced in MC are consistent with the differences observed between a Cherenkov imaging camera and in-vivo dose measurement with OSLD on patient skin. According to MC studies of a flat rectangular PVC board, the difference between Cherenkov and dose is mostly due to the spoiler. This is confirmed by observing consistent dose and Cherenkov distributions in clinical measurements on a PVC board without the spoiler. The accumulated dose and Cherenkov distributions for each patient are obtained by projecting the MC output of the 6 postures of the TSET treatment together onto a finite element model of the patient.

Cherenkov imaging for Total Skin Electron Therapy - an evaluation of dose uniformity.

Total Skin Electron Therapy (TSET) utilizes high-energy electrons to treat cancers on the entire body surface. The otherwise invisible radiation beam can be observed via the optical Cherenkov photons emitted from interaction between the high-energy electron beam and tissue. Cherenkov emission can be used to evaluate the dose uniformity on the surface of the patient in real-time using a time-gated intensified camera system. Each patient was monitored during TSET by in-vivo detectors (IVD) as well as Scintillators. Patients undergoing TSET in various conditions (whole body and half body) were imaged and analyzed. A rigorous methodology for converting Cherenkov intensity to surface dose as products of correction factors, including camera vignette correction factor, incident radiation correction factor, and tissue optical properties correction factor. A comprehensive study has been carried out by inspecting various positions on the patients such as vertex, chest, perineum, shins, and foot relative to the umbilicus point (the prescription point).

Reactive oxygen species explicit dosimetry to predict local tumor growth for Photofrin-mediated photodynamic therapy.

Although photodynamic therapy (PDT) is an established modality for cancer treatment, current dosimetric quantities, such as light fluence and PDT dose, do not account for the differences in PDT oxygen consumption for different fluence rates (ϕ). A macroscopic model was adopted to calculate reactive oxygen species concentration ([ROS]rx) to predict Photofrin-PDT outcome in mice bearing radiation-induced fibrosarcoma (RIF) tumors. Singlet oxygen is the primary cytotoxic species for ROS, which is responsible for cell death in type II PDT, although other type I ROS is included in the parameters used in our model. Using a combination of fluences (50-250 J∕cm2) and ϕ (75 or 150 mW∕cm2), tumor regrowth rate, k, was determined for each condition by fitting the tumor volume versus time to V0 *exp(k*t). Treatment was delivered with a collimated laser beam of 1 cm diameter at 630 nm. Explicit dosimetry of light fluence rate on tissue surface, tissue oxygen concentration, tissue optical properties, and Photofrin concentration were performed. Light fluence rate at 3 mm depth (ϕ 3mm) was determined for the treatment volume based on Monte-Carlo simulations and measured tissue optical properties. Initial tissue oxygenation [3 O 2]0 was measured by an Oxylite oxygen probe before PDT and used to calculate [ROS]rx,calc. This value was compared to [ROS]rx,meas as calculated with the entire tissue oxygen spectrum [3 O 2](t), measured over the duration of light delivery for PDT. Cure index, CI = 1-k/kctr , for tumor growth up to 14 days after PDT was predicted by four dose metrics: light fluence, PDT dose, and [ROS]rx,calc, and [ROS]rx,meas. PDT dose was defined as the product of the time-integral of photosensitizer concentration and ϕ at a 3 mm tumor depth. These studies show that [ROS]rx,meas best correlates with CI and is an effective dosimetric quantity that can predict treatment outcome.

Evaluation of Light Fluence Distribution Using an IR Navigation System for HPPH-mediated Pleural Photodynamic Therapy (pPDT).

Uniform light fluence distribution for patients undergoing photodynamic therapy (PDT) is critical to ensure predictable PDT outcomes. However, current practice when delivering intrapleural PDT uses a point source to deliver light that is monitored by seven isotropic detectors placed within the pleural cavity to assess its uniformity. We have developed a real-time infrared (IR) tracking camera to follow the movement of the light point source and the surface contour of the treatment area. The calculated light fluence rates were matched with isotropic detectors using a two-correction factor method and an empirical model that includes both direct and scattered light components. Our clinical trial demonstrated that we can successfully implement the IR navigation system in 75% (15/20) of the patients. Data were successfully analyzed in 80% (12/15) patients because detector locations were not available for three patients. We conclude that it is feasible to use an IR camera-based system to track the motion of the light source during PDT and demonstrate its use to quantify the uniformity of light distribution, which deviated by a standard deviation of 18% from the prescribed light dose. The navigation system will fail when insufficient percentage of light source positions is obtained (<30%) during PDT.

Validation of tissue optical properties measurement using diffuse reflectance spectroscopy (DRS).

The effectiveness of photodynamic treatment depends on several factors including an accurate knowledge of optical properties of the tissue to be treated. Transmittance and diffuse reflectance spectroscopic techniques are commonly used to determine tissue optical properties. Although transmittance spectroscopy technique is accurate in determining tissue optical properties, it is only valid in an infinite medium and can only be used for interstitial measurements. Diffuse reflectance spectroscopy, on the other hand, is easily adapted to most tissue geometries including skin measurements that involve semi-infinte medium. However, the accuracy of the measured optical properties can be affected by uncertainty in the measurements themselves and/or due to the uncertainty in the fitting algorithm. In this study, we evaluate the accuracy of optical properties determination using diffuse reflectance spectroscopy implemented using a contact probe setup. We characterized the error of the optical properties fitted using two fitting algorithms, a wavelength wise fitting algorithm and a full reflectance spectral fitting algorithm. By conducting systematic investigation of the measurements and fitting algorithm of DRS, we gained an understanding of the uncertainties in the measured optical properties and outlined improvement measures to minimize these errors.

A summary of light dose distribution using an IR navigation system for Photofrin-mediated Pleural PDT.

Uniform delivery of light fluence is an important goal for photodynamic therapy. We present summary results for an infrared (IR) navigation system to deliver light dose uniformly during intracavitory PDT by tracking the movement of the light source and providing real-time feedback on the light fluence rate on the entire cavity surface area. In the current intrapleural PDT protocol, 8 detectors placed in selected locations in the pleural cavity monitor the light doses. To improve the delivery of light dose uniformity, an IR camera system is used to track the motion of the light source as well as the surface contour of the pleural cavity. A MATLAB-based GUI program is developed to display the light dose in real-time during PDT to guide the PDT treatment delivery to improve the uniformity of the light dose. A dualcorrection algorithm is used to improve the agreement between calculations and in-situ measurements. A comprehensive analysis of the distribution of light fluence during PDT is presented in both phantom conditions and in clinical cases.

An improved analytic function for predicting light fluence rate in circular fields on a semi-infinite geometry.

Accurate determination of in-vivo light fluence rate is critical for preclinical and clinical studies involving photodynamic therapy (PDT). This study compares the longitudinal light fluence distribution inside biological tissue in the central axis of a 1 cm diameter circular uniform light field for a range of in-vivo tissue optical properties (absorption coefficients (µa) between 0.01 and 1 cm-1 and reduced scattering coefficients (µs') between 2 and 40 cm-1). This was done using Monte-Carlo simulations for a semi-infinite turbid medium in an air-tissue interface. The end goal is to develop an analytical expression that would fit the results from the Monte Carlo simulation for both the 1 cm diameter circular beam and the broad beam. Each of these parameters is expressed as a function of tissue optical properties. These results can then be compared against the existing expressions in the literature for broad beam for analysis in both accuracy and applicable range. Using the 6-parameter model, the range and accuracy for light transport through biological tissue is improved and may be used in the future as a guide in PDT for light fluence distribution for known tissue optical properties.