Conference Presentation Video

Conference Presentation Video

Abstract

Abstract.
We measured the optical properties of freshly excised kidneys with renal parenchymal tumors to assess the feasibility of photodynamic therapy (PDT) in these patients. Kidneys were collected from 16 patients during surgical nephrectomies. Spatially resolved, white light, steady-state diffuse reflectance measurements were performed on normal and neoplastic tissue identified by a pathologist. Reflectance data were fit using a radiative transport model to obtain absorption (μa) and transport scattering coefficients (μs′), which define a characteristic light propagation distance, δ. Monte Carlo (MC) simulations of light propagation from cylindrical diffusing fibers were run using the optical properties extracted from each of the kidneys. Interpretable spectra were obtained from 14 kidneys. Optical properties of human renal cancers exhibit significant inter-lesion heterogeneity. For all diagnoses, however, there is a trend toward increased light penetration at longer wavelengths. For renal cell carcinomas (RCC), mean values of δ increase from 1.28 to 2.78 mm as the PDT treatment wavelength is increased from 630 to 780 nm. MC simulations of light propagation from interstitial optical fibers show that fluence distribution in tumors is significantly improved at 780 versus 630 nm. Our results support the feasibility of PDT in selected renal cancer patients, especially with photosensitizers activated at longer wavelengths.

Figures in this Article

Conference Presentation Video

Introduction

Conference Presentation Video

Renal parenchymal tumors, most frequently renal cell carcinomas (RCC), are the sixth most common malignancy diagnosed in men and the eighth most common in women in the United States. In 2012, it is predicted that there will be more than 64,000 new cases of renal cancer in the United States, with roughly 14,000 deaths occurring from the disease.1 Traditionally, these tumors were detected by the classical triad of gross hematuria, flank pain, and an abdominal mass,2 but currently fewer than 10% of patients have these symptoms or signs, and more than half the tumors are detected as incidental findings on abdominal imaging performed for unrelated symptoms, such as abdominal pain. Most incidentally found renal tumors are fairly small (less than 4 cm), and they generally have an excellent prognosis with minimally invasive procedures. Traditionally, nonmetastatic RCCs have been treated with curative intent by radical nephrectomy, a major surgery with concomitant postoperative complications and frequently accompanied by worsening overall kidney function. This procedure can now frequently be performed laparoscopically using robotic technology, reducing morbidity over an open surgical approach without compromising oncologic outcome in appropriately selected cases.3 Furthermore, there is compelling evidence that, for appropriately selected patients, nephron-sparing surgery, in which a partial rather than a total nephrectomy is performed, is as effective as radical nephrectomy.4 These nephron-sparing treatments, originally performed for obligate reasons (solitary kidney, severe renal impairment, bilateral tumors, etc.), have become increasingly popular and now represent standard therapy for tumors that are anatomically amenable to such an approach. Again, based on technical and anatomical considerations, such therapies can often be done laparoscopically with as much success as open surgery in appropriately selected cases.

Conference Presentation Video

Over the past five years, there has been increasing interest in further reducing morbidity by the use of minimally invasive techniques,5 particularly radiologically driven ones. These have primarily focused on thermotherapies, including radiofrequency ablation (RFA)6 and cryotherapy.7 These techniques, though still considered experimental by some, are becoming increasingly mainstream. In short- and medium-range follow-up, outcomes approach those for extirpative surgeries.8 Patient and lesion selection appears to be directly proportional to outcomes; smaller, more peripheral lesions tend to have better outcomes.

Conference Presentation Video

Nephron-sparing minimally invasive treatments are particularly helpful in patients with hereditary causes of RCC, such as von Hippel-Lindau and Birt-Hogg-Dubé syndromes, as these patients tend to have multiple and recurrent renal tumors.9 However, many of these tumors arise in anatomical locations that would prevent them from being amenable to percutaneous RFA or cryotherapeutic approaches. These lesions include centrally located tumors near the major renal blood vessels, in which a “heat sink” effect significantly reduces the efficacy of thermotherapies, primarily by moderating temperatures. Although there have been various attempts to make tumors more sensitive to those therapies by using pretreatment with targeted antiangiogenic therapies, the efficacy of such approaches remains unknown and certainly may add complexity and expense to existing therapies. Furthermore, roughly 20% to 30% of small renal space occupying lesions are indeterminant cysts (Bozniak 2 s and 3 s), which are not ideal thermotherapeutic approaches and can pose increasing challenges for laparoscopic partial nephrectomy.10 Thus, approximately 40% of newly diagnosed small renal lesions, which represent as many as two-thirds of the newly diagnosed renal tumors, are not always amenable to current minimally invasive nephron-sparing interventions.

