Salivary Gland Tissue Chip
Our previously reported salivary gland tissue chip17 was leveraged for high-content drug screening to identify novel radioprotective compounds. The chip platform consists seeding primary salivary gland cell clusters 20-100 µm (Fig. 1A), suspended in a poly(ethylene glycol) (PEG) hydrogel precursor and MMP-degradable peptide crosslinker solution together with the photoinitiator LAP (Fig. 1B), into an array of near-spherical microbubble (MB) cavities (Fig. 1C) formed in poly(dimethylsiloxane) (PDMS). Each chip, containing 40-50 MBs, is affixed within wells of a 96-well plate. In situ polymerization of the hydrogel was achieved using long-wave, low intensity UV light. Over time, the cell clusters aggregate and proliferate to form SGm (Fig. 1D).
High-throughput methods to assess drug radioprotection
In prior studies17, immunohistochemical (IHC) staining was used to quantify the number of γH2AX and 53BP1 puncta within nuclei, which are sensitive markers for double-stranded DNA breaks50 and direct measures of radiation damage. IHC staining is, however, a laborious time and resource-consuming process requiring retrieval of the SGm from the MB array chip, tissue sectioning, staining, and imaging. While IHC enables sensitive analyses of radiation damage, it completely abrogates the goal of in situ high-content screening for which the tissue chip was developed. As mentioned above, several assays commonly used to assess radiation-induced cellular damage were tested in the tissue chip format at various time points post-radiation (Table S.1). Based on signal-to-noise ratio and reproducibility, the glutathione 18 and senescence19 assays were selected for further development. The glutathione assay was used to measure glutathione levels at various time points post-radiation to determine the optimal time for detecting differences between 0 Gy and 15 Gy. Based on the data, the greatest signal separation was measured at 4 days post-irradiation (Fig. S1). This time point is similar to previous reports on decreases in glutathione post-radiation18 and was used for all experiments moving forward.
Next, WR-1065 was tested to identify effective concentration(s) that prevented changes in glutathione levels post-radiation. A SGm tissue chip was cultured for 4 days to allow spheres to form17 . The chip was then treated with WR-1065 30 min before and during radiation, followed by drug wash out with media 30 min post-radiation (Fig. 2A). This dosing scheme is consistent with the use of Amifostine clinically14 and with our previous work17. The glutathione assay was performed 4 days post-radiation (Fig. 2A). Example images show high levels of glutathione at 0 Gy (Fig. 2B) that is decreased by 15 Gy radiation (Fig. 2C) and maintained with 4 mM WR-1065 (Fig. 2D). Quantification shows that 0.1 mM and 0.4 mM WR-1065 were ineffective at preventing radiation damage. In contrast, 1 mM and 4 mM WR-1065 provided significant protection (Fig. 2E). The 4 mM dose corresponds with our previous work on DNA damage markers γH2AX and 53BP117 and values from literature51 and establishes the range of effective concentrations for WR-1065 treatment in vitro. Moreover, it is clinically relevant as 15-30 minutes before radiation, patients are administered Amifostine intravenously at 200 mg/m2 15, 40. Assuming an average adult body surface area of 17,000 cm2 and a blood volume of 1.35 gallons (5.1 liters)52, 53, Amifostine is administered at 300 µM.
Since WR-1065 is an antioxidant and mediates radioprotective effects through free radical scavenging and induction of superoxide dismutase expression54, we tested other antioxidants implicated as radioprotective (Tempol, Edaravone, N-acetylcysteine). Tempol exhibited excellent radioprotection at both 1 and 4 mM (Fig. 2F) which is consistent with prior studies20, 49, 55. Edaravone showed complete radioprotection at 1 mM but not at 4 mM (Fig. 2G). Edaravone maintained ~43% of glutathione levels at 0.1 mM (Fig. S2A), suggesting that the optimal concentration range for Edaravone might be lower than WR-1065, which is consistent with literature indicating radio-protective dose ranges are 0.1-1 mM23, 56, 57. For N-acetylcysteine (NAC), no improvement in glutathione levels was observed at 1 or 4 mM compared to untreated SGm (Fig. 2H); however,glutathione levels were rescued by 56% when treated with 10 mM NAC (Fig. S2B), also consistent with literature22.
Drugs reported with non-antioxidant radioprotective mechanisms were also tested using the glutathione assay. Rapamycin is an mTOR inhibitor reported to restore salivary flow rate post-irradiation in swine24. Ex-Rad reduces p53-dependent and independent apoptosis25. Palifermin (keratinocyte growth factor) has been reported to stimulate salivary gland stem/progenitor cell expansion post-radiation27. Using drug concentrations based on literature26, 28, our data shows some protection against glutathione level changes resulting from radiation with 50 µM rapamycin (58%, Fig. 2I) and 100 ng/mL Palifermin (47%, Fig. 2J), but no protection from 50 µM Ex-Rad (Fig. 2K).
For senescence, a protocol similar to the glutathione assay was followed to determine the optimal time point post-radiation for detecting a change in senescence between 0 Gy and 15 Gy, as measured by senescence-associated β-galactosidase activity. A time point of 5 days post-radiation was optimal (Fig. S3), similar to a previous study19. WR-1065 radioprotection was tested by adding drug to the chips 30 min before radiation followed by wash out 30 min post-radiation (Fig. 3A). An expected increase in senescence was detected for SGm exposed to 15 Gy compared to 0 Gy (Figs. 3B,C) , which was restored to levels equivalent to 0 Gy with the addition of WR-1065 (Fig. 3D). Quantification shows that both 1 mM and 4 mM WR-1065 treatment resulted in complete radioprotection (Fig. 3E), consistent with the glutathione assay results.
Next, the level of senescence was measured using the other reported radioprotective drugs. Tempol (Fig. 3F) and Edaravone (Fig. 3G) reduced senescence at both 1 mM and 4 mM by 75% and 91% for Tempol, 94% and 113% for Edaravone, at 1 and 4 mM, respectively, versus untreated, irradiated controls. Edaravone also reduced senescence by 96% at 0.1 mM versus untreated, irradiated controls (Fig. S4A). For NAC, 83% and 57% reduction in senescence was observed with 1 mM and 4 mM (Fig. 3H) and 82% at 10 mM (Fig. S4B) versus untreated, irradiated controls. For drugs with non-antioxidant mechanisms, rapamycin showed complete protection (Fig. 3I), whereas Palifermin provided 64% protection (Fig. 3J) and Ex-Rad conferred only 45% radioprotection (Fig. 3K).
A summary of results from the glutathione and senescence assays shows similar trends for radioprotection (Table S.5). The few differences may result from different mechanisms of action and/or assay targets. Logically, the glutathione assay may be more sensitive to antioxidant function (NAC, Tempol, Edaravone), while the senescence assay may be more appropriate for drugs such as rapamycin, which has anti-senescence properties58. These results highlight the trade-offs in developing screening assays and point to the benefit of screening with two assays. Although the assays developed are indirect measures of radiation-induced DNA damage, they nonetheless were validated to detect radiation induced cell-damage and drug radioprotection. More importantly, these assays can be used for in situ high-content drug screening with multiple replicates (40-50) per test and enhanced throughput compared to IHC staining for γH2AX.
Drug library screening identified several promising radioprotective compounds
The glutathione and senescence assays were used to screen a library of FDA-approved drugs (Selleck Chemicals) at 100 µM. Drugs were first screened using the glutathione assay according to the timeline in Fig. 2A. Any compound resulting in statistically equivalent glutathione levels compared to the 0 Gy control was considered a hit (Fig. 4, orange circles). Hits with the glutathione assay were then tested uisng the senescence assay and considered a “double hit” if senescence levels were statistically equivalent to levels at 0 Gy (Fig. 4, blue circles). A list of the 438 drugs screened and relevant statistics are shown in Table S.3. Overall, 438 drugs from the library were tested with a hit rate of 5.7%, for a total of 25 double hits (Fig. 4B) listed in Table 1. While this hit rate is higher than many other drug screening reports (0.1 – 0.3%)59, 60, this may be due to the high statistical rigor afforded by the tissue chip format. Additionally, phenotypic screens generally have higher hit rates than target-based screens and they maybe a more successful strategy for discovery of first-in-class medicines (> 1%)61-63.
Identification of potential radioprotective drug mechanisms
Of the 25 potential radioprotective compounds, 20 have known interactions with proteins involved in calcium signaling identified within the BioAssay database in PubChem (Fig. 5A). These compounds may impact secretory signaling in the salivary gland, which can be radioprotective64, 65. While degranulation may not be key to radioprotection, secretory stimulation may play a role in proliferation and survival of the secretory cells64. Similarly, using the Drug Set Enrichment Analysis (DSEA) tool to identify pathways, the Reactome analysis related to secretion appear to be upregulated by many of the identified drugs, supporting this potential mechanism (Fig. 5B,C, and Fig. S5). Interestingly, only 9 of the 25 compounds have known antioxidant properties, and 12 are anti-inflammatory. This is critical data indicating that alternate mechanisms of radioprotection may be achievable and represented in the identified hits. A reduction in pathway activity related to cell adhesion, cell-cell and cell-matrix interactions also represents a potential area of exploration. Multiple studies have identified changes in integrin expression, cell adhesion and matrix interactions upon radiation exposure, which may be linked to cell response1, 66-69. Manipulation of these mechanisms in the salivary gland may convey radioprotection.
Hit down-selection using drug promiscuity data and EC50 values
A systematic approach was used to down selection the 25 double hits for in vivo testing. Since drugs within the library are FDA-approved, considerable information on their pharmacology in mice and humans is readily available through resources such as PubChem. Within PubChem, the BioAssay database was created by the National Institute of Health (NIH) as an open repository containing results of small molecule screening data38. We used the BioAssay data to analyze drug promiscuity, which refers to the ability of a drug to bind multiple molecular targets with distinct pharmacological outcomes, often causing unwanted side effects108. The drugs exhibiting bioactivity in a large number of assays were deprioritized. Data for each double hit was obtained from the database and promiscuity was calculated as the percent of assays reported as “active” (Table 1). Drugs with high activity (>10%) were excluded from further testing. Additionally, Etidronate, Melatonin, and Albendazole were excluded due to poor bioavailability109-111, and Eplerenone was excluded due to solubility concerns112.
The remaining 13 drugs were tested for glutathione levels post-irradiation in dose-limiting experiments (1 -100 µM) to identify effective concentrations. Many drugs were only effective at 100 µM and were excluded due to a lack of potency. Forthe remaining seven drugs (Phenylbutazone, Meropenem, Diethylstilbestrol, Prazosin, Enoxacin, Glipizide, and Doripenem) dose responses and heat maps summarizing the dose-response results for the glutathione and senescence assays are presented in Fig. 6 and Fig. 7, respectively. Dose-response curves for the glutathione assay show that Phenylbutazone (Fig. 6B) and Meropenem (Fig. 6C) exhibited radioprotection equivalent to 0 Gy over a concentration range 0.1 – 100 µM and Diethylstilbestrol (Fig. 6D) was radioprotective at 10-100 µM. Prazosin (Fig. 6E), Enoxacin (Fig. 6F), Glipizide (Fig. 6G), and Doripenem (Fig. 6H) showed protection only at higher concentrations.
The radioprotection trends based on the senescence assay (Fig. 7) differed somewhat, with Phenylbutazone (Fig. 7B) and Meropenem (Fig. 7C) showing only partial protection and Diethylstilbestrol (Fig. 7D) exhibiting radioprotection equivalent to 0 Gy between 50-100 µM concentrations whereas Prazosin (Fig. 7E) showed complete protection between 0.1-100 µM. Glipizide (Fig. 7G) and Doripenem (Fig. 7H) also showed variable protection.
EC50 values extrapolated from dose response curves are shown in Table 2. Phenylbutazone showed the most promising results, with low EC50 values for both the glutathione (0.08 µM) and senescence (0.05 µM) assays. Phenylbutazone is a non-steroidal anti-inflammatory drug (NSAID) that inhibits cyclooxygenases (COX-1 and COX-2), enzymes that produce prostaglandins113. Prostaglandins, specifically PGE2 signaling, have been shown to increase in irradiated salivary glands, and mitigation of salivary gland damage was achieved through treatment with the anti-inflammatory drug Indomethacin1, 114. Indomethacin also showed radioprotection in our drug screen but was ineffective at concentrations lower than 100 µM. Indomethacin only blocks COX-1, underscoring the greater efficacy of Phenylbutazone. Phenylbutazone was originally developed for chronic pain for conditions such as arthritis but has since been restricted to treating ankylosing spondylitis due to induction of rare but severe blood disorders, including anemia and leukopenia113. However, doses ranged from 300-1000 mg, generating a plasma concentration of 30-50 µg/mL113. In contrast, the EC50 of 0.08 µM for protecting against radiation-induced glutathione changes established in this study equates to a 26 ng/mL concentration. Thus, the risk of severe adverse effects may be greatly diminished for doses necessary for radioprotection. Additionally, Phenylbutazone has excellent bioavailability (up to 90%) 113, 115 and long half-life (50-105 hrs)113, which may enable dose de-escalation, further decreasing risks.
Enoxacin is an antibacterial agent used for treating urinary tract infections116 that was previously identified as radioprotective117. Using a high-throughput screening method with the viability of lymphocytes as the primary readout, two classes of antibiotics (tetracyclines and fluoroquinolones) were identified as robust radioprotectors, including Enoxacin117. This observation corroborates our drug screening results, in which several antibiotics were identified as double hits, including Enoxacin, Meropenem, Doripenem hydrate, Rifampin, and Rifabutin. The Enoxacin EC50 of 2.4 µM for the glutathione assay is similar to the 13 µM EC50 reported for viability of lymphocyte cells117. Notably, five other fluoroquinolones reported as radioprotectors (Levofloxacin, Gatifloxacin, Ofloxacin, Moxifloxacin, and Norfloxacin)117 were not hits in our drug screen (Table S3). These disparities may be related to differences in cell type (salivary gland versus lymphocyte) or readouts (glutathione/senescence versus viability).
Based upon drug down selection data and measured EC50 values, Phenylbutazone, Enoxacin, and Doripenem hydrate were selected for in vivo validation in mice with γH2AX foci per nucleus IHC staining as an outcome measure consistent with prior work showing correlation with the development of xerostomia17, 41, 42, 46. Vehicle-treated SGm exposed to 15 Gy radiation exhibited a 3.3-fold increase in the number of γH2AX foci per nucleus, indicating a significant increase in double-stranded DNA breaks due to radiation exposure (Figs. 8A-D, I). Treatment with WR-1065 resulted in a 0.5-fold reduction in γH2AX foci per nucleus compared to 15 Gy controls (Figs. 8D, E, I). No significant differences existed between 0 Gy controls and WR-1065 treated SGm exposed to 15 Gy (Figs. 8C, E, I). These results are similar to prior studies utilizing WR-1065 in vitro and in vivo via retrograde ductal injection17, 42 . Treatment with the test compounds, Phenylbutazone, and Enoxacin resulted in 0.4- and 0.5-fold reduction in γH2AX foci per nucleus compared to 15 Gy controls, respectively (Figs. 8E-G, I). Results observed after Phenylbutazone and Enoxacin treatment were not significantly different from 0 Gy controls (Figs. 8C, F, G, I). Treatment with Doripenem hydrate did not reduce γH2AX foci per nucleus relative to 15 Gy controls and showed a 2.7-fold increase compared to 0 Gy controls (Figs. 8C, D, H, I).