Treatment planning for photodynamic therapy of abscess cavities using patient-specific optical properties measured prior to illumination

Photodynamic therapy (PDT) is an effective antimicrobial therapy that we used to treat human abscess cavities in a Phase 1 clinical trial. This trial included pre-PDT measurements of abscess optical properties, which affect light dose (light fluence) at the abscess wall and PDT response. This study simulated PDT treatment planning for 13 subjects that received optical spectroscopy prior to clinical PDT, to determine the impact of measured optical properties on ability to achieve fluence rate targets in 95% of the abscess wall. Retrospective treatment plans were evaluated for 3 conditions: (1) clinically delivered laser power and assumed, homogeneous optical properties, (2) clinically delivered laser power and measured, homogeneous optical properties, and (3) with patient-specific treatment planning using measured, homogeneous optical properties. Treatment plans modified delivered laser power, intra-cavity Intralipid (scatterer) concentration, and laser fiber type. Using flat-cleaved laser fibers, the proportion of subjects achieving 95% abscess wall coverage decreased significantly relative to assumed optical properties when using measured values for 4 mW cm-2(92% versus 38%,p= 0.01) and 20 mW cm-2(62% versus 15%,p= 0.04) thresholds. When measured optical properties were incorporated into treatment planning, the 4 mW cm-2target was achieved for all cases. After treatment planning, optimal Intralipid concentration across subjects was 0.14 ± 0.09%, whereas 1% was used clinically. Required laser power to achieve the 4 mW cm-2target was significantly correlated with measured abscess wall absorption (ρ= 0.7,p= 0.008), but not abscess surface area (ρ= 0.2,p= 0.53). When using spherical diffuser fibers for illumination, both optimal Intralipid concentration (p= 0.0005) and required laser power (p= 0.0002) decreased compared to flat cleaved fibers. At 0% Intralipid concentration, the 4 mW cm-2target could only be achieved for 69% of subjects for flat-cleaved fibers, compared to 100% for spherical diffusers. Based on large inter-subject variations in optical properties, individualized treatment planning is essential for abscess photodynamic therapy. (Clinical Trial Registration: The parent clinical trial from which these data were acquired is registered on ClinicalTrials.gov as 'Safety and Feasibility Study of Methylene Blue Photodynamic Therapy to Sterilize Deep Tissue Abscess Cavities,' with ClinicalTrials.gov identifier NCT02240498).

Hybrid heterogeneous phantoms for biomedical applications: a demonstration to dosimetry validation

Phantoms simultaneously mimicking anatomical and optical properties of real tissues can play a pivotal role for improving dosimetry algorithms. The aim of the paper is to design and develop a hybrid phantom model that builds up on the strengths of solid and liquid phantoms for mimicking various anatomical structures for prostate cancer photodynamic therapy (PDT) dosimetry validation. The model comprises of a photosensitizer-embedded gelatin lesion within a liquid Intralipid prostate shape that is surrounded by a solid silicone outer shell. The hybrid phantom was well characterized for optical properties. The final assembled phantom was also evaluated for fluorescence tomographic reconstruction in conjunction with SpectraCure's IDOSE software. The developed model can lead to advancements in dosimetric evaluations. This would improve PDT outlook as a clinical treatment modality and boost phantom based standardization of biophotonic devices globally.

Treatment planning for photodynamic therapy of abscess cavities using patient-specific optical properties measured prior to illumination

Background

Photodynamic therapy (PDT) is an effective antimicrobial therapy that we used to treat human abscess cavities in a recently completed Phase 1 clinical trial. This trial included pre-PDT measurements of abscess optical properties, which affect the expected light dose to the abscess wall and eventual PDT response.

Purpose

The objective of this study was to simulate PDT treatment planning for the 13 subjects that received optical spectroscopy prior to clinical abscess PDT. Our goal was to determine the impact of these measured optical properties on our ability to achieve fluence rate targets in 95% of the abscess wall.

Methods

During a Phase 1 clinical trial, 13 subjects received diffuse reflectance spectroscopy prior to PDT in order to determine the optical properties of their abscess wall. Retrospective treatment plans seeking to achieve fluence rate targets in 95% of the abscess wall were evaluated for all subjects for 3 conditions: (1) at the laser power delivered clinically with assumed optical properties, (2) at the laser power delivered clinically with measured optical properties, and (3) with patient-specific treatment planning using these measured optical properties. Factors modified in treatment planning included delivered laser power and intra-cavity Intralipid (scatterer) concentration. The effects of laser fiber type were also simulated.

Results

Using a flat-cleaved laser fiber, the proportion of subjects that achieved 95% abscess wall coverage decreased significantly when incorporating measured optical properties for both the 4 mW/cm 2 (92% vs. 38%, p=0.01) and 20 mW/cm 2 (62% vs. 15%, p=0.04) fluence rate thresholds. However, when measured optical properties were incorporated into treatment planning, a fluence rate of 4 mW/cm 2 was achieved in 95% of the abscess wall for all cases. In treatment planning, the optimal Intralipid concentration across subjects was found to be 0.14 ± 0.09% and the optimal laser power varied from that delivered clinically but with no clear trend (p=0.79). The required laser power to achieve 4 mW/cm 2 in 95% of the abscess wall was significantly correlated with measured µ a at the abscess wall (ρ=0.7, p=0.008), but not abscess surface area (ρ=0.2, p=0.53). When using spherical diffuser fibers as the illumination source, the optimal intralipid concentration decreased to 0.028 ± 0.026% (p=0.0005), and the required laser power decreased also (p=0.0002), compared to flat cleaved fibers. If the intra-cavity lipid emulsion (Intralipid) was replaced with a non-scattering fluid, all subjects could achieve the 4 mW/cm 2 fluence rate threshold in 95% of the abscess wall using a spherical diffuser, while only 69% of subjects could reach the same criterion using a flat cleaved fiber.

Conclusions

The range of optical properties measured in human abscesses reduced coverage of the abscess wall at desirable fluence rates. Patient-specific treatment planning including these measured optical properties could bring the coverage back to desirable levels by altering the Intralipid concentration and delivered optical power. These results motivate a future Phase 2 clinical trial to directly compare the efficacy of patient-specific-treatment planning with fixed doses of Intralipid and light.Clinical Trial Registration: The parent clinical trial from which these data were acquired is registered on ClinicalTrials.gov as "Safety and Feasibility Study of Methylene Blue Photodynamic Therapy to Sterilize Deep Tissue Abscess Cavities," with ClinicalTrials.gov identifier NCT02240498 .

