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

Haixia Qiu; Michele M. Kim; Rozhin Penjweini; Timothy C. Zhu

doi:10.1117/1.JBO.21.8.088002

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

Haixia Qiu, Michele M. Kim, Rozhin Penjweini, Timothy C. Zhu

Author Affiliations +

Journal of Biomedical Optics, Vol. 21, Issue 8, 088002 (August 2016). https://doi.org/10.1117/1.JBO.21.8.088002

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 ( [ O 2 1 ] 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 / cm 2 ) and ϕ (50, 75, and 150 mW / cm 2 ), tumor regrowth rate, k , was determined for each condition. A tumor cure index, CI = 1 − k / k control , was calculated based on the k between PDT-treated groups and that of the control,

1. Introduction

As an established methodology for cancer treatment, photodynamic therapy (PDT) has been approved by the US Food and Drug Administration for the treatment of some cancers and precancers, such as Photofrin-PDT for esophageal cancer, nonsmall cell lung cancer, Barrett’s esophagus,¹^–³ and 5-aminolevulinic acid (ALA)-PDT for actinic keratosis.⁴ Furthermore, off-label PDT applications as well as clinical trials have been performed.⁵^–⁷ PDT is a photochemical process combining the interactions among photosensitizer, light, and oxygen. During this process, the ground-state photosensitizer ( $S_{0}$ ) is excited by light and converted to an excited state ( $S_{1}$ ). The $S_{1}$ state is a short-lived state and the photosensitizer will lose its energy and return to the ground state by either emitting a fluorescent photon or decaying to the ground state nonradiatively, or it will undergo a process known as intersystem crossing and converted to a triplet state ( $T_{1}$ ). While in the $T_{1}$ state, the photosensitizer can be involved in two types of processes called type I and type II. In type I photodynamic interactions, the excited triplet can transfer an electron to ground-state molecular oxygen ( ${}^{3}{O_{2}}$ ) to form a charged superoxide anion, which in turn, form reactive oxygen species (ROS) via secondary chemical reactions. ROS interact with substrates to cause biological damage.⁸^–¹⁰ The photosensitizer triplet may react directly with a biological substrate under hypoxic conditions. In type II interactions, the excited triplet will transfer energy to ${}^{3}{O_{2}}$ to produce excited-state singlet oxygen ( ${}^{1}{O_{2}}$ ), the major cytotoxic species causing biological damage.⁸ Although the precise mechanisms of PDT are not yet fully understood, it is generally accepted that ${}^{1}{O_{2}}$ is primarily responsible for cell death and PDT outcome for type II photosensitizers.¹¹ Although the clinical use and indications of PDT are increasing slowly, an ideal dosimetric predictor for PDT efficacy is still elusive, which to some extent, has hindered the clinical application of PDT. In most clinical protocols, dosimetry is still based on the administered photosensitizer dosage and the delivered light fluence under a range of fluence rates ( $φ$ ). This method of clinical dosimetry does not take into account the absorption and scattering of light or the uptake of the photosensitizer in the target tissue.¹² This might be responsible for the unpredictable response and variability between clinical dosimetry and PDT outcome.¹³^,¹⁴

To find a better solution for PDT dosimetry, explicit PDT dosimetry has been introduced. Unlike using the given photosensitizer dosage and light fluence as PDT metrics, in explicit PDT dosimetry, the light dose absorbed by the photosensitizer in tissue is measured and defined as PDT dose, which is determined by the product of photosensitizer concentration and light fluence. PDT dose may serve as a good predictor of outcome in ${}^{3}{O_{2}}$ enriched conditions.¹² However, in hypoxic conditions, it becomes less accurate for predicting the PDT efficacy since it does not consider the rate of PDT consumption of ${}^{3}{O_{2}}$ for different $φ$ .¹⁵^,¹⁶ With integration of the variation of ${}^{3}{O_{2}}$ concentration ( $[{}^{3}{O_{2}}]$ ) into explicit dosimetry, a so-called singlet oxygen explicit dosimetry method based on mathematical modeling to calculate the ${}^{1}{O_{2}}$ has made it possible to predict the PDT outcome.¹⁷ For this purpose, an empirical four-parameter macroscopic ${}^{1}{O_{2}}$ model was developed, and the reacted singlet oxygen ( ${[{}^{1}{O_{2}}]}_{rx}$ ) was proposed as a dosimetry quantity for the accumulated concentration of ${}^{1}{O_{2}}$ that is responsible for cellular killing and clinical outcomes.¹⁸ In this study, this model was adopted to evaluate the potential of using calculated ${[{}^{1}{O_{2}}]}_{rx}$ as a dosimetric predictor for Photofrin-mediated PDT on tumor growth in a mouse tumor model.

2. Materials and Methods

2.1.

Mouse Tumor Model

Radiation-induced fibrosarcoma (RIF) cells of logarithmic growth phase ( $30 μ l$ of $1 \times 10^{7} cells / ml$ ) were injected subcutaneously in the right shoulders of 6- to 8-week-old female C3H mice (NCI-Frederick, Frederick, Maryland). The fur of the tumor region was clipped prior to cell inoculation. PDT was carried out around 5 to 10 days after inoculation when tumors reached 3 to 4 mm in length and in height. A few days before PDT, the tumor and the surrounding area was depilated with Nair (Church & Dwight Co., Inc., Ewing, New Jersey). All procedures were approved by the University of Pennsylvania Institutional Animal Care and Use Committee. Animal husbandry was provided by the University of Pennsylvania Laboratory Animal Resources in association with Assessment and Accreditation of Laboratory Animal Care.

2.2.

Photodynamic Therapy Protocol

Photofrin (Pinnacle Biologics, Chicago, Illinois) at a dosage of $5 mg / kg$ was injected through the mouse tail vein as described previously.¹⁹^,²⁰ At a 24-h drug–light interval, superficial irradiation of the tumor was performed with a 630-nm laser (Biolitec AG., A-1030, Vienna). A microlens fiber was coupled to the laser to irradiate the tumor uniformly. Animals were assigned to four light dose groups, and each group was comprised of 2 to 3 subgroups with different $φ$ . There were a total of 11 treatment groups: $50 J / {cm}^{2}$ at 50, 75, and $150 mW / {cm}^{2}$ , $135 J / {cm}^{2}$ at 50, 75, and $150 mW / {cm}^{2}$ , $200 J / {cm}^{2}$ at 50 and $150 mW / {cm}^{2}$ , and $250 J / {cm}^{2}$ at 50, 75, and $150 mW / {cm}^{2}$ . Tumor-bearing mice that received neither light irradiation nor Photofrin were used as controls.

