Translator Disclaimer
7 January 2020 Reactive oxygen species explicit dosimetry to predict tumor growth for benzoporphyrin derivative-mediated vascular photodynamic therapy
Author Affiliations +

Photodynamic therapy (PDT) is a well-established treatment modality for cancer and other malignant diseases; however, quantities such as light fluence and PDT dose do not fully account for all of the dynamic interactions between the key components involved. In particular, fluence rate (ϕ) effects, which impact the photochemical oxygen consumption rate, are not accounted for. In this preclinical study, reacted reactive oxygen species ([ROS]rx) was investigated as a dosimetric quantity for PDT outcome. The ability of [ROS]rx to predict the cure index (CI) of tumor growth, CI  =  1  −  k  /  kctr, where k and kctr are the growth rate of tumor under PDT study and the control tumor without PDT, respectively, for benzoporphyrin derivative (BPD)-mediated PDT, was examined. Mice bearing radiation-induced fibrosarcoma (RIF) tumors were treated with different in-air fluences (Φ  =  22.5 to 166.7  J  /  cm2) and in-air fluence rates (ϕair  =  75 to 250   mW  /  cm2) with a BPD dose of 1  mg  /  kg and a drug-light interval (DLI) of 15 min. Treatment was delivered with a collimated laser beam of 1-cm-diameter at 690 nm. Explicit measurements of in-air light fluence rate, tissue oxygen concentration, and BPD concentration were used to calculate for [ROS]rx. Light fluence rate at 3-mm depth (ϕ3  mm), determined based on Monte-Carlo simulations, was used in the calculation of [ROS]rx at the base of tumor. CI was used as an endpoint for three dose metrics: light fluence, PDT dose, and [ROS]rx. PDT dose was defined as the product of the time-integral of photosensitizer concentration and ϕ3  mm. Preliminary studies show that [ROS]rx best correlates with CI and is an effective dosimetric quantity that can predict treatment outcome. The threshold dose for [ROS]rx for vascular BPD-mediated PDT using DLI of 15 min is determined to be 0.26 mM and is about 3.8 times smaller than the corresponding value for conventional BPD-mediated PDT using DLI of 3 h.



Photodynamic therapy (PDT) is a treatment modality for cancer and other localized diseases. PDT is not only “dynamic” but also multifaceted.1,2 PDT incorporates light, photosensitizer, and oxygen to create reactive oxygen species (ROS) to kill cells. Unlike radiotherapy and chemotherapy, PDT causes fewer side effects, as it does not involve ionizing radiation and can be well-localized.3,4 PDT is uniquely advantageous compared to other treatment modalities, as it is also associated with fast postoperative recovery and better cosmetic outcome. However, widespread use of PDT has been stymied due to the difficulty in accurately quantifying the dose. Furthermore, assessment of PDT efficacy is difficult due to the lack of a well-defined dose metric that accurately predicts biological response. ROS are acceptable to be the cytotoxic agents causing therapeutic outcome in PDT. Direct detection of ROS can provide the most accurate quantity to guide treatments and predict treatment outcomes. However, in vivo detection of ROS during clinical PDT is very challenging due to its weak signal and short lifetime. To overcome this, a macroscopic reactive oxygen species explicit dosimetry (ROSED) model was recently developed to calculate for the accumulated reacted ROS concentration ([ROS]rx) that is predictive of PDT treatment outcome.58

Benzoporphyrin derivative monoacid ring A (BPD-MA, trademark Visudyne®) is a commonly used photosensitizer that has been approved by the U.S. Food and Drug Administration for the treatment of wet age-related macular degeneration.9 Using a macroscopic ROSED model of light fluence (rate), BPD drug concentration, and tissue oxygen concentration ([O23]), [ROS]rx can be determined to evaluate its effectiveness as a dosimetrical predictor for BPD-mediated PDT outcome. Vascular-targeted PDT can be achieved using a short (15 min) drug-light interval (DLI), defined as the time interval between the PS drug injection and the start of PDT treatment.10,11 By inducing vascular shutdown, nutrient supply and removal of metabolic waste is halted, which results in radiation-induced fibrosarcoma (RIF) tumor cell death. This is beneficial as targeting tumor vasculature is easier to access, more efficient in cancer cell killing, and has a lower likelihood of developing drug resistance.

