Monte Carlo modelling of daylight activated photodynamic therapy

Simulating photodynamic therapy for the treatment of glioblastoma using Monte Carlo radiative transport

Significance

Glioblastoma (GBM) is a rare but deadly form of brain tumor with a low median survival rate of 14.6 months, due to its resistance to treatment. An independent simulation of the INtraoperative photoDYnamic therapy for GliOblastoma (INDYGO) trial, a clinical trial aiming to treat the GBM resection cavity with photodynamic therapy (PDT) via a laser coupled balloon device, is demonstrated.

Aim

To develop a framework providing increased understanding for the PDT treatment, its parameters, and their impact on the clinical outcome.

Approach

We use Monte Carlo radiative transport techniques within a computational brain model containing a GBM to simulate light path and PDT effects. Treatment parameters (laser power, photosensitizer concentration, and irradiation time) are considered, as well as PDT's impact on brain tissue temperature.

Results

The simulation suggests that 39% of post-resection GBM cells are killed at the end of treatment when using the standard INDYGO trial protocol (light fluence = 200    J / cm 2 at balloon wall) and assuming an initial photosensitizer concentration of 5    μ M . Increases in treatment time and light power (light fluence = 400    J / cm 2 at balloon wall) result in further cell kill but increase brain cell temperature, which potentially affects treatment safety. Increasing the p hotosensitizer concentration produces the most significant increase in cell kill, with 61% of GBM cells killed when doubling concentration to 10    μ M and keeping the treatment time and power the same. According to these simulations, the standard trial protocol is reasonably well optimized with improvements in cell kill difficult to achieve without potentially dangerous increases in temperature. To improve treatment outcome, focus should be placed on improving the photosensitizer.

Conclusions

With further development and optimization, the simulation could have potential clinical benefit and be used to help plan and optimize intraoperative PDT treatment for GBM.

A singlet state oxygen generation model based on the Monte Carlo method of visible antibacterial blue light inactivation

Visible antibacterial blue light (VABL) has received much attention recently as a nondestructive inactivation approach. However, due to the sparse distribution of bacteria, the light energy evaluation method used in existing studies is inaccurate. Thus, the sensitivity of microorganisms to VABL in different experiments cannot be compared. In this paper, a Monte Carlo-based photon transport model with the optimized scattering phase function was constructed. The model calculated the spatial light energy distribution and the temporal distribution of cumulative singlet state oxygen (CSO) under various cell and medium parameters. The simulation results show that when the cells are sparsely distributed, <30% of light energy from the light source is absorbed by microbes and participates in photochemical reactions. The CSO produced increases with cell density and cell size. Little light energy is available, and thus, the concentration of CSO produced is insufficient to inactivate microbes at deeper depths. As the light intensity and inactivation time increased, the production of singlet state oxygen tended to level off. The model proposed here can quantify the generation of singlet state oxygen and provide a more accurate light energy guide for the VABL inactivation process.

Meshless Monte Carlo radiation transfer method for curved geometries using signed distance functions

SIGNIFICANCE

Monte Carlo radiation transfer (MCRT) is the gold standard for modeling light transport in turbid media. Typical MCRT models use voxels or meshes to approximate experimental geometry. A voxel-based geometry does not allow for the precise modeling of smooth curved surfaces, such as may be found in biological systems or food and drink packaging. Mesh-based geometry allows arbitrary complex shapes with smooth curved surfaces to be modeled. However, mesh-based models also suffer from issues such as the computational cost of generating meshes and inaccuracies in how meshes handle reflections and refractions.

AIM

We present our algorithm, which we term signedMCRT (sMCRT), a geometry-based method that uses signed distance functions (SDF) to represent the geometry of the model. SDFs are capable of modeling smooth curved surfaces precisely while also modeling complex geometries.

APPROACH

We show that using SDFs to represent the problem's geometry is more precise than voxel and mesh-based methods.

RESULTS

sMCRT is validated against theoretical expressions, and voxel and mesh-based MCRT codes. We show that sMCRT can precisely model arbitrary complex geometries such as microvascular vessel network using SDFs. In comparison with the current state-of-the-art in MCRT methods specifically for curved surfaces, sMCRT is more precise for cases where the geometry can be defined using combinations of shapes.

CONCLUSIONS

