Rozhin Penjweini

Study of tissue oxygen supply rate in a macroscopic photodynamic therapy singlet oxygen model.

Abstract. An appropriate expression for the oxygen supply rate (Γs) is required for the macroscopic modeling of the complex mechanisms of photodynamic therapy (PDT). It is unrealistic to model the actual heterogeneous tumor microvascular networks coupled with the PDT processes because of the large computational requirement. In this study, a theoretical microscopic model based on uniformly distributed Krogh cylinders is used to calculate Γs=g (1−[O32]/[O32]0) that can replace the complex modeling of blood vasculature while maintaining a reasonable resemblance to reality; g is the maximum oxygen supply rate and [O32]/[O32]0 is the volume-average tissue oxygen concentration normalized to its value prior to PDT. The model incorporates kinetic equations of oxygen diffusion and convection within capillaries and oxygen saturation from oxyhemoglobin. Oxygen supply to the tissue is via diffusion from the uniformly distributed blood vessels. Oxygen can also diffuse along the radius and the longitudinal axis of the cylinder within tissue. The relations of Γs to [3O2]/[3O2]0 are examined for a biologically reasonable range of the physiological parameters for the microvasculature and several light fluence rates (ϕ). The results show a linear relationship between Γs and [3O2]/[3O2]0, independent of ϕ and photochemical parameters; the obtained g ranges from 0.4 to 1390  μM/s.

Explicit dosimetry for 2-(1-hexyloxyethyl)-2-devinyl pyropheophorbide-a-mediated photodynamic therapy: macroscopic singlet oxygen modeling

Abstract. Type II photodynamic therapy (PDT) is based on the photochemical reactions mediated through an interaction between a photosensitizer, ground-state oxygen ([O32]), and light excitation at an appropriate wavelength, which results in production of reactive singlet oxygen ([O12]rx). We use an empirical macroscopic model based on four photochemical parameters for the calculation of [O12]rx threshold concentration ([O12]rx,sh) causing tissue necrosis in tumors after PDT. For this reason, 2-(1-hexyloxyethyl)-2-devinyl pyropheophorbide-a (HPPH)-mediated PDT was performed interstitially on mice with radiation-induced fibrosarcoma (RIF) tumors. A linear light source at 665 nm with total energy released per unit length of 12 to 100  J/cm and source power per unit length (LS) of 12 to 150  mW/cm was used to induce different radii of necrosis. Then the amount of [O12]rx calculated by the macroscopic model incorporating explicit PDT dosimetry of light fluence distribution, tissue optical properties, and HPPH concentration was correlated to the necrotic radius to obtain the model parameters and [O12]rx,sh. We provide evidence that [O12]rx is a better dosimetric quantity for predicting the treatment outcome than PDT dose, which is proportional to the time integral of the products of the photosensitizer concentration and light fluence rate.

Erratum: Macroscopic singlet oxygen modeling for dosimetry of Photofrin-mediated photodynamic therapy: an in-vivo study

Abstract. Although photodynamic therapy (PDT) is an established modality for cancer treatment, current dosimetric quantities, such as light fluence and PDT dose, do not account for the differences in PDT oxygen consumption for different fluence rates (φ). A macroscopic model was adopted to evaluate using calculated reacted singlet oxygen concentration ([O21]rx) to predict Photofrin-PDT outcome in mice bearing radiation-induced fibrosarcoma tumors, as singlet oxygen is the primary cytotoxic species responsible for cell death in type II PDT. Using a combination of fluences (50, 135, 200, and 250  J/cm2) and φ (50, 75, and 150  mW/cm2), tumor regrowth rate, k, was determined for each condition. A tumor cure index, CI=1−k/kcontrol, was calculated based on the k between PDT-treated groups and that of the control, kcontrol. The measured Photofrin concentration and light dose for each mouse were used to calculate PDT dose and [O21]rx, while mean optical properties (μa=0.9  cm−1, μs′=8.4  cm−1) were used to calculate φ for all mice. CI was correlated to the fluence, PDT dose, and [O21]rx with R2=0.35, 0.79, and 0.93, respectively. These results suggest that [O21]rx serves as a better dosimetric quantity for predicting PDT outcome.

Intracellular dynamics and fate of polystyrene nanoparticles in A549 Lung epithelial cells monitored by image (cross-) correlation spectroscopy and single particle tracking.

Novel insights in nanoparticle (NP) uptake routes of cells, their intracellular trafficking and subcellular targeting can be obtained through the investigation of their temporal and spatial behavior. In this work, we present the application of image (cross-) correlation spectroscopy (IC(C)S) and single particle tracking (SPT) to monitor the intracellular dynamics of polystyrene (PS) NPs in the human lung carcinoma A549 cell line. The ensemble kinetic behavior of NPs inside the cell was characterized by temporal and spatiotemporal image correlation spectroscopy (TICS and STICS). Moreover, a more direct interpretation of the diffusion and flow detected in the NP motion was obtained by SPT by monitoring individual NPs. Both techniques demonstrate that the PS NP transport in A549 cells is mainly dependent on microtubule-assisted transport. By applying spatiotemporal image cross-correlation spectroscopy (STICCS), the correlated motions of NPs with the early endosomes, late endosomes and lysosomes are identified. PS NPs were equally distributed among the endolysosomal compartment during the time interval of the experiments. The cotransport of the NPs with the lysosomes is significantly larger compared to the other cell organelles. In the present study we show that the complementarity of ICS-based techniques and SPT enables a consistent elaborate model of the complex behavior of NPs inside biological systems.