2.3.

Photodynamic Therapy Efficiency Assessment

Mouse tumors were measured daily using a sliding caliper for up to 14 days post-PDT. The tumor volume ( $V$ ) was calculated by the formula $V = π / 6 \times a^{2} \times b$ ; where $a$ and $b$ refer to the width and length of the tumor, respectively.²¹ Since initial tumor volumes at the time of treatment were not identical among mice, daily tracked tumor volumes were scaled according to the initial volume; the treatment day volume is normalized to be $\sim 12 {mm}^{3}$ for all mice. Tumors were selected for various treatment groups randomly. By fitting the tumor volumes to an exponential growth equation, $V = A e^{k d}$ the tumor regrowth rate ( $k$ ) was obtained, where $A$ is the amplitude and $d$ represents the number of the days after PDT. $k$ was not allowed to be negative in the fit. The tumor cure index (CI) was calculated by the following expression:²²^–²⁴

Eq. (1)

CI = 1 - \frac{k}{k_{control}},

where

k_{control}

is the tumor regrowth rate for the control mice without PDT.

2.4.

Determination of the Light Distribution in Tumors with Different Optical Properties

In-vivo measurements of the tissue optical properties were performed in another set of mice administered with $5 mg / kg$ Photofrin with a two-catheter method described in detail elsewhere.¹⁸ Two parallel catheters were inserted into the tumor with one centrally positioned and the other a fixed distance away. A 2-mm point source fiber was coupled to the 630-nm diode laser, and light fluence profiles were obtained along the length of the second catheter using an isotropic detector. Profiles were then fit to a diffusion equation to extract the absorption coefficient ( $μ_{a}$ ) and reduced scattering coefficient ( $μ_{s}^{'}$ ) of the tumor tissue for the mice.²⁵ An analytical functional fitting to Monte Carlo (MC) simulation was used to calculate the spatial distribution of $φ$ in the tumors with different $μ_{a}$ and $μ_{s}^{'}$ using a method that has been published previously.²⁶ Figure 1(a) shows the spatial distribution of the ratio of fluence rate and in-air fluence rate ( $φ / φ_{air}$ ). The deviation of $φ$ due to the effect of various $μ_{a}$ and $μ_{s}^{'}$ has been shown in Fig. 1(b). In-air fluence rate was measured on the surface of the tumors for each individual mouse that was used for this study. This independent verification indicated that calculation using the mean optical properties is accurate to within 10%.

Fig. 1

MC simulation of the spatial distribution of light fluence rate ( $φ$ ) in RIF tumors with different optical properties ( $μ_{a}$ and $μ_{s}^{'}$ ). (a) The ratio of light fluence rate and in-air fluence rate ( $φ / φ_{air}$ ) versus tumor depth and (b) the ratio of light fluence rate for each condition and the mean fluence rate ( $φ / φ_{Mean}$ ) versus tumor depth. $φ_{Mean}$ was calculated from the mean optical properties of ${\bar{μ}}_{a} = 0.9 {cm}^{- 1}$ and ${\bar{μ}}_{s}^{'} = 8.4 {cm}^{- 1}$ .

2.5.

Measurement of the Photofrin Concentration in Tumors

A custom-designed multifiber probe connected to a CCD camera (InSpectrum 150, Roper Scientific, Princeton, New Jersey) was used to collect surface fluorescence of Photofrin at the tumor site. A 405-nm laser was used as an excitation light to obtain the fluorescence signals. The fluorescence measurement was conducted immediately before and after PDT. The raw fluorescence spectrum was fit to the basis spectrum of Photofrin and autofluorescence in the absence of the photosensitizer using a singular value decomposition (SVD) method described by Finlay et al.²⁷^,²⁸ The attenuation of Photofrin fluorescence signal due to the light absorption and scattering by tissue was corrected by applying an empirical correction factor (CF)

Eq. (2)

CF = \frac{C_{01} + C_{02} μ_{s}^{'}}{μ_{s}^{'}} \cdot \exp [(b_{1} + b_{2} μ_{s}^{'}) . \sqrt{3 μ_{a} μ_{s}^{'}}],

where constants

C_{01} = 3.1

,

C_{02} = 0.26

,

b_{1} = 0.13

, and

b_{2} = - 0.00012

were determined from fitting the fluorescence SVD for phantoms with different

μ_{a}

and

μ_{s}^{'}

. To obtain the values of the parameters in the expression for CF, a set of liquid phantoms were prepared with various optical properties (

μ_{a} = 0

to

0.9 {cm}^{- 1}

,

μ_{s}^{'} = 5

to

15 {cm}^{- 1}

) and a constant Photofrin concentration of

2 mg / kg

(

3.3 μ M

) using Intralipid (Fresenius Kabi, Uppsala, Sweden) as the scatterer and ink (Parker^® Quink^®) as the absorber. Optical properties of each phantom were measured and fluorescence spectra were obtained. The corrected SVD (

{SVD}_{cor} = SVD \times CF

) with the appropriate CF should have the same magnitude for phantoms with the same Photofrin concentration. Fitting the raw SVD to this enabled for fitting and determining the parameters necessary in CF. Furthermore, a calibration curve was determined between

{SVD}_{cor}

and various Photofrin concentrations ranging from 0 to

7 mg / kg

(0 to

11.7 μ M

). Raw fluorescence SVD measurements performed in vivo could then be corrected using Eq. (6), and absolute Photofrin concentration could be determined using the calibration curve from

{SVD}_{cor}

values in vivo [see Figs. 2(a) and 2(b)].

Fig. 2

(a) Fluorescence singular decomposition (SVD) amplitude for phantom experiments with different optical properties but the same Photofrin concentration. (b) Photofrin concentration (in $mg / kg$ ) versus corrected fluorescence SVD ( ${SVD}_{cor}$ ). (c) The in-vivo measured photosensitizer concentration using the contact probe (based on fluorescence spectroscopy) versus ex-vivo measured Photofrin concentration.

Using a method described by Penjweini et al.,¹⁶ ex-vivo measurements of the Photofrin concentration were performed in another set of mice administered with the photosensitizer at the same concentration as the PDT-treated mice. After the 24-h drug–light interval, interstitial measurements were performed using the same fluorescence method as that of the PDT-treated mice. The tumors were then immediately excised and frozen. Homogenized solutions of the tumors were prepared using Soluble (PerkinElmer, Waltham, Massachusetts) and fluorescence of the homogenized sample was measured by a spectrofluorometer (FluoroMax-3; Jobin Yvon, Inc.) with an excitation at 405 nm and an emission range from 590 to 740 nm with an emission maximum at 630 nm. The photosensitizer concentration in the tissue was calculated based on the change in fluorescence peak magnitude resulting from the addition of a known amount of Photofrin to each sample after its initial reading. As shown in Fig. 2(c), the ex-vivo measurements of Photofrin concentration were compared to those obtained in vivo using the interstitial method to evaluate the in-vivo acquired Photofrin concentrations; the linear fit, $y = 0.98 x$ , to the data (shown as a solid line) with the fitting goodness of $R^{2} = 0.99$ shows their reasonable agreement, validating the in-vivo measurements. The dashed line represents the line for $y = x$ , if the two measurements were completely in agreement.

For actual application of CF [Eq. (2)], the mean tissue optical properties ( $μ_{a} = 0.9 {cm}^{- 1}$ , $μ_{s}^{'} = 8.4 {cm}^{- 1}$ ) were used for all mice. This may introduce an additional maximum CF variation of $\pm 4 %$ due to the variation of tissue optical properties as shown in Fig. 1.

2.6.

Macroscopic Singlet Oxygen Model for Photodynamic Therapy

An empirical macroscopic ${}^{1}{O_{2}}$ model derived from reaction rate equations for a type II PDT mechanism was adopted in this study.¹⁸ The complete set of the reaction rate equations and their derivations have been described in detail elsewhere.¹⁸^,²⁹^–³¹ In this model, the reaction rate equations are simplified as the following, which incorporates a set of PDT kinetic equations:

Eq. (3)

\frac{d [S_{0}]}{d t} + (ξ σ \frac{φ ([S_{0}] + δ) [{}^{3}{O_{2}}]}{[{}^{3}{O_{2}}] + β}) [S_{0}] = 0,

Eq. (4)

\frac{d [{}^{3}{O_{2}}]}{d t} + (ξ \frac{φ [S_{0}]}{[{}^{3}{O_{2}}] + β}) [{}^{3}{O_{2}}] - g (1 - \frac{[{}^{3}{O_{2}}]}{[{}^{3}{O_{2}}] (t = 0)}) = 0,

Eq. (5)

\frac{d {[{}^{1}{O_{2}}]}_{rx}}{d t} + (ξ \frac{φ [S_{0}] [{}^{3}{O_{2}}]}{[{}^{3}{O_{2}}] + β}) = 0,

where

S

represents the light source power.

ξ

,

σ

,

δ

,

β

, and

g

are specific PDT photochemical parameters with definitions and magnitudes listed in Table 1.

Table 1

Model parameters used in the macroscopic kinetics equations for Photofrin.

Parameter	Definition	Value	References
$ϵ$ ( ${cm}^{- 1} μ M^{- 1}$ )	Photofrin extinction coefficient	0.0035	18
$ξ$ ( ${cm}^{2} s^{- 1} {mW}^{- 1}$ )	Specific oxygen consumption rate	$3.7 \times 10^{- 3}$	18,32,33
$σ$ ( $μ M^{- 1}$ )	Specific photobleaching ratio	$7.6 \times 10^{- 5}$	18,32
$β$ ( $μ M$ )	Oxygen quenching threshold concentration	11.9	18,32
$δ$ ( $μ M$ )	Low-concentration correction	33	18,34
$g$ ( $μ M / s$ )	Macroscopic oxygen maximum perfusion rate	0.76	18
${[{}^{3}{O_{2}}]}_{0}$ ( $μ M$ )	Initial ground-state oxygen concentration	40	35,36

The calculation of the equations was performed using MATLAB^® R2015b (Natick, Massachusetts).

2.7.

Statistical Analysis

The tumor regrowth rate ( $k$ ) and CI of each condition have been expressed as the $mean \pm standard deviation$ . Differences in the tumor regrowth rate and CI between each PDT group and control group were analyzed by Mann–Whitney test. Analyses were carried out using SPSS 22.0 software. For all tests, $p < 0.05$ level (95% confidence level) was considered statistically significant.

3. Results

3.1.

Tumor Regrowth Curves

To compare the regrowth rate between different tumors, volumes were normalized so that the initial volumes on day 0 were matched to be the same among all tumors, $\sim 12 {mm}^{3}$ . Figure 3 shows the normalized tumor volume versus time (in days) for the 11 treatment groups and the control group along with the fits to the data with an exponential growth equation. These exponential fits to the data determine the value of $k$ for each treatment group of mice. The statistical analyses showed a reduction of tumor regrowth rate for all treated groups compared to the control (all with $p < 0.05$ ). Some mice within the same group had different CIs. For the group with a total fluence of $135 J / {cm}^{2}$ and $φ = 50 mW / {cm}^{2}$ , one out of four (25%) mice showed a complete response (no tumor regrowth for 14 days post-PDT). This is reflected in Table 2 under the column labeled CI variation count. At the same total fluence, all mice (100%) for the group treated at $φ = 75 mW / {cm}^{2}$ and no mice (0%) for the group treated at $φ = 150 mW / {cm}^{2}$ exhibited a complete response. At $200 J / {cm}^{2}$ , no tumors showed a complete response with either 50 or $150 mW / {cm}^{2}$ . At $250 J / {cm}^{2}$ , three out of five (60%) tumors showed a complete response at $50 mW / {cm}^{2}$ , two out of four (50%) had a complete response at $75 mW / {cm}^{2}$ , and two out of nine (22.2%) had a complete response with $150 mW / {cm}^{2}$ (Table 2).

Fig. 3

Exponential fitting of the change in normalized tumor volume over time after PDT. Before PDT, all the tumor volumes in each group had no significant difference between each other ( $p > 0.05$ ). After PDT, comparing with the control group, each PDT group could effectively inhibit the regrowth of tumor ( $p < 0.05$ ).

Table 2

The drug concentration, tumor regrowth rate, CI, and [O21]rx of each PDT group.

Groups	Fluence rate (mW/cm2)	Fluence (J/cm2)	PDT dose (μM J/cm2)	Mean Photofrin (μM)	Photofrin variationc (μM)	Mean [O21]rxb (mM)	[O21]rx variationc (mM)	k (days−1)	CI variation countd	CI
1	50	50	$126.5 \pm 27.6$	$3.5 \pm 0.9$	—	$0.35 \pm 0.06$	—	$0.31 \pm 0.02$	$0 / 4$	$0.24 \pm 0.05$
2	75	50	$235.5 \pm 55.8$	$6.0 \pm 1.3$	—	$0.40 \pm 0.03$	—	$0.30 \pm 0.07$	$0 / 5$	$0.28 \pm 0.17$
3	150	50	$199.9 \pm 112.8$	$4.9 \pm 2.6$	—	$0.24 \pm 0.02$	—	$0.32 \pm 0.01$	$0 / 4$	$0.23 \pm 0.02$
4	50	135	$229.7 \pm 38.8$	$3.0 \pm 1.0$	$2.5 / 4.3$ a	$0.66 \pm 0.19$	$0.57 / 0.91$ e	$0.19 \pm 0.13$	$1 / 4$	$0.54 \pm 0.21$
5	75	135	$661.9 \pm 327.1$	$7.2 \pm 3.2$	—	$1.05 \pm 0.21$	—	0	$5 / 5$	1
6	150	135	$322.3 \pm 42.0$	$3.5 \pm 0.4$	—	$0.53 \pm 0.03$	—	$0.23 \pm 0.04$	$0 / 4$	$0.45 \pm 0.09$
7	50	200	$204.8 \pm 37.1$	$2.2 \pm 1.1$	—	$0.63 \pm 0.29$	—	$0.21 \pm 0.01$	$0 / 5$	$0.49 \pm 0.02$
8	150	200	$417.8 \pm 137.9$	$3.8 \pm 1.0$	—	$0.76 \pm 0.09$	—	$0.17 \pm 0.02$	$0 / 5$	$0.59 \pm 0.04$
9	50	250	$362.1 \pm 124.0$	$3.4 \pm 1.1$	$3.2 / 3.7$	$1.03 \pm 0.35$	$0.97 / 1.09$	$0.21 \pm 0.18$	$3 / 5$	$0.77 \pm 0.31$
10	75	250	$263.0 \pm 51.9$	$2.4 \pm 0.4$	$2.1 / 2.7$	$0.74 \pm 0.11$	$0.65 / 0.81$	$0.18 \pm 0.11$	$2 / 4$	$0.57 \pm 0.31$
11	150	250	$393.0 \pm 188.2$	$3.1 \pm 1.8$	$2.5 / 5.3$ a	$0.81 \pm 0.27$	$0.69 / 1.04$ e	$0.20 \pm 0.10$	$2 / 9$	$0.52 \pm 0.23$
Control	0	0	0	0	—	0	—	$0.42 \pm 0.07$	—	0

Note: In group 4, 9, 10 and 11, cured and not cured tumors coexisted in the same PDT group.

^aPhotofrin concentration in tumors with complete and partial response is significantly different, p<0.05.

^bThe average value of reacted singlet oxygen concentration at 3-mm depth of the tumor among mice used.

^cVariation is shown for (tumors with complete response)/(tumors with partial response).

^dCI variation is shown as (tumors with a complete response)/(total number of tumors in the group).

^e[O21]rx concentration in tumors with complete and partial response is significantly different, p<0.05.

3.2.

Variation in Photofrin Concentration

With the same administered Photofrin dose of $5 mg / kg$ , the uptake of photosensitizer in tumors varied from 0.95 to $10.53 μ M$ ; the mean concentration was $3.93 \pm 2.15 μ M$ ; and the median value was $3.54 μ M$ . As shown in Table 2, the mean uptake of Photofrin in tumors of each treatment group ranged from 2 to $7 μ M$ . Treatment groups with mixed tumor control results (complete and partial responses) reflected this difference in Photofrin concentration as well. To evaluate the in-vivo measurements, ex-vivo measurements of the Photofrin concentration were performed in another set of mice administered the same amounts of the photosensitizer. As shown in Fig. 2(c), the in-vivo measured concentrations have been compared with those obtained ex-vivo; each individual data point represents the average value of three measurements with the standard error of the mean. A linear fit of $y = 0.98 x$ (solid black line) with $R^{2} = 0.99$ shows a good correlation of the in-vivo and ex-vivo data.

3.3.

Variation in Calculated [¹O₂]_rx

Figure 4 shows that the spatial distribution of calculated ${[{}^{1}{O_{2}}]}_{rx}$ along the tumor depth. With increasing tumor depth, ${[{}^{1}{O_{2}}]}_{rx}$ is decreased due to the distribution of treatment light inside the tumor. Figure 5 shows the impact of fluence on mean ${[{}^{1}{O_{2}}]}_{rx}$ , which is the mean reacted singlet oxygen produced at 3-mm depth inside the tumor based on a previous study.¹⁶ Although the amount of ${[{}^{1}{O_{2}}]}_{rx}$ is increased with an increase in fluence, the generation of ${[{}^{1}{O_{2}}]}_{rx}$ varied with $φ$ . For example, at the lowest fluence of $50 J / {cm}^{2}$ , the ${[{}^{1}{O_{2}}]}_{rx}$ calculated for groups treated with 50, 75, and $150 mW / {cm}^{2}$ was 0.35, 0.40, and 0.24, respectively. In all 11 PDT treatment groups, the highest amount of ${[{}^{1}{O_{2}}]}_{rx}$ was found in the group treated with $75 mW / {cm}^{2}$ to a fluence of $135 J / {cm}^{2}$ , while the lowest ${[{}^{1}{O_{2}}]}_{rx}$ was found in the group treated with $150 mW / {cm}^{2}$ to a fluence of $50 J / {cm}^{2}$ . Table 2 shows the variation of ${[{}^{1}{O_{2}}]}_{rx}$ within PDT treatment groups with a mix of tumors that exhibited complete and partial treatment response (group 4, 9, 10, and 11). The mean calculated ${[{}^{1}{O_{2}}]}_{rx}$ between the tumors that showed complete response and partial response was significantly different in groups 4 and 11 ( $p < 0.05$ ), whereas there was no significant difference ( $p > 0.05$ ) in groups 9 and 10.

Fig. 4

Spatial distribution of reacted singlet oxygen ( ${[{}^{1}{O_{2}}]}_{rx}$ ) in tumors of each PDT group. Here, the Photofrin concentration was the mean calculated concentration of each PDT condition.

Fig. 5

The impact of PDT fluence on reacted singlet oxygen generation at a 3-mm tumor depth. The black solid line shows the best-fit to the data using functional form, $y = 0.003 x$ with a goodness of fit of $R^{2} = 0.57$ . The gray area shows the upper and lower bounds of the fit with their 95% confidence level.

3.4.

Correlation Between Three Photodynamic Therapy Metrics and Cure Index

The relationship between CI and three PDT dose metrics, i.e., fluence, PDT dose, and ${[{}^{1}{O_{2}}]}_{rx}$ , was further investigated. Figures 6(a)–6(c) show their correlations. The solid lines show the best linear fit to the data using functional forms $y = 0.003 x$ for (a), $y = 0.001 x$ for (b), and $y = 0.795 x$ for (c) with a goodness of fit of $R^{2} = 0.35$ , 0.79, and 0.93 for (a), (b), and (c), respectively. The gray area shows the upper and lower bounds of the fit with their 95% confidence levels. As shown in Fig. 6(a), the relationship between fluence and CI was poor with $R^{2} = 0.35$ . At the same fluence, CI varied greatly with $φ$ .

Fig. 6

Linear fitting of tumor CI versus (a) fluence, (b) PDT dose, and (c) mean ${[{}^{1}{O_{2}}]}_{rx}$ at a 3-mm tumor depth in each group. The gray area in each figure shows the upper and lower bounds of the fit with 95% confidence level. The solid lines show the best-fit to the data using functional forms, $y = 0.003 x$ for (a), $y = 0.001 x$ for (b), and $y = 0.795 x$ for (c) with goodness of fit of $R^{2} = 0.35$ , 0.79, and 0.93 for (a), (b), and (c), respectively.

The relationship of CI and ${[{}^{1}{O_{2}}]}_{rx}$ could be fit with a linear function of the form $y = 0.795 x$ with a goodness of fit of $R^{2} = 0.93$ [see Fig. 6(c)]. With increasing amounts of ${[{}^{1}{O_{2}}]}_{rx}$ , CI was also improved. Among the 11 PDT treatment conditions, $135 J / {cm}^{2}$ at $75 mW / {cm}^{2}$ resulted in the highest calculated ${[{}^{1}{O_{2}}]}_{rx}$ concentration as well as the best control of tumor regrowth. However, there was one PDT treatment condition ( $150 mW / {cm}^{2}$ , $250 J / {cm}^{2}$ ) that did not follow this trend. While the amount of ${[{}^{1}{O_{2}}]}_{rx}$ of this condition was higher than that of groups treated with $135 J / {cm}^{2}$ at $50 mW / {cm}^{2}$ , $250 J / {cm}^{2}$ at $75 mW / {cm}^{2}$ , and $200 J / {cm}^{2}$ at $150 mW / {cm}^{2}$ , the CI was lower. Details of the results are listed in Table 2. Upon further investigation of the tumor CI and the mean calculated ${[{}^{1}{O_{2}}]}_{rx}$ for individual tumors that are a part of treatment groups where tumors exhibited both complete and partial response from treatment (Table 2), the results show that when the calculated ${[{}^{1}{O_{2}}]}_{rx}$ was larger than 1 mM, tumors resulted in a complete response.

4. Discussion

The three elements of PDT (light, drug, and ${}^{3}{O_{2}}$ ) are dynamically changing and interacting with each other during PDT, which makes PDT dosimetry very complex and challenging. As the primary cause of cell damage during PDT, ${}^{1}{O_{2}}$ has gained special attention as the dosimetric index both by direct experimental methods as well as indirect mathematical modeling methods, but it is still difficult for direct measurements in in-vivo PDT studies.³⁷^,³⁸ Here, the predictive potential of calculated ${[{}^{1}{O_{2}}]}_{rx}$ using an indirect method for PDT outcome was evaluated.

Based on four different PDT light doses, plus a series of $φ$ , the effectiveness of Photofrin-PDT on tumor regrowth inhibition was investigated. Although the results showed that all of these PDT conditions were effective in controlling the regrowth of tumors, CI varied with changes in total fluence and $φ$ . Furthermore, unlike treatment with the monomeric photosensitizers, benzoporphyrin derivative monoacid ring A- or 2-(1-Hexyloxyethyl)-2-devinylpyropheophorbide-mediated PDT,²²^,²³ Photofrin-mediated PDT-treated tumors exhibited both complete and partial response in the same treatment group at certain PDT conditions.

In this study, tumors with complete response and partial response coexisted in the group treated with $135 J / {cm}^{2}$ at $50 mW / {cm}^{2}$ (group 4) as well as $250 J / {cm}^{2}$ at 50, 75, and $150 mW / {cm}^{2}$ (groups 9, 10, and 11). Varying results for the same PDT condition were also reported by Sitnik and Henderson³⁹ and Wang et al.¹⁹ for a wide range of $φ$ from 10 to $200 mW / {cm}^{2}$ . Furthermore, in this study, the linear fitting of fluence with CI showed that the relationship between them was very poor with $R^{2} = 0.35$ . All these data suggest the inaccuracy of the fluence alone for the prediction of PDT outcome.

PDT dose is a metric that incorporates the light dose absorbed by the photosensitizer in tumors.¹² As indicated in this study, PDT dose was more correlated to CI than fluence with a goodness of fit of $R^{2} = 0.79$ . However, CI was not always improved with increasing PDT dose. For instance, a PDT dose of $204.84 μ M J / {cm}^{2}$ resulted a higher CI ( $0.49 \pm 0.02$ ) than that of 235.47 and $322.21 μ M J / {cm}^{2}$ , which resulted in a CI of $0.28 \pm 0.17$ and $0.45 \pm 0.09$ , respectively. A PDT dose of $362.13 μ M J / {cm}^{2}$ also resulted in a higher CI ( $0.77 \pm 0.31$ ) than that of 417.76 and $392.98 μ M J / {cm}^{2}$ , which resulted in a CI of $0.59 \pm 0.31$ and $0.52 \pm 0.23$ , respectively. To some extent, the differences of the in-vivo photosensitizer concentration from the administrated dose could explain the different predictive ability between PDT dose and fluence. As demonstrated in this study, at a given Photofrin dose of $5 mg / kg$ , the uptake of Photofrin varied from tumor to tumor. In PDT conditions of $135 J / {cm}^{2}$ at $50 mW / {cm}^{2}$ (group 4) and $250 J / {cm}^{2}$ at $150 mW / {cm}^{2}$ (group 11), the in-vivo Photofrin concentration within the same group between tumors with complete and partial response was significantly different ( $p < 0.05$ ). The tendency to have large variations in Photofrin uptake by the tumors is in agreement with the other studies.¹⁸^,⁴⁰

The calculated ${[{}^{1}{O_{2}}]}_{rx}$ did not follow the trend of increasing light dose or $φ$ . Among the 11 PDT treatment groups with light ranging from 50 to $250 J / {cm}^{2}$ , the highest amount of ${[{}^{1}{O_{2}}]}_{rx}$ was found with the group treated to $135 J / {cm}^{2}$ with $75 mW / {cm}^{2}$ . This is consistent with the $φ$ dependence of PDT oxygen consumption. Higher $φ$ tends to create hypoxic environments that limit the production of ${}^{1}{O_{2}}$ .⁴¹ Furthermore, this treatment group had the highest measured in-vivo Photofrin concentration at $7.2 \pm 3.2 μ M$ . These results emphasize the importance of explicit dosimetry of light, ${}^{3}{O_{2}}$ , and photosensitizer concentration to utilize ${[{}^{1}{O_{2}}]}_{rx}$ as a predictive dosimetric quantity for outcome.

The potential of ${[{}^{1}{O_{2}}]}_{rx}$ as a predictor for PDT efficiency was evaluated by analyzing the correlation between ${[{}^{1}{O_{2}}]}_{rx}$ and CI. In most PDT conditions, the amount of ${[{}^{1}{O_{2}}]}_{rx}$ calculated was consistent with CI; CI improved with increasing ${[{}^{1}{O_{2}}]}_{rx}$ . The correlation of ${[{}^{1}{O_{2}}]}_{rx}$ with CI in treatment groups 4, 9,10, and 11 was further analyzed due to the mixed response from tumors. The results indicated that when ${[{}^{1}{O_{2}}]}_{rx}$ was more than 0.81 mM, tumors exhibited a complete response in all PDT conditions except with a fluence of $250 J / {cm}^{2}$ at $50 mW / {cm}^{2}$ . In the same condition, when the mean ${[{}^{1}{O_{2}}]}_{rx}$ concentration was of 0.97 mM, only partial tumor control was achieved. All these results indicated that the threshold dose of ${[{}^{1}{O_{2}}]}_{rx}$ to induce tumor control was more than 1 mM. Previous studies demonstrated the threshold dose of ${[{}^{1}{O_{2}}]}_{rx}$ to produce necrosis in the same RIF tumor model was about $0.74 \pm 0.25 mM$ for Photofrin-PDT (Table 1).⁴² As for inducing a tumor cure for PDT, it is reasonable to believe that the mean values of ${[{}^{1}{O_{2}}]}_{rx}$ are larger than that necessary to induce tumor necrosis. However, further studies with a wide range of PDT fluences and a larger sample size are needed to validate these results.

Among the three metrics, i.e., fluence, PDT dose, and ${[{}^{1}{O_{2}}]}_{rx}$ , ${[{}^{1}{O_{2}}]}_{rx}$ had the best correlation with CI with a goodness of fit of $R^{2} = 0.93$ . These results indicated that ${[{}^{1}{O_{2}}]}_{rx}$ better predicts PDT outcome than fluence or PDT dose.

5. Conclusions

In this empirical macroscopic ${}^{1}{O_{2}}$ model, a set of equations were used to simplify the energy transfer processes in type II PDT, and calculated ${[{}^{1}{O_{2}}]}_{rx}$ was proposed as a dosimetric quantity. These equations can be solved by inputting parameters including the light fluence (rate), tissue optical properties, photosensitizer concentration, and photochemical parameters. The predictive ability of calculated ${[{}^{1}{O_{2}}]}_{rx}$ , as a dosimetric quantity, was evaluated by analyzing the relationship between ${[{}^{1}{O_{2}}]}_{rx}$ and CI due to PDT. The relationship between CI with two other commonly used metrics, fluence, and PDT dose, was also investigated. Preliminary results showed that compared with fluence or PDT dose, ${[{}^{1}{O_{2}}]}_{rx}$ was the most accurate metric and could serve as a better dosimetric quantity. Further studies are required to validate these results and establish the threshold dose of ${[{}^{1}{O_{2}}]}_{rx}$ to induce a complete response in this RIF tumor model. Furthermore, other factors, such as the photosensitizer photobleaching and the ${}^{3}{O_{2}}$ variation during PDT should be considered to improve the predictive ability of this macroscopic model.

Acknowledgments

We would like to thank Dr. Jarod C. Finlay for his helpful discussions regarding the theory and analysis and Dr. Theresa M. Busch for helping us with the mouse study protocols. This work was supported by grants from the National Institutes of Health Nos. R01 CA154562 and P01 CA87971.

References

1.

T. J. Dougherty et al., “Photodynamic therapy,” J. Natl. Cancer Inst., 90 (12), 889 –905 (1998). http://dx.doi.org/10.1093/jnci/90.12.889 Google Scholar

2.

S. A. Gross and H. C. Wolfsen, “The role of photodynamic therapy in the esophagus,” Gastrointest. Endoscopy Clin. N. Am., 20 (1), 35 –53 (2010). http://dx.doi.org/10.1016/j.giec.2009.07.008 Google Scholar

3.

L. Corti et al., “Long-term survival of patients treated with photodynamic therapy for carcinoma in situ and early non-small-cell lung carcinoma,” Lasers Surg. Med., 39 (5), 394 –402 (2007). http://dx.doi.org/10.1002/(ISSN)1096-9101 Google Scholar

4.

I. J. Vergilis-Kalner and J. L. Cohen, “The use of photodynamic therapy as chemoprevention for the treatment of actinic keratoses and reduction in the number of non-melanoma skin cancers,” J. Drugs Dermatol., 12 (10), 1085 –1086 (2013). Google Scholar

5.

P. K. Lee and A. Kloser, “Current methods for photodynamic therapy in the US: comparison of MAL/PDT and ALA/PDT,” J. Drugs Dermatol., 12 (8), 925 –930 (2013). Google Scholar

6.

M. J. Bader et al., “Photodynamic therapy of bladder cancer—a phase I study using hexaminolevulinate (HAL),” Urol. Oncol. Semin. Orig. Invest., 31 (7), 1178 –1183 (2013). http://dx.doi.org/10.1016/j.urolonc.2012.02.007 Google Scholar

7.

C. B. Simone and K. A. Cengel, “Photodynamic therapy for lung cancer and malignant pleural mesothelioma,” Semin. Oncol., 41 (6), 820 –830 (2014). http://dx.doi.org/10.1053/j.seminoncol.2014.09.017 Google Scholar

8.

T. C. Zhu et al., “Singlet oxygen threshold doses for PDT,” Photonics Lasers Med., 4 (1), 59 –71 (2015). http://dx.doi.org/10.1515/plm-2014-0037 Google Scholar

9.

A. Greer, “Christopher Foote’s discovery of the role of singlet oxygen [

{}^{1}{O_{2}}

(

{{}^{1}Δ}_{g}

)] in photosensitized oxidation reactions,” Acc. Chem. Res., 39 (11), 797 –804 (2006). http://dx.doi.org/10.1021/ar050191g Google Scholar

10.

C. S. Foote, “Definition of type I and type II photosensitized oxidation,” Photochem. Photobiol., 54 (5), 659 (1991). http://dx.doi.org/10.1111/php.1991.54.issue-5 Google Scholar

11.

A. Juarranz et al., “Photodynamic therapy of cancer. Basic principles and applications,” Clin. Transl. Oncol., 10 (3), 148 –154 (2008). http://dx.doi.org/10.1007/s12094-008-0172-2 Google Scholar

12.

B. C. Wilson, M. S. Patterson and L. Lilge, “Implicit and explicit dosimetry in photodynamic therapy: a new paradigm,” Laser Med. Sci., 12 (3), 182 –199 (1997). http://dx.doi.org/10.1007/BF02765099 Google Scholar

13.

J. D. Vollet-Filho et al., “Non-homogeneous liver distribution of photosensitizer and its consequence for photodynamic therapy outcome,” Photodiagn. Photodyn. Ther., 7 (3), 189 –200 (2010). http://dx.doi.org/10.1016/j.pdpdt.2010.07.002 Google Scholar

14.

T. G. Vulcan et al., “Comparison between isotropic and nonisotropic dosimetry systems during intraperitoneal photodynamic therapy,” Lasers Surg. Med., 26 (3), 292 –301 (2000). http://dx.doi.org/10.1002/(ISSN)1096-9101 Google Scholar

15.

T. M. Sitnik, J. A. Hampton and B. W. Henderson, “Reduction of tumour oxygenation during and after photodynamic therapy in vivo: effects of fluence rate,” Br. J. Cancer, 77 (9), 1386 –1394 (1998). http://dx.doi.org/10.1038/bjc.1998.231 Google Scholar

16.

R. Penjweini et al., “Explicit dosimetry for 2-(1-hexyloxyethyl)-2-devinyl pyropheophorbide-a-mediated photodynamic therapy: macroscopic singlet oxygen modeling,” J. Biomed. Opt., 20 (12), 128003 (2015). http://dx.doi.org/10.1117/1.JBO.20.12.128003 Google Scholar

17.

T. C. Zhu, “Dosimetry in pleural photodynamic therapy,” J. Natl. Compr. Cancer Network, 10 (Suppl. 2), S60 –S64 (2012). Google Scholar

18.

K. K. Wang et al., “Explicit dosimetry for photodynamic therapy: macroscopic singlet oxygen modeling,” J. Biophotonics, 3 (5–6), 304 –318 (2010). http://dx.doi.org/10.1002/jbio.v3:5/6 Google Scholar

19.

H. W. Wang et al., “Effect of photosensitizer dose on fluence rate responses to photodynamic therapy,” Photochem. Photobiol., 83 (5), 1040 –1048 (2007). http://dx.doi.org/10.1111/php.2007.83.issue-5 Google Scholar

20.

R. C. Mesquita et al., “Tumor blood flow differs between mouse strains: consequences for vasoresponse to photodynamic therapy,” PLoS One, 7 (5), e37322 (2012). http://dx.doi.org/10.1371/journal.pone.0037322 POLNCL 1932-6203 Google Scholar

21.

T. M. Busch et al., “Fluence rate-dependent intratumor heterogeneity in physiologic and cytotoxic responses to Photofrin photodynamic therapy,” Photochem. Photobiol. Sci., 8 (12), 1683 –1693 (2009). http://dx.doi.org/10.1039/b9pp00004f Google Scholar

22.

M. M. Kim, R. Penjweini and T. C. Zhu, “In vivo outcome study of BPD-mediated PDT using a macroscopic singlet oxygen model,” Proc. SPIE, 9308 93080A (2015). http://dx.doi.org/10.1117/12.2077803 Google Scholar

23.

R. Penjweini, M. M. Kim and T. C. Zhu, “In-vivo outcome study of HPPH mediated PDT using singlet oxygen explicit dosimetry (SOED),” Proc. SPIE, 9308 93080N (2015). http://dx.doi.org/10.1117/12.2076441 PSISDG 0277-786X Google Scholar

24.

H. Qiu et al., “Dosimetry study of PHOTOFRIN-mediated photodynamic therapy in a mouse tumor model,” Proc. SPIE, 9694 96940T (2016). http://dx.doi.org/10.1117/12.2211169 PSISDG 0277-786X Google Scholar

25.

A. Dimofte, J. C. Finlay and T. C. Zhu, “A method for determination of the absorption and scattering properties interstitially in turbid media,” Phys. Med. Biol., 50 (10), 2291 –2311 (2005). http://dx.doi.org/10.1088/0031-9155/50/10/008 Google Scholar

26.

T. C. Zhu, A. Lu and Y. H. Ong, “An improved analytic function for predicting light fluence rate in circular fields on a semi-infinite geometry,” Proc. SPIE, 9706 97061D (2016). http://dx.doi.org/10.1117/12.2213052 PSISDG 0277-786X Google Scholar

27.

J. C. Finlay, T. C. Zhu and A. Dimofte, “Diffuse reflectance spectra measured in vivo in human tissues during Photofrin-mediated pleural photodynamic therapy,” Proc. SPIE, 6139 61390O (2006). http://dx.doi.org/10.1117/12.647016 PSISDG 0277-786X Google Scholar

28.

J. C. Finlay et al., “Porphyrin bleaching and PDT-induced spectral changes are irradiance dependent in ALA-sensitized normal rat skin in vivo,” Photochem. Photobiol., 73 (1), 54 –63 (2001). http://dx.doi.org/10.1562/0031-8655(2001)073<0054:PBAPIS>2.0.CO;2 Google Scholar

29.

K. K. Wang et al., “Optimization of physiological parameter for macroscopic modeling of reacted singlet oxygen concentration in an in-vivo model,” Proc. SPIE, 7164 71640O (2009). http://dx.doi.org/10.1117/12.809024 PSISDG 0277-786X Google Scholar

30.

T. C. Zhu et al., “Macroscopic modeling of the singlet oxygen production during PDT,” Proc. SPIE, 6427 642708 (2007). http://dx.doi.org/10.1117/12.701387 PSISDG 0277-786X Google Scholar

31.

T. C. Zhu, B. Liu and R. Penjweini, “Study of tissue oxygen supply rate in a macroscopic photodynamic therapy singlet oxygen model,” J. Biomed. Opt., 20 (3), 038001 (2015). http://dx.doi.org/10.1117/1.JBO.20.3.038001 Google Scholar

32.

I. Georgakoudi, M. G. Nichols and T. H. Foster, “The mechanism of Photofrin photobleaching and its consequences for photodynamic dosimetry,” Photochem. Photobiol., 65 (1), 135 –144 (1997). http://dx.doi.org/10.1111/php.1997.65.issue-1 Google Scholar

33.

M. G. Nichols and T. H. Foster, “Oxygen diffusion and reaction kinetics in the photodynamic therapy of multicell tumour spheroids,” Phys. Med. Biol., 39 (12), 2161 –2181 (1994). http://dx.doi.org/10.1088/0031-9155/39/12/003 Google Scholar

34.

J. S. Dysart, G. Singh and M. S. Patterson, “Calculation of singlet oxygen dose from photosensitizer fluorescence and photobleaching during mTHPC photodynamic therapy of MLL cells,” Photochem. Photobiol., 81 (1), 196 –205 (2005). http://dx.doi.org/10.1562/2004-07-23-RA-244.1 Google Scholar

35.

J. P. Whiteley, D. J. Gavaghan and C. E. Hahn, “Mathematical modelling of oxygen transport to tissue,” J. Math. Biol., 44 (6), 503 –522 (2002). http://dx.doi.org/10.1007/s002850200135 Google Scholar

36.

A. Carreau et al., “Why is the partial oxygen pressure of human tissues a crucial parameter? Small molecules and hypoxia,” J. Cell. Mol. Med., 15 (6), 1239 –1253 (2011). http://dx.doi.org/10.1111/jcmm.2011.15.issue-6 Google Scholar

37.

M. Niedre, M. S. Patterson and B. C. Wilson, “Direct near-infrared luminescence detection of singlet oxygen generated by photodynamic therapy in cells in vitro and tissues in vivo,” Photochem. Photobiol., 75 (4), 382 –391 (2002). http://dx.doi.org/10.1562/0031-8655(2002)0750382DNILDO2.0.CO2 Google Scholar

38.

K. K. Wang, S. Mitra and T. H. Foster, “A comprehensive mathematical model of microscopic dose deposition in photodynamic therapy,” Med. Phys., 34 (1), 282 –293 (2007). http://dx.doi.org/10.1118/1.2401041 Google Scholar

39.

T. M. Sitnik and B. W. Henderson, “The effect of fluence rate on tumor and normal tissue responses to photodynamic therapy,” Photochem. Photobiol., 67 (4), 462 –466 (1998). http://dx.doi.org/10.1111/php.1998.67.issue-4 Google Scholar

40.

M. Korbelik, “Cellular delivery and retention of Photofrin: II. The effects of human versus mouse and bovine serum,” Photochem. Photobiol., 56 (3), 391 –397 (1992). http://dx.doi.org/10.1111/php.1992.56.issue-3 Google Scholar

41.

T. H. Foster et al., “Oxygen consumption and diffusion effects in photodynamic therapy,” Radiat. Res., 126 (3), 296 –303 (1991). http://dx.doi.org/10.2307/3577919 Google Scholar

42.

T. C. Zhu et al., “Comparison of singlet oxygen threshold dose for PDT,” Proc. SPIE, 8931 89310I (2014). http://dx.doi.org/10.1117/12.2039719 PSISDG 0277-786X Google Scholar

Biography

Haixia Qiu received her MD in 2006 in laser medicine from the Medical College of Chinese People's Liberation Army. She is currently a visiting scholar in the Department of Radiation Oncology at the University of Pennsylvania. Her current research interests include PDT dosimetry, the underlying mechanism of PDT and the clinical application of PDT.

Michele M. Kim received her BA and MS degrees in physics from the University of Pennsylvania in 2012 and is currently a PhD candidate at the University of Pennsylvania in physics while pursuing a certificate in medical physics. Her research topics include preclinical and clinical PDT dosimetry.

Rozhin Penjweini received her PhD degree in 2012 in physics from the University of Vienna. She is currently a postdoctoral researcher in the Department of Radiation Oncology at the University of Pennsylvania. Her current research interest is in vivo explicit singlet oxygen dosimetry for PDT. She also has practical experience in various fluorescence microscopy techniques for studying the structure, transport, and stability of nanomedicines in living cells.

Timothy C. Zhu received his PhD degree in 1991 in physics from Brown University. He is currently a professor in the Department of Radiation Oncology at the University of Pennsylvania. His current research interests include in vivo explicit PDT dosimetry, singlet oxygen explicit dosimetry, integrated system for interstitial and intracavitory PDT, diffuse optical tomography, and external beam radiation transport.

Citation Download Citation

Haixia Qiu, Michele M. Kim, Rozhin Penjweini, and Timothy C. Zhu "Macroscopic singlet oxygen modeling for dosimetry of Photofrin-mediated photodynamic therapy: an in-vivo study," Journal of Biomedical Optics 21(8), 088002 (23 August 2016). https://doi.org/10.1117/1.JBO.21.8.088002

Published: 23 August 2016

Access the abstract

JOURNAL ARTICLE
9 PAGES

DOWNLOAD PAPER SAVE TO MY LIBRARY

GET CITATION

CITATIONS

Cited by 42 scholarly publications.

Explore citations on Lens.org

RIGHTS & PERMISSIONS

Get copyright permission Get copyright permission on Copyright Marketplace

KEYWORDS

Photodynamic therapy

Tumors

Oxygen

In vivo imaging

Luminescence

Optical properties

Tissue optics

1.

Introduction

2.

Materials and Methods