Local monitoring of photosensitizer transient states provides feedback for enhanced efficiency and targeting selectivity in photodynamic therapy

Photodynamic therapy (PDT) fundamentally relies on local generation of PDT precursor states in added photosensitizers (PS), particularly triplet and photo-radical states. Monitoring these states in situ can provide important feedback but is difficult in practice. The states are strongly influenced by local oxygenation, pH and redox conditions, often varying significantly at PDT treatment sites. To overcome this problem, we followed local PDT precursor state populations of PS compounds, via their fluorescence intensity response to systematically varied excitation light modulation. Thereby, we could demonstrate local monitoring of PDT precursor states of methylene blue (MB) and IRdye700DX (IR700), and determined their transitions rates under different oxygenation, pH and redox conditions. By fiber-optics, using one fiber for both excitation and fluorescence detection, the triplet and photo-radical state kinetics of locally applied MB and IR700 could then be monitored in a tissue sample. Finally, potassium iodide and ascorbate were added as possible PDT adjuvants, enhancing intersystem crossing and photoreduction, respectively, and their effects on the PDT precursor states of MB and IR700 could be locally monitored. Taken together, the presented procedure overcomes current methodological limitations and can offer feedback, guiding both excitation and PDT adjuvant application, and thereby more efficient and targeted PDT treatments.

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.

J, M.D.P. DE. C, RAJENDRAM S. Photodynamic Therapy : An Overview and Insights into a Prospective Mainstream Anticancer Therapy

Photodynamic therapy (PDT) procedure has minimum invasiveness in contrast to conventional anticancer surgical procedures. Although clinically approved a few decades ago, it is not commonly used due to its poor efficacy, mainly due to poor light penetration into deeper tissues. PDT uses a photosensitizer (PS), which is photoactivated on illumination by light of appropriate wavelength and oxygen in the tissue, leading to a series of photochemical reactions producing reactive oxygen species (ROS) triggering various mechanisms resulting in lethal effects on tumor cells. This review looks into the fundamental aspects of PDT, such as photochemistry, photobiological effects, and the current clinical applications in the light of improving PDT to become a mainstream therapeutic procedure against a broad spectrum of cancers and malignant lesions. The side effects of PDT, both early and late-onset, are elaborated on in detail to highlight the available options to minimize side effects without compromising therapeutic efficacy. This paper summarizes the benefits, drawbacks, and limitations of photodynamic therapy along with the recent attempts to achieve improved therapeutic efficacy via monitoring various cellular and molecular processes through fluorescent imagery aided by suitable biomarkers, prospective nanotechnology-based targeted delivery methods, the use of scintillating nanoparticles to deliver light to remote locations and also combining PDT with conventional anticancer therapies have opened up new dimensions for PDT in treating cancers. This review inquires and critically analyses prospective avenues in which a breakthrough would finally enable PDT to be integrated into mainstream anticancer therapy.

Scalable and accessible personalized photodynamic therapy optimization with FullMonte and PDT-SPACE

SIGNIFICANCE

Open-source software packages have been extensively used in the past three decades in medical imaging and diagnostics, aiming to study the feasibility of the application ex vivo. Unfortunately, most of the existing open-source tools require some software engineering background to install the prerequisite libraries, choose a suitable computational platform, and combine several software tools to address different applications.

AIM

To facilitate the use of open-source software in medical applications, enabling computational studies of treatment outcomes prior to the complex in-vivo setting.

APPROACH

FullMonteWeb, an open-source, user-friendly web-based software with a graphical user interface for interstitial photodynamic therapy (iPDT) modeling, visualization, and optimization, is introduced. The software can perform Monte Carlo simulations of light propagation in biological tissues, along with iPDT plan optimization. FullMonteWeb installs and runs the required software and libraries on Amazon Web Services (AWS), allowing scalable computing without complex set up.

RESULTS

FullMonteWeb allows simulation of large and small problems on the most appropriate compute hardware, enabling cost improvements of 10  ×   versus always running on a single platform. Case studies in optical property estimation and diffuser placement optimization highlight FullMonteWeb's versatility.

CONCLUSION

The FullMonte open source suite enables easier and more cost-effective in-silico studies for iPDT.

In Vivo Models for Studying Interstitial Photodynamic Therapy of Locally Advanced Cancer

Interstitial photodynamic therapy (I-PDT) is a promising therapy considered for patients with locally advanced cancer. In I-PDT, laser fibers are inserted into the tumor for effective illumination and activation of the photosensitizer in a large tumor. The intratumoral light irradiance and fluence are critical parameters that affect the response to I-PDT. In vivo animal models are required to conduct light dose studies, to define optimal irradiance and fluence for I-PDT. Here we describe two animal models with locally advanced tumors that can be used to evaluate the response to I-PDT. One model is the C3H mouse bearing large subcutaneous SCCVII carcinoma (400-600 mm³). Using this murine model, multiple light regimens with one or two optical fibers with cylindrical diffuser ends (cylindrical diffuser fiber, CDF) can be used to study tumor response to I-PDT. However, tissue heating may occur when 630 nm therapeutic light is delivered through CDF at an intensity ≥60 mW/cm and energy ≥100 J/cm. These thermal effects can impact tumor response while treating locally advanced mice tumors. Magnetic resonance imaging and thermometry can be used to study these thermal effects. A larger animal model, New Zealand White rabbit with VX2 carcinoma (~5000 mm³) implanted in either the sternomastoid (neck implantation model) or the biceps femoris muscle (thigh implantation model), can be used to study I-PDT with image-based pretreatment planning using computed tomography. In the VX2 model, the light delivery can include the use of multiple laser fibers to test light dosimetry and delivery that are relevant for clinical use of I-PDT.

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.

Perspectives on interstitial photodynamic therapy for malignant tumors

SIGNIFICANCE

Despite remarkable advances in the core modalities used in combating cancer, malignant diseases remain the second largest cause of death globally. Interstitial photodynamic therapy (IPDT) has emerged as an alternative approach for the treatment of solid tumors.

AIM

The aim of our study is to outline the advancements in IPDT in recent years and provide our vision for the inclusion of IPDT in standard-of-care (SoC) treatment guidelines of specific malignant diseases.

APPROACH

First, the SoC treatment for solid tumors is described, and the attractive properties of IPDT are presented. Second, the application of IPDT for selected types of tumors is discussed. Finally, future opportunities are considered.