Conference Presentation Video

Photodynamic therapy (PDT) using laser excitation of a photosensitizer to generate reactive oxygen species in tumors has been approved by health regulatory agencies worldwide and has been used in a wide variety of cancers.11,12 Though it has not been used in clinical scenarios for renal parenchymal tumors, there is little reason to expect, particularly based upon the highly vascular nature of these tumors, that they would not be amenable—indeed ideal—for PDT approaches combined with state-of-the-art imaging and percutaneous applications.

Conference Presentation Video

The purpose of this study was to obtain optical property information—the absorption and scattering coefficients—from lesions in freshly excised kidneys. This data may then be incorporated into treatment planning models for interstitial PDT to assess the feasibility of PDT to treat renal cancers.

Conference Presentation Video

Research Population Recruitment and Tissue Handling

Conference Presentation Video

Method of subject identification

Conference Presentation Video

Subjects who were patients of attending physicians in the Department of Urology at the University of Rochester Medical Center and were scheduled to receive surgery for kidney cancer were eligible to participate in this study.

Conference Presentation Video

Process of consent

Conference Presentation Video

The study, risks, and benefits of participation were presented to subjects. Consent forms approved by the University of Rochester’s Research Subjects Review Board (RSRB) were presented to patients at the time of their office visit.

Conference Presentation Video

Inclusion/exclusion criteria

Conference Presentation Video

Subjects enrolled in this study were males and females over the age of 18 who were undergoing surgery at Strong Memorial Hospital for renal cancer. Pediatric patients (<18years of age) were not enrolled. Samples from patients with known infectious diseases such as AIDS, hepatitis, and tuberculosis were not collected. Optical property measurements were performed on 16 kidneys collected over two years, with interpretable spectra obtained from 14 of these subjects.

Conference Presentation Video

Sample collection

Conference Presentation Video

Kidneys were collected during surgical nephrectomies. There was no deviation from routine clinical surgical practice. Specimens collected during surgery were handled under the supervision of a pathologist prior to, during, and after the optical property measurements. If the surgeon decided that any excised kidney required immediate evaluation by pathology, that specimen was not included in this study.

Conference Presentation Video

Nephrectomy specimens were brought directly from the operating room to the laboratory for optical measurements. A pathologist performed the initial dissection of the specimen as per standard surgical pathology grossing procedures on a bench adjacent to the diffuse reflectance spectroscopy setup. This ensured a clean, flat-cut surface against which the spectroscopy probe could be positioned with good contact. By visual inspection, the pathologist then identified normal and neoplastic tissue on which optical measurements would be performed. A representative dissected kidney is shown in Fig. 1. The measurements were taken using a custom diffuse reflectance spectroscopy system. The spectroscopy probe, shown schematically in Fig. 2, made gentle contact with the surface of the tissue being interrogated. The nephrectomy specimen was then immediately taken to the surgical pathology laboratory, where it was accessioned and processed as per routine practice.

Conference Presentation Video

Schematic of the optical probe used in contact with the tissue surface for diffuse reflectance measurements. The source fiber launched broadband light into the tissue, and 10 detection fibers, located at distances of 1 to 10 mm from the source fiber, collected the diffusely reflected signal.

Conference Presentation Video

Sample storage and confidentiality

Conference Presentation Video

The kidney samples were assigned a code number at the time of collection. We created a database of the tumor and normal kidney optical properties, medical record numbers, and the pathological diagnosis.

Conference Presentation Video

Diffuse Reflectance Spectroscopy

Conference Presentation Video

