Search

We have developed a comprehensive theoretical model for rigorously describing the spatial and temporal dynamics of oxygen (3O2) consumption and transport and microscopic photodynamic dose deposition during photodynamic therapy (PDT) in vivo. Previously published models have been improved by considering perfused vessels as a time-dependent 3O2 source and linking the 3O2 concentration in the vessel to that within the tissue through the Hill equation. The time-dependent photochemical 3O2 consumption rate incorporates sensitizer photobleaching effects and an experimentally determined initially nonuniform photosensitizer distribution. The axial transport of 3O2 is provided for in the capillaries and in the surrounding tissue. A self-sensitized singlet oxygen (1O2)-mediated bleaching mechanism and the measured, initially nonuniform distribution of mesotetrahydroxyphenyl chlorin at 3 h after intravascular administration were used to demonstrate the capabilities of the model. Time-evolved distributions of 3O2 concentration were obtained by numerically solving two-dimensional diffusion-with-reaction equations both in the capillary and the adjacent tissue. Using experimentally established physiological and photophysical parameters, the mathematical model allows computation of the dynamic variation of hemoglobin-3O2 saturation (SO2) within the vessels, irreversible sensitizer degradation due to photobleaching, and the microscopic distributions of 3O2, sensitizer concentration, and 1O2 dose deposition under various irradiation conditions. The simulations reveal severe axial gradients in 3O2 and in photodynamic dose deposition in response to a wide range of clinically relevant treatment parameters. Thus, unlike former Krogh cylinder-based models, which assume a constant 3O2 concentration at the vessel, this new model identifies conditions in which 3O2 depletion and minimal deposition of reacting 1O2 exist near the end of axial segments of vessels and shows that treatment-limiting 3O2 depletion is induced at fluence rates as low as 10 mW cm(-2). These calculations also demonstrate that intercapillary heterogeneity of photosensitizer contributes significantly to the distribution of photodynamic dose. This more rigorous mathematical model enables comparison with experimentally observable, volume-averaged quantities such as SO2 and the loss of sensitizer fluorescence through bleaching that have not been included in previous analyses. Further, it establishes some of the intrinsic limitations of such measurements. Specifically, our simulations demonstrate that tissue measurements of SO2 and of photobleaching are necessarily insensitive to microscopic heterogeneity of photodynamic dose deposition and are sensitive to intercapillary spacing. Because prior knowledge of intercapillary distances in tumors is generally unavailable, these measurements must be interpreted with caution. We anticipate that this model will make useful dosimetry predictions that should inform optimal treatment conditions and improve current clinical protocols.

We present two forward-adjoint models for recovering intrinsic fluorescence spectra and hemoglobin oxygen saturation of turbid samples. The first fits measured diffuse reflectance spectra to obtain the absorption and scattering spectra of the medium, and these are then used to correct distortions imposed on the fluorescence spectrum by absorption and scattering. The second fits only the measured fluorescence spectrum to determine simultaneously the amplitudes of absorption and fluorescence basis spectra and scattering parameters. Both methods are validated with Monte Carlo simulations and experimentally in scattering phantoms containing nicotinamide adenine dinucleotide and human erythrocytes. Preliminary measurements from murine tumors in vivo are presented.

We present a method for the analysis of steady state diffuse reflectance spectra obtained from vascularized tissue or from tissue simulating phantoms at a single, short source-detector separation. This method uses reasonable assumptions about the structure of the reduced scattering spectrum and basis absorption spectra for oxy- and deoxyhemoglobin, which dominate tissue absorption in the visible region of the spectrum. Using a hybrid P3-diffusion description of light propagation, described originally by Hull and Foster [J. Opt. Soc. Am. A 18, 584-599 (2001)] and suitable for short (approximately 1 mm) source-detector separations and optical properties of tissue at visible wavelengths, we create a forward model of the diffuse reflectance with four free parameters. We demonstrate that this model is able to recover accurately the hemoglobin concentrations and scattering properties from synthetic data generated by Monte Carlo simulation and from reflectance spectra acquired from tissue-simulating phantoms containing intact human erythrocytes. We show also that the method is capable of monitoring carbogen-induced changes in murine tumor oxygenation in vivo. The successful implementation of single, short detector separations enables the measurement of intratumor heterogeneities in hemoglobin oxygen saturation and responses to carbogen using a simple fiber-based probe design.

We report results of experiments that evaluated the influence of oxygenation on the penetration of red light in tissue, with particular emphasis on 630 and 650 nm laser wavelengths commonly used in photodynamic therapy (PDT) of solid tumours. Direct measurements in tissue-simulating phantoms comprised of intact human erythrocytes suspended in a scattering emulsion demonstrated significant enhancements in fluence rate at depths of 0.5-2 cm from the irradiated surface when the cells were fully oxygenated versus fully deoxygenated. The 630 and 650 nm fluence rates at depth in the homogeneous phantoms continued to increase when examined over a continuous range of oxygen partial pressures from 0 to 160 Torr. When considered as a function of haemoglobin oxygen saturation, the largest increases in fluence rate were observed as the saturation increased beyond 70%. Dramatic increases in optical fluence rate were measured at the base of 1-cm-thick subcutaneous EMT6 mammary carcinomas in vivo when the tumour-bearing mouse was subjected to carbogen through a nose cone. These results indicate that improved tumour oxygenation is important in PDT not only for the maintenance of the oxygen-dependent photochemistry but, through the effects reported here, may also enable more efficient treatment of thicker lesions.

We present direct experimental evidence of the fluence-rate-dependent, radiation-induced variations in intratumor oxygen partial pressure (pO(2)) in HT29 human colon adenocarcinoma xenografts subjected to meta-tetra(hydroxyphenyl)chlorin (mTHPC)-based photodynamic therapy (PDT). The data establish a correlation between tumor oxygenation and treatment outcome. Tumor-bearing mice were injected with 0.3 mg/kg photosensitizer and subjected 72 h later to a 12 J/cm(2) red light dose administered at fluence rates of 5, 30, 90 and 160 mW/cm(2). A significant decrease in mean and median pO(2) was registered at approximately half of the total radiation fluence was delivered in tumors treated at rates of 160 and 90 mW/cm(2). Conversely, with the two lower fluence rates, intratumor pO(2) was maintained at levels comparable to those measured before illumination. Tumor oxygenation values registered shortly after every treatment protocol were at least equal to baseline levels, thus excluding the possibility of significant acute vessel damage during illumination. The tumor regrowth profile correlated with the pO(2) values monitored during irradiation. Tumors treated with fluence rates of 5 and 30 mW/cm(2) exhibited significantly longer tumor quadrupling times than those treated at 160 and 90 mW/cm(2). Improved tumor destruction could be expected by reducing the rate and the extent of oxygen depletion during meta-tetra(hydroxyphenyl)chlorin photodynamic therapy using low fluence rates.