We believe that SDF-based MCRT models are a complementary method to voxel and mesh models in terms of being able to model complex geometries and accurately treat curved surfaces, with a focus on precise simulation of reflections and refractions. sMCRT is publicly available at https://github.com/lewisfish/signedMCRT.

An open-label prospective study to assess short incubation time white LED light photodynamic therapy in the treatment of superficial basal cell carcinoma

Artificial white LED light photodynamic therapy (awl-PDT) is an effective, pain-free treatment for actinic keratosis. The efficacy of awl-PDT in the treatment of superficial basal cell carcinoma (sBCC) has not been assessed. Patients with histologically confirmed sBCC underwent two treatments of awl-PDT 1 week apart. Lesions were incubated with methyl 5-aminolaevulinic acid for 30 min and then illuminated using the Maquet Power LED 500 theatre light (405-800nm, 140 000 lux) to deliver an equivalent red light dose of 75 J/cm² at a rate of 55 mW/cm² . Pain was measured using a visual analogue scale during treatment. Clinical response was assessed at day 28. Follow-up continued 3 months for 1 year. Cosmetic outcome was assessed at 3 months and 1 year. Twenty-eight patients with 36 lesions and a mean age of 63.64 (SD 2.62) were recruited. The median lesion size was 15 mm (IQR 8.75). The response rate at day 28 was 100%. Recurrence rates were 3/36 (8.3%) at 3 months, 6/36 (16.7%) at 6 months, 10/36 (27.8%) at 9 months and 11/36 (30.6%) at 1 year. Median pain scores were 0/100 (IQR 0) and 0/100 (IQR 5) during treatments one and two, respectively. Cosmetic outcome was excellent or good in the majority of cases. Although initially effective for sBCC at 28 days, 30.6% of lesions recurred 1 year after awl-PDT. Pain scores were negligible, and the cosmetic outcome was favourable. Further head-to-head studies with optimised protocols are required to determine if awl-PDT has a role in the treatment of sBCC.

Monte Carlo simulations of photodynamic therapy in human blood model

This study aims to simulate a therapeutic plan for a normal human blood model under various patho-physiological conditions, such as the development of leukemia/blood diseases, by means of Monte Carlo multilayered simulation. The photosensitizing compound selectively accumulates in the target cells. A superficial treatment of a blood sample was performed at different ratios of oxygen saturation ([Formula: see text]) under the concentration ([Formula: see text] = 30 µM) effect of merocyanine 540 (MC540) in the blood irradiation. This was done under the application of visible light of wavelength ~ [Formula: see text] at an exposure time ~ 60 s. The dose of photodynamic therapy (PDT) was evaluated for the biological damage, leading to necrosis and blood damage during the treatment. In addition, the effect of PDT treatment response in the blood is related to hemoglobin oxygen saturation, resulting in an excellent relationship between the changes caused by the treatment in the blood at a peculiar oxygen saturation rate (for the highest response: [Formula: see text] 50%) and a light dose (LD) of 3.83 [Formula: see text] above the minimal toxicity of normal tissues. The photodynamic dose is related to the depth of necrosis and the time of treatment for the achievement of the LD delivery at the PDT of blood.

Development of a Predictive Monte Carlo Radiative Transfer Model for Ablative Fractional Skin Lasers

It is possible to enhance topical drug delivery by pretreatment of the skin with ablative fractional lasers (AFLs). However, the parameters to use for a given AFL to achieve the desired depth of ablation or the desired therapeutic or cosmetic outcome are hard to predict. This leaves open the real possibility of overapplication or underapplication of laser energy to the skin. In this study, we developed a numerical model consisting of a Monte Carlo radiative transfer (MCRT) code coupled to a heat transfer and tissue damage algorithm. The simulation is designed to predict the depth effects of AFL on the skin, verified with in vitro experiments in porcine skin via optical coherence tomography (OCT) imaging. Ex vivo porcine skin is irradiated with increasing energies (50-400 mJ/pixel) from a CO₂ AFL. The depth of microscopic treatment zones is measured and compared with our numerical model. The data from the OCT images and MCRT model complement each other well. Nonablative thermal effects on surrounding tissue are also discussed. This model, therefore, provides an initial step toward a predictive determination of the effects of AFL on the skin. Lasers Surg. Med. © 2020 The Authors. Lasers in Surgery and Medicine published by Wiley Periodicals LLC.

