A dynamic model for ALA-PDT of skin: simulation of temporal and spatial distributions of ground-state oxygen, photosensitizer and singlet oxygen

Singlet oxygen (¹O₂) direct dosimetry and photosensitizer fluorescence photobleaching are being investigated and applied as dosimetric tools during 5-aminolevulinic acid (ALA)-induced protophorphyrin IX (PpIX) photodynamic therapy (PDT) of normal skin and skin cancers. The correlations of photosensitizer fluorescence and singlet oxygen luminescence (SOL) emission signals to ¹O2 distribution and cumulative ¹O₂dose are difficult to interpret because of the temporal and spatial variations of three essential components (light fluence rate, photosensitizer concentration and oxygen concentration) in PDT. A one-dimensional model is proposed in this paper to simulate the dynamic process of ALA-PDT of normal human skin in order to investigate the time-resolved evolution of PpIX, ground-state oxygen (³O₂and ¹O₂ distributions. The model incorporates a simplified three-layer semi-infinite skin tissue, Monte Carlo simulations of excitation light fluence and both PpIX fluorescence and SOL emission signals reaching the skin surface, ¹O₂-mediated photobleaching mechanism for updating PpIX, ³O₂ and ¹O₂ distributions after the delivery of each light dose increment, ground-state oxygen supply by diffusion from the atmosphere and perfusion from blood vessels, a cumulative ¹O₂-dependent threshold vascular response, and the initial non-uniform distribution of PpIX. The PpIX fluorescence simulated using this model is compared with clinical data reported by Cottrell et al (2008 Clin. Cancer Res. 14 4475-83) for a range of irradiances (10-150 mW cm⁻²). Except for the vascular response, one set of parameters is used to fit data at all irradiances. The time-resolved depth-dependent distributions of PpIX, ³O₂ and ¹O₂ at representative irradiances are presented and discussed in this paper, as well as the PDT-induced vascular response at different depths. Tissue hypoxia and shutdown of oxygen supply occur in the upper dermis, where PpIX is also preserved at the end of treatment.

[1]  D. Kessel,et al.  DIRECT DETECTION OF SINGLET OXYGEN SENSITIZED BY HAEMATOPORPHYRIN AND RELATED COMPOUNDS , 1986, Photochemistry and photobiology.

[2]  Elena Salomatina,et al.  Optical properties of normal and cancerous human skin in the visible and near-infrared spectral range. , 2006, Journal of biomedical optics.

[3]  Ross W. Boyle,et al.  Photodynamic Therapy and the Development of Metal-Based Photosensitisers , 2008, Metal-based drugs.

[4]  Arjen Amelink,et al.  Monitoring ALA‐induced PpIX Photodynamic Therapy in the Rat Esophagus Using Fluorescence and Reflectance Spectroscopy , 2008, Photochemistry and photobiology.

[5]  Thomas H Foster,et al.  Irradiance-Dependent Photobleaching and Pain in δ-Aminolevulinic Acid-Photodynamic Therapy of Superficial Basal Cell Carcinomas , 2008, Clinical Cancer Research.

[6]  Michael S Patterson,et al.  Photobleaching kinetics, photoproduct formation, and dose estimation during ALA induced PpIX PDT of MLL cells under well oxygenated and hypoxic conditions , 2006, Photochemical & photobiological sciences : Official journal of the European Photochemistry Association and the European Society for Photobiology.

[7]  B. Wilson,et al.  The physics, biophysics and technology of photodynamic therapy , 2008, Physics in medicine and biology.

[8]  Henricus J C M Sterenborg,et al.  Protoporphyrin IX Fluorescence Photobleaching and the Response of Rat Barrett's Esophagus Following 5-aminolevulinic Acid Photodynamic Therapy , 2006, Photochemistry and photobiology.

[9]  Henricus J C M Sterenborg,et al.  Quantitative Model Calculation of the Time-dependent Protoporphyrin IX Concentration in Normal Human Epidermis After Delivery of ALA by Passive Topical Application or Iontophoresis¶ , 2002, Photochemistry and photobiology.

[10]  K. Stamnes,et al.  Choice of Optimal Wavelength for PDT: The Significance of Oxygen Depletion , 2005, Photochemistry and photobiology.

[11]  Ken Kang-Hsin Wang,et al.  Simulations of measured photobleaching kinetics in human basal cell carcinomas suggest blood flow reductions during ALA‐PDT , 2009, Lasers in surgery and medicine.

[12]  Brian C. Wilson,et al.  A Monte Carlo model of detected singlet oxygen luminescence and photosensitizer fluorescence during ALA-PDT of skin , 2009, World Congress of the International Photodynamic Association.

[13]  B. Wilson,et al.  Comparison of the In Vivo Photodynamic Threshold Dose for Photofrin, Mono‐ and Tetrasulfonated Aluminum Phthalocyanine Using a Rat Liver Model , 1998, Photochemistry and photobiology.

[14]  Ken Kang-Hsin Wang,et al.  A comprehensive mathematical model of microscopic dose deposition in photodynamic therapy. , 2006, Medical physics.

[15]  B. Wilson,et al.  Forward-adjoint fluorescence model: Monte Carlo integration and experimental validation. , 1997, Applied optics.

[16]  T. Foster,et al.  Singlet Oxygen‐Versus Nonsinglet Oxygen‐Mediated Mechanisms of Sensitizer Photobleaching and Their Effects on Photodynamic Dosimetry , 1998, Photochemistry and photobiology.

[17]  P. R. Bevington,et al.  Data Reduction and Error Analysis for the Physical Sciences , 1969 .

[18]  P. Vaupel.  Effect of percentual water content in tissues and liquids on the diffusion coefficients of O2, CO2, N2, and H2 , 1976, Pflügers Archiv.

[19]  L Wang,et al.  MCML--Monte Carlo modeling of light transport in multi-layered tissues. , 1995, Computer methods and programs in biomedicine.

[20]  L. Wang,et al.  Optimal beam size for light delivery to absorption-enhanced tumors buried in biological tissues and effect of multiple-beam delivery: a Monte Carlo study. , 1997, Applied optics.

[21]  J. Hudson,et al.  An analysis of factors affecting tissue oxygen tension , 1964, Proceedings of the Royal Society of London. Series B. Biological Sciences.

[22]  Jarod C Finlay,et al.  Photobleaching kinetics of Photofrin in vivo and in multicell tumour spheroids indicate two simultaneous bleaching mechanisms. , 2004, Physics in medicine and biology.

[23]  D. Kessel,et al.  Tumor blood-flow changes following protoporphyrin IX-based photodynamic therapy in mice and humans. , 1999, Journal of photochemistry and photobiology. B, Biology.

[24]  D. Lübbers,et al.  Initial Results of Intracutaneous Measurements of PO 2 Profiles , 1987 .

[25]  Thomas H Foster,et al.  Photophysical Parameters, Photosensitizer Retention and Tissue Optical Properties Completely Account for the Higher Photodynamic Efficacy of meso-Tetra-Hydroxyphenyl-Chlorin vs Photofrin¶ , 2005, Photochemistry and photobiology.

[26]  A. Hill,et al.  The possible effects of the aggregation of the molecules of haemoglobin on its dissociation curves , 1910 .

[27]  Wen Wang.  Oxygen Partial Pressure in Outer Layers of Skin: Simulation Using Three‐Dimensional Multilayered Models , 2005, Microcirculation.

[28]  Michael S. Patterson,et al.  Calculation of Singlet Oxygen Dose from Photosensitizer Fluorescence and Photobleaching During mTHPC Photodynamic Therapy of MLL Cells¶ , 2005 .

[29]  D. Lübbers,et al.  The cutaneous uptake of atmospheric oxygen contributes significantly to the oxygen supply of human dermis and epidermis , 2002, The Journal of physiology.

[30]  Tayyaba Hasan,et al.  In-vivo singlet oxygen dosimetry of clinical 5-aminolevulinic acid photodynamic therapy. , 2008, Journal of biomedical optics.

[31]  Jarod C Finlay,et al.  Explicit dosimetry for photodynamic therapy: macroscopic singlet oxygen modeling , 2010, Journal of biophotonics.