Explicit macroscopic singlet oxygen modeling for benzoporphyrin derivative monoacid ring A (BPD)-mediated photodynamic therapy.

Photodynamic therapy (PDT) is an effective non-ionizing treatment modality that is currently being used for various malignant and non-malignant diseases. In type II PDT with photosensitizers such as benzoporphyrin monoacid ring A (BPD), cell death is based on the creation of singlet oxygen (1O2). With a previously proposed empirical five-parameter macroscopic model, the threshold dose of singlet oxygen ([1O2]rx,sh]) to cause tissue necrosis in tumors treated with PDT was determined along with a range of the magnitude of the relevant photochemical parameters: the photochemical oxygen consumption rate per light fluence rate and photosensitizer concentration (ξ), the probability ratio of 1O2 to react with ground state photosensitizer compared to a cellular target (σ), the ratio of the monomolecular decay rate of the triplet state photosensitizer (β), the low photosensitizer concentration correction factor (δ), and the macroscopic maximum oxygen supply rate (g). Mice bearing radiation-induced fibrosarcoma (RIF) tumors were treated interstitially with a linear light source at 690nm with total energy released per unit length of 22.5-135J/cm and source power per unit length of 12-150mW/cm to induce different radii of necrosis. A fitting algorithm was developed to determine the photochemical parameters by minimizing the error function involving the range between the calculated reacted singlet oxygen ([1O2]rx) at necrosis radius and the [1O2]rx,sh. [1O2]rx was calculated based on explicit dosimetry of the light fluence distribution, the tissue optical properties, and the BPD concentration. The initial ground state oxygen concentration ([3O2]0) was set to be 40μM in this study. The photochemical parameters were found to be ξ=(55±40)×10-3cm2mW-1s-1, σ=(1.8±3)×10-5μM-1, and g=1.7±0.7μMs-1. We have taken the literature values for δ=33μM, and β=11.9μM. [1O2]rx has shown promise to be a more effective dosimetry quantity for predicting necrosis than either light dose or PDT dose, where the latter is simplistically a temporal integral of the products of the photosensitizer concentration and light fluence rate.

In vivo outcome study of BPD-mediated PDT using a macroscopic singlet oxygen model

Macroscopic modeling of the apparent reacted singlet oxygen concentration ([1O2]rx) for use with photodynamic therapy (PDT) has been developed and studied for benzoporphryin derivative monoacid ring A (BPD), a common photosensitizer. The four photophysical parameters (ξ, σ, β, δ) and threshold singlet oxygen dose ([1O2]rx, sh) have been investigated and determined using the RIF model of murine fibrosarcomas and interstitial treatment delivery. These parameters are examined and verified further by monitoring tumor growth post-PDT. BPD was administered at 1 mg/kg, and mice were treated 3 hours later with fluence rates ranging between 75 – 150 mW/cm2 and total fluences of 100 – 350 J/cm2. Treatment was delivered superficially using a collimated beam. Changes in tumor volume were tracked following treatment. The tumor growth rate was fitted for each treatment condition group and compared using dose metrics including total light dose, PDT dose, and reacted singlet oxygen. Initial data showing the correlation between outcomes and various dose metrics indicate that reacted singlet oxygen serves as a good dosimetric quantity for predicting PDT outcome.

A Comparison of Singlet Oxygen Explicit Dosimetry (SOED) and Singlet Oxygen Luminescence Dosimetry (SOLD) for Photofrin-Mediated Photodynamic Therapy

Accurate photodynamic therapy (PDT) dosimetry is critical for the use of PDT in the treatment of malignant and nonmalignant localized diseases. A singlet oxygen explicit dosimetry (SOED) model has been developed for in vivo purposes. It involves the measurement of the key components in PDT—light fluence (rate), photosensitizer concentration, and ground-state oxygen concentration ([3O2])—to calculate the amount of reacted singlet oxygen ([1O2]rx), the main cytotoxic component in type II PDT. Experiments were performed in phantoms with the photosensitizer Photofrin and in solution using phosphorescence-based singlet oxygen luminescence dosimetry (SOLD) to validate the SOED model. Oxygen concentration and photosensitizer photobleaching versus time were measured during PDT, along with direct SOLD measurements of singlet oxygen and triplet state lifetime (τΔ and τt), for various photosensitizer concentrations to determine necessary photophysical parameters. SOLD-determined cumulative [1O2]rx was compared to SOED-calculated [1O2]rx for various photosensitizer concentrations to show a clear correlation between the two methods. This illustrates that explicit dosimetry can be used when phosphorescence-based dosimetry is not feasible. Using SOED modeling, we have also shown evidence that SOLD-measured [1O2]rx using a 523 nm pulsed laser can be used to correlate to singlet oxygen generated by a 630 nm laser during a clinical malignant pleural mesothelioma (MPM) PDT protocol by using a conversion formula.