RESULTS

Strong research efforts in academic, clinical, and industrial settings have led to significant improvements in the current implementation of IPDT, and these studies have demonstrated the unique advantages of this modality for the treatment of solid tumors. It is envisioned that further randomized prospective clinical trials and treatment optimization will enable a wide acceptance of IPDT in the clinical community and inclusion in SoC guidelines for well-defined clinical indications.

CONCLUSIONS

The minimally invasive nature of this treatment modality combined with the relatively mild side effects makes IPDT a compelling alternative option for treatment in a number of clinical applications. The adaptability of this technique provides many opportunities to both optimize and personalize the treatment.

Optimizing Interstitial Photodynamic Therapy Planning With Reinforcement Learning-Based Diffuser Placement

Interstitial photodynamic therapy (iPDT) has shown promising results recently as a minimally invasive stand-alone or intra-operative cancer treatment. The development of non-toxic photosensitizing drugs with improved target selectivity has increased its efficacy. However, personalized treatment planning that determines the number of photon emitters, their positions and their input powers while taking into account tissue anatomy and treatment response is still lacking to further improve outcomes.

OBJECTIVE

To develop new algorithms that generate high-quality plans by optimizing over the light source positions, along with their powers, to minimize the damage to organs-at-risk while eradicating the tumor. The optimization algorithms should also accurately model the physics of light propagation through the use of Monte-Carlo simulators.

METHODS

We use simulated-annealing as a baseline algorithm to place the sources. We propose different source perturbations that are likely to provide better outcomes and study their impact. To minimize the number of moves attempted (and effectively runtime) without degrading result quality, we use a reinforcement learning-based method to decide which perturbation strategy to perform in each iteration. We simulate our algorithm on virtual brain tumors modeling real glioblastoma multiforme cases, assuming a 5-ALA PpIX induced photosensitizer that is activated at [Formula: see text] wavelength.

RESULTS

The algorithm generates plans that achieve an average of 46% less damage to organs-as-risk compared to the manual placement used in current clinical studies.

SIGNIFICANCE

Having a general and high-quality planning system makes iPDT more effective and applicable to a wider variety of oncological indications. This paves the way for more clinical trials.

Photodynamic therapy dosimetry using multiexcitation multiemission wavelength: toward real-time prediction of treatment outcome

Evaluating the optical properties of biological tissues is needed to achieve accurate dosimetry during photodynamic therapy (PDT). Currently, accurate assessment of the photosensitizer (PS) concentration by fluorescence measurements during PDT is typically hindered by the lack of information about tissue optical properties. In the present work, a hand-held fiber-optic probe instrument monitoring fluorescence and reflectance is used for assessing blood volume, reduced scattering coefficient, and PS concentration facilitating accurate dosimetry for PDT. System validation was carried out on tissue phantoms using nonlinear least squares support machine regression analysis. It showed a high correlation coefficient (>0.99) in the prediction of the PS concentration upon a large variety of phantom optical properties. In vivo measurements were conducted in a PDT chlorine e6 dose escalating trial involving 36 male Swiss mice with Ehrlich solid tumors in which fluences of 5, 15, and 40  J cm  -  2 were delivered at two fluence rates (100 and 40  mW cm  -  2). Remarkably, quantitative measurement of fluorophore concentration was achieved in the in vivo experiment. Diffuse reflectance spectroscopy (DRS) system was also used to independently measure the physiological properties of the target tissues for result comparisons. Then, blood volume and scattering coefficient measured by the fiber-optic probe system were compared with the corresponding result measured by DRS and showed agreement. Additionally, tumor hemoglobin oxygen saturation was measured using the DRS system. Overall, the system is capable of assessing the implicit photodynamic dose to predict the PDT outcome.

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.

Optimizing interstitial photodynamic therapy with custom cylindrical diffusers

Interstitial photodynamic therapy (iPDT) has shown promise recently as a minimally invasive cancer treatment, partially due to the development of non-toxic photosensitizers in the absence of activation light. However, a major challenge in iPDT is the pre-treatment planning process that specifies the number of diffusers needed, along with their positions and allocated powers, to confine the light distribution to the target volume as much as possible. In this work, a new power allocation algorithm for cylindrical light diffusers including those that can produce customized longitudinal (tailored) emission profiles is introduced. The proposed formulation is convex to guarantee the minimum over-dose possible on the surrounding organs-at-risk. The impact of varying the diffuser lengths and penetration angles on the quality of the plan is evaluated. The results of this study are demonstrated for different photosensitizers activated at different wavelengths and simulated on virtual tumors modeling virtual glioblastoma multiforme cases. Results show that manufacturable cylindrical diffusers with tailored emission profiles can significantly outperform those with conventional flat profiles with an average damage reduction on white matter of 15% to 55% and on gray matter of 23% to 58%.

Solving analytically the simplified spherical harmonics equations in cylindrical turbid media

A methodology is presented for analytically solving simplified spherical harmonics equations (SP_N) in a finite homogeneous absorbing and scattering cylindrical medium. The SP_N equations are a reliable approximation to the radiative transfer equation for describing light propagation inside turbid media. The equations consist of a set of coupled partial differential equations (PDEs). The analytical solution developed here is for a steady-state isotropic point source located at an arbitrary point inside a cylindrical turbid medium. Partial-reflection boundary conditions are considered, as they realistically model the refractive index mismatch between a turbid medium and its surroundings (air), as occurs in practice in biomedical optics. The eigen method is used to decouple the set of SP_N PDEs. The methodology is applied to the SP₃, which has proved to be sufficiently accurate in practice, but it is readily generalizable to higher orders. The solution is compared with the analytical solution of the diffusion equation as well as to gold standard Monte Carlo simulations for validation, against which it shows good agreement. This work is important, as it provides an additional tool for validating numerical solutions of SP_N equations for curved geometries, namely, cylindrical shapes, which are often used in practice.

Automatic interstitial photodynamic therapy planning via convex optimization