Measuring Daylight: A Review of Dosimetry in Daylight Photodynamic Therapy

Successful daylight photodynamic therapy (DPDT) relies on the interaction of light, photosensitisers and oxygen. Therefore, the 'dose' of light that a patient receives during treatment is a clinically relevant quantity, with a minimum dose for effective treatment recommended in the literature. However, there are many different light measurement methods used in the published literature, which may lead to confusion surrounding reliable and traceable dose measurement in DPDT, and what the most appropriate method of light measurement in DPDT might be. Furthermore, for the majority of practitioners who do not carry out any formal dosimetry and for the patients receiving DPDT, building confidence in the evidence supporting this important treatment option is of key importance. This review seeks to clarify the methodology of DPDT and discusses the literature relating to DPDT dosimetry.

Monte Carlo Simulations of Heat Deposition During Photothermal Skin Cancer Therapy Using Nanoparticles

Photothermal therapy using nanoparticles is a promising new approach for the treatment of cancer. The principle is to utilise plasmonic nanoparticle light interaction for efficient heat conversion. However, there are many hurdles to overcome before it can be accepted in clinical practice. One issue is a current poor characterization of the thermal dose that is distributed over the tumour region and the surrounding normal tissue. Here, we use Monte Carlo simulations of photon radiative transfer through tissue and subsequent heat diffusion calculations, to model the spatial thermal dose in a skin cancer model. We validate our heat rise simulations against experimental data from the literature and estimate the concentration of nanorods in the tumor that are associated with the heat rise. We use the cumulative equivalent minutes at 43 °C (CEM43) metric to analyse the percentage cell kill across the tumour and the surrounding normal tissue. Overall, we show that computer simulations of photothermal therapy are an invaluable tool to fully characterize thermal dose within tumour and normal tissue.

Modeling PpIX effective light fluence at depths into the skin for PDT dose comparison

BACKGROUND

Daylight-activated PDT has seen increased support in recent years as a treatment method for actinic keratosis and other non-melanoma skin cancers. The inherent variability observed in broad-spectrum light used in this methodology makes it difficult to plan and monitor light dose, or compare to lamp light doses.

METHODS

The present study expands on the commonly used PpIX-weighted effective surface irradiance metric by introducing a Monte Carlo method for estimating effective fluence rates into depths of the skin. The fluence rates are compared between multiple broadband and narrowband sources that have been reported in previous studies, and an effective total fluence for various treatment times is reported. A dynamic estimate of PpIX concentration produced during pro-drug incubation and treatment is used with the fluence estimates to calculate a photodynamic dose.

RESULTS

Even when there is up to a 5x reduction between the effective surface irradiance of the broadband light sources, the effective fluence below 250 μm depth is predicted to be relatively equivalent. An effective threshold fluence value (0. 70Jeff/cm2) is introduced based on a meta-analysis of previously published ALA-PpIX induced cell death. This was combined with a threshold PpIX concentration (50 nM) to define a threshold photodynamic dose of 0.035 u M Jeff/cm2.

CONCLUSIONS

The threshold was used to generate lookup tables to prescribe minimal treatment times to achieve depth-dependent cytotoxic effect based on incubation times and irradiance values for each light source.

Quantifying Direct DNA Damage in the Basal Layer of Skin Exposed to UV Radiation from Sunbeds

Nonmelanoma and melanoma skin cancers are attributable to DNA damage caused by ultraviolet (UV) radiation exposure. One DNA photoproduct, the cyclobutane pyrimidine dimer (CPD), is believed to lead to DNA mutations caused by UV radiation. Using radiative transfer simulations, we compare the number of CPDs directly induced by UV irradiation from artificial and natural UV sources (a standard sunbed and the midday summer Mediterranean sun) for skin types I and II on the Fitzpatrick scale. We use Monte Carlo radiative transfer (MCRT) modeling to track the progression of UV photons through a multilayered three dimensional (3D) grid that simulates the upper layers of the skin. By recording the energy deposited in the DNA-containing cells of the basal layer, the number of CPDs formed can be quantified. The aim of this work was to compare the number of CPDs formed in the basal layer of the skin and by implication the risk of developing cancer, as a consequence of irradiation by artificial and natural sources. Our simulations show that the number of CPDs formed per second during sunbed irradiation is almost three times that formed during solar irradiation.

