Photodynamic therapy (PDT) is an emerging cancer treatment that uses photosensitizers activated with visible wavelength light. The photoactivated photosensitizer (PS) excites oxygen to the state via the type-II photochemical pathway.1 Previous studies have provided strong evidence that singlet oxygen is the active species in cancer cell or endothelial cell necrosis.2, 3, 4, 5, 6 Thus, a device that facilitates the measurement of the singlet oxygen in vivo could provide a crucial parameter in PDT dosimetry and the potential of improved and even individualized therapeutic design.
There has been considerable interest in developing a sensor for singlet oxygen that could be used as a real-time dosimeter during PDT treatments. 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 Correlations of the singlet oxygen produced with treatment efficacy could be one important use of such a sensor. Some researchers have attempted to develop dosimeters based on the prompt, radiatively allowed fluorescence of the PS in the tumor, but simple PS fluorescence detection may not accurately reflect the complex and dynamic interactions that lead to the therapeutic effect.3 Pogue 11 used small electrodes to measure total oxygen content in tumors during PDT in animal studies and demonstrated deoxygenation during treatments. Other researchers have shown that oxygenation of tissue can enhance PDT efficiency.12, 13, 14
Some of the earlier work used sensitive but slow (bandwidths on the order of ) detectors for the singlet oxygen emission. While these proved useful for in vitro studies, this approach was unsuccessful in vivo due to the slow detector response. Recently, a high quantum efficiency, high speed , low dark current photomultiplier tube (PMT) became commercially available, and this greatly enhanced detection of the weak singlet oxygen emission near . Using this detection method, we have developed a sensitive, diode laser-based monitor for singlet oxygen produced by PDT.15, 16, 17 In Fig. 1 , we illustrate our detection strategy that includes both temporal discrimination and spectral filtering to maximize the sensitivity of our system to singlet oxygen emission. The pioneering work of Wilson, Patterson, and Lilge8 has also described a singlet oxygen monitor based on similar detection methods. However, there are significant differences between the two approaches. Wilson’s group used a Q-switched, frequency-doubled, Nd:YAG laser (pulse lengths ) in contrast to our relatively long pulses. Diode laser pulses do not produce significant energy compression, and the peak power of the diode laser is much less than that of Q-switched devices. In addition, our system is fiber coupled; the diode laser beam delivery and the singlet oxygen collection are both facilitated by fiber optic cables, where Wilson, Patterson, and Lilge used free-space optics. More recently, Maisch 18 used a Q-switched Nd:YAG laser and the PMT approach to detect singlet oxygen produced by PDT in bacteria. Yamamoto 19 have used a Q-switched Nd:YAG laser-pumped dye laser to produce and monitor singlet oxygen both in vitro and in vivo.
Previously, we described an early version of our singlet oxygen monitor and demonstrated that we could measure the production of singlet oxygen in a variety of singlet oxygen quenching media including water, methanol, acetone, and protein-rich aqueous solutions.15, 16, 17 Here we report the detection of singlet oxygen on tumor-laden rats with the photosensitizers chlorin e6 and 5-aminolevulinic acid-induced protoporphyrin (ALA-induced Pp IX). We observed a positive correlation between the observed singlet oxygen signal and the post-treatment regression of tumors on several rats.
Apparatus and Measurement of Singlet Oxygen
We used time-resolved emission measurements and optical filtering to enhance the sensitivity for detecting the singlet oxygen, as shown in Fig. 1. The temporal filtering method is outlined in Fig. 1a, and the spectral filtering to further isolate the singlet oxygen emission (near ) from long wavelength PS fluorescence and/or phosphorescence is shown in Fig. 1b. It is well established that singlet oxygen is severely quenched in aqueous media. Indeed, the lifetime is reduced from about in the gas phase to approximately in aqueous media. In biological media, however, the lifetime is further reduced to around .13, 14 This is one of the major challenges of detecting singlet oxygen both in vitro and in vivo. The high bandwidth, near infrared (IR) PMT and a temporally gated detection system help isolate the longer lived singlet oxygen emission from the much faster PS fluorescence. Two diode lasers were used in the singlet oxygen sensor: 1. , pulse width, repetition rate for , and 2. , pulse width, repetition rate for ALA.
The singlet oxygen detection system includes: 1. the diode laser module; 2. optical filters/PMT detection system; and 3. the data acquisition system with a photon counting board. Both the pulsed diode laser light delivery and singlet oxygen emission collection are fiber optic-coupled via a bifurcated probe, and this facilitates application of the sensor to in vivo studies. The excitation light was delivered to the site being examined by one leg of the fiber optic, and the near-IR emission was collected by the other leg. The detected emission included the singlet oxygen signal and other possible near-IR radiation such as PS fluorescence, phosphorescence, and tissue autofluorescence. At the distal end of the detection fiber, the near-IR radiation was collimated and sent through a series of three narrow bandpass filters ( bandwidths and transmission center wavelengths of 1.22, 1.27, and ) contained in a sliding filter holder configuration. At each bandpass filter wavelength, the signal was accumulated over 30,000 laser pulses, typically with operation. The filters provided spectral discrimination of the singlet oxygen emission from the interferences mentioned before. The emission feature of singlet oxygen is centered in the wavelength region of . The two filters centered at 1.22 and provided measurements of the background emission that did not contain singlet oxygen emission. The optical radiation transmitted by the filters was detected by a thermoelectrically cooled PMT (Hamamatsu model H9170-45, NJ, USA) run in a photon counting mode. The output current pulses from the PMT were amplified with a high bandwidth amplifier and fed to the data acquisition system. We combined the acquisition system with a fast photon counting board (Becker and Hickl model MSA-300, MA, USA) that can handle up to a operation rate for our system configuration.
For our initial in vivo measurements with , we used only one filter centered at a wavelength of . The background signal subtraction was done with a signal measured at an off-tumor site distinct from the tumor site. The off-tumor site retains much less PS, so the detected signal is predominantly due to scattered light from the irradiated tissue. However, we found large variations of background signal due to differing skin optical characteristics at different locations within the same animal. With only one optical filter at , these variations severely affected the background subtraction method used to extract the singlet oxygen emission signal, especially for weaker singlet oxygen signals. Thus, to minimize signal variations due to skin optical characteristics and to enhance the detection sensitivity, we introduced the three filter method. This allowed measurement of both the singlet oxygen and background emissions over the same area of skin.
Our system produces the detected singlet oxygen emission by the type-II process. The pulsed diode laser in our sensor is distinct from the cw PDT treatment laser. The pulsed diode laser produces singlet oxygen from ground state oxygen within the tissue that is available at various times before, during, and after the treatments by the cw PDT laser. The singlet oxygen that we measure is proportional to the product of the available ground state oxygen concentration multiplied by the PS concentration (i.e., ). In our system, this measurement is done at various times before and during the PDT treatments. The cw treatment laser is shut off for a short period , and the pulsed diode laser is used to produce singlet oxygen. Thus, our monitor measures the capacity of the irradiated tissue to produce singlet oxygen, and the measured signal is assumed to be proportional to the singlet oxygen produced when the PDT treatment laser is on. In future work, we plan to pulse the treatment laser to provide singlet oxygen emission signals directly during PDT treatment.
Protocol for Animal Study
We measured singlet oxygen production in tumor-laden rats during PDT treatment using the singlet oxygen monitor at the Wellman Laboratories of Photomedicine at Massachusetts General Hospital (MGH). The animal model used the R3327-MatLyLu prostate cancer cell line that somewhat follows the human disease pattern and metastasizes into lymph nodes and lungs. Since the R3327-MatLyLu cells are syngeneic, there is no need for immune suppression in the rats, which enhances the model’s clinical relevance. This study was performed under MGH protocol and was reviewed and approved by the MGH Subcommittee on Research Animal Care (SRAC)-OLW Assurance (title of protocol: Measurement of Singlet Oxygen Formation During Photodynamic Therapy Using a Diode Laser-Based System: Rats Model).
Tumors were induced by subcutaneous injection of a suspension of R3327-MatLyLu cells into the flanks of two-month-old male Copenhagen rats of 150- to weight. Once the tumors had reached treatment size, the animals received an injection of either intravenous or intraperitoneal ALA. After a waiting period of that allows for accumulation of PS in the tumor, light exposure was done in the anesthetized animals through the shaved skin.20, 21 Treatment irradiations of the tumors were performed using cw diode laser sources (HPD, Incorporated) with a wavelength that matched the absorption profile of the PS; diode laser for , and diode laser for ALA. The incident power density of the treatment laser was .
The singlet oxygen probe was placed above the skin of the subject animals. The singlet oxygen emissions were recorded at several times: 1. before infusion of a PS, 2. immediately before starting the therapeutic light exposure but after the PS had spread systemically for a specific PS incubation time period, 3. at several separate times during the therapeutic irradiation, and 4. at the end of the light exposure. Each measurement contained 30,000 pulses from the diode laser source. The measurement required interrupting the therapeutic light exposure, but only for approximately . Using the singlet oxygen monitor in this way does not influence PDT treatment protocols.
We recorded tumor growth after PDT treatment with a control group and different ALA dosage groups: 100 and body weight of ALA, and 50 and of total light dosage. For the tumor regression study, implant sizes were repeatedly assessed by caliper measurement. Tumor volumes were calculated using the formula for the ellipsoid . Tumors were observed while growing to treatment dimensions of about , which requires approximately . The animals were sacrificed at from the day of PDT treatment.
Results and Discussion
Kinetic Model for Singlet Oxygen
We interpret the observed signals of the singlet oxygen emission using a kinetic model described earlier17 and summarized next. Assuming that the singlet oxygen is produced via the type-II mechanism and treating the excitation pulse as an “instantaneous” source, leads to a time-dependent solution for the singlet oxygen concentration in after the termination of the excitation pulse, given by:is the number of photons per in the excitation pulse incident on the sample, is the absorption cross section of the PS in , is the concentration of PS ground state, is the quantum yield of singlet oxygen, and and are lifetimes of PS triplet state and singlet oxygen, respectively.
This is identical to the result discussed in Refs. 8, 9, 10, where they used a Q-switched, frequency-doubled, Nd:YAG laser or an optical parametric oscillator as the excitation source. Both produce very short pulses , and Eq. 1 is an adequate description. However, our diode laser approach uses much longer pulses , and the instantaneous excitation model is no longer valid while the diode laser is on. In our case, the diode laser intensity is constant over the duration of the pulse, and the kinetics of singlet oxygen production and destruction during the laser pulse become non-negligible. To gain a better understanding of this long pulse limit, we solved this kinetic model numerically using MathCAD. An example of this model study is shown in Fig. 2 . During the diode laser pulse, the population of the PS singlet state quickly reaches a steady state value [Fig. 2a]. The PS triplet state population grows and populates the oxygen singlet state via energy transfer collisions during the laser pulse [Fig. 2b]. The total temporal profile [Fig. 2c] is the sum of the emissions from both the PS and the singlet oxygen. We readily observed these processes in several in vitro studies, as shown in Fig. 2d.17 Some data observed in the present in vivo studies also displayed both the growth of the singlet oxygen signal in the early part of the diode laser pulse and its decay subsequent to the laser pulse.
Figures 2e and 2f show comparisons of this model to our data from aqueous solutions with two different excitation pulse widths: 1 and . (Note that we collect data for prior to the initiation of the diode laser pulse to provide a “zero” level for the subsequent data reduction.) Many photosensitizers produce some prompt PS singlet state emission even at the singlet oxygen emission wavelength of . This PS fluorescence decays in a few nanoseconds after the diode laser pulse is terminated, as shown in both our model and experimental measurements. Using PDT treatment parameters such as PS concentration, laser intensity, PS triplet lifetime, singlet oxygen lifetime, and energy transfer rates, the model provides a guide for the interpretation of the singlet oxygen production in the solution phase. The model also implies that there may be an optimum diode laser pulse width to maximize the singlet oxygen produced per pulse, and this may enable one to develop optimized conditions of PDT treatments.
Singlet Oxygen Detection in Tumor-Laden Rats
Our first experiments to investigate the detection sensitivity for singlet oxygen in vivo used and only the bandpass filter for spectral discrimination of the singlet oxygen emission. We recorded the temporal evolutions of the emission at and integrated the emission signal counts over the time interval of after the termination of the diode laser pulse. The singlet oxygen and background signals were measured at two distinct locations on each rat. The background emission signals were recorded at a location that did not include tumor tissue. The singlet oxygen emission plus any background signals were measured on the tumor site. The signal of singlet oxygen produced from the tumor site was calculated by subtracting the background emission (off-tumor site) from the emission observed on the tumor site. This measurement (in photoelectron counts) provided a relative value for the singlet oxygen concentration produced by the pulsed diode laser. For the in vivo studies, we used two diode laser sources: a continuous wave (cw) treatment laser and a pulsed laser for the singlet oxygen measurements. We controlled the total PDT treatment light fluence of the cw laser, and made the singlet oxygen measurements during short periods of time when the cw laser was off.
Figures 3a, 3b, 3c show typical temporal profiles (background subtracted) in photoelectron counts at immediately after a light treatment with a diode laser for three tumor-laden rats. One had no , one was infused to a concentration of body weight, and the third was infused to a concentration of body weight. The data indicate a correlation between the temporal profiles of singlet oxygen emission and the initial concentrations of the infused . Figure 3d summarizes the observed singlet oxygen production with different PDT light irradiations for three tumor-laden rats with . These data represent the summation of photoelectron counts (after background subtraction) at . For these data, the singlet oxygen was measured at two locations on the tumor in each rat, and the data show the two-site average signal of the singlet oxygen production as a function of total PDT treatment light fluence. We observed photoelectrons that we interpret as being due to singlet oxygen emission in the control group that had no . This is most likely due to singlet oxygen production from naturally existing porphyrins in the skin, as reported by other groups.10 After of PDT treatment with the cw treatment laser, there was a distinct increase in the singlet oxygen signal. Subsequent to further treatment, the data show a reduction of the singlet oxygen signal as the total light dose was increased. The causes of this diminution are likely a combination of tumor hypoxia and photobleaching.
While the background subtraction method was adequate for the photosensitizer that produced relatively large singlet oxygen emissions, other photosensitizers, e.g., ALA, produce considerably smaller singlet oxygen signals as we observed in previous in vitro studies.15, 16, 17 Consequently, to enhance the sensitivity of the singlet oxygen monitor, we added two optical filters that pass near-IR radiation at 1.22 and (shorter and longer wavelengths than the singlet oxygen emission wavelength). All three filters were contained in a sliding filter holder to allow rapid and reproducible insertion of each filter into the optical detection path. The observed signals with the bandpass filters centered at 1.22 and were used to improve the background subtraction routine and provide more accurate determination of the singlet oxygen signal. The average value of the signals at 1.22 and was subtracted from the signal recorded with the bandpass filter to obtain the singlet oxygen signal. These measurements were made sequentially but at the same location on the rats. The background signals may contain photosensitizer fluorescence, phosphorescence, and autofluorescence of any of the optical components, including the collection fiber optics.
In the ALA study, we used this three-filter detection method and a diode laser for the PDT treatment light source. The singlet oxygen production was measured with a pulsed diode laser 1. before ALA injection, 2. just before PDT treatment after of ALA incubation time period, 3. between the PDT treatment, and 4. at the end of the PDT treatment. The tumor size was recorded for following PDT treatment.
Figures 4a and 4b show fits to data for a single animal but are typical for the ALA study. Figure 4a shows the averages of 30,000 diode laser pulses prior to injection with ALA, each at the three wavelengths transmitted by the narrow bandpass filters. These data demonstrate that optical signals are observed at all three wavelengths, probably due to broadband PS fluorescence and phosphorescence. Since we extract the singlet oxygen signal using the difference between the signal at and the average of the two signals at 1.22 and , Fig. 4a implies essentially no detectable singlet oxygen signal. In contrast, Fig. 4b shows similar data, but after injection with ALA and after 50J of PDT treatment. The signal at is clearly larger than at the other two, off-band wavelengths and persists for a longer time.
During the PDT treatments, the optical fiber probe of the singlet oxygen monitor was moved to allow unimpeded access of the tumor by the cw treatment laser beam that covered the entire area of the tumor site ( diam). Subsequent to the delivery of a selected treatment fluence, the cw laser was shut off, and the fiber probe of the singlet oxygen monitor was repositioned from the tumor site. The diameter of the monitor diode laser beam on the tumor site was . This removal and repositioning of the monitor beam before and after a treatment caused variations in the observed signal levels, presumably due to nonuniformities in the optical properties of the tissue (including the skin) of the subject animal. However, we measure the difference between the average of the two, out-of-band (with respect to singlet oxygen) emissions at 1.22 and and the in-band signal at to determine the relative singlet oxygen concentration. This variation is clearly shown in Figs. 4a and 4b. In Fig. 4a (no PS), we observed signals at all three wavelengths that were nearly identical and larger than those of Fig. 4b (with PS). However, the signal detected at is larger than the signals at 1.22 and in Fig. 4b. This is consistent with the production of singlet oxygen in the tumor.
The singlet oxygen signal was monitored in each animal, and Figs. 4c and 4d show the singlet oxygen production (the measured counts at minus the counts averaged at 1.22 and ) for two rats as a function of PDT light dosage. The data shown in Fig. 4d indicate a negative value for the singlet oxygen concentration for the no ALA case. This is a result of the background subtraction routine described before and represents the systematic uncertainty for the weak signals observed in the absence of a PS. The singlet oxygen signal was clearly higher after the initial PDT treatment. This may be due to increased oxygenation or blood flow to the treated area. Indeed, several groups have reported enhanced blood flow and increased ground state oxygenation early in the PDT treatment process.13, 14 Cottrell have also recently observed a formation of a photoproduct during ALA photobleaching, and the photoproduct is a PDT active species producing singlet oxygen.22 They reported that this process increased the concentration of the photoproduct until about of fluence had been delivered. Note that some researchers23, 24 have attempted to develop PDT dosimeters based on the fluorescence intensity of the PS in the tumor, but photobleaching of the PS precludes this as an accurate method. The singlet oxygen production depends on all the parameters that affect PDT treatment: [PS], light dosage, and .
Tumor Regression Study
Following the PDT treatment with different ALA dosages and light dosages, the tumor size was measured each day for . In Fig. 5 , we present results of tumor growth and regression following the treatment. Figure 5a illustrates the tumor regression measured for three rats as a function of the singlet oxygen measured during the PDT treatment. One was in a control group that had no PDT treatment light or ALA. One rat produced a relatively small singlet oxygen signal of 155 counts and showed steady tumor growth similar to the control. The rat that produced the most singlet oxygen of 921 counts (approximately six times that of the one with 155 counts singlet oxygen) showed a tumor regression of two days after the treatment. Tumor growth then began again. The treatments were not intended to be curative.
In Fig. 5b, we plot the tumor regression as a function of the singlet oxygen measured during the PDT light treatment in a group of seven rats. This plot shows a positive correlation between singlet oxygen production and tumor volume regression. The error bars represent the estimated statistical uncertainties in measuring tumor size and the standard deviations of the number of singlet oxygen photoelectrons detected. The data indicate larger systematic uncertainties in singlet oxygen production and tumor regression for different animals. Even with the rather large uncertainties and our very limited study, we observed a clear trend of larger tumor volume removal with more singlet oxygen production during PDT light treatment. While these are preliminary data from a limited number of animals, there does appear to be a positive correlation of the volume of tumor regression and the amount of singlet oxygen produced. An extensive animal study is planned to further investigate this relationship with the goal of developing a statistically significant database.
We described a novel, fiber-coupled, diode laser-based system for measuring the production of singlet oxygen during PDT, and successfully applied this device to measure PDT produced singlet oxygen in vivo. Using tumor-laden rats, we observed singlet oxygen produced by PDT treatments with both and ALA. We also observed tumor regression that correlates with the measured singlet oxygen produced. Although these initial results display only moderate signal-to-noise ratios, there is a clear relationship between singlet oxygen production and tumor regression. This could potentially provide the critical parameter in PDT dosimetry and facilitate individual therapeutic design. We plan to enhance the sensitivity of the monitor and will complete a much more extensive animal study to further test the capability of monitoring singlet oxygen in vivo.
This work was supported under a grant from the National Cancer Institute of the National Institutes of Health (grant number 5 R44 CA096243-03). We are grateful for this support. The authors would like to thank B. W. Pogue (Dartmouth College) and B. C. Wilson (University of Toronto, Canada) for helpful discussions.