Finding a high-quality treatment plan is an essential, yet difficult, stage of Photodynamic therapy (PDT) as it will determine the therapeutic efficacy in eradicating malignant tumors. A high-quality plan is patient-specific, and provides clinicians with the number of fiber-based spherical diffusers, their powers, and their interstitial locations to deliver the required light dose to destroy the tumor while minimizing the damage to surrounding healthy tissues. In this work, we propose a general convex light source power allocation algorithm that, given light source locations, guarantees optimality of the resulting solution in minimizing the over/under-dosage of volumes of interest. Furthermore, we provide an efficient framework for source selection with concomitant power reallocation to achieve treatment plans with a clinically feasible number of sources and comparable quality. We demonstrate our algorithms on virtual test cases that model glioblastoma multiforme tumors, and evaluate the performance of four different photosensitizers with different activation wavelengths and specific tissue uptake ratios. Results show an average reduction of the damage to organs-at-risk (OAR) by 29% to 31% with comparable runtime to existing power allocation techniques.

Reconstruction of a Deformed Tumor Based on Fiducial Marker Registration: A Computational Feasibility Study

Interstitial photodynamic therapy has shown promising results in the treatment of locally advanced head and neck cancer. In this therapy, systemic administration of a light-sensitive drug is followed by insertion of multiple laser fibers to illuminate the tumor and its margins. Image-based pretreatment planning is employed in order to deliver a sufficient light dose to the complex locally advanced head-and-neck cancer anatomy, in order to meet clinical requirements. Unfortunately, the tumor may deform between pretreatment imaging for the purpose of planning and intraoperative imaging when the plan is executed. Tumor deformation may result from the mechanical forces applied by the light fibers and variation of the patient’s posture. Pretreatment planning is frequently done with the assistance of computed tomography or magnetic resonance imaging in an outpatient suite, while treatment monitoring and control typically uses ultrasound imaging due to considerations of costs and availability in the operation room. This article presents a computational method designed to bridge the gap between the 2 imaging events by taking a tumor geometry, reconstructed during preplanning, and by following the displacement of fiducial markers, which are initially placed during the preplanning procedure. The deformed tumor shape is predicted by solving an inverse problem, seeking for the forces that would have resulted in the corresponding fiducial marker displacements. The computational method is studied on spheres of variable sizes and demonstrated on computed tomography reconstructed locally advanced head and neck cancer model. Results of this study demonstrate an average error of less than 1 mm in predicting the deformed tumor shape, where 1 mm is typically the order of uncertainty in distance measurements using magnetic resonance imaging or computed tomography imaging and high-quality ultrasound imaging. This study further demonstrates that the deformed shape can be calculated in a few seconds, making the proposed method clinically relevant.

Advances in antimicrobial photodynamic inactivation at the nanoscale

The alarming worldwide increase in antibiotic resistance amongst microbial pathogens necessitates a search for new antimicrobial techniques, which will not be affected by, or indeed cause resistance themselves. Light-mediated photoinactivation is one such technique that takes advantage of the whole spectrum of light to destroy a broad spectrum of pathogens. Many of these photoinactivation techniques rely on the participation of a diverse range of nanoparticles and nanostructures that have dimensions very similar to the wavelength of light. Photodynamic inactivation relies on the photochemical production of singlet oxygen from photosensitizing dyes (type II pathway) that can benefit remarkably from formulation in nanoparticle-based drug delivery vehicles. Fullerenes are a closed-cage carbon allotrope nanoparticle with a high absorption coefficient and triplet yield. Their photochemistry is highly dependent on microenvironment, and can be type II in organic solvents and type I (hydroxyl radicals) in a biological milieu. Titanium dioxide nanoparticles act as a large band-gap semiconductor that can carry out photo-induced electron transfer under ultraviolet A light and can also produce reactive oxygen species that kill microbial cells. We discuss some recent studies in which quite remarkable potentiation of microbial killing (up to six logs) can be obtained by the addition of simple inorganic salts such as the non-toxic sodium/potassium iodide, bromide, nitrite, and even the toxic sodium azide. Interesting mechanistic insights were obtained to explain this increased killing.

Endobronchial ultrasound-guidance for interstitial photodynamic therapy of locally advanced lung cancer-a new interventional concept

Recent advances in interventional pulmonology led to a significant expansion of the diagnostic and therapeutic role of endobronchial ultrasound. In this paper, we describe a new concept for using endobronchial ultrasound to guide interstitial photodynamic therapy (PDT). For this purpose, we conducted in vitro and in vivo experiments using a phantom and animal models, respectively. A new 0.5 mm optical fiber, with cylindrical diffuser end, was used to deliver the therapeutic light through the 21-gauge endobronchial ultrasound needle. The animal experiments were performed under real-time ultrasonography guidance in mice and rabbits' tumor models. Safe and effective fiber placements and tumor illumination was accomplished. In addition, computer simulation of light propagation suggests that locally advanced lung cancer tumor can be illuminated. This study demonstrates the potential feasibility of this new therapeutic modality approach, justifying further investigation in the treatment of locally advanced lung cancers.

Vascular targeted photodynamic therapy with TOOKAD® Soluble (WST11) in localized prostate cancer: efficiency of automatic pre-treatment planning

Vascular targeted photodynamic therapy (VTP) with WST11 is a novel non-thermal focal treatment for localized prostate cancer that has shown favorable and early efficacy results in previously published studies. In this work, we investigate the efficiency of automatic dosimetric treatment planning. An action model established in a previous study was used in an image-guided optimization scheme to define the personalized optimal light dose for each patient. The calculated light dose is expressed as the number of optical cylindrical fibers to be used, their positions according to an external insertion grid, and the lengths of their diffuser parts. Evaluation of the method was carried out on data collected from 17 patients enrolled in two multi-centric clinical trials. The protocol consisted of comparing the method-simulated necrosis to the result observed on day 7 MR enhanced images. The method performances showed that the final result can be estimated with an accuracy of 10%, corresponding to a margin of 3 mm. In addition, this process was compatible with clinical conditions in terms of calculation times. The overall process took less than 10 min. Different aspects of the VTP procedure were already defined and optimized. Personalized treatment planning definition remained as an issue needing further investigation. The method proposed herein completes the standardization of VTP and opens new pathways for the clinical development of the technique.

Surface markers for guiding cylindrical diffuser fiber insertion in interstitial photodynamic therapy of head and neck cancer

BACKGROUND AND OBJECTIVES

Image-based treatment planning can be used to compute the delivered light dose during interstitial photodynamic therapy (I-PDT) of locally advanced head and neck squamous cell carcinoma (LA-HNSCC). The objectives of this work were to evaluate the use of surface fiducial markers and flexible adhesive grids in guiding interstitial placement of laser fibers, and to quantify the impact of discrepancies in fiber location on the expected light dose volume histograms (DVHs).