Comparison of Blue and White Lamp Light with Sunlight for Daylight-Mediated, 5-ALA Photodynamic Therapy, in vivo

Daylight-mediated photodynamic therapy (d-PDT) as a treatment for actinic keratosis (AK) is an increasingly common technique due to a significant reduction in pain, leading to better patient tolerability. While past studies have looked at different light sources and delivery methods, this study strives to provide equivalent PpIX-weighted light doses with the hypothesis that artificial light sources could be equally as effective as natural sunlight if their PpIX-weighted fluences were equalized. Normal mouse skin was used as the model to compare blue LED light, metal halide white light and natural sunlight, with minimal incubation time between topical ALA application and the onset of light delivery. A total PpIX-weighted fluence of 20 J_eff cm^-2 was delivered over 2 h, and the efficacy of response was quantified using three acute bioassays for PDT damage: PpIX photobleaching, Stat3 crosslinking and quantitative histopathology. These bioassays indicated blue light was slightly inferior to both sunlight and white light, but that the latter two were not significantly different. The results suggest that metal halide white light could be a reasonable alternative to daylight PDT, which should allow a more controlled treatment that is independent of weather and yet should have similar response rates with limited pain during treatment.

Quantifying the radiant exposure and effective dose in patients treated for actinic keratoses with topical photodynamic therapy using daylight and LED white light

Daylight photodynamic therapy (dl-PDT) is as effective as conventional PDT (c-PDT) for treating actinic keratoses but has the advantage of reducing patient discomfort significantly. Topical dl-PDT and white light-PDT (wl-PDT) differ from c-PDT by way of light sources and methodology. We measured the variables associated with light dose delivery to skin surface and influence of geometry using a radiometer, a spectral radiometer and an illuminance meter. The associated errors of the measurement methods were assessed. The spectral and spatial distribution of the radiant energy from the LED white light source was evaluated in order to define the maximum treatment area, setup and treatment protocol for wl-PDT. We compared the data with two red LED light sources we use for c-PDT. The calculated effective light dose is the product of the normalised absorption spectrum of the photosensitizer, protoporphyrin IX (PpIX), the irradiance spectrum and the treatment time. The effective light dose from daylight ranged from 3  ±  0.4 to 44  ±  6 J cm^-2due to varying weather conditions. The effective light dose for wl-PDT was reproducible for treatments but it varied across the treatment area between 4  ±  0.1 J cm^-2 at the edge and 9  ±  0.1 J cm^-2 centrally. The effective light dose for the red waveband (615-645 nm) was 0.42  ±  0.05 J cm^-2 on a clear day, 0.05  ±  0.01 J cm^-2 on an overcast day and 0.9  ±  0.01 J cm^-2 using the white light. This compares with 0.95  ±  0.01 and 0.84  ±  0.01 J cm^-2 for c-PDT devices. Estimated errors associated with indirect determination of daylight effective light dose were very significant, particularly for effective light doses less than 5 J cm^-2 (up to 83% for irradiance calculations). The primary source of error is in establishment of the relationship between irradiance or illuminance and effective dose. Use of the O'Mahoney model is recommended using a calibrated logging illuminance meter with the detector in the plane of the treatment area.

Monte Carlo modelling of photodynamic therapy treatments comparing clustered three dimensional tumour structures with homogeneous tissue structures

We explore the effects of three dimensional (3D) tumour structures on depth dependent fluence rates, photodynamic doses (PDD) and fluorescence images through Monte Carlo radiation transfer modelling of photodynamic therapy. The aim with this work was to compare the commonly used uniform tumour densities with non-uniform densities to determine the importance of including 3D models in theoretical investigations. It was found that fractal 3D models resulted in deeper penetration on average of therapeutic radiation and higher PDD. An increase in effective treatment depth of 1 mm was observed for one of the investigated fractal structures, when comparing to the equivalent smooth model. Wide field fluorescence images were simulated, revealing information about the relationship between tumour structure and the appearance of the fluorescence intensity. Our models indicate that the 3D tumour structure strongly affects the spatial distribution of therapeutic light, the PDD and the wide field appearance of surface fluorescence images.