This study, to our knowledge, is the first study to investigate the relationship between various dose metrics (fluence, PDT dose, and [ROS]rx) and cure index (CI) at 14 days in an in vivo mouse model for BPD-mediated vascular PDT. ROSED was performed to evaluate the treatment outcomes of BPD-mediated vascular PDT in mice bearing RIF tumors. The major photochemical parameters in the macroscopic ROSED model have been investigated and determined for the photosensitizer BPD-MA for DLI of 3 h12,13 and are found to be similar for DLI of 15 min.14


Materials and Methods


Tumor Model

RIF cells were cultured and 30  μl were injected at 1×107  cells/ml in the right shoulders of 6- to 8-week-old female C3H mice (NCI-Frederick, Frederic), as described previously.5,1517 The resulting RIF tumors are subcutaneous. Animals were under the care of the University of Pennsylvania Laboratory Animal Resources. All studies were approved by the University of Pennsylvania Institutional Animal Care and Use Committee. Tumors were treated when they were 3 to 5 mm in diameter. The fur of the tumor region was clipped prior to cell inoculation, and the treatment area was depilated with Nair (Church and Dwight Co., Inc., Ewing, New Jersey) at least 24 h before measurements. Mice were provided a chlorophyll-free (alfalfa-free) rodent diet (Harlan Laboratories Inc., Indianapolis, Indiana, United States) starting for at least 10 days prior to treatment to eliminate the fluorescence signal from chlorophyll-breakdown products, which have a similar emission range to the BPD fluorescence spectra used to determine the concentration of BPD in the tumor. During the whole treatment, mice were kept under anesthesia on a heat pad at 38°C [see Fig. 1(a)].

Fig. 1

Experiment setup with the (a) multifiber contact spectroscopy probe, tissue pO2 was recorded before and during PDT treatment. (b) A handheld broadband reflectance spectroscopy contact probe was used to measure the optical properties and drug concentration before and after PDT. (c) Oxylite prooxygen monitor with a fluorescence-based bare-fiber oxygen probe (Oxford Optronix, Oxford, UK)



PDT Treatment Conditions

Treatment delivery was done using an optical fiber with a microlens attachment coupled to a diode laser. A 690-nm laser (B&W Tek Inc., Newark, Delaware) was used for PDT after a 15-min DLI. The in-air fluence rate (ϕair) is defined as the calculated irradiance determined by the laser power divided by the treatment area (1-cm diameter spot size). The in-air fluence was calculated by multiplying the in-air fluence rate by the treatment time. RIF tumor-bearing mice with no photosensitizer and no light excitation were used as controls (n=5). The photochemical parameters are summarized in Table 1.

Table 1

Photochemical parameters for BPD based on Refs. 18 and 19.

Photochemical parameterDefinitionValueReferences
ε (cm1μM1)Photosensitizer extinction coefficient0.078318 and 19
δ (μM)Low-concentration correction3318 and 19
β (μM)Oxygen quenching threshold concentration11.918 and 19
σ (μM1)Specific photobleaching ratio(1.8±0.3)×10518 and 19
ξ (cm2mW1s1)Specific oxygen consumption rate(55±15)×10318 and 19
g (μMs1)Macroscopic maximum oxygen supply rate1.7±0.418 and 19
[ROS]rx,sh (mM)Singlet oxygen threshold dose for tumor regrowth0.26±0.05This study
[D]0 (μMJ/cm2)PDT dose, the product of the drug concentration, and light fluence7.5±1.0This study


Photodynamic Therapy Treatment Protocol

BPD (trademark Visudyne®) at a dosage of 1  mg/kg was injected through the mouse tail vein as described previously.16,19 At a 15-min DLI, superficial irradiation of the tumor was performed with a 690-nm diode laser with a maximum output power of 8 W (B&W Tek Inc., Newark, Delaware). A microlens fiber was coupled to the laser to irradiate the tumor uniformly [see Fig. 1(a)]. Mice were treated with in-air fluence rates (ϕair) of 75 to 250  mW/cm2 and total in-air fluences of 22.5 to 166.7  J/cm2 to induce different PDT outcomes and to assess the reciprocity between BPD concentration and light dose. The “in-air fluence rate” is defined as the calculated irradiance determined by laser power divided by the treatment area. The “in-air fluence” was calculated by multiplying the “in-air fluence rate” by the treatment time. Animals were assigned to six light dose groups (see Table 2). RIF tumor-bearing mice that received neither light irradiation nor BPD were used as controls.

Table 2