White light, steady-state diffuse reflectance measurements were performed at each of the normal and neoplastic regions of the specimen using techniques that have been developed and optimized in the laboratory over the past decade. Our instrumentation and method of data reduction for performing quantitative diffuse reflectance spectroscopy have been described in detail elsewhere.13,14 Broadband light from a xenon lamp was delivered to the surface of the tissue via an optical fiber with a 400-μm core diameter, which is terminated in a probe of our own design and construction (Fig. 2). The probe also holds a linear array of 20 detection fibers (200-μm core diameter) positioned at regular 1-mm spacings at distances of 1 to 20 mm from the source fiber. In the studies reported here, adequate signal-to-noise was obtained from the first 10 of these fibers, corresponding to source-detector separations ranging from 1 to 10 mm. We confirmed that measurements from these 10 fibers were sufficient by systematically deleting data from the most remote fiber and refitting the remaining points until a change in the returned optical properties was observed. For all kidneys analyzed, five to eight source-detector separations were adequate for accurate recovery of optical properties.

Conference Presentation Video

The flat, polished surface of the probe made good contact with the clean, flat-cut surface of the tissue under investigation, allowing white light from the source fiber to be injected into the specimen surface. The distal ends of the detection fibers were gathered into a bundle that terminated in a ferrule. At the end of the ferrule, the fibers were arranged in a row and positioned in the focal plane of a 0.275-m imaging spectrograph (SpectraPro-275, Princeton Instruments/Acton, Acton, MA). The fibers were imaged on to the surface of a thermoelectrically cooled, 512×512, 16-bit CCD camera (Pixis512, Princeton Instruments/Acton, Trenton, NJ). The detector fiber position was imaged along one axis of the CCD array, and the spectral information over the range of 630 to 800 nm was recorded along the other. Thus, at each wavelength in the data set, the experiment captured a spatially resolved attenuation measurement. Typical signal integration times for measurements described here were approximately 5 s. The raw data were corrected for system throughput and fiber cross-talk and then processed through a Savitzky-Golay smoothing filter. From these attenuation measurements, the absorption (μa) and transport scattering coefficients (μs′) of the tissue were extracted.

Conference Presentation Video

Extraction of Optical Properties

Conference Presentation Video

A broadly accepted theoretical description of light propagation in biological tissues is the Boltzmann transport equation.15 This equation treats light propagation as energy transport, in which absorption and scattering are the dominant attenuation mechanisms. We have previously demonstrated the accuracy of an approximation to this expression known as the P3 approximation, which expands the Boltzmann transport equation in Legendre polynomials up to third order.16 The diffuse reflectance predicted by this model explicitly depends on the optical properties of the sample and can therefore be used to extract these properties from measured data. The P3, rather than the diffusion approximation, was used here because of its validity over a wider range of possible optical properties.

Conference Presentation Video

When light is introduced into tissue, it is attenuated by wavelength-dependent scattering and absorption. In a diffuse reflectance measurement, a fraction of the multiply scattered light is emitted from the tissue surface, where it is collected by one of the probe fibers, as illustrated schematically in Fig. 2. With increasing distance from the source fiber, the amplitude of the reflectance signal is diminished, because light exiting the tissue at larger source-detector separations has traversed a greater path length. These spatially resolved diffuse reflectance data may be fit using an appropriate theoretical model to obtain μa and μs′ without assumptions regarding the composition of the tissue. Here, the P3 model of light transport was fit to the measured reflectance using a constrained nonlinear optimization routine found in the MATLAB optimization toolbox (Mathworks, Natick, MA). The fitting function, χ2, to be minimized was the sum of the squared error between the measured dataset at each wavelength and the P3 model, scaled by the measured diffuse reflectance at each source-detector separation, as shown in the equation Display Formula

χ2=∑j=110{[D(r=rj)−R(r=rj,μa,μs′)]/D(r=rj)}2,(1)

where D(rj) are radially resolved measurements and R(rj,μa,μs′) are reflectance values at positions rj computed using the P3 approximation with test optical properties. All 10 source-detector separations, rj, were used simultaneously to compute χ2 for each set of test optical properties. Constrained minimization was performed using the interior point algorithm included in MATLAB’s fmincon function, with μa and μs′ extracted as free parameters. The values of μa and μs′ were both constrained to be positive in the fit, and the value of the scattering anisotropy (g) was assumed to be 0.85.16,17 The transport scattering coefficient, μs′, was constrained to be monotonically decreasing with increasing wavelength, as this has been demonstrated experimentally.18,19 Values of γ=1.85 and δ=2.57 were assumed in the P3 approximation, using a Henyey-Greenstein phase function with g=0.85. The values of μa and μs′ corresponding to a minimized χ2 were stored in an array. The radially resolved measurement for each wavelength in the dataset was analyzed using an identical procedure.

