1.IntroductionThe lungs are a significant reservoir of metastatic burden, with isolated pulmonary metastases frequently affecting patients with sarcomas and colorectal carcinomas.1 Approximately, 20% of sarcoma patients and 10% to 25% of colorectal cancer patients will develop isolated pulmonary metastases.2,3 To date, surgical metastasectomy is the mainstay of treatment, yet only a small subset (5% to 15%) of patients with isolated disease qualify for surgical resection. Because surgery can only remove visible nodules, micrometastatic disease, which makes up a significant proportion of the metastatic burden, is inadequately targeted, driving recurrence rates. Chemotherapy is commonly administered as an adjuvant treatment but fails to control or eradicate the occult metastatic burden. Several factors have been postulated to explain the poor response to chemotherapy, including aggressive tumor biology, the development of chemotherapy-resistant clones of tumor cells, and inadequate chemotherapy doses limited by systemic toxicities.4 Even in the setting of optimal treatment, long-term survival is limited, with a 5-year overall survival of 20% to 40% for sarcoma and 54% to 68% for colorectal cancer patients.2,3,5,6 Accordingly, alternative treatments are needed to effectively eradicate micrometastatic disease across large tissue volumes without the possibility of selecting therapy resistant clones, leading to treatment failure. 1.1.Photodynamic Therapy and DosimetryPhotodynamic therapy (PDT) is a promising cancer treatment modality that is less susceptible to the genetic expression of tumors and can avoid treatment resistance. PDT is based on the administration of a photosensitive drug which, upon absorption of light of a specific wavelength, generates cytotoxic reactive oxygen species, particularly singlet oxygen, leading to cell death when the photosensitizer-light product exceeds a particular minimum or threshold. PDT was used clinically as early as the 1980s, albeit some limitations have prevented it from gaining widespread use in oncology.7–9 Initial photosensitizers have poor tumor selectivity and long retention times in proliferating tissues, rendering patients sensitive to sunlight for several days to weeks. Given the high attenuation of the photons activating the photosensitizers, particularly by hemoglobin and cytochromes, and the increased path length due to high photon scattering by the lung’s connective tissue and frequent air-tissue interfaces, light penetration in the lung is poor. Hence, PDT was used only for small, early-stage tumors, accessible by bronchoscope or as a palliative treatment to relieve obstruction.10–12 Newer photosensitizers and photosensitizing precursors such as Chlorins, marketed as Temoporfin or Foscan in Europe, and 5-aminolevulinic acid (5-ALA), and others,10,13,14 respectively, improve tumor selectivity versus normal cells, while reducing the time for skin photosensitivity. Selective uptake and retention of the photosensitizer in tumors results in higher singlet oxygen concentration in the tumor relative to surrounding healthy tissue, eradicating it while limiting collateral damage. Another development is advances in treatment planning and dosimetry. Akin to radiation oncology, titrating treatment effects to correctly predict the target destruction while minimizing toxicity to healthy tissue is essential. For effective treatment in well-oxygenated tissues such as the lung, including their micrometastases, the PDT dose must exceed a critical threshold in the tumor tissue while not exceeding it in the normal lung.15–17 The threshold value considers a tissue’s ability to mitigate the cytotoxic burden. It is given by its photosensitizer concentration [PS], molar extinction coefficient, , and the optical fluence or photon density, , according to .18,19 The fluence must be adjusted to meet the required threshold conditions throughout the clinical target volume for a given (PS) in the tumor and the normal host tissue. 1.2.Lung PerfusionIsolated lung perfusion techniques were conceived to deliver treatment to the lung without risk of systemic toxicity.20–22 Based on the Toronto “ex vivo lung perfusion,” a version of isolated “in vivo lung perfusion” (IVLP) was developed employing the identical circuit and perfusion technique to provide a repeatable and non-injurious lung platform, whereby it causes no long-term reduction in lung function. The safety of this technique was demonstrated in pre-clinical acute and survival porcine models for the delivery of chemotherapy.23–25 IVLP drains blood from the lung, resulting in minimal cellular and hemoglobin materials in the lung vasculature after replacement with a low cellular perfusion solution, STEEN (XViVO, Göteborg, Sweden), overcoming one of the light penetration limitations. As shown in our previous work, the low-cellular perfusate results in a significant reduction in the effective attenuation coefficient driven mostly by a decrease in the absorption coefficient. This translated to a 3.3-mm increase in the light penetration depth, equivalent to a three-fold increase in the potential treatment volumes when using a 2-cm long cylindrical diffuser.26 IVLP also provides surgical access to the chest cavity for light delivery and allows control of the lung treatment environment (e.g., temperature, pH, drug concentrations, cointerventions). 1.3.Rationale and AimsTwo key parameters necessary for lung PDT treatment planning, especially in the context of lung perfusion, have been previously identified.26 The first parameter is the optical properties of the lung during both STEEN and blood perfusion, including the effective attenuation, absorption, and reduced scattering coefficients at three wavelengths relevant to the activation spectra of different photosensitizers. The second parameter is the photodynamic threshold of 5-ALA-induced protoporphyrin IX (PpIX) or Chlorin e6-mediated PDT of normal porcine lung tissue.26 This study aims to develop a protocol for whole-lung perfusion-assisted PDT with an adequate safety profile and high therapeutic potential. This requires developing and evaluating a light delivery system for homogenous illumination of the entire lung to the minimum required fluence () across the treatment planning volume, here ideally an entire lung. 2.Materials and Methods2.1.AnimalsAll animals received humane care as per the Principles of Laboratory Animal Care formulated by the National Society for Medical Research and the Guide for the Care of Laboratory Animals. The study protocol was approved by the Animal Care Research Committee at the Toronto General Hospital Research Institute. A total of 10 Yorkshire male pigs weighing 35 to 40 kg were used in this study. 2.2.Testing the Performance of the Light-Emitting Diode Light Delivery SystemFor this experiment, three Yorkshire male pigs weighing 35 to 40 kg were used and underwent a simplified version of a previously described porcine IVLP survival model.25 Briefly, a left thoracotomy was performed with dissection of the left pulmonary artery, left pulmonary vein, and left atrial cuff. After IV administration of 5000 IU of heparin, the left pulmonary artery and atrial cuff were clamped. Following cannulation of the pulmonary artery, a venotomy was created in the left atrium between the left upper and lower pulmonary veins, and the lung was flushed anterograde through the pulmonary artery cannula with 0.5 to 0.75 L of Perfadex, low potassium dextran solution (XVIVO, Goteborg, Sweden). After the lung flush, up to four custom-built, 5-cm diameter and 1.5-cm thick cylindrical, light-emitting diode (LED) discs, emitting 630 or 660 nm light were placed into the thorax in transparent ultrasound covers (CIVCO, Coralville, Iowa) to maintain sterility [Fig. 1(a)]. Each LED’s emission was quantified at five positions before placement using an externally calibrated fiberoptic probe (IP85 isotropic probe, Écublens, Switzerland) with coupled power meters (PM130D, Thorlabs, Thorlabs Canada ULC, Saint-Laurent, QC, Canada). The inter- and intra-source light fluence rate was and at 630 and 660 nm, respectively. Two light sources were placed on each of the upper and lower lobes of the lung. Three calibrated isotropic optical fiber detectors (IP85, Medlight, Switzerland) for local fluence rate quantification were placed endobronchially using a bronchoscope (Olympus Medical Canada, Richmond Hill, Ontario, Canada) to ensure the detectors were guided into separate lung segments [Fig. 1(b)]. These fluence-rate detectors were placed deep within the lung past the third order bronchi, and four isotropic fibers were placed between the light source and the lung surface to monitor the sources’ irradiance. The fluence-rate detectors were connected to an in-house built eight-channel LabVIEW-controlled photodiode detector system.27 The fluence rate sensor responsivity was corrected for the change of the detectors’ surrounding medium refractive index as required.28 A CT scan was performed to provide the spatial relationship between isotropic fibers and light sources [Fig. 1(c)]. As the light sources have metal components, 3D-printed identically sized hollow phantoms replaced the sources for imaging. After the experiment, animals were exsanguinated under anesthesia. To create an in silico model visualizing the fluence-rate distribution in the lung, based on the CT scan, the thorax was segmented into regions comprising light sources, optical fiber detectors, the bronchial tree, lung parenchyma, and bone. These regions were rendered into volumes represented tetrahedrally, thus creating the in silico model to visualize the fluence-rate distribution in the lung. The bony thorax had an insignificant effect on the fluence rate distribution and was excluded to reduce the model’s file size. The fluence rate distribution inside the lung was modeled by Monte-Carlo using FullMonte using the tissue optical properties for the various optical components as listed in Table 1.26,31 The validity of the simulation model has been previously described in various contexts, including the lung.31–34 Table 1Optical properties used for the FullMonte-based light propagation simulations.26,29,30
Simulations were performed for all sources independently and with all sources together. The fluence rate measured by the isotropic detector fibers was plotted as a function of the distance to each source and compared with the in silico simulation results to confirm the optical properties and the resulting fluence rate gradient. For experiment 1, the simulated and the measured fluence rates were also plotted against each other for comparison. Fluence rate histograms of the pig lung volume were created to estimate the effective lung volume targetable by this illumination approach. 2.3.Photodynamic Therapy Protocol Safety and Feasibility AssessmentSeven Yorkshire male pigs weighing 35 to 40 kg were used, four were administered 5-ALA (Sigma-Aldrich, Mississauga, Ontario, Canada), and three were administered Chlorin e6 (Synverdis GmbH, Heidelberg, Germany). The study used an accelerated dose-escalation study design, whereby the dose of PDT was doubled with each sequential case if no PDT-related toxicity was observed during the previous case.35 If the starting dose was not tolerated, the subsequent case used half the previous photosensitizer dose. For the first two cases with 5-ALA and the first with Chlorin e6, different radiant exposures () were delivered to the upper and lower lobes by adjusting the irradiation time, providing additional data points for the escalation study. Cases to confirm a maximally tolerated dose level were performed with only one radiant exposure for the upper and lower lobes, respectively, to achieve accurate fluence-volume-histograms. Lung physiology was assessed during the baseline procedure, at reperfusion, and 72 h post-PDT. Gas exchange was evaluated through arterial blood gas analysis. Airway dynamics were assessed through the measurement of dynamic and static compliance. Pigs underwent a simplified lung perfusion procedure, as described in Sec. 2.2. Once the lung flush had concluded, the PDT procedure was performed. Photosensitizers were given intravenously apart from the third case with Chlorin e6. Following 5-ALA administration and a 3- to 4-h interval, light exposure was started. The starting 5-ALA dose was . The three Chlorin e6 cases were administered and had drug light intervals of 1 h, 15 min, and 0 min, respectively. For the third chlorin case with a 0-min drug-light interval, a Chlorin e6 loading dose of Perfadex flush, equivalent to , was administered to achieve the highest tissue concentration over the entire light illumination time, as a worst-case scenario. Because the half-life of Chlorin e6 in plasma is , a continuous flush of of Chlorin e6 was maintained during the PDT illumination time.36 For the four cases administered 5-ALA and the first case using Chlorin e6, two LED sources, wrapped in transparent ultrasound covers, were placed on the upper and lower lobes, respectively. Isotropic detectors were secured atop the light source to monitor the irradiance. Additional isotropic detectors were placed on the lung surface and within the lung under bronchoscopy guidance, as described above in Sec. 2.2. 5-ALA-induced PpIX and Chlorin e6 were activated at 630 nm delivering 6 to and 660 nm delivering 7.2 to , respectively. The second and third Chlorin e6 cases utilized a fiberoptic-coupled modulight laser (665/808 OEM; Modulight Inc., Tampere, Finland) emitting 665 nm radiation with the goal to achieve a significantly higher administered PDT dose. The optical fiber was equipped with a microlens illuminating a 3.5-cm diameter spot on the lung surface. Irradiance was delivering a total fluence of 240 to . At the conclusion of PDT, the pigs were recovered and treated over a 72-h survival period. Pigs underwent a final procedure where endpoint lung toxicity assessments were made prior to euthanasia by exsanguination under anesthesia. Lung toxicity assessment included physiological, gross, CT, and histological examination. A CT scan was performed at 72 h immediately before the final assessment. The lungs were harvested while inflated after sacrificing for assessment of gross appearance. Tissue biopsies were collected directly before and after PDT from the lung periphery to quantify the photosensitizer concentration. At 72 h, lung biopsies were taken from the upper and lower lobes centered around areas directly abutting the LED sources. Tissue samples were fixed in 10% phosphate-buffered formalin, embedded in paraffin, sectioned, and stained in hematoxylin and eosin for assessment by light microscopy using a pathologic acute lung injury assessment system described previously described.37 Toxicity was defined as a clinically significant impairment of lung function with a (P/F) ratio , a reduction of compliance in the treated lung below and/or signs of severe acute lung injury on gross, CT, and histological examination. 2.4.Photosensitizer Threshold AssessmentThe Chlorin e6 and ALA-induced PpIX concentrations were determined from pre- and post-PDT biopsies by tissue solubilization and fluorescence spectroscopy as published.38 The depth of necrosis was measured on H&E slides for all cases with visible necrosis. To quantify the PDT threshold value of both photosensitizers for normal lung tissue, the fluence, , at the necrotic boundary, as simulated using the previously determined tissue optical properties and the measured photosensitizer concentration from the tissue biopsies, was utilized. The molar extinction coefficients of PpIX ( at 635 nm) and Chlorin e6 ( at 660 nm) were taken from the Oregon Medical Laser Center website.39 The derived ALA-induced PpIX PDT threshold values were compared against our previously determined threshold values.26 3.Results3.1.Testing the Light-Emitting Diode Light Delivery System PerformanceFigure 2 shows segmentation examples of the bronchial tree, lung parenchyma, sources, detectors, and surrounding bony ribcage used in Monte Carlo photon propagation simulations. Figure 3(a) shows the measured fluence rate from sequential illumination by a single source versus a sensor’s relative distance to that source’s center position for the third case of ALA-induced PpIX. This slope provides the experimentally calculated effective attenuation coefficient compared with the lung tissue optical properties, including absorption and scattering coefficients, used in the Monte Carlo simulations. The experimentally derived effective attenuation coefficient was whereas simulated cases used an effective attenuation coefficient of for lung tissue. Figure 3(b) shows the comparison of the simulated and measured fluence rates for one experiment. Overall, the predicted fluence rate as modeled in the simulation, based on our previously determined optical properties26 and compared with the actual measurements from the detectors, shows a strong predictive value of the simulation model toward the fluence rate throughout the clinical target volume. To determine the attainable fluence rate distribution in the lung based on three 5-cm diameter sources operated simultaneously, dose-volume histograms based on the simulations model demonstrating the light fluence throughout the entire organ (Fig. 4). This figure shows the fluence rate for a lung irradiated with three LED sources emitting 665 nm photons with a total power of 2.43 W. The variation of the fluence rate was within two orders of magnitude over of the lung volume despite covering less than 15% of the lung surface by the three light sources. The dose-volume histograms from simulations and our direct measurements consistently demonstrated fluence rates and across and 80% of a single lung, respectively, comprising the clinical target volume. 3.2.Photodynamic Therapy Protocol Safety and Feasibility AssessmentAll animals tolerated the procedure well, with no adverse events from the surgical procedure. All cases demonstrated preserved lung function and physiology at the end of the PDT treatment with a stable P/F ratio above 300 mmHg, preserved lung compliance, and no signs of significant lung injury on gross or histological appearance. All animals survived to the 72 h endpoint in good clinical condition. 3.2.1.ALA dose de-escalationAn overview of the cases administered ALA is displayed in Fig. 5. Lung physiological assessments for the ALA cases are shown in Fig. 6. The first case using an ALA dose of demonstrated severe lung injury with a P/F ratio of 97, dynamic and static compliance of , signs of necrosis and severe inflammation on both gross and histological examinations. Injury in the upper lobe, irradiated at , was greater than for the lower lobe, irradiated at , as shown by a P/F ratio of 46 versus 189 from the left upper and lower pulmonary veins, respectively. The ALA dose was halved for case 2 and yielded a milder but still clinically significant injury to the normal lung with a P/F ratio of 440, dynamic and static compliance of , and a shallower necrosis and lung injury pattern on gross, CT, and histological examination. There was no noticeable difference in the injury extent between the upper and lower lobes irradiated at 12 and , respectively. Case 3 was repeated at but with a single light dose of to confirm toxicity at this optical dose, given case 2 had borderline clinically relevant lung injury. This case showed a more moderate lung injury with a P/F ratio of 248, dynamic and static compliance of , and a more profound degree of necrosis and inflammation on gross, CT, and histological appearance. The ALA dose was lowered to for case 4. This case demonstrated no lung injury with a P/F ratio of 456, dynamic and static compliance of 13 and , respectively, and no evidence of necrosis or inflammation on gross or histological examination. 3.2.2.Chlorin e6 dose escalationFigure 7 provides an overview of the cases administered Chlorin e6 and Fig. 8 shows lung physiological assessments for the Chlorin e6 cases. The first Chlorin e6 case delivered PDT comprised and a radiant exposure of 12 and for the upper and lower lobes, respectively, resulting in no evidence of necrosis or inflammation on gross, CT, or histological examination. The P/F ratio remained normal at 593, and dynamic and static compliance were 13 and , respectively. Given the lack of demonstratable PDT effect and our previous study showing no Chlorin e6 uptake in healthy tissue after a 1-h drug-light interval,26 an increase in drug dose far beyond clinical protocols may be required to produce a PDT effect in normal tissue. The subsequent two cases were performed to elucidate the minimum PDT threshold or maximally tolerated dose in normal lung tissue rather than a strict dose-escalation design. Case 2 was performed with a radiant exposure of via irradiation of a 5- to 6-cm diameter spot by collimated laser light to cause local necrosis. A shorter drug-light interval of 15 min was employed to retain some photosensitizer in the normal lung tissue. Again, no evidence of necrosis or inflammation on gross, CT, or histological examination was noted, with a P/F ratio of 459 and dynamic and static compliance of 13 and , respectively. To further minimize the drug-light interval, the final case was performed with Chlorin e6 administered through the lung flush continuously during the delivery of radiant exposure, whereby the photosensitizer acted on the tissue and the vasculature. This resulted in deep local lung injury over the irradiated area with a P/F ratio of 331 and dynamic and static compliance of . Necrosis extended up to 11 mm into the irradiated lung. 3.2.3.Photosensitizer threshold assessmentHistologic assessment of the depth of necrosis (Fig. 9) was used to establish the photodynamic threshold values for both photosensitizers, listed in Table 2. The necrosis threshold values for ALA-induced PpIX based on all data sets was . Given that the Chlorin e6 lung tissue concentrations were below the detection limit for the first two cases, it was impossible to calculate the photodynamic threshold of Chlorin e6 following i.v. administration. Based on the third case of Chlorin e6, the minimum threshold can be estimated to be at least when using Chlorin e6 as a vascular photosensitizer. Table 2Case-based factors measured for photodynamic threshold determination, including light dose delivered, tissue photosensitizer concentration, and depth of necrosis.
For ALA cases 1 and 2, the two depths of necrosis presented correspond to the two light doses administered. 4.DiscussionThis technique has several important implications when considering PDT treatments in the lung. The delivery of a selective PDT protocol requires accurate dosimetry that can be used to predict a desired treatment effect. Tissue response to PDT can be described by a threshold concept, which is defined as the optical fluence and tissue photosensitizer concentration required to result in PDT-induced necrosis of the target tissue.19 Tumor tissues often demonstrate a higher PDT threshold than cancer cells due to differences in their resilience from oxidative stress. Conversely, many photosensitizers work by exploiting increased accumulation and retention of photosensitizers in cancer cells compared with normal cells.13 These phenomena need to be balanced in favor of the host tissue to allow for a potential therapeutic window for PDT treatment above the threshold of the cancer cells but below that of normal tissue across the entire lung volume. The creation of a whole-lung PDT procedure targeting lung metastases has not been previously reported. The challenge is providing a sufficient fluence rate across the entire organ while maintaining sufficient selectivity of the PDT effect toward the metastases to spare healthy tissue. Herein, a path toward overcoming these challenges and delivering whole-lung PDT for diffuse metastatic disease is presented. A light delivery system for improved illumination of the lung with increased treatment volumes comprised of multiple cylindrical LED light sources at the desired light wavelengths, here 630 or 660 nm, is presented. The light source aperture area comprises up to and can deliver up to 3.235 W optical power. Our previous work compared the optical properties of the lung during normal blood perfusion conditions with perfusion using an acellular perfusion solution. Acellular perfusion resulted in significant reduction in the tissue absorption coefficient resulting in a 3.3-mm increased penetration depth of light and accordingly a three-fold increase in effective treatment volumes.26 Using these previously determined optical properties for the lung tissue as input parameters for Monte Carlo simulations of the fluence rate throughout the lung, based on the source positions obtained from CT imaging. The source intensities are subject to assumption of homogeneous, Lambertian emission. The lung tissue was assumed to be a turbid homogenous medium beyond the 8th generation airways as resolving for smaller airways has no significant effect on fluence rate calculations.40 Up to the 8th generation, the initial bronchial tree was segmented as empty air spaces. Comparison of the simulated fluence rate and effective attenuation with direct measurements taken during the cases exhibited high concordance, as shown in Fig. 3, validating this technique to generate dose-volume histograms. While there is a significant increase in the effective treatment volume of the lung, a further expansion of the light-emitting surface area is required, as at present of the surface area is directly illuminated. Other groups have previously developed flexible radiation-emitting devices that conform to complex surfaces, such as the Freiburg Flap41 or light-emitting textiles.42,43 The former, initially developed for brachytherapy application, can be outfitted with optical fibers in place of the radiation sources, to provide reliable light intra-operatively to the lung that is amenable to dosimetry planning.44 However, both techniques emit light toward the lung and the pleural cavity, and additional reflective materials must be placed between the light source and the pleural wall. A limitation of the LED device presented here may be the low attainable fluence rates of 1 to , which are significantly below 50 to typically employed in PDT protocols, leading to long irradiation times in the operating room, and increased power output without causing thermal heating is necessary. The required fluence rate or photodynamic threshold for PDT of occult micrometastatic disease is currently unknown. Previous clinical protocols have delivered high PDT radian exposures () to small target areas containing gross tumors. The feasibility and safety profile of whole-lung PDT using large-area light delivery devices placed in the thoracic cavity were demonstrated. Photosensitizers are usually given intravenously at an appropriate drug-light interval. Here, a modified lung flush procedure was employed to drain the organ of blood and allow improved light penetration, resulting in a 3.3-mm increased light penetration. Examining the cases administered 5-ALA, the maximally tolerated dose for normal lung tissue was of 5-ALA and a radiant exposure of . ALA-induced PPIX-mediated PDT demonstrated a dose-response relationship dependent on both, the photosensitizer concentration and delivered radiant exposure. The calculated threshold value for a normal porcine lung is , representing the upper limit of the therapeutic window, not to be exceeded throughout the clinical target volume. Regarding the Chlorin e6 cases, no PDT effects or photosensitizer uptake was demonstrated in the lung tissue at clinically used drug-light intervals, even after administering highly elevated radiant exposures. In an additional case to estimate the minimum threshold dose, Chlorin e6 was employed as a vascular-acting photosensitizer by delivering it via the lung flush solution while simultaneously illuminating select lung surface areas. As anticipated, this case demonstrated a vigorous PDT response with deep local injury to the lung and a minimum threshold value of lower than other PDT thresholds reported in the literature for cellular-acting photosensitizers. This protocol does not correspond to the usual drug-light interval used for Chlorin e6, which is aimed at targeting the tumor cellular compartment. The low threshold values could be explained as cellular and vascular PDT are performed. The clinical protocol involves a typical 3- to 4-h drug-light interval, based on the longer retention time of Chlorin e6 in tumor cells while clearing mostly from normal cells. Hence, the calculated PDT threshold values represent the upper therapeutic window limit for normal lung tissues. Previous work has demonstrated the selective uptake of Chlorin e6 up to nine-fold greater in tumor cells compared with normal cells, including sarcoma, small cell and non-small cell lung cancers, and nasopharyngeal carcinoma, in various small animal tumor and xenograft models.45–47 Similarly, chlorin e6 delivered by nano-emulsion has demonstrated significantly greater accumulation in lung metastases compared with normal lung tissue in a mouse tumor model.48 The significance of our results is that if it demonstrates similar selective uptake in tumor versus healthy lung tissues, Chlorin e6 could be an ideal photosensitizer for whole-lung PDT. Rapid clearing from normal lung tissue permits delivery of high irradiances and radiant exposure doses to maximize metastatic tumor destruction over an entire lung while avoiding normal tissue toxicity. However, exploiting this therapeutic window requires dosimetry that allows safe and accurate targeting of the lung while sparing other organs or structures at risk within the high fluence field. A light delivery system similar to that demonstrated here, in combination with the patient-specific simulations, can provide the required dosimetric control. Simulations based on CT images providing dose-volume histograms pre-procedure can guide the required number, area, position, and irradiances of sources to achieve the desired PDT treatment effect. This study has limitations that must be considered when extrapolating the results presented herein. From the light delivery system and simulation perspective, a porcine model was used with adolescent, healthy pigs, which may not accurately represent a human clinical setting. We previously reported the optical properties of human lungs;26 however, the validity of these properties for human lungs with different pathologies and baseline characteristics is unknown, and hence, the accuracy of the simulations may be limited. Regarding the potential selectivity and efficacy of our perfusion-assisted PDT protocol, this study did not include a tumor model. Therefore, we cannot confirm that the photosensitizers would effectively kill tumor cells and provide a targetable therapeutic index without significant toxicity. Selectivity is extrapolated from established threshold values for ALA-induced PpIX and Chlorin e6 for various cancers, which are lower than those we have reported here in normal lung tissue.19,49,50 Hence, in vivo studies with metastatic lung tumors are required to demonstrate the efficacy of this protocol. In addition, while this study has demonstrated the safety of this perfusion-assisted PDT protocol for ALA-induced PpIX and Chlorin e6 during a subacute timeframe, long-term safety was not examined. However, PDT-induced necrosis occurs within a 72-h timeframe and is not known to provide late side effects in other settings.13 Finally, using a large animal model, this study necessitated a low sample size, and thus, the results are susceptible to significant case-by-case variation. However, as a proof-of-concept study demonstrating a novel potential treatment protocol, determining exact safety parameters was not the objective of this work. 5.ConclusionThis study demonstrated a whole-lung perfusion-assisted PDT protocol. A light delivery system for improved PDT illumination homogeneity of the lung was presented. A simulation-based approach to precise PDT dose treatment planning was presented. The use of the light delivery system was demonstrated, identifying the maximally tolerated dose and feasibility of a perfusion-assisted whole-lung PDT protocol with 5-ALA and Chlorin e6. This protocol demonstrated therapeutic potential for both photosensitizers, with Chlorin e6, demonstrating minimal uptake in healthy lung tissue, with the possibility of sufficiently high therapeutic selectivity to target multifocal lung cancer. This work will inform further study confirming the efficacy of this protocol in tumor models and providing the basis for future clinical translation. Furthermore, it outlines an important approach to dosimetry in the lungs and other solid organs that can be used for personalized PDT protocols. DisclosuresThe authors have no relevant financial interests in the paper or other potential conflicts of interest to disclose. Code and Data AvailabilitySoftware code for Fullmonte is available at https://gitlab.com/FullMonte/FullMonteSW/-/wikis/home. AcknowledgmentsThe study was funded through grants from the Ontario Research Fund by the Ministry of Economic Development and Trade (Grant No. ORF 08-23) and the Ontario Ministry of Health and Long-Term Care and the Canadian Cancer Society. The study was further funded through the Toronto General and Western Hospital Foundation. Dr. Ramadan was supported by the Surgeon Scientist Training Program at the University of Toronto, the Thoracic Surgery Foundation Resident Research Fellowship, and the Hold’em For Life Oncology Clinician Scientist Award. Dr. Lilge was supported by the Princess Margaret Cancer Foundation. The authors acknowledge the gift of Chlorin e6 from Dr. Dirk Huettenberger and Synverdis GmbH. ReferencesU. Pastorino, P. McCormack and R. Ginsberg,
“A new staging proposal for pulmonary metastases: results of analysis of 5206 cases of resected pulmonary metastases,”
Chest Surg. Clin. North Amer., 8 197
(1998).
Google Scholar
K. G. Billingsley et al.,
“Pulmonary metastases from soft tissue sarcoma,”
Ann. Surg., 229
(5), 602 https://doi.org/10.1097/00000658-199905000-00002
(1999).
Google Scholar
E. Mitry et al.,
“Epidemiology, management and prognosis of colorectal cancer with lung metastases: a 30-year population-based study,”
Gut, 59
(10), 1383
–1388 https://doi.org/10.1136/gut.2010.211557 GUTTAK 0017-5749
(2010).
Google Scholar
C. M. Worsley, E. S. Mayne and R. B. Veale,
“Clone wars: the evolution of therapeutic resistance in cancer,”
Evol. Med. Public Health, 2016
(1), 180
–181 https://doi.org/10.1093/emph/eow015
(2016).
Google Scholar
M. H. Ihn et al.,
“Curative resection for metachronous pulmonary metastases from colorectal cancer: analysis of survival rates and prognostic factors,”
Cancer Res. Treat., 49
(1), 104
–115 https://doi.org/10.4143/crt.2015.367
(2017).
Google Scholar
S. Yokoyama et al.,
“Survival after initial lung metastasectomy for metastatic colorectal cancer in the modern chemotherapeutic era,”
BMC Surg., 17
(1), 54 https://doi.org/10.1186/s12893-017-0252-8
(2017).
Google Scholar
D. E. Dolmans, D. Fukumura and R. K. Jain,
“Photodynamic therapy for cancer,”
Nat. Rev. Cancer, 3
(5), 380
–387 https://doi.org/10.1038/nrc1071 NRCAC4 1474-175X
(2003).
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 https://doi.org/10.1016/S1572-1000(05)00007-4
(2004).
Google Scholar
Y. Hayata et al.,
“Hematoporphyrin derivative and laser photoradiation in the treatment of lung cancer,”
Chest, 81
(3), 269
–277 https://doi.org/10.1378/chest.81.3.269 CHETBF 0012-3692
(1982).
Google Scholar
N. Ikeda, J. Usuda and S. Maehara,
“Photodynamic therapy for central-type early-stage lung cancer,”
Gen. Thorac. Cardiovasc. Surg., 68
(7), 679
–683 https://doi.org/10.1007/s11748-019-01240-1
(2020).
Google Scholar
C. Lin et al.,
“Analysis of the short-term effect of photodynamic therapy on primary bronchial lung cancer,”
Lasers Med. Sci., 36
(4), 753
–761 https://doi.org/10.1007/s10103-020-03080-5
(2020).
Google Scholar
R. Jayadevappa et al.,
“Outcomes of patients with advanced non-small cell lung cancer and airway obstruction treated with photodynamic therapy and non-photodynamic therapy ablation modalities,”
J. Thorac. Dis., 11
(10), 4389
–4399 https://doi.org/10.21037/jtd.2019.04.60
(2019).
Google Scholar
C. A. Robertson, D. H. Evans and H. Abrahamse,
“Photodynamic therapy (PDT): a short review on cellular mechanisms and cancer research applications for PDT,”
J. Photochem. Photobiol. B: Biol., 96
(1), 1
–8 https://doi.org/10.1016/j.jphotobiol.2009.04.001 JPPBEG 1011-1344
(2009).
Google Scholar
G. Shafirstein et al.,
“Photodynamic therapy of non-small cell lung cancer: narrative review and future directions,”
Ann. Amer. Thorac. Soc., 13
(2), 265
–275 https://doi.org/10.1513/AnnalsATS.201509-650FR
(2015).
Google Scholar
B. C. Wilson and M. S. Patterson,
“The physics, biophysics and technology of photodynamic therapy,”
Phys. Med. Biol., 53
(9), R61
–R109 https://doi.org/10.1088/0031-9155/53/9/R01 PHMBA7 0031-9155
(2008).
Google Scholar
T. J. Farrell et al.,
“Comparison of the in vivo photodynamic threshold dose for photofrin, mono- and tetrasulfonated aluminum phthalocyanine using a rat liver model,”
Photochem. Photobiol., 68
(3), 394
–399 https://doi.org/10.1111/j.1751-1097.1998.tb09698.x PHCBAP 0031-8655
(1998).
Google Scholar
L. Lilge et al.,
“The sensitivity of normal brain and intracranially implanted VX2 tumour to interstitial photodynamic therapy,”
Br. J. Cancer, 73
(3), 332
–343 https://doi.org/10.1038/bjc.1996.58 BJCAAI 0007-0920
(1996).
Google Scholar
M. S. Patterson, B. C. Wilson and R. Graff,
“In vivo tests of the concept of photodynamic threshold dose in normal rat liver photosensitized by aluminum chlorosulphonated phthalocyanine,”
Photochem. Photobiol., 51
(3), 343
–349 https://doi.org/10.1111/j.1751-1097.1990.tb01720.x PHCBAP 0031-8655
(1990).
Google Scholar
L. D. Lilge et al.,
“Determination of the photodynamic threshold for normal rabbit brain and for intracranially implanted VX2 tumors,”
Proc. SPIE, 1882 60
–73 https://doi.org/10.1117/12.147692
(1993).
Google Scholar
M. R. Johnston et al.,
“Isolated lung perfusion with adriamycin. A preclinical study,”
Cancer, 52
(3), 404
–409 https://doi.org/10.1002/1097-0142(19830801)52:3<404::AID-CNCR2820520304>3.0.CO;2-J CANCAR 0008-543X
(1983).
Google Scholar
M. R. Johnston et al.,
“Isolated total lung perfusion as a means to deliver organ-specific chemotherapy: long-term studies in animals,”
Surgery, 98
(1), 35
–44 SURGAZ 0039-6060
(1985).
Google Scholar
M. R. Johnston, R. F. Minchen and C. A. Dawson,
“Lung perfusion with chemotherapy in patients with unresectable metastatic sarcoma to the lung or diffuse bronchioloalveolar carcinoma,”
J. Thorac. Cardiovasc. Surg., 110
(2), 368
–373 https://doi.org/10.1016/S0022-5223(95)70232-6 JTCSAQ 0022-5223
(1995).
Google Scholar
P. R. dos Santos et al.,
“Modified in vivo lung perfusion allows for prolonged perfusion without acute lung injury,”
J. Thorac. Cardiovasc. Surg., 147
(2), 774
–782 https://doi.org/10.1016/j.jtcvs.2013.10.009 JTCSAQ 0022-5223
(2014).
Google Scholar
P. Reck dos Santos et al.,
“Modified in vivo lung perfusion for local chemotherapy: a preclinical study with doxorubicin,”
Ann. Thorac. Surg., 101
(6), 2132
–2140 https://doi.org/10.1016/j.athoracsur.2015.12.043
(2016).
Google Scholar
K. Ramadan et al.,
“A model to assess acute and delayed lung toxicity of oxaliplatin during in vivo lung perfusion,”
J. Thorac. Cardiovasc. Surg., 161
(5), 1626
–1635 https://doi.org/10.1016/j.jtcvs.2020.03.033 JTCSAQ 0022-5223
(2021).
Google Scholar
K. T. Ramadan et al.,
“Determination of optical properties and photodynamic threshold of lung tissue for treatment planning of in vivo lung perfusion assisted photodynamic therapy,”
Photodiagn. Photodyn. Ther., 35 102353 https://doi.org/10.1016/j.pdpdt.2021.102353
(2021).
Google Scholar
R. A. Weersink et al.,
“Techniques for delivery and monitoring of TOOKAD (WST09)-mediated photodynamic therapy of the prostate: clinical experience and practicalities,”
J. Photochem. Photobiol. B: Biol., 79
(3), 211
–222 https://doi.org/10.1016/j.jphotobiol.2005.01.008 JPPBEG 1011-1344
(2005).
Google Scholar
J. P. Marijnissen and W. M. Star,
“Calibration of isotropic light dosimetry probes based on scattering bulbs in clear media,”
Phys. Med. Biol., 41
(7), 1191
–1208 https://doi.org/10.1088/0031-9155/41/7/008 PHMBA7 0031-9155
(1996).
Google Scholar
J. L. Sandell and T. C. Zhu,
“A review of in-vivo optical properties of human tissues and its impact on PDT,”
J. Biophotonics, 4
(11–12), 773
–787 https://doi.org/10.1002/jbio.201100062
(2011).
Google Scholar
A. N. Bashkatov, E. A. Genina and V. V. Tuchin,
“Optical properties of skin, subcutaneous, and muscle tissues: a review,”
J. Innov. Opt. Health Sci., 4
(1), 9
–38 https://doi.org/10.1142/S1793545811001319
(2011).
Google Scholar
J. Cassidy et al.,
“High-performance, robustly verified Monte Carlo simulation with FullMonte,”
J. Biomed. Opt., 23
(8), 085001 https://doi.org/10.1117/1.JBO.23.8.085001 JBOPFO 1083-3668
(2018).
Google Scholar
B. C. Wilson and G. Adam,
“A Monte Carlo model for the absorption and flux distributions of light in tissue,”
Med. Phys., 10
(6), 824
–830 https://doi.org/10.1118/1.595361 MPHYA6 0094-2405
(1983).
Google Scholar
L. Wang, S. Ren and X. Chen,
“Comparative evaluations of the Monte Carlo-based light propagation simulation packages for optical imaging,”
J. Innov. Opt. Health Sci., 11
(1), 1750017 https://doi.org/10.1142/S1793545817500171
(2017).
Google Scholar
C. McFadden et al.,
“Transbronchial light illumination for peripheral lung cancer: a numerical feasibility study,”
in Clin. and Preclin. Opt. Diagn. II,
11079_21
(2019). https://doi.org/10.1117/12.2527570 Google Scholar
C. Le Tourneau, J. J. Lee and L. L. Siu,
“Dose escalation methods in phase I cancer clinical trials,”
J. Natl. Cancer Inst., 101
(10), 708
–720 https://doi.org/10.1093/jnci/djp079 JNCIEQ
(2009).
Google Scholar
X. Li et al.,
“Quantification of a novel photosensitizer of Chlorin e6-C15-monomethyl ester in beagle dog plasma using HPLC: application to pharmacokinetic studies,”
Molecules, 22
(5), 693 https://doi.org/10.3390/molecules22050693
(2017).
Google Scholar
H. S. Ginsberg et al.,
“Role of early genes in pathogenesis of adenovirus pneumonia,”
Proc. Natl. Acad. Sci. U. S. A., 87
(16), 6191
–6195 https://doi.org/10.1073/pnas.87.16.6191
(1990).
Google Scholar
L. Lilge, C. O’Carroll and B. C. Wilson,
“A solubilization technique for photosensitizer quantification in ex vivo tissue samples,”
J. Photochem. Photobiol. B: Biol., 39
(3), 229
–235 https://doi.org/10.1016/S1011-1344(97)00010-9 JPPBEG 1011-1344
(1997).
Google Scholar
C. McFadden,
“Photodynamic Therapy Dose Planning for Peripheral Lung Cancer,”
https://hdl.handle.net/1807/103179
().
Google Scholar
J. M. Park et al.,
“Evaluation of treatment plans using various treatment techniques for the radiotherapy of cutaneous Kaposi’s sarcoma developed on the skin of feet,”
J. Appl. Clin. Med. Phys., 15
(6), 173
–187 https://doi.org/10.1120/jacmp.v15i6.4970
(2014).
Google Scholar
S. Mordon et al.,
“Light emitting fabric technologies for photodynamic therapy,”
Photodiagn. Photodyn. Ther., 12
(1), 1
–8 https://doi.org/10.1016/j.pdpdt.2014.11.002
(2015).
Google Scholar
A. George and P. S. Shrivastav,
“Photodynamic therapy with light emitting fabrics: a review,”
Arch. Dermatol. Res., 314
(10), 929
–936 https://doi.org/10.1007/s00403-021-02301-3
(2022).
Google Scholar
S. Chamberlain et al.,
“An optical surface applicator for intraoperative photodynamic therapy,”
Lasers Surg. Med., 52
(6), 523
–529 https://doi.org/10.1002/lsm.23168 LSMEDI 0196-8092
(2020).
Google Scholar
B. Ramaswamy et al.,
“Photodynamic diagnosis of a human nasopharyngeal carcinoma xenograft model using the novel Chlorin e6 photosensitizer Fotolon,”
Int. J. Oncol., 26
(6), 1501
–1506 IJONES 1019-6439
(2005).
Google Scholar
W. W. Chin, P. W. Heng and M. Olivo,
“Chlorin e6 – polyvinylpyrrolidone mediated photosensitization is effective against human non-small cell lung carcinoma compared to small cell lung carcinoma xenografts,”
BMC Pharmacol., 7 15 https://doi.org/10.1186/1471-2210-7-15
(2007).
Google Scholar
S. Liao et al.,
“Antitumor effect of photodynamic therapy/sonodynamic therapy/sono-photodynamic therapy of Chlorin e6 and other applications,”
Mol. Pharm., 20
(2), 875
–885 https://doi.org/10.1021/acs.molpharmaceut.2c00824
(2023).
Google Scholar
C. Park et al.,
“Chlorin e6-Loaded PEG-PCL nanoemulsion for photodynamic therapy and in vivo drug delivery,”
Int. J. Mol. Sci., 20
(16), 3958 https://doi.org/10.3390/ijms20163958 1422-0067
(2019).
Google Scholar
J. van den Boogert et al.,
“5-Aminolaevulinic acid-induced protoporphyrin IX accumulation in tissues: pharmacokinetics after oral or intravenous administration,”
J. Photochem. Photobiol. B: Biol., 44
(1), 29
–38 https://doi.org/10.1016/S1011-1344(98)00102-X JPPBEG 1011-1344
(1998).
Google Scholar
R. C. M. C. Ferraz et al.,
“Determination of threshold dose of photodynamic therapy to measure superficial necrosis,”
Photomed. Laser Surg., 27
(1), 93
–99 https://doi.org/10.1089/pho.2007.2207
(2009).
Google Scholar
BiographyKhaled Ramadan is a general surgeon, currently performing additional training in thoracic surgery at the University of Toronto. He received his BHSc from McMaster University in 2011 and his MD and PhD from the University of Toronto in 2015 and 2023, respectively. His research interests include lung perfusion techniques in thoracic oncology and PDT applications in thoracic surgery. Tina Saeidi is currently a 4th year PhD student in the Department of Medical Biophysics at the University of Toronto. Tina received her MSc degree in laser and photonics from Ruhr-Universität Bochum, Germany, in 2020. Her current project focuses on the optimization of PDT treatment planning by predicting the spatial distribution of the photosensitizer. Vanderlei Bagnato received his PhD from the Massachusetts Institute of Technology in 1987 and his double major in material science engineering from the Universidade Federal de São Carlos, Brasil, in 1981, and in physics from the Universidade de São Paulo, Brasil, in 1981. He received the title “Livre–Docente” from the University of São Paulo in 1989. He is a full-time professor at the University of São Paulo in 1993. Recently, he also joined the Texas A&M University in biomedical engineering. His activities are based on laser cooling, trapping of neutral atoms, and applications of optics and laser to health science. He has published over 800 papers in international journals with more than 25,000 citations. During his academic life, he has supervised more than 130 graduate students. He is a member of the following academies: National Academy of Science NAS (USA), Pontifical Academy of Science (Vaticano), TWAS, Brazilian Academy of Science (ABC), and Latin America Academy of Science (ACAL). Lothar Lilge is a professor in medical biophysics at the University of Toronto. He received his MSc in 1987 from the University of Frankfurt, Germany, and his PhD from the University of Muenster, Germany, in 1992. He joined the Faculty at the University of Toronto in 1997 and currently holds a full professorship. He is working on optical spectroscopy for breast cancer risk assessment and PDT. The latter includes the assessment of novel photosensitizers, PDT dose monitoring techniques, and treatment planning. |