In-air light fluence, in-air light fluence rate ϕair, photosensitizer concentrations pre-PDT [BPD]pre, PDT dose at 3-mm depth, reacted oxygen species concentration [ROS]rx at 3-mm depth, tumor regrowth rate k, and CI for each PDT treatment group. Number of mice per group is shown in the second column.

Group ## miceIn-air fluence (J/cm2)ϕair (mW/cm2)Time (s)a[BPD]pre (μM)[O23] (μM)Fluenceb (J/cm2)PDT dosec (μM J/cm2)[ROS]rxd (mM)k (1/days)CIe variance countCI=1−k/kctr

aTotal treatment time.

bLight fluence is determined at 3-mm depth based on measured mean in-air light fluence rate at surface and ϕ/ϕair, Eq. (6), using the mean tissue optical properties: μa=6.9  cm−1, μs′=11  cm−1.18 The actual measured in-air fluence rate at surface can be ±10% different from the nominal values and they are accounted in the fluence (J/cm2) at 3 mm.

cPDT dose is defined as the time integral of the product of the photosensitizer concentration and the light fluence rate (ϕ), [S0](t) is calculated using the parameters in Table 1.

d[ROS]rx is calculated using Eq. (5) using the parameters listed in Table 1.

eCI variance count gives number of mice with no tumor regrowth after PDT, out of the total number of mice in the group.


Oxygen Measurements

The in vivo tissue oxygen partial pressure pO2 was measured during PDT treatment using a phosphorescence-based O23 probe (OxyLite Pro, Oxford Optronix, Oxford, United Kingdom). A bare-fiber-type probe (NX-BF/O/E, Oxford Optronix, Oxford, United Kingdom) was placed inside the tumor at a 3-mm depth from the treatment surface. The O23 concentration ([O23]) was calculated by multiplying the measured pO2 with the O23 solubility in tissue, which is 1.295  μM/mmHg.20,21 Measured [O23]0 and [O23](t) was used to calculate for [ROS]rx using the macroscopic ROSED model.20,22


Measurement of BPD Concentration

Following the DLI of 15 min, measurements of light fluence rate, photosensitizer concentration, and [O23] were performed. BPD fluorescence spectra were obtained using a custom-made multifiber contact probe before and after PDT.1 The probe was connected to a 405-nm laser (Power Technology Inc., Little Rock, Arkansas) for the fluorescence excitation of BPD and a multichannel CCD spectrograph (InSpectrum, Princeton Instruments, Trenton, New Jersey) for the collection of the fluorescence spectra. The in vivo photosensitizer concentration was obtained by comparing the measured BPD spectra with those of phantoms with known photosensitizer concentrations. The attenuation of the fluorescence signal due to the light absorption and scattering by tissues was corrected by applying an empirical correction factor described elsewhere.19 The accuracy of in vivo measurements was validated by ex vivo measurements in separate mice.


Ex Vivo Validation of BPD Concentration

In vivo fluorescence measurements of the photosensitizer concentration as described above were performed for all tumors before PDT. To evaluate the accuracy of the in vivo fluorescence measurements, ex vivo measurements of the BPD concentration were performed in separate set of mice and compared with the BPD concentration determined from in vivo measurements. All five mice were administered BPD at different concentrations between (0.25 to 1.25  mg/kg). In vivo fluorescence measurements were taken from each mouse at 15 min and 3 h after BPD administration. After fluorescence measurements at 3-h time point were taken, mice were euthanized and the tumors were excised, protected from light, and stored at 80°C. For ex vivo analyses, homogenized solutions of the tumors were prepared using Solvable (PerkinElmer, Waltham, Massachusetts). The fluorescence of the homogenized sample was measured by a spectrofluorometer (FluoroMax-3; Jobin Yvon, Inc.) with an excitation wavelength of 405 nm and an emission range from 630 to 750 nm with an emission maximum at 667 nm. The photosensitizer concentration in the tissue was calculated based on the change in fluorescence resulting from the addition of a known amount of BPD to each sample after its initial reading. The in vivo measurements were correlated to ex vivo data using a linear fit to examine their agreement based on the goodness of the fit (R2) (see Fig. 2). The ex vivo measurements of BPD concentration were compared to those obtained in vivo using the contact probe method to evaluate for the accuracy of the in vivo acquired BPD concentrations. The linear fit to 3 h DLI results (shown as a solid line) shows close agreement between the in vivo and ex vivo BPD concentrations, y=1.041x, with a fitting goodness of R2=0.9858; similarly, linear fit to 15 min DLI results (shown as a dash-dot line) shows reasonable agreement between the in vivo and ex vivo BPD concentrations, y=0.89x, with a fitting goodness of R2=0.9847. The dashed line represents the line for y=x, if the two measurements were completely in agreement.