Conference Presentation Video

Monte Carlo Model of Interstitial PDT

Conference Presentation Video

PDT of renal tumors would be performed by introducing an optical fiber or fibers directly into the tumor under image guidance. Simulations of light propagation in tissue from cylindrical diffusing tip fibers were performed using Monte Carlo (MC) modeling. The particular code used was developed in our laboratory and incorporates a rigorous model of cylindrical diffusing fibers implemented on a graphics processing unit.20 Briefly, the model uses a realistic description of cylindrical diffusing fibers to launch photon packets into the surrounding tissue, with standard MC techniques being used to control photon propagation.21 Interactions between photon packets and the diffuser after launch are also considered.

Conference Presentation Video

MC simulations were used to model two scenarios. The first consisted of an ellipsoidal tumor embedded in a normal kidney. The tumor had equatorial radii of 1.1 cm and a polar radius of 1.3 cm, with a semi-infinite region of normal kidney tissue surrounding it. This tumor size was determined by averaging measured tumor dimensions from clinical images and assuming an ellipsoidal tumor shape. In this case, the optical properties used were those extracted from P3 fits to measured diffuse reflectance in both normal kidney and tumor tissue, with each being assigned to the appropriate region of the sample. Based on the treatment wavelengths used for currently available and anticipated PDT photosensitizers, simulations were run at 630, 652, 670, 692, and 780 nm for each kidney sample.

Conference Presentation Video

The second model consisted of a homogeneous sample of either normal kidney or tumor tissue. Optical properties were set uniformly to match those extracted from P3 fits to measured diffuse reflectance in normal or tumor tissue. Simulations of both normal and tumor tissue were run at 50 discrete wavelengths from 630 to 800 nm for each kidney sample. This allowed for modeling of a characteristic light propagation distance at each of the measured wavelengths.

Conference Presentation Video

For all MC simulations, 1,000,000 photons were launched from a 2-cm-long diffuser. In all cases, the tissue sample was represented with a 100×100×100 grid of cuboid voxels, each with dimensions of 0.026 cm. A Henyey-Greenstein phase function with g=0.85 was used for all simulations.

Conference Presentation Video

Diffuse Reflectance Spectroscopy

Conference Presentation Video

We analyzed data from 14 kidneys. The results obtained from two kidneys were not interpretable due to low signal levels. Figure 3 shows representative diffuse reflectance spectra acquired from the RCC region of a freshly excised human kidney. The spectra were corrected for measured system throughput, background signal, and cross-talk in the CCD. Individual traces correspond to different source-detector separations on the diffuse reflectance probe, which range from approximately 1 to 10 mm.

Conference Presentation Video

Measured diffuse reflectance spectra at the first five source-detector separations collected from a kidney region identified as RCC. Spectra were corrected for measured optical system throughput, background, and cross-talk in the CCD. Individual traces correspond to different source-detector separations on the diffuse reflectance probe, with their locations relative to the source fiber given in the legend.

Conference Presentation Video

Kidney Optical Properties

Conference Presentation Video

The optical properties of tumor and normal kidney tissue were extracted at each wavelength by fitting a P3 model of light transport to the measured data using a constrained nonlinear optimization routine as described in Sec. 2. An example of such a fit is shown in Fig. 4. The individual data points are corrected, normalized diffuse reflectance measurements made at 780 nm in a representative RCC. Each point was acquired at a different source-detector separation [r (mm)] using the probe as described above. The fit shown in Fig. 4 yielded optical properties of μa=0.072mm−1 and μs′=1.40mm−1.

Conference Presentation Video

The individual data points are corrected, normalized diffuse reflectance measurements made at 780 nm in a representative RCC. Each point was acquired at a different source-detector separation [r (mm)]. Error bars are the standard deviations of three repeated measurements of the same sample. The fit of the P3 model to these data yielded optical properties of μa=0.072mm−1 and μs′=1.40mm−1.

Conference Presentation Video