METHODS

Seven gel-based phantoms were made to mimic geometries of LA-HNSCC. Clinical flexible grids and fiducial markers were used to guide the insertion of optically transparent catheters, which are used to place cylindrical diffuser fibers within the phantoms. A computed tomography (CT) was used to image the markers and phantoms before and after catheter insertion and to determine the difference between the planned and actual location of the catheters. A finite element method was utilized to compute the light DVHs. Statistical analysis was employed to evaluate the accuracy of fiber placement and to investigate the correlation between the location of the fibers and the calculated DVHs.

RESULTS

There was a statistically significant difference (P = 0.018) between all seven phantoms in terms of the mean displacement. There was also statistically significant correlation between DVHs and depth of insertion (P = 0.0027), but not with the lateral displacement (P = 0.3043). The maximum difference between actual and planned DVH was related to the number of fibers (P = 0.0025) and the treatment time.

CONCLUSIONS

Surface markers and a flexible grid can be used to assist in the administration of a prescribed DVH within 15% of the target dose provided that the treatment fibers are placed within 1.3 cm of the planned depth of insertion in anatomies mimicking LA-HNSCC. The results suggest that the number of cylindrical diffuser fibers and treatment time can impact the delivered DVHs. Lasers Surg. Med. 49:599-608, 2017. © 2017 Wiley Periodicals, Inc.

Realtime Tracking of the Photobleaching Trajectory During Photodynamic Therapy

OBJECTIVE

Photodynamic therapy (PDT) is an alternative treatment for cancer, which involves the administration of a photosensitizing agent that is activated by light at a specific wavelength. This illumination causes after a sequence of photoreactions, the production of reactive oxygen species responsible for the death of the tumor cells but also the degradation of the photosensitizing agent, which then loose the fluorescence properties. The phenomenon is commonly known as the photobleaching process and can be considered as a therapy efficiency indicator.

METHODS

This paper presents the design and validation of a real-time controller able to track a preset photobleaching trajectory by modulating the light impulses width during the treatment sessions.

RESULTS

This innovative solution was validated by in vivo experiments that have shown a significantly improvement of reproducibility of the interindividual photobleaching kinetic.

CONCLUSION

We believe that this approach could lead to personalized PDT modalities.

SIGNIFICANCE

This work may open new perspectives in the control and optimization of photodynamic treatments.

Focal therapy as primary treatment for localized prostate cancer: definition, needs and future

Focal therapy (FT) may offer a promising treatment option in the field of low to intermediate risk localized prostate cancer. The aim of this concept is to combine minimal morbidity with cancer control as well as maintain the possibility of retreatment. Recent advances in MRI and targeted biopsy has improved the diagnostic pathway of prostate cancer and increased the interest in FT. However, before implementation of FT in routine clinical practice, several challenges are still to overcome including patient selection, treatment planning, post-therapy monitoring and definition of oncologic outcome surrogates. In this article, relevant questions regarding the key steps of FT are critically discussed and the main available energy modalities are analyzed taking into account their advantages and unmet needs.

Photodynamic Therapy of Non-Small Cell Lung Cancer. Narrative Review and Future Directions

Photodynamic therapy (PDT) is an established treatment modality for non-small cell lung cancer. Phototoxicity, the primary adverse event, is expected to be minimized with the introduction of new photosensitizers that have shown promising results in phase I and II clinical studies. Early-stage and superficial endobronchial lesions less than 1 cm in thickness can be effectively treated with external light sources. Thicker lesions and peripheral lesions may be amenable to interstitial PDT, where the light is delivered intratumorally. The addition of PDT to standard-of-care surgery and chemotherapy can improve survival and outcomes in patients with pleural disease. Intraoperative PDT has shown promise in the treatment of non-small cell lung cancer with pleural spread. Recent preclinical and clinical data suggest that PDT can increase antitumor immunity. Crosslinking of signal transducer and activator of transcription-3 molecules is a reliable biomarker to quantify the photoreaction induced by PDT. Randomized studies are required to test the prognosis value of this biomarker, obtain approval for the new photosensitizers, and test the potential efficacy of interstitial and intraoperative PDT in the treatment of patients with non-small cell lung cancer.

Revisiting photodynamic therapy dosimetry: reductionist & surrogate approaches to facilitate clinical success

Photodynamic therapy (PDT) can be a highly complex treatment, with many parameters influencing treatment efficacy. The extent to which dosimetry is used to monitor and standardize treatment delivery varies widely, ranging from measurement of a single surrogate marker to comprehensive approaches that aim to measure or estimate as many relevant parameters as possible. Today, most clinical PDT treatments are still administered with little more than application of a prescribed drug dose and timed light delivery, and thus the role of patient-specific dosimetry has not reached widespread clinical adoption. This disconnect is at least partly due to the inherent conflict between the need to measure and understand multiple parameters in vivo in order to optimize treatment, and the need for expedience in the clinic and in the regulatory and commercialization process. Thus, a methodical approach to selecting primary dosimetry metrics is required at each stage of translation of a treatment procedure, moving from complex measurements to understand PDT mechanisms in pre-clinical and early phase I trials, towards the identification and application of essential dose-limiting and/or surrogate measurements in phase II/III trials. If successful, identifying the essential and/or reliable surrogate dosimetry measurements should help facilitate increased adoption of clinical PDT. In this paper, examples of essential dosimetry points and surrogate dosimetry tools that may be implemented in phase II/III trials are discussed. For example, the treatment efficacy as limited by light penetration in interstitial PDT may be predicted by the amount of contrast uptake in CT, and so this could be utilized as a surrogate dosimetry measurement to prescribe light doses based upon pre-treatment contrast. Success of clinical ALA-based skin lesion treatment is predicted almost uniquely by the explicit or implicit measurements of photosensitizer and photobleaching, yet the individualization of treatment based upon each patients measured bleaching needs to be attempted. In the case of ALA, lack of PpIX is more likely an indicator that alternative PpIX production methods must be implemented. Parsimonious dosimetry, using surrogate measurements that are clinically acceptable, might strategically help to advance PDT in a medical world that is increasingly cost and time sensitive. Careful attention to methodologies that can identify and advance the most critical dosimetric measurements, either direct or surrogate, are needed to ensure successful incorporation of PDT into niche clinical procedures.

Online dosimetry for temoporfin-mediated interstitial photodynamic therapy using the canine prostate as model