Fig. 2

Comparison between in vivo and ex vivo measured BPD concentrations. Solid line is a linear fit for ex vivo versus in-vivo measured BPD concentration at 3 h. Dash-dot line is a linear fit to 15 min results. Ex vivo measurements were made at 3 h only since BPD concentration at 15 min is not detectable ex vivo. Dashed line is for y=x.



Tumor Regrowth Rate Analysis

Tumor volumes were tracked daily, for 14 days, after PDT. Width (a) and length (b) were measured with slide calibers, and tumor volumes (V) were calculated using V=π×a2×b/6.23 Tumor regrowth factor (k) was calculated by the best exponential fit [with a form f(d)=Aekd] to the measured volumes over the days (d). CI was calculated for each treatment group as

Eq. (1)

where k is the tumor regrowth factor for each group and kctr is the regrowth factor for the control group, which received no injection of BPD and light illumination.


Reactive Oxygen Species Explicit Dosimetry

PDT process can be described by a set of kinetic equations that can be simplified to describe the creation of [ROS]rx.8,24,25 These equations are dependent on the temporal and spatial distribution of ϕ, photosensitizer concentration ([S0]), ground state oxygen concentration ([O23]), oxygen supply rate (g), and the photosensitizer-specific reaction-rate parameters (δ, β, σ, and ξ). The relevant equations are

Eq. (2)


Eq. (3)


Eq. (4)


The details of the five parameters involved in the kinetic equations can be found elsewhere (see Table 1).18,19 ξ is the photochemical oxygen consumption rate per light fluence rate and photosensitizer concentration under ample O23 supply. σ is the probability ratio of an ROS molecule to react with a ground state photosensitizer compared to the ROS molecule reacting with a cellular target. β represents the ratio of the monomolecular decay rate of the triplet state photosensitizer to the bimolecular rate of the triplet photosensitizer quenching by O23. δ is the low-concentration correction factor, and g is the maximum macroscopic oxygen perfusion rate. [ROS]rx was calculated by integrating the term of the right-hand side of Eq. (4) over the time course of PDT treatment using the measured ϕ, [S0], and [O23]:

Eq. (5)

[ROS]rx is determined at 3-mm depth using the calculated ϕ and measured [O23] concentration at 3-mm depth to ensure that its minimum value covers the maximum extent of RIF tumors used in this study. [S0] is assumed to be uniform throughout the tumor. If one uses 1 or 2 mm instead of 3 mm, the value of [ROS]rx will increase, thus the resulting threshold [ROS]rx value but the general curve shape [Fig. 6(c)] will not change.



BPD-mediated PDT with different in-air fluences and different ϕair, and different exposure times was performed in mouse models bearing RIF tumors. Light fluence rate, photosensitizer concentration, and tissue oxygenation were measured to calculate PDT dose and [ROS]rx. Table 2 summarizes all treatment conditions as well as the measured and calculated quantities using the photochemical parameters summarized in Table 1.

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, 12  mm3. Figure 3 shows the normalized tumor volume versus time (in days) for the 10 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 #4 with a total fluence of 30  J/cm2 and ϕ=75  mW/cm2, one out of three (33.3%) mice showed a complete response (no tumor regrowth at 14 days post-PDT). This is reflected in Table 2 under the column labeled CI variation count.

Fig. 3

Tumor volumes over days after V=V0ekd PDT treatment. Solid lines are the exponential fit to the data with a functional form of ekd, where d is days ater PDT treatment. The resulting tumor regrowth rates k and its uncentrainty δk are listed in Table 2. The legend for each group lists: in-air fluence rate, in-air fluence, treatment time (in s), BPD concentration (in μM), tissue oxygen concentration (in μM), tumor regrowth rate, and R2 of the fitting the the exponential equation.


Figure 4 shows the temporal dependence of photosensitizer uptake.

Fig. 4

The temporal changes of BPD concentration versus time at 3-mm depth for the treatment conditions. The lines represent linear fits to BPD concentrations during treatment. The average PS concentration [PS] is given in the figure legend for each condition. The uncertainty of [PS] is listed in Table 2.


