PDT Study Using a Model Incorporating Initial Oxygen Concentration and Blood Flow Increase

R. Penjweini[1], T. C. Zhu[1],
[1] Department of Radiation Oncology, University of Pennsylvania, Philadelphia, PA, USA
Published in 2015

Introduction: Type II photodynamic therapy (PDT) is an experimental modality for cancer treatment based on the combined action of a photosensitizing drug (photosensitizer), a special wavelength of light and singlet oxygen (1O2) generation; the cell killing is caused by the reaction of cellular acceptors with 1O2. A mathematical model has been previously developed to incorporate the macroscopic kinetic equations for 1O2 generation, photosensitizers in ground and triplet states, oxygen, and tissue acceptors along with the diffusion equation for the light transport in tissue. Initial oxygen concentration ([3O2]0) and the changes of the blood flow during PDT might affect the treatment efficacy. In this study, the effects of [3O2]0, as well as the blood flow increase during PDT, on the magnitude of 1O2 generation is studied.

Methods: By simplifying and combining the energy transfer processes in PDT, a set of equations are produced, which describes the creation of 1O2. These equations are dependent on various parameters such as the light source (LS), optical absorption and scattering coefficients (µa and µ's) and photochemical parameters of the photosensitizer. In this model, the spatial distribution of light fluence (ϕ) in the tumor is calculated via Eq. (1), based on the diffusion approximation. Spatial and temporal photosensitizer distribution (S0), 3O2 and 1O2 concentrations are obtained by solving a set of coupled time-dependent differential equations [1].

The symbol Γs = g (1-[3O2]/[3O2]0), where g is the maximum oxygen supply rate and δ is the correction parameter for low photosensitizer concentration; [1O2]rx is defined as the 1O2 effectively leading cell death. The initial conditions are: ϕ= 0, [1O2]rx = 0, [3O2] = [3O2]0, and [S0] = [S0]0 = 5 M.

Application of COMSOL Multiphysics® software: The forward calculation of the macroscopic kinetic equations was done in COMSOL Multiphysics® software for the modeling of ³O₂ and 1O2 generation. The finite-element calculation was implemented within COMSOL software by varying the input parameters, such as [³O₂]₀ and ϕ. Based on our mice study, [³O₂]₀ was varying between 5–60 μM. The blood flow was considered by making g to be time dependent based on the published results [2-3].

Results: Using the obtained results, one can correlate the tissue oxygenation with the macroscopic quantity measured in the blood vessel. Moreover, the effects of the blood flow changes during PDT on 1O2 generation, [³O₂] as well as photosensitizer photobleaching is obtained. This model also shows that [³O₂]₀ has no effect on the long term 1O2 generation.

Conclusion: Several improvements of the model formulations have been made to consider [³O₂]₀ for both tumor and normal tissues during PDT. In the well-oxygenated tissue, the exact [3O₂] will have little influence on treatment efficacy. When [3O₂] becomes limiting, small changes in ϕ or [3O₂]0 have large effects.