By analyzing spatially resolved reflectance at each measured wavelength, absorption and scattering spectra over the measured spectral range were determined for each kidney sample. Examples from two lesion types are shown in Fig. 5. In these two cases, μa was lower in the tumor than in the normal kidney tissue. However, this trend was not reproducible across all of the kidney samples (data not shown). The spectra exhibit decreased absorption and scattering with increasing wavelength, as expected. The absorption feature near 759 nm corresponds to a peak characteristic of deoxyhemoglobin.

Conference Presentation Video

Values of (a) μa and (b) μs′ extracted from fitting the P3 model to measured diffuse reflectance spectra from a kidney diagnosed with RCC. In this case, both μa and μs′ are lower for tumor tissue than for kidney parenchyma. However, this was not reproducible for all RCC samples. (c) μa and (d) μs′ extracted from fitting the P3 model to measured diffuse reflectance spectra from a urothelial carcinoma and adjacent normal parenchyma. In this case, μs′ is higher in the tumor.

Conference Presentation Video

The recovered values of μa ranged from 0.016 to 0.14mm−1 for RCC conventional; 0.003 to 0.22mm−1 for oncocytoma; 0.02 to 0.062mm−1 for urothelial carcinoma; 0.14 to 0.31mm−1 for RCC, papillary type 1; 0.5 to 0.75mm−1 for RCC, papillary multifocal; and 0.017 to 0.13mm−1 for kidney parenchyma. The recovered values of μs′ ranged from 0.48 to 1.92mm−1 for RCC conventional; 1.23 to 2.5mm−1 for oncocytoma; 1.06 to 1.33mm−1 for urothelial carcinoma; 0.25 to 0.3mm−1 for RCC, papillary type 1; 1.32 to 1.83mm−1 for RCC, papillary multifocal; and 0.11 to 5.51mm−1 for kidney parenchyma. These recovered optical properties yielded values of the transport albedo [a′≡μs′/(μa+μs′)] that ranged from 0.49 to 0.99. The recovery of transport albedos lower than 0.98 indicates that the P3 approximation was an appropriate choice over the diffusion approximation.16

Conference Presentation Video

Monte Carlo Modeling

Conference Presentation Video

Monte Carlo simulations of light propagation from cylindrical diffusing fibers were run using the optical properties extracted from each of the kidney samples. A representative fluence distribution generated by the delivery of 100Jcm−1 of 630-nm light from a 2-cm-long diffuser is shown in Fig. 6. This wavelength corresponds to the long wavelength absorption maximum of the FDA-approved photosensitizer Photofrin. The figure shows fluence deposited in the surrounding tissue using optical properties that were derived from a kidney diagnosed with urothelial carcinoma. In order to examine this fluence distribution more carefully, cuts were made through the fluence distribution perpendicular to the diffuser axis. The results of this are shown in Fig. 7(a) for a 2-cm-long diffuser delivering 100Jcm−1, with each contour line corresponding to a fluence increment of 50Jcm−2. The outer black line in the figure represents the boundary between tumor and normal tissue, as described previously. The inner black line represents the boundary between tumor tissue and the diffuser.

Conference Presentation Video

3D rendering of the simulated fluence distribution around an optical fiber with the proximal end of a 2-cm-long diffusing segment of the fiber positioned at z=0.3cm above the boundary separating the tumor and normal tissue regions. Tumor optical properties were set to μa=0.062mm−1, μs′=1.34mm−1, and g=0.85, which correspond to values extracted from an urothelial carcinoma at 630 nm, as shown in Fig. 5(c) and 5(d). Voxel size was 0.026×0.026×0.026cm3.

Conference Presentation Video

Contour plots of cuts through simulated fluence distributions around a 2-cm-long diffusing fiber at (a) 630 nm and (b) 780 nm, with a total fluence of 100 J cm−1 delivered. Each contour line corresponds to a fluence increment of 50 J cm−2. The simulated volume consists of an ellipsoidal tumor region, with equatorial radii of 1.1 cm and a polar radius of 1.3 cm, surrounded by a semi-infinite region of normal tissue. The boundary between tumor and normal tissue is shown as the outer black line; the boundary with the diffuser is shown as the inner black line. The optical properties correspond to values extracted from the tumor and normal regions of a kidney diagnosed with urothelial carcinoma, as shown in Fig. 5(c) and 5(d).

Conference Presentation Video