In-vivo outcome study of HPPH mediated PDT using singlet oxygen explicit dosimetry (SOED)

Type II photodynamic therapy (PDT) is based on the use of photochemical reactions mediated through an interaction between a tumor-selective photosensitizer, photoexcitation with a specific wavelength of light, and production of reactive singlet oxygen. However, the medical application of this technique has been limited due to inaccurate PDT dosimetric methods. The goal of this study is to examine the relationship between outcome (in terms of tumor growth rate) and calculated reacted singlet oxygen concentration [1O2]rx after HPPH-mediated PDT to compare with other PDT dose metrics, such as PDT dose or total light fluence. Mice with radiation-induced fibrosarcoma (RIF) tumors were treated with different light fluence and fluence rate conditions. Explicit measurements of photosensitizer drug concentration and tissue optical properties via fluorescence and absorption measurement with a contact probe before and after PDT were taken to then quantify total light fluence, PDT dose, and [1O2]rx based on a macroscopic model of singlet oxygen. In addition, photobleaching of photosenitizer were measured during PDT as a second check of the model. Changes in tumor volume were tracked following treatment and compared to the three calculated dose metrics. The correlations between total light fluence, PDT dose, reacted [1O2]rx and tumor growth demonstrate that [1O2]rx serves as a better dosimetric quantity for predicting treatment outcome and a clinically relevant tumor growth endpoint.

Evaluation of singlet oxygen explicit dosimetry for predicting treatment outcomes of benzoporphyrin derivative monoacid ring A-mediated photodynamic therapy

Abstract. Existing dosimetric quantities do not fully account for the dynamic interactions between the key components of photodynamic therapy (PDT) or the varying PDT oxygen consumption rates for different fluence rates. Using a macroscopic model, reacted singlet oxygen ([O21]rx) was calculated and evaluated for its effectiveness as a dosimetric metric for PDT outcome. Mice bearing radiation-induced fibrosarcoma tumors were treated with benzoporphyrin derivative monoacid ring A (BPD) at a drug-light interval of 3 h with various in-air fluences (30 to 350  J/cm2) and in-air fluence rates (50 to 150  mW/cm2). Explicit measurements of BPD concentration and tissue optical properties were performed and used to calculate [O21]rx, photobleaching ratio, and PDT dose. For four mice, in situ measurements of O23 and BPD concentration were monitored in real time and used to validate the in-vivo photochemical parameters. Changes in tumor volume following treatment were used to determine the cure index, CI=1−k/kctr, where k and kctr are the tumor regrowth rates with PDT and without PDT, respectively. The correlation between CI and the dose metrics showed that the calculated [O21]rx at 3 mm is an effective dosimetric quantity for predicting treatment outcome and a clinically relevant tumor regrowth endpoint.

Investigating the impact of oxygen concentration and blood flow variation on photodynamic therapy

Type II photodynamic therapy (PDT) is used for cancer treatment based on the combined action of a photosensitizer, a special wavelength of light, oxygen (3O2) and generation of singlet oxygen (1O2). Intra-patient and inter-patient variability of oxygen concentration ([3O2]) before and after the treatment as well as photosensitizer concentration and hemodynamic parameters such as blood flow during PDT has been reported. Simulation of these variations is valuable, as it would be a means for the rapid assessment of treatment effect. A mathematical model has been previously developed to incorporate the diffusion equation for light transport in tissue and the macroscopic kinetic equations for simulation of [3O2], photosensitizers in ground and triplet states and concentration of the reacted singlet oxygen ([1O₂]rx) during PDT. In this study, the finite-element based calculation of the macroscopic kinetic equations is done for 2-(1- Hexyloxyethyl)-2-devinyl pyropheophorbide (HPPH)-mediated PDT by incorporating the information of the photosensitizer photochemical parameters as well as the tissue optical properties, photosensitizer concentration, initial oxygen concentration ([3O2]0), blood flow changes and Φ that have been measured in mice bearing radiation-induced fibrosarcoma (RIF) tumors. Then, [1O2]rx calculated by using the measured [3O2] during the PDT is compared with [1O2]rx calculated based on the simulated [3O₂]; both calculations showed a reasonably good agreement. Moreover, the impacts of the blood flow changes and [3O2]0 on [1O2]rx have been investigated, which showed no pronounced effect of the blood flow changes on the long-term 1O2 generation. When [3O2]0 becomes limiting, small changes in [3O₂] have large effects on [1O2]rx.