Online light dosimetry with real-time feedback was applied for temoporfin-mediated interstitial photodynamic therapy (PDT) of dog prostate. The aim was to investigate the performance of online dosimetry by studying the correlation between light dose plans and the tissue response, i.e., extent of induced tissue necrosis and damage to surrounding organs at risk. Light-dose planning software provided dose plans, including light source positions and light doses, based on ultrasound images. A laser instrument provided therapeutic light and dosimetric measurements. The procedure was designed to closely emulate the procedure for whole-prostate PDT in humans with prostate cancer. Nine healthy dogs were subjected to the procedure according to a light-dose escalation plan. About 0.15 mg/kg temoporfin was administered 72 h before the procedure. The results of the procedure were assessed by magnetic resonance imaging, and gross pathology and histopathology of excised tissue. Light dose planning and online dosimetry clearly resulted in more focused effect and less damage to surrounding tissue than interstitial PDT without dosimetry. A light energy dose-response relationship was established where the threshold dose to induce prostate gland necrosis was estimated from 20 to 30  J/cm2.

An image guided treatment platform for prostate cancer photodynamic therapy

This study describes a multimodality images based platform to drive photodynamic therapies of prostate cancer using WST 11 TOOKAD Soluble drug. The platform integrates a pre-treatment planning tool based on magnetic resonance imaging and a per-treatment guidance tool based on transrectal ultrasound images. Evaluation of the platform on clinical data showed that prediction of the therapy outcome was possible with an accuracy of 90 %.

A new finite element approach for near real-time simulation of light propagation in locally advanced head and neck tumors

Background and Objectives

Several clinical studies suggest that interstitial photodynamic therapy (I‐PDT) may benefit patients with locally advanced head and neck cancer (LAHNC). For I‐PDT, the therapeutic light is delivered through optical fibers inserted into the target tumor. The complex anatomy of the head and neck requires careful planning of fiber insertions. Often the fibers' location and tumor optical properties may vary from the original plan therefore pretreatment planning needs near real‐time updating to account for any changes. The purpose of this work was to develop a finite element analysis (FEA) approach for near real‐time simulation of light propagation in LAHNC.

Methods

Our previously developed FEA for modeling light propagation in skin tissue was modified to simulate light propagation from interstitial optical fibers. The modified model was validated by comparing the calculations with measurements in a phantom mimicking tumor optical properties. We investigated the impact of mesh element size and growth rate on the computation time, and defined optimal settings for the FEA. We demonstrated how the optimized FEA can be used for simulating light propagation in two cases of LAHNC amenable to I‐PDT, as proof‐of‐concept.

Results

The modified FEA was in agreement with the measurements (P = 0.0271). The optimal maximum mesh size and growth rate were 0.005–0.02 m and 2–2.5 m/m, respectively. Using these settings the computation time for simulating light propagation in LAHNC was reduced from 25.9 to 3.7 minutes in one case, and 10.1 to 4 minutes in another case. There were minor differences (1.62%, 1.13%) between the radiant exposures calculated with either mesh in both cases.

Conclusions

Our FEA approach can be used to model light propagation from diffused optical fibers in complex heterogeneous geometries representing LAHNC. There is a range of maximum element size (MES) and maximum element growth rate (MEGR) that can be used to minimize the computation time of the FEA to 4 minutes. Lasers Surg. Med. 47:60–67, 2015. © 2015 The Authors. Lasers in Surgery and Medicine Published by Wiley Periodicals, Inc.

Comparison of flat cleaved and cylindrical diffusing fibers as treatment sources for interstitial photodynamic therapy

PURPOSE

For interstitial photodynamic therapy (iPDT) of bulky tumors, careful treatment planning is required in order to ensure that a therapeutic dose is delivered to the tumor, while minimizing damage to surrounding normal tissue. In clinical contexts, iPDT has typically been performed with either flat cleaved or cylindrical diffusing optical fibers as light sources. Here, the authors directly compare these two source geometries in terms of the number of fibers and duration of treatment required to deliver a prescribed light dose to a tumor volume.

METHODS

Treatment planning software for iPDT was developed based on graphics processing unit enhanced Monte Carlo simulations. This software was used to optimize the number of fibers, total energy delivered by each fiber, and the position of individual fibers in order to deliver a target light dose (D90) to 90% of the tumor volume. Treatment plans were developed using both flat cleaved and cylindrical diffusing fibers, based on tissue volumes derived from CT data from a head and neck cancer patient. Plans were created for four cases: fixed energy per fiber, fixed number of fibers, and in cases where both or neither of these factors were fixed.

RESULTS

When the number of source fibers was fixed at eight, treatment plans based on flat cleaved fibers required each to deliver 7180-8080 J in order to deposit 90 J/cm(2) in 90% of the tumor volume. For diffusers, each fiber was required to deliver 2270-2350 J (333-1178 J/cm) in order to achieve this same result. For the case of fibers delivering a fixed 900 J, 13 diffusers or 19 flat cleaved fibers at a spacing of 1 cm were required to deliver the desired dose. With energy per fiber fixed at 2400 J and the number of fibers fixed at eight, diffuser fibers delivered the desired dose to 93% of the tumor volume, while flat cleaved fibers delivered this dose to 79%. With both energy and number of fibers allowed to vary, six diffusers delivering 3485-3600 J were required, compared to ten flat cleaved fibers delivering 2780-3600 J.

CONCLUSIONS

For the same number of fibers, cylindrical diffusers allow for a shorter treatment duration compared to flat cleaved fibers. For the same energy delivered per fiber, diffusers allow for the insertion of fewer fibers in order to deliver the same light dose to a target volume.

Luminescence and photosensitivity of gadolinium labeled hematoporphyrin monomethyl ether

Photodynamic therapy for deep-lying lesions needs an appropriate imaging modality, precise evaluation of tissue oxygen and an effective photosensitizer. Gadolinium based metalloporphyrins Gd(III)-HMME is proposed in this study as a potential multifunctional theranostic agent, as photosensitizer, ratiometric oxygen sensor and MRI contrast agent. The time resolved spectroscopy revealed the luminescence peak of Gd(III)-HMME at 710 and 779 nm with a lifetime of 64 μs in oxygen-free methanol to be phosphorescent. This phosphorescence is strongly dependent on dissolved oxygen concentration. Its intensity in oxygen saturated methanol solution is 21% of that in deoxygenated solution. The singlet oxygen quantum yields ΦΔ of HMME and Gd(III)-HMME in air saturated methanol solution were determined to be 0.79 and 0.40 respectively using comparative spectra method. These phenomena indicate that the oxygen sensibility and production of singlet oxygen of Gd(III)-HMME can fulfill the requirement of PDT treatment.