Figure 5 reproduces in the human kidney the well-known observation that absorption and scattering in tissue both decrease with increasing wavelength in this spectral range. It is therefore appealing to consider the use of photosensitizers with longer wavelength absorption in interstitial PDT of the kidney. This potential is shown in Fig. 7, where fluence deposition is modeled at 630 and 780 nm. At the longer wavelength, the penetration of light is greatly improved compared to 630 nm. For a delivered fluence of 100Jcm−1, ∼20% of the tumor volume received a light dose of at least 50Jcm−2 at 630 nm. This increased to approximately 50% of the tumor volume at 780 nm. Additionally, nearly 30% of the tumor received a light dose of 100Jcm−2 or more at 780 nm, compared to ∼11% at 630 nm. This trend was reproducible across all samples, with the percentage of the tumor receiving a prescribed light dose increasing by an average of 2.2-fold as a result of changing the wavelength from 630 to 780 nm. Increasing the fluence delivered by the cylindrical diffuser to 600Jcm−1 extended the 50Jcm−2 contour to the simulated boundary of the tumor (not shown).

Conference Presentation Video

In order to describe a characteristic distance of light propagation from a cylindrical diffusing fiber, MC simulations were run using homogeneous optical properties in a semi-infinite tissue volume surrounding a 2-cm-long diffuser. A zeroth order modified Bessel function of the second kind, K0, is a useful analytic approximation to light emitted from an infinitely long cylinder. Thus, a function of the form Display Formula

ϕ(r)=K0(r/δ)(2)

was fit to radially averaged cuts through simulated fluence distributions at the axial midpoint of the diffuser. Here, ϕ(r) is the fluence at radial position r, K0 is a zeroth order modified Bessel function of the second kind, and δ is a characteristic light propagation distance, which we define as the distance from the diffuser over which light is attenuated to ∼40% of its value at the fiber surface. A representative fit to a simulation run with μa=0.017mm−1 and μs′=1.55mm−1 is shown in Fig. 8. The optical properties used correspond to values extracted from a typical RCC at 780 nm. The fit shown is typical and of good quality. The δ value extracted from this fit was 3.43 mm. The results of this are shown in Fig. 9 for optical properties from a typical kidney diagnosed with RCC. An increased penetration depth with increasing wavelength was found for all measured tumor samples. The mean penetration depths for the five tumor types in this study and for the normal kidneys are collected in Table 1 for five wavelengths selected for relevance to clinical PDT.

Conference Presentation Video

Results of fitting a zeroth order modified Bessel function of the second kind, K0 [Eq. (2)], to a cut through the simulated fluence around a 2-cm-long diffuser. The fluence plotted is the radially averaged simulated fluence at the axial midpoint of the diffuser. Optical properties were set to μa=0.017mm−1, μs′=1.55mm−1, and g=0.85, which correspond to values extracted from a representative RCC at 780 nm. From such fits, the characteristic propagation distance, δ, is obtained.

Conference Presentation Video

Discussion

Conference Presentation Video

Renal parenchymal tumors, most frequently RCC, are increasingly being detected earlier as incidental radiological findings. These tumors are often detected when they are small enough such that radical nephrectomy, either by open surgery or laparoscopy, is not necessary. Patients who have only one kidney, have bilateral disease, or are too frail to undergo surgery are not normally treated by nephrectomy. Minimally invasive, kidney-sparing therapies based on percutaneous approaches such as RFA and cryotherapy have their own limitations, particularly in well perfused organs like the kidney. There is therefore a compelling rationale for considering image-guided, interstitial PDT in this context, and this motivated the study of the optical properties of RCC, other less common lesions, and normal kidney parenchyma.

Conference Presentation Video

The light dose delivered by a cylindrical diffusing fiber in interstitial PDT is generally specified in joules per unit length of diffuser, Jcm−1. This fluence is provided by delivering light to the fiber at a given fluence rate, in mWcm−1, for a prescribed treatment time. These fluence rates22,23 are typically in the range of 100 to 400mWcm−1, with total fluences of 150 to 300Jcm−1. The fluence delivered by the single diffuser in the simulations reported here ranged from 100 to 600Jcm−1, which, at a typical fluence rate of 150mWcm−1, corresponds to treatment times of 11 to 66 min. If deemed clinically necessary, these times could be reduced by increasing the fluence rate or inserting additional fibers.