Measured [O23] was used to refine the photochemical parameters previously determined for the reactive oxygen species explicit dosimetry model used to calculate [ROS]rx. Measured data of [O23] are shown with symbols in Fig. 5. There is no significant change in [O23] during the treatment.

Fig. 5

The temporal dependence of [O23] concentration for different treatment conditions. The concentration of [O23] and its uncertainty are listed in Table 2.


Fluence, PDT dose, and calculated [ROS]rx at 3 mm were compared as dosimetric quantities to correlate with the treatment outcome of BPD-mediated PDT for RIF tumors on a mouse model. The outcome was evaluated by the calculation of CI. No tumor regrowth up to 14 days after treatment resulted in a CI of 1. PDT dose is calculated using the product of PS uptake and measured light fluence rate at 3 mm. We used Eq. (4) and the photophysiological parameters shown in Table 2 to calculate [ROS]rx. The goodness of the fit and the corresponding upper and lower bounds of the fit with 95% confidence interval (gray area) to the fluence, PDT dose, and the calculated [ROS]rx are presented in Fig. 6. Figure 6(a) shows that, while fluence correlates sigmoidal with the PDT outcome, it exhibits large uncertainties as defined by the large bounds of the gray area, as well as by the low value of R2=0.6616. As shown in Fig. 6(b), PDT dose allows for reduced subject variation and improved predictive efficacy as compared to fluence alone. PDT dose showed a better correlation with CI with a higher value of R2=0.9331 and a narrower band of gray area as it accounts for both light dose and tissue [BPD] levels. However, PDT dose overestimates [ROS]rx in the presence of hypoxia as it does not account for the oxygen dependence of ROS quantum yield. The goodness of fit R2=0.9911 and the narrowest gray area in Fig. 6(c) shows that the measured [ROS]rx correlates the best with CI. [ROS]rx accounts for the key quantities of light fluence, photosensitizer concentration, and tissue oxygen level, respectively.

Fig. 6

CI plotted against (a) fluence at a 3-mm tumor depth, (b) calculated PDT dose at 3-mm depth, and (c) mean reacted oxygen species at 3-mm depth ([ROS]rx) calculated using Eqs. (3)–(5) and the parameters summarized in Table 1. The solid lines show the best-fit to the data with functional forms CI=1/(1+35.99e0.09476x), 1/(1+1647e0.9867x), and 1/(1+152.7e19.12x) with R2=0.6616, 0.9331, and 0.9918 for (a), (b), and (c), respectively. The gray region indicates the upper and lower bounds of the fit with 95% confidence level. The gray curves are obtained from an MC simulation of the sigmoid model (see text for details).




As shown in Fig. 2, the accuracy of the in vivo BPD concentration determined by our fluorescence spectroscopy is validated by comparing it with ex vivo measurements. In vivo fluorescence measurements were taken at 15 min and 3 h after BPD administration on the same tumors, and ex vivo fluorescence measurements were performed only once, at 3-h incubation time. The agreement between in vivo and ex vivo BPD concentration at 3 h is within 2%. Comparing the 15-min DLI in vivo BPD concentration with the 3-h DLI ex vivo BPD concentration, the in vivo measurements were about 11% lower, indicating a lower tissue uptake at shorter DLI.

Compared to control mice, all treated mice with total fluences larger than 30  J/cm2 had significant control of the tumor regrowth after PDT (see Fig. 3). The regrowth rate decreased when in-air light fluence increased (e.g., comparing the group of mice treated to 166.7  J/cm2 and that of 30  J/cm2). However, in-air light fluence alone is not a good predictor of the tumor regrowth because of significant variations of either BPD in vivo concentration or oxygen concentration ([O23]). As a result, we have split each (in-air fluence rate and in-air fluence) group further into subgroups depending on either BPD concentration or [O23], whenever there are a significant difference.

Based on a previous study,26 an empirical six-parameter fitting equation was used to fit the Monte-Carlo (MC) simulated data for a 1-cm diameter field, with μa=0.69  cm1, and μs=11  cm1.19 The equation is of the following form:19

Eq. (6)

where the parameters λ1 (16.23  cm1), λ2 (5.58  cm1), λ3 (9.72  cm1), b (0.37), C2 (5.08), and C3 (0.03) are functions of μa(0.69  cm1) and μs (11  cm1) and details of each can be found elsewhere.26 INV=[SSD/(SSD+d)]2, where the source-to-surface distance (SSD)=9.34  cm based on the measurement of light fluence rate in water for the same collimated beam as a function of depth.