Monitoring oxygen concentration during photodynamic therapy using prompt photosensitizer fluorescence

A novel technique is described that uses either time-resolved or steady state prompt photosensitizer fluorescence to measure local oxygen concentration. Solution experiments conducted with Al(III) phthalocyanine chloride tetrasulfonic acid confirmed that the steady state fluorescence signal is dependent on the oxygen concentration and fluence rate. A relationship between prompt sensitizer fluorescence and sensitizer triplet quenching efficiency is derived which does not require knowledge of the Stern-Volmer constant. Similar relationships are also derived for sensitizer delayed fluorescence and phosphorescence. An explicit photodynamic therapy (PDT) dose metric that incorporates light dosimetry, sensitizer dosimetry, and triplet quenching efficiency is introduced. All components of this metric can be determined by optical measurements.

Protoporphyrin IX fluorescence and photobleaching during interstitial photodynamic therapy of malignant gliomas for early treatment prognosis

BACKGROUND AND OBJECTIVE

Interstitial photodynamic therapy (iPDT) of non-resectable recurrent glioblastoma using 5-aminolevulinic acid (ALA)-induced protoporphyrin IX (PpIX) has shown a promising outcome. It remained unclear, however, to what extent inter- and intra-tumoural differences of PpIX concentrations influence the efficacy of iPDT. In the current pilot study, we analysed PpIX concentrations quantitatively and assessed PpIX induced fluorescence and photobleaching intraoperatively.

MATERIALS AND METHODS

Five patients harbouring non-resectable glioblastomas were included. ALA (20 or 30 mg/kg body weight) was given 5-8 hours before treatment. Stereotactic biopsies were taken throughout the tumour volume for both histological analysis and determination of PpIX concentrations, which were measured by chemical extraction. Cylindrical light diffusors were stereotactically implanted. Prior to and after irradiation, fluorescence measurements were performed. Outcome measurement was based on clinical and neuro-radiological follow up.

RESULTS

In three patients, a strong PpIX fluorescence was seen before treatment, which was completely photobleached after iPDT. High concentrations of PpIX could be detected in viable tumour parts of these patients (mean PpIX uptake per tumour: 1.4-3.0 µM). In the other two patients, however, no or only low PpIX uptake (0-0.6 µM) could be detected. The patients with strong PpIX uptake showed treatment response and long-term clinical stabilisation (no progression in 29, 30 and 36 months), early treatment failure was seen in the remaining two patients (death after 3 and 9 months).

CONCLUSIONS

Intra-tumoural PpIX concentrations exhibited pronounced inter- and intra-tumoural variations in glioblastoma, which are directly linked to variable degrees of fluorescence intensity. High intra-tumoural PpIX concentrations with strong fluorescence intensity and complete photobleaching after iPDT seem to be associated with favourable outcome. Real-time monitoring of PpIX fluorescence intensity and photobleaching turned out to be feasible and safe and might be employed for early treatment prognosis of iPDT.

Separation of absorption and scattering properties of turbid media using relative spectrally resolved cw radiance measurements

We present a new method for extracting the effective attenuation coefficient and the diffusion coefficient from relative spectrally resolved cw radiance measurements using the diffusion approximation. The method is validated on both simulated and experimental radiance data sets using Intralipid-1% as a test platform. The effective attenuation coefficient is determined from a simple algebraic expression constructed from a ratio of two radiance measurements at two different source-detector separations and the same 90° angle. The diffusion coefficient is determined from another ratio constructed from two radiance measurements at two angles (0° and 180°) and the same source-detector separation. The conditions of the validity of the method as well as possible practical applications are discussed.

Optical property measurements establish the feasibility of photodynamic therapy as a minimally invasive intervention for tumors of the kidney

