Photosensitizer absorption coefficient modeling and necrosis prediction during photodynamic therapy

The development of accurate predictive models for Photodynamic Therapy (PDT) has emerged as a valuable tool to adjust the current therapy dosimetry to get an optimal treatment response, and definitely to establish new personal protocols. Several attempts have been made in this way, although the influence of the photosensitizer depletion on the optical parameters has not been taken into account so far. We present a first approach to predict the spatio-temporal variation of the photosensitizer absorption coefficient during PDT applied to dermatological diseases, taking into account the photobleaching of a topical photosensitizer. This permits us to obtain the photons density absorbed by the photosensitizer molecules as the treatment progresses and to determine necrosis maps to estimate the short term therapeutic effects in the target tissue. The model presented also takes into account an inhomogeneous initial photosensitizer distribution, light propagation in biological media and the evolution of the molecular concentrations of different components involved in the photochemical reactions. The obtained results allow to investigate how the photosensitizer depletion during the photochemical reactions affects light absorption by the photosensitizer molecules as the optical radiation propagates through the target tissue, and estimate the necrotic tumor area progression under different treatment conditions.