Conference Presentation Video

The characteristic propagation distance, δ, which we have defined, corresponds to the distance in the tissue at which the deposited fluence falls to ∼40% of the maximum value immediately adjacent to the fiber. The propagation distances shown in Table 1 translate to diameters around a single diffuser of ∼2.6mm at 630 nm to 5.6 mm at 780 nm. It is important to emphasize that this does not define the zone that can be treated effectively. As illustrated in Fig. 7, the area that receives a specified, clinically relevant fluence is typically larger than the characteristic propagation distance. Also, as suggested above, multiple diffuser fibers can be used to treat a larger volume. This is commonly done in PDT treatment of the prostate,24 where fibers are inserted on a grid specified by a template, with an approximate spacing of 5 to 7 mm between diffusers. The size and arrangement of this grid can be tailored to fit the tumor of an individual patient. Treatment plans designed to optimize diffuser placement and the total fluence delivered may be developed based on estimates of optical properties or on their direct interstitial measurement. The heterogeneity in tumor optical properties reported in this study provides strong rationale for direct pretreatment optical property determinations using interstitial optical fiber-based spectroscopy. PDT may also be informed by radiation therapy treatment plans, which specify therapeutic doses to target volumes and identify sensitive surrounding healthy tissue.

Conference Presentation Video

The wavelengths used for MC simulation were chosen to correspond to the absorption maxima of several clinically relevant PDT photosensitizers. The wavelength of 630 nm corresponds to Photofrin, a porphyrin photosensitizer, which is currently approved to treat a number of diseases, including esophageal, endobronchial, and bladder cancers and Barrett’s esophagus with high grade dysplasia.12,25 Foscan, also known as Temoporfin or mTHPC, has an absorption maximum of 652 nm and is approved in Europe for recurrent, inoperable head and neck cancers, and it is being investigated in ongoing clinical trials for the treatment of other neoplastic conditions.25,26 Simulations performed at 670 nm were chosen to correspond to the absorption of photosensitizers such as HPPH or methylene blue. HPPH has been used in the treatment of thoracic malignancies,27,28 Barrett’s esophagus,29 and early cancers and precancers of the oral cavity.30 Methylene blue is FDA-approved as a diagnostic dye and is the only known treatment for methemoglobinemia. Methylene blue has also been used successfully as a photosensitizer in a number of PDT studies.31–33 Benzoporphyrin derivative (BPD), also known as verteporfin, has an absorption maximum of 692 nm and is approved worldwide for the treatment of age-related macular degeneration.12 The wavelength of 780 nm was chosen to correspond to the absorption of some emerging photosensitizers, such as TOOKAD34 and others still in preclinical development,35 which are designed to take advantage of increased penetration of light in this spectral range. Based on our findings, these next-generation sensitizers with high absorption at 780 nm and beyond would be favorable candidates for clinical trials of interstitial PDT of renal cancers.

Conference Presentation Video

Acknowledgments

Conference Presentation Video

The authors thank Susan Schoen, Elizabeth M. Smith, Maureen Kiernan, and Kay Rust in the Department of Urology for help with patient recruiting and tissue accrual. We thank Laurie Baxter from the Department of Pathology for coming to our laboratory to section the kidneys. We are grateful for JoAnne McNamara’s help with Research Subject Review Board protocols and consent forms. This work was supported by NIH Grant CA68409 awarded by the National Cancer Institute.

Figures

Schematic of the optical probe used in contact with the tissue surface for diffuse reflectance measurements. The source fiber launched broadband light into the tissue, and 10 detection fibers, located at distances of 1 to 10 mm from the source fiber, collected the diffusely reflected signal.

The individual data points are corrected, normalized diffuse reflectance measurements made at 780 nm in a representative RCC. Each point was acquired at a different source-detector separation [r (mm)]. Error bars are the standard deviations of three repeated measurements of the same sample. The fit of the P3 model to these data yielded optical properties of μa=0.072mm−1 and μs′=1.40mm−1.