Figure 4 shows no variation of PS concentration during PDT for vascular BPD-mediated PDT. This is understandable because most of the BPD drug is located in the vasculature and PDT is obviously insufficient to deplete the BPD concentration [S0]. Similarly, Fig. 5 shows that insignificant variation of [O23] concentration is observed for each group. Comparing the measured [S0] and [O23] with Eqs. (2) and (3), respectively, we come to the conclusion that d[S0]/dt=d[O23]/dt=0 for the vascular PDT condition. This is completely understandable because of the ample oxygen supply in the blood vessel and presumably ample BPD concentration in the vasculature.

Fluence, PDT dose at 3 mm, and [ROS]rx at 3 mm were compared as dosimetric quantities to estimate the outcome of BPD-mediated vascular PDT for RIF tumors on a mouse model. Outcome was evaluated by the calculation of CI. No tumor regrowth up to 14 days after treatment resulted in a CI of 1. The goodness of the fit and the corresponding upper and lower bounds of the fits (gray area) to the fluence, PDT dose, and mean [ROS]rx are presented in Fig. 6. MC simulation was used to produce the gray area due to the uncertainty in the parameter a, b for the sigmoid curve y=1/[1+a×exp(bx)]. The means and standard deviations of the simulation parameters a and b were obtained using the global optimization toolbox of MATLAB® (cftool.m). In the MC simulation, we selected 1000 parameter pairs (a,b) within the standard deviation (δa,δb) using a random number generator with normal distributions, the resulting calculated y=1/[1+a×exp(bx)] for all (a,b) for a particular dosimetric metric (fluence, PDT dose, or [ROS]rx) was used to generate a cumulative probability distribution for each of the y values between [0, 1]. We then found corresponding x values for the 2.5% tiles and 97.5% tiles of the cumulative probability distribution and they form the two bounds of the gray zone, the left and right bounds were joined to form the uncertainty (gray) areas. Figure 6(a) shows that, while fluence correlates following a sigmoid curve with PDT outcome, it exhibits large uncertainties as defined by the large bounds of the gray area as well as by the low value of R2=0.66. As shown in Fig. 6(b), PDT dose allows for reduced subject variation and improved predictive efficacy as compared to fluence. PDT dose showed a better correlation with CI with a higher value of R2=0.93 and a narrower band of gray area as it accounts for both light dose and tissue [BPD] levels. However, PDT dose overestimates [ROS]rx in the presence of hypoxia as it does not account for the oxygen dependence of ROS (mostly O12) quantum yield. The goodness of fit R2=0.99 and the narrowest gray area in Fig. 6(c) shows that the mean [ROS]rx correlates the best with CI among the three.