We measured the optical properties of freshly excised kidneys with renal parenchymal tumors to assess the feasibility of photodynamic therapy (PDT) in these patients. Kidneys were collected from 16 patients during surgical nephrectomies. Spatially resolved, white light, steady-state diffuse reflectance measurements were performed on normal and neoplastic tissue identified by a pathologist. Reflectance data were fit using a radiative transport model to obtain absorption (μa) and transport scattering coefficients (μs'), which define a characteristic light propagation distance, δ. Monte Carlo (MC) simulations of light propagation from cylindrical diffusing fibers were run using the optical properties extracted from each of the kidneys. Interpretable spectra were obtained from 14 kidneys. Optical properties of human renal cancers exhibit significant inter-lesion heterogeneity. For all diagnoses, however, there is a trend toward increased light penetration at longer wavelengths. For renal cell carcinomas (RCC), mean values of δ increase from 1.28 to 2.78 mm as the PDT treatment wavelength is increased from 630 to 780 nm. MC simulations of light propagation from interstitial optical fibers show that fluence distribution in tumors is significantly improved at 780 versus 630 nm. Our results support the feasibility of PDT in selected renal cancer patients, especially with photosensitizers activated at longer wavelengths.

Experimental spectro-angular mapping of light distribution in turbid media

We present a new approach to the analysis of radiance in turbid media. The approach combines data from spectral, angular and spatial domains in a form of spectro-angular maps. Mapping provides a unique way to visualize details of light distribution in turbid media and allows tracking changes with distance. Information content of experimental spectro-angular maps is verified by a direct comparison with simulated data when an analytical solution of the radiative transfer equation is used. The findings deepen our understanding of the light distribution in a homogenous turbid medium and provide a first step toward applying the spectro-angular mapping as a diagnostic tool for tissue characterization.

Temoporfin (Foscan®, 5,10,15,20-tetra(m-hydroxyphenyl)chlorin)--a second-generation photosensitizer

This review traces the development and study of the second-generation photosensitizer 5,10,15,20-tetra(m-hydroxyphenyl)chlorin through to its acceptance and clinical use in modern photodynamic (cancer) therapy. The literature has been covered up to early 2011.

A model to estimate the outcome of prostate cancer photodynamic therapy with TOOKAD Soluble WST11

Interstitial photodynamic therapy is becoming an interesting modality to treat some early stage prostate cancers. A light-sensitive drug is injected to the patient and activated by light using optical fibres inserted inside the prostate. In this work, we were interested in the characterization of the light action model for the WST11 (Tookad® Soluble) drug. A retrospective analysis was performed on results from 28 patients enrolled in phase I and II trials with the WST11 drug. A drug dose of 4 mg/kg patient, dose light of 200 J cm(-1) and wavelength of 753 nm were used. Correlation between the illuminated volume and the obtained necrosis, measured at day 7 MR images, was clearly established. This result suggests that photodynamic therapy planning is possible based on this model.

Role of optical spectroscopy using endogenous contrasts in clinical cancer diagnosis

Optical spectroscopy has been intensively studied for cancer management in the past two decades. This review paper first introduces the background of optical spectroscopy for cancer management, which includes the advantages of optical techniques compared to other established techniques, the principle of optical spectroscopy and the typical setup of instrumentation. Then the recent progress in optical spectroscopy for cancer diagnosis in the following organs is reviewed: the brain, breast, cervix, lung, stomach, colon, prostate and the skin. Reviewed papers were selected from the PubMed database with keywords combining the terms of individual optical spectroscopy techniques and cancers. The primary focus is on the in vivo applications of optical spectroscopy in clinical studies. Ex vivo studies are also included for some organs to highlight special applications or when there are few in vivo results in the literature. Practical considerations of applying optical spectroscopy in clinical settings such as the speed, cost, complexity of operation, accuracy and clinical value are discussed. A few commercially available clinical instruments that are based on optical spectroscopy techniques are presented. Finally several technical challenges and standard issues are discussed and firm conclusions are made.

System for interstitial photodynamic therapy with online dosimetry: first clinical experiences of prostate cancer

The first results from a clinical study for Temoporfin-mediated photodynamic therapy (PDT) of low-grade (T1c) primary prostate cancer using online dosimetry are presented. Dosimetric feedback in real time was applied, for the first time to our knowledge, in interstitial photodynamic therapy. The dosimetry software IDOSE provided dose plans, including optical fiber positions and light doses based on 3-D tissue models generated from ultrasound images. Tissue optical property measurements were obtained using the same fibers used for light delivery. Measurements were taken before, during, and after the treatment session. On the basis of these real-time measured optical properties, the light-dose plan was recalculated. The aim of the treatment was to ablate the entire prostate while minimizing exposure to surrounding organs. The results indicate that online dosimetry based on real-time tissue optical property measurements enabled the light dose to be adapted and optimized. However, histopathological analysis of tissue biopsies taken six months post-PDT treatment showed there were still residual viable cancer cells present in the prostate tissue sections. The authors propose that the incomplete treatment of the prostate tissue could be due to a too low light threshold dose, which was set to 5 J∕cm2.

How tissue optics affect dosimetry of photodynamic therapy

We describe three lessons learned about how tissue optics affect the dosimetry of red to near-infrared treatment light during PDT, based on working with Dr. Tayyaba Hasan. Lesson 1-The optical fluence rate φ near the tissue surface exceeds the delivered irradiance (E). A broad beam penetrates into tissue to a depth (z) as φ=Eke(-μz), with an attenuation constant μ and a backscatter term k. In tissues, k is typically in the range 3-5, and 1∕μ equals δ, the 1∕e optical penetration depth. Lesson 2-Edge losses at the periphery of a uniform treatment beam extend about 3δ from the beam edge. If the beam diameter exceeds 6δ, then there is a central zone of uniform fluence rate in the tissue. Lesson 3-The depth of treatment is linearly proportional to δ (and the melanin content of pigmented epidermis in skin) while proportional to the logarithm of all other factors, such as irradiance, exposure time, or the photosensitizer properties (concentration, extinction coefficient, quantum yield for oxidizing species). The lessons illustrate how tissue optics play a dominant role in specifying the treatment zone during PDT.

Next-generation acceleration and code optimization for light transport in turbid media using GPUs

A highly optimized Monte Carlo (MC) code package for simulating light transport is developed on the latest graphics processing unit (GPU) built for general-purpose computing from NVIDIA - the Fermi GPU. In biomedical optics, the MC method is the gold standard approach for simulating light transport in biological tissue, both due to its accuracy and its flexibility in modelling realistic, heterogeneous tissue geometry in 3-D. However, the widespread use of MC simulations in inverse problems, such as treatment planning for PDT, is limited by their long computation time. Despite its parallel nature, optimizing MC code on the GPU has been shown to be a challenge, particularly when the sharing of simulation result matrices among many parallel threads demands the frequent use of atomic instructions to access the slow GPU global memory. This paper proposes an optimization scheme that utilizes the fast shared memory to resolve the performance bottleneck caused by atomic access, and discusses numerous other optimization techniques needed to harness the full potential of the GPU. Using these techniques, a widely accepted MC code package in biophotonics, called MCML, was successfully accelerated on a Fermi GPU by approximately 600x compared to a state-of-the-art Intel Core i7 CPU. A skin model consisting of 7 layers was used as the standard simulation geometry. To demonstrate the possibility of GPU cluster computing, the same GPU code was executed on four GPUs, showing a linear improvement in performance with an increasing number of GPUs. The GPU-based MCML code package, named GPU-MCML, is compatible with a wide range of graphics cards and is released as an open-source software in two versions: an optimized version tuned for high performance and a simplified version for beginners (http://code.google.com/p/gpumcml).

Photodynamic therapy: superficial and interstitial illumination

Photodynamic therapy (PDT) is reviewed using the treatment of skin tumors as an example of superficial lesions and prostate cancer as an example of deep-lying lesions requiring interstitial intervention. These two applications are among the most commonly studied in oncological PDT, and illustrate well the different challenges facing the two modalities of PDT-superficial and interstitial. They thus serve as good examples to illustrate the entire field of PDT in oncology. PDT is discussed based on the Lund University group's over 20 yr of experience in the field. In particular, the interplay between optical diagnostics and dosimetry and the delivery of the therapeutic light dose are highlighted. An interactive multiple-fiber interstitial procedure to deliver the required therapeutic dose based on the assessment of light fluence rate and sensitizer concentration and oxygen level throughout the tumor is presented.

Photodynamic therapy for focal ablation of the prostate

Although in early stages of clinical development, photodynamic therapy (PDT) shows promise in delivering focal treatment of both primary and post-radiotherapy prostate cancer. This article will review the mechanism of action of PDT, previous research using PDT for treating prostate cancer including the development of newer vascular-acting photosensitizers, and the potential advantages and disadvantages of PDT in delivering focal therapy.