Contour plots of cuts through simulated fluence distributions around a 2-cm-long diffusing fiber at (a) 630 nm and (b) 780 nm, with a total fluence of 100 J cm−1 delivered. Each contour line corresponds to a fluence increment of 50 J cm−2. The simulated volume consists of an ellipsoidal tumor region, with equatorial radii of 1.1 cm and a polar radius of 1.3 cm, surrounded by a semi-infinite region of normal tissue. The boundary between tumor and normal tissue is shown as the outer black line; the boundary with the diffuser is shown as the inner black line. The optical properties correspond to values extracted from the tumor and normal regions of a kidney diagnosed with urothelial carcinoma, as shown in Fig. 5(c) and 5(d).

Measured diffuse reflectance spectra at the first five source-detector separations collected from a kidney region identified as RCC. Spectra were corrected for measured optical system throughput, background, and cross-talk in the CCD. Individual traces correspond to different source-detector separations on the diffuse reflectance probe, with their locations relative to the source fiber given in the legend.

Values of (a) μa and (b) μs′ extracted from fitting the P3 model to measured diffuse reflectance spectra from a kidney diagnosed with RCC. In this case, both μa and μs′ are lower for tumor tissue than for kidney parenchyma. However, this was not reproducible for all RCC samples. (c) μa and (d) μs′ extracted from fitting the P3 model to measured diffuse reflectance spectra from a urothelial carcinoma and adjacent normal parenchyma. In this case, μs′ is higher in the tumor.

3D rendering of the simulated fluence distribution around an optical fiber with the proximal end of a 2-cm-long diffusing segment of the fiber positioned at z=0.3cm above the boundary separating the tumor and normal tissue regions. Tumor optical properties were set to μa=0.062mm−1, μs′=1.34mm−1, and g=0.85, which correspond to values extracted from an urothelial carcinoma at 630 nm, as shown in Fig. 5(c) and 5(d). Voxel size was 0.026×0.026×0.026cm3.

Results of fitting a zeroth order modified Bessel function of the second kind, K0 [Eq. (2)], to a cut through the simulated fluence around a 2-cm-long diffuser. The fluence plotted is the radially averaged simulated fluence at the axial midpoint of the diffuser. Optical properties were set to μa=0.017mm−1, μs′=1.55mm−1, and g=0.85, which correspond to values extracted from a representative RCC at 780 nm. From such fits, the characteristic propagation distance, δ, is obtained.

Return to: Optical property measurements establish the feasibility of photodynamic therapy as a minimally invasive intervention for tumors of the kidney

Copyright in the material you requested is held by SPIE (unless otherwise noted).
This email ability is provided as a courtesy, and by using it you agree that you are requesting the material
solely for personal, non-commercial use, and that it is subject to SPIE'sTerms of Use.
The information provided in order to email this topic will not be used to send unsolicited email, nor will it be
furnished to third parties. Please refer to SPIE'sPrivacy Policy for further information.

Shibboleth is an access management service that provides single sign-on to protected resources. It replaces the multiple user names and passwords necessary to access subscription-based content with a single user name and password that can be entered once per session. It operates independently of a user's location or IP address. If your institution uses Shibboleth authentication, please contact your site administrator to receive your user name and password.

Shibboleth is an access management service that provides single sign-on to protected resources. It replaces the multiple user names and passwords necessary to access subscription-based content with a single user name and password that can be entered once per session. It operates independently of a user's location or IP address. If your institution uses Shibboleth authentication, please contact your site administrator to receive your user name and password.

Shibboleth is an access management service that provides single sign-on to protected resources. It replaces the multiple user names and passwords necessary to access subscription-based content with a single user name and password that can be entered once per session. It operates independently of a user's location or IP address. If your institution uses Shibboleth authentication, please contact your site administrator to receive your user name and password.

Shibboleth is an access management service that provides single sign-on to protected resources. It replaces the multiple user names and passwords necessary to access subscription-based content with a single user name and password that can be entered once per session. It operates independently of a user's location or IP address. If your institution uses Shibboleth authentication, please contact your site administrator to receive your user name and password.

Access to SPIE eBooks is limited to subscribing institutions and is not available
as part of a personal subscription. Print or electronic versions of individual SPIE
books may be purchased via SPIE.org.

Shibboleth is access management services that provide single sign-on to protected resources. It replace the multiple user names and passwords necessary to access subscription-based content with a single user name and password that can be entered once per session. It operates independently of a user's location or IP address. If your institution uses Shibboleth authentication, please contact your site administrator to receive your user name and password.