Based on the findings of this study, PDT dose and [ROS]rx exhibit threshold dose behavior as they can be fitted by a sigmoid function {S(x)=1/(1+e[(xx0)/w0]}, where x0=7.5  μMJ/cm2 with uncertainty w0=1.0  μMJ/cm2 and x0=0.26  mM with uncertainty w0=0.05 for PDT dose and [ROS]rx, respectively. For PDT dose, x0 can be converted to the absorbed dose by BPD by multiplying the extinction coefficient (ε=0.0783  μM1cm1), resulting in 0.59  J/cm3, which corresponds to (2.0±0.2)×1018  photons/cm3 (by dividing the energy per photon hc/λ=2.88×1019  J for λ=690  nm). The mean PDT dose threshold for BPD at DLI 15 min (7.5  μMJ/cm2) is 7.7 times lower than those reported for BPD at DLI 3 h (58  μMJ/cm2).19 The mean [ROS]rx threshold concentration of x0=0.26±0.05  mM for DLI of 15 min is 3.8 times lower than to the published result for BPD (0.98 mM) for DLI of 3 h.19 The decrease for PDT dose can be explained by a decrease of BPD uptake between the tissue and the vasculature, i.e., our result implies that BPD concentration in the vessel is 7.7 times higher than those in the tissue. This is consistent with the published literature where a simulation using a diffusing model showed much lower concentration of BPD in the vessel than that in the surrounding tissue.27 The decrease for [ROS]rx threshold dose for vascular BPD PDT is due to the fact that both BPD concentration and [O23] in the vessel are substantially higher than those in the surrounding tissue. However, to estimate the resulting [O23] difference between the vessel and the surrounding tissue requires an estimate of the light fluence inside the vessel versus the light fluence rate inside tissue. Some literature reports the former to be lower than the latter by 20%.28



The response of mouse RIF tumors to PDT depends on the tissue oxygenation, photosensitizer uptake, total energy delivered, and the ϕ, in which the treatment is delivered. An accurate dosimetry quantity for the evaluation of the treatment outcome should account for all of these parameters. This study evaluated the efficacy and outcomes of different PDT treatments and how fluence, PDT dose, and [ROS]rx compare as dosimetric quantities. The correlation between CI and [ROS]rx suggests that [ROS]rx at 3 mm is the best quantity to predict the treatment outcome for a clinically relevant tumor regrowth endpoint. PDT dose is a better dosimetry quantity when compared to fluence but is worse than [ROS]rx as it does not account for the consumption of [O23] for different ϕ. For BPD in RIF tumors, our measurements show constant temporal dependence of in vivo oxygen concentration during PDT, which cannot be well modeled by our macroscopic model; thus it is necessary to make [O23] measurements during PDT to obtain [ROS]rx. We find the threshold value of [ROS]rx for BPD-mediated vascular PDT at DLI 15 min to be 3.8 times smaller than the corresponding value for BPD-mediated PDT at DLI 3 h. This is being reported for the first time.


There is no conflict of interest to declare for any authors.


The authors would like to thank Dr. Theresa M. Busch, Min Yuan, Joann Miller, and Shirron Carter for their advice concerning the mouse studies and protocols. This work was supported by grants from the U.S. National Institutes of Health (NIH) (Nos. R44 CA183236 and P01 CA87971).



R. Penjweini et al., “Evaluation of the 2-(1-hexyloxyethyl)-2-devinyl pyropheophorbide (HPPH) mediated photodynamic therapy by macroscopic singlet oxygen modeling,” J. Biophotonics, 9 1344 –1354 (2016). Google Scholar


H. Qiu et al., “Macroscopic singlet oxygen modeling for dosimetry of Photofrin-mediated photodynamic therapy: an in-vivo study,” J. Biomed. Opt., 21 (8), 088002 (2016). JBOPFO 1083-3668 Google Scholar


P. Agostinis et al., “Photodynamic therapy of cancer: an update,” CA Cancer J. Clin., 61 (4), 250 –281 (2011). CAMCAM 0007-9235 Google Scholar


A. P. Castano, T. N. Demidova and M. R. Hamblin, “Mechanisms in photodynamic therapy: part one—photosensitizers, photochemistry and cellular localization,” Photodiagn. Photodyn. Ther., 1 (4), 279 –293 (2004). Google Scholar


K. K. Wang et al., “Explicit dosimetry for photodynamic therapy: macroscopic singlet oxygen modeling,” J. Biophotonics, 3 (5–6), 304 –318 (2010). Google Scholar


X. Liang, K. K. Wang and T. C. Zhu, “Singlet oxygen dosimetry modeling for photodynamic therapy,” Proc. SPIE, 8210 82100T (2012). Google Scholar


D. D. McMillan et al., “Parameter determination for singlet oxygen modeling of BPD-mediated PDT,” Proc. SPIE, 8568 856810 (2013). Google Scholar


T. C. Zhu et al., “Macroscopic modeling of the singlet oxygen production during PDT,” Proc. SPIE, 6427 642708 (2007). Google Scholar


R. F. Spaide, J. Sorenson and L. Maranan, “Photodynamic therapy of subfoveal choroidal neovascularization in age-related macular degeneration with Verteporfin: one-year results of 2 randomized clinical trials—TAP report 1,” Arch. Ophthalmol., 117 (10), 1329 –1345 (1999). AROPAW 0003-9950 Google Scholar


B. Chen et al., “Tumor vascular permeabilization by vascular-targeting photosensitization: effects, mechanism, and therapeutic implications,” Clin. Cancer Res., 12 (3), 917 –923 (2006). Google Scholar


T. Sheng et al., “Reactive oxygen species explicit dosimetry to predict tumor growth for Photofrin-mediated photodynamic therapy,” Proc. SPIE, 10860 108600V (2019). Google Scholar


M. M. Kim et al., “Explicit macroscopic singlet oxygen modeling for benzoporphyrin derivative monoacid ring A (BPD)-mediated photodynamic therapy,” J. Photochem. Photobiol. B., 164 314 –322 (2016). JPPBEG 1011-1344 Google Scholar


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). Google Scholar


M. M. Kim et al., “Parameter determination for BPD-mediated vascular PDT,” Proc. SPIE, 8931 89311D (2014). Google Scholar


T. C. Zhu et al., “In-vivo singlet oxygen threshold doses for PDT,” Photonics Lasers Med., 4 (1), 59 –71 (2015). Google Scholar


B. Liu et al., “Comparison of PDT parameters for RIF and H460 tumor models during HPPH-mediated PDT,” Proc. SPIE, 8931 89311C (2014). Google Scholar


T. C. Zhu et al., “Comparison of singlet oxygen threshold dose for PDT,” Proc. SPIE, 8931 89310I (2014). Google Scholar


M. M. Kim et al., “On the in-vivo photocheical rate parameters for PDT reactive oxygen species modeling,” Phys. Med. Biol., 62 (5), R1 –R48 (2017). PHMBA7 0031-9155 Google Scholar


M. M. Kim, R. Penjweini and T. C. Zhu, “Evaluation of singlet oxygen explicit dosimetry for predicting treatment outcomes of benzoporphyrin derivative monoacid ring A-mediated photodynamic therapy,” J. Biomed. Opt., 22 (2), 028002 (2017). JBOPFO 1083-3668 Google Scholar


R. Penjweini et al., “Explicit dosimetry for 2-(1-hexyloxyethyl)-2-devinyl pyropheophorbide-a (HPPH) mediated photodynamic therapy: macroscopic singlet oxygen modeling,” J. Biomed. Opt., 20 (12), 128003 (2015). JBOPFO 1083-3668 Google Scholar


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). JBOPFO 1083-3668 Google Scholar


H. Qiu et al., “Dosimetry study of PHOTOFRIN-mediated photodynamic therapy in a mouse tumor model,” Proc. SPIE, 9694 96940T (2016). Google Scholar


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). PPSHCB 1474-905X Google Scholar


X.-H. Hu et al., “Modeling of a type II Photofrin-mediated photodynamic therapy process in a heterogeneous tissue phantom,” Photochem. Photobiol., 81 1460 –1468 (2005). PHCBAP 0031-8655 Google Scholar


M. M. Kim et al., “Singlet oxygen explicit dosimetry to predict long-term local tumor control for BPD-mediated photodynamic therapy,” Proc. SPIE, 10047 100470X (2018). Google Scholar


Y. H. Ong and T. C. Zhu, “An analytic function for predicting light fluence rate of circular fields on a semi-infinite turbid medium,” Opt. Express, 24 (23), 26261 –26281 (2016). OPEXFF 1094-4087 Google Scholar


X. Zhou et al., “Analysis of effective molecular diffusion rates for Verteporfin in subcutaneous versus orthotopic duning prostate tumors,” Photochem. Photobiol., 79 323 –331 (2004). PHCBAP 0031-8655 Google Scholar


E. P. M. LaRochelle et al., “Modeling PpIX effective light fluence at depths into the skin for PDT dose comparison,” Photodiagn. Photodyn. Ther., 25 425 –435 (2019). Google Scholar


Tianqi Sheng received her PhD in optics from Jilin University, China, in 2013 and she was a postdoctoral researcher at the Department of Radiation Oncology, University of Pennsylvania during this study. Her research interest includes preclinical photodynamic therapy dosimetry.

Yihong Ong received his PhD in biomedical engineering in 2015 from Nanyang Technological University, Singapore. Currently, he is a postdoctoral researcher in the Department of Radiation Oncology at the University of Pennsylvania. His current research interests include development of reactive oxygen species explicit dosimetry (ROSED) and diffuse optics for blood flow measurement during PDT.

Wensheng Guo received his PhD in 1998 in Biostatistics from University of Michigan. Currently, he is a professor in the Department of Biostatistics and Epidemiology at the University of Pennsylvania. His current research interest includes spline smoothing, time series analysis, functional models, and longitudinal data analysis.

Timothy C. Zhu received his PhD 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 explicit PDT dosimetry, ROSED, integrated system for interstitial and intracavitory PDT, diffuse optical tomography, in vivo dosimetry, and external beam radiation transport.

© The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Tianqi Sheng, Yi Hong Ong, Wensheng Guo, and Timothy C. Zhu "Reactive oxygen species explicit dosimetry to predict tumor growth for benzoporphyrin derivative-mediated vascular photodynamic therapy," Journal of Biomedical Optics 25(6), 063805 (7 January 2020).
Received: 18 September 2019; Accepted: 3 December 2019; Published: 7 January 2020

Back to Top