Abstract
Small-molecule inhibitors of PARP are thought to mediate their antitumor effects as catalytic inhibitors that block repair of DNA single-strand breaks (SSB). However, the mechanism of action of PARP inhibitors with regard to their effects in cancer cells is not fully understood. In this study, we show that PARP inhibitors trap the PARP1 and PARP2 enzymes at damaged DNA. Trapped PARP–DNA complexes were more cytotoxic than unrepaired SSBs caused by PARP inactivation, arguing that PARP inhibitors act in part as poisons that trap PARP enzyme on DNA. Moreover, the potency in trapping PARP differed markedly among inhibitors with niraparib (MK-4827) > olaparib (AZD-2281) >> veliparib (ABT-888), a pattern not correlated with the catalytic inhibitory properties for each drug. We also analyzed repair pathways for PARP–DNA complexes using 30 genetically altered avian DT40 cell lines with preestablished deletions in specific DNA repair genes. This analysis revealed that, in addition to homologous recombination, postreplication repair, the Fanconi anemia pathway, polymerase β, and FEN1 are critical for repairing trapped PARP–DNA complexes. In summary, our study provides a new mechanistic foundation for the rational application of PARP inhibitors in cancer therapy. Cancer Res; 72(21); 5588–99. ©2012 AACR.
Introduction
PARP1 is an abundant nuclear protein and the founding member of the PARP family (1–4). It binds damaged DNA through its N-terminal zinc finger motifs, which activates its catalytic C-terminal domain to hydrolyze NAD+ and produce linear and branched PAR chains that can extend over hundreds of ADP-ribose units (1–4; see Fig. 7B). PARP2, although less abundant and contributing 5% to 10% of the total PARP activity (1, 5), is essential for the viability of PARP1 knockout mice (6). The rapid binding of PARP1 to DNA strand breaks is critical for the resealing of DNA single-strand breaks (SSB) during base excision repair (BER; refs. 7–9) and for the repair of topoisomerase I cleavage complexes (10) and DNA double-strand breaks (DSB; refs. 11, 12). AutoPARylation of PARP1 and PARP2, and PARylation of chromatin proteins promote the recruitment of DNA repair factors (1, 13, 14). Because of the high negative charge of PAR polymers, extensive autoPARylation of PARP1 and PARP2 leads to their dissociation from DNA, which is required for DNA repair completion (15).
A large number of SSBs are generated endogenously as well as BER intermediates (16). When SSBs are encountered by replication forks, they generate DSBs that need to be repaired by homologous recombination (HR; refs. 17, 18). Accordingly, PARP inhibition results in the accumulation of recombinogenic substrates marked by RAD51 and γ-H2AX nuclear foci (19, 20). This explains why cancer cells deficient in BRCA1 and BRCA2 are selectively hypersensitive to PARP inhibitors (21–23). Recent studies have also suggested that disruptions of Fanconi anemia (FA)- or other HR-related genes also sensitize cells to PARP inhibitors (24, 25), providing a strong rationale for the development of PARP inhibitors for cancer treatment (4, 26, 27).
Olaparib (AZD2281), veliparib (ABT-888), and niraparib (MK-4827) are PARP inhibitors in advance clinical trials. They are highly efficacious PARP catalytic inhibitors with IC50 values reaching the low nanomolar range (4, 27–29). Their structure includes a nicotinamide moiety that competes with NAD+ (Supplementary Fig. S1 and see Fig. 7B). Assuming that the synthetic lethality of PARP inhibitors in the context of HR or BRCA deficiencies is solely explained by catalytic inhibition, the drugs should have a similar effect as PARP deletion. However, wild-type cells are more sensitive to a PARP inhibitor combined with the alkylating agent methylmethane sulfonate (MMS) than Parp1−/− mouse cells (30–32). Furthermore, PARP inhibition delays SSB repair to a greater extent than PARP depletion (33). To explain these results, a PARP1-trapping model has been proposed (18, 32). This model is based on the idea that PARP1 is trapped on DNA by PARP inhibitors, and PARP1–DNA complexes can interfere with DNA replication. Recently, increased PARP1 association with DNA in alkylation damage by a PARP inhibitor was reported using chromatin immunoprecipitation. However, accumulation of PARP1 was not detected in chromatin-bound fraction at the protein level (32).
Our study shows that, in addition to catalytic inhibition, PARP inhibitors induce cytotoxic PARP–DNA complexes, and we provide methods to detect such complexes. Secondly, we show that clinically relevant PARP inhibitors differ markedly in their potency to trap PARP–DNA complexes. Finally, we reveal previously unknown repair factors/pathways for the PARP–DNA complexes, and we discuss their importance in the context of personalized medicine.
Materials and Methods
Cell lines
The DT40 cell lines used in this study were obtained from the Laboratory of Radiation Genetics Graduate School of Medicine in Kyoto University (Kyoto, Japan). The summary of the DT40 mutant cell lines is described in Supplementary Table S1. All human cancer cell lines were obtained from the National Cancer Institute Developmental Therapeutics Program (Frederick, MD).
Immunoblotting
To prepare whole cell lysates, cells were lysed with CelLytic M Lysis Reagent (C2978, Sigma-Aldrich). After thorough mixing and incubation at 4°C for 30 minutes, lysate were centrifuged at 15,000 × g at 4°C for 10 minutes, and supernatants were collected. To prepare subcellular fraction of nuclear soluble and chromatin-bound fraction, 10 million DT40 cells with 10 mL medium in 15 mL tube or semiconfluent human cells with 10 mL medium in 10 cm dish were treated with indicated drugs for 30 minutes or 4 hours, respectively, and then cells were collected. For the fractionation, we used a Subcellular Protein Fractionation Kit from Thermo Scientific (78840) following the manufacturer's instructions. Immunoblotting was carried out using standard procedures. Densitometric analyses of immunoblots were carried out using Image J software (NIH). Further explanation of fractionation methods is available in the Supplementary Experimental Procedures.
Antibodies
Rabbit polyclonal anti-PARP1 antibody (sc-7150), mouse monoclonal anti–proliferating cell nuclear antigen (PCNA) antibody (sc-56), and mouse monoclonal anti-FANCD2 antibody (sc-20022) were from Santa Cruz Biotechnology. Mouse monoclonal anti–γ-H2AX antibody (05-636) and rabbit polyclonal anti-histone H3 antibody (07-690) were from Upstate Biotechnology. Mouse monoclonal anti-actin antibody (A3853) and rabbit anti–γ-tubulin antibody (T3559) were from Sigma. Mouse monoclonal anti-Top1 antibody (#556597) was from BD Biosciences. Rabbit polyclonal anti-PAR polymer antibody (#4336-BPC-100) was from Trevigen. Rabbit polyclonal anti-PARP2 antibody (ab93416) was from Abcam. Secondary antibodies were horseradish peroxidase–conjugated antibodies to mouse or rabbit IgG (GE Healthcare).
Immunoassay for PAR substrate
The validated chemiluminescent immunoassay for PAR using commercially available reagents is described previously (34). The detailed laboratory procedures can be viewed from NCI website at http://dctd.cancer.gov/ResearchResources/biomarkers/PolyAdenosylRibose.htm
Measurement of cellular sensitivity to DNA-damaging agents
To measure the sensitivity of cells to drugs, cells were continuously exposed to various concentrations of the drugs for 72 hours in triplicate. For DT40 cells, 200 cells were seeded into 384-well white plate (#6007680 PerkinElmer Life Sciences) in 40 μL of medium per well. For human cells, 1,500 cells of DU145 or 500 cells of SF295 were seeded in 96-well white plate (#6005680 PerkinElmer Life Sciences) in 100 μL of medium per well. Cell survival was determined using the ATPlite 1-step kit (PerkinElmer). Briefly, 20 μL or 50 μL ATPlite solution for 384-well or 96-well plate, respectively, was added to each well. After 5 minutes, luminescence was measured by EnVision 2104 Multilabel Reader (PerkinElmer). The ATP level in untreated cells was defined as 100%. Percent viability of treated cells was defined as treated cells/untreated cells × 100.
Evaluation of the relative cellular sensitivity across a panel of DT40 mutant cells
The cellular sensitivity assay using a 384-well plate allowed us to examine 6 different cell lines at the same time, while we always included wild-type cells in each plate for normalization. To evaluate the relative cellular sensitivity of each mutant to wild-type cells, we drew the single logarithm sensitivity curve by setting the survival of untreated cells as 100%. The concentration that causes 90% reduction of ATP activity (inhibition concentration 90%; IC90) by olaparib or niraparib, and the concentration that causes 50% reduction (IC50) by veliparib, were determined manually from the sensitivity curves. The IC90 and IC50 of each mutant were divided by the IC90 and IC50 of wild-type cells on the same plate, respectively. Finally, the quotient was converted into logarithmic scale (base 2).
siRNA transfection
Gene-specific siRNAs (mix of 4 sequences) for PARP1 (L-006656-00-0005), PARP2 (L-010127-02-0005), and negative control siRNA (D-001810-10) were products of Dharmacon. Ten nanomolar of each siRNA was transfected to DU145 and SF295 cells with Lipofectamine RNAiMAX Reagent (13778, Invitorogen) according to the manufacturer's instructions. Culture medium was changed 6 to 8 hours after the transfection. The transfected cells were trypsinized and seeded on a 96-well plate 48 hours after the transfection. Drugs were added 12 hours after the seeding, and cells were incubated for another 72 hours before ATP activity was measured. To confirm the downregulation of each protein, the transfected cells were collected for the Western blot analysis 72 hours after the transfection.
Drug preparation
Olaparib, veliparib, and niraparib were obtained from the Drug Synthesis and Chemistry Branch, National Cancer Institute (Bethesda, MD). 4-Amino-1,8-naphthalimide (4-AN) was from Sigma-Aldrich (A0966). Drug stock solutions were made in dimethyl sulfoxide (DMSO) at 10 μmol/L. The stock solutions were stored at −20°C in dark and diluted in the culture medium immediately before use. One percent or 10% MMS was prepared fresh each time from 99% MMS (129925, Sigma-Aldrich) in PBS, and then diluted in culture medium to final concentration.
Fluorescence anisotropy DNA-binding assay
The fluorescence anisotropy experiment was carried out with a 30 bp duplex labeled with 5′-Alexa Fluor488. The deoxyoligonucleotide (sequence: 5′-Alexa Fluor 488-ACCCTGCTGTGGGCdUGGAGAACAAGGTGAT) was annealed to its complementary DNA strand in buffer containing 50 mmol/L potassium acetate, 20 mmol/L tris-acetate, 10 mmol/L magnesium acetate, and 1 mmol/L dithiothreitol, pH 7.9. All oligonucleotides were purchased from Integrated DNA Technologies. Uracil-DNA glycosylase and APE1 (New England Biolabs) were added to the annealed DNA sample and incubated at 37°C for 1 hour. The resulting DNA construct contains a DNA nick and a 5′-deoxyribose-5-phosphate (5′-dRP) at the nicked site. The completion of digestion is verified by denaturing PAGE. For fluorescence anisotropy measurements, 250 nmol/L recombinant PARP1 (a kind gift from Dr. Valerie Schreiber, University of Strasbourg, France), 1 nmol/L DNA construct, and increasing concentration of PARP inhibitors were combined and incubated for 30 minutes in buffer containing 50 mmol/L Tris–HCl (pH 8.0), 4 mmol/L MgCl2, 100 mmol/L NaCl, and 100 ng/μL BSA. A total of 1 mM NAD+ was added to the samples to initiate the experiment, and the fluorescence anisotropy values were measured at indicated time using an EnVision 2104 Multilabel Reader equipped with a FITC FP Label (PerkinElmer). The control samples lack NAD+ or PARP inhibitor. The fluorescence anisotropy values reported here were average of 2 independent experiments (each experiment was measured in duplicates).
Results
PARP inhibition by olaparib is more cytotoxic than genetic depletion of PARP
To compare the effect of PARP deletion and PARP inhibition, we took advantage of the fact that avian cells genetically lack PARP2 (35). Therefore, PARP1−/− avian B-lymphoblast DT40 cells are equivalent to PARP1 and PARP2 double knockout (35). Moreover, the catalytic domains and DNA-binding zinc finger motifs are highly conserved between human and chicken PARP1 (Supplementary Fig. S2). Because PARP1 and PAR levels were undetectable in PARP1−/− cells (Fig. 1A–C), and the PAR level of wild-type cells treated with 0.1 μmol/L olaparib was as low as in PARP1−/− cells (Fig. 1C and 4B), we can compare PARP deletion and catalytic PARP inhibition using PARP1−/− cells and wild-type DT40 cells treated with olaparib. We tested the sensitivity of wild-type, PARP1−/−, and BRCA2tr/− (BRCA2 truncated mutant; ref. 36) DT40 cells to olaparib (Fig. 1D). As expected, BRCA2tr/− cells showed greater sensitivity than wild-type. PARP1−/− cells were highly resistant to olaparib, showing the drug selectivity for PARP1. Cell-cycle analyses revealed that olaparib had no impact on the cell cycle of PARP1−/− cells, whereas it induced G2 accumulation in wild-type and massive apoptosis in BRCA2tr/− cells (Fig. 1E). γ-H2AX levels were also significantly increased by olaparib in wild-type cells but not in PARP1−/− cells (Fig. 1F). These results show that PARP1 is required for the cytotoxicity of olaparib and that PARP inhibition is not equivalent to PARP deletion.
PARP1 is required for cell killing by olaparib. A, Western blot analysis of whole cell lysates prepared from wild-type and PARP1−/− DT40 cells. Blots were probed with the indicated antibodies. B, nuclear localization of PARP1 in wild-type DT40 cells and loss of signal in PARP1−/− DT40 cells. Immunofluorescence microscopy was conducted with anti-PARP1 antibody and nuclei were stained with 4′,6-diamidino-2-phenylindole (DAPI). C, PAR level in wild-type and PARP1−/− DT40 cells measured by ELISA. Cells were treated with the indicated olaparib concentrations for 2 hours. D, survival curves after continuous olaparib treatment for 72 hours. Cellular ATP activity was used to measure cell viability. The survival of untreated cells was set as 100%. Error bars represent SD (n = 3). E, cell-cycle analyses after olaparib treatment. The indicated cell lines were treated with 10 μmol/L olaparib for 24 hours; after which, cells were analyzed by flow cytometry analysis. The left and right peaks indicate G1 and late S and G2 populations, respectively. PI, propidium iodide. F, DNA damage measured by immunofluorescence microscopy with anti–γ-H2AX antibody after 2 hours treatments with or without 10 μmol/L olaparib. Left, representative images and right plot, a quantitation of the number of γ-H2AX foci per individual cells. The number of cells analyzed (N) and the average number of γ-H2AX foci (Ave.) are shown above the plot.
PARP1 is required for cell killing by olaparib. A, Western blot analysis of whole cell lysates prepared from wild-type and PARP1−/− DT40 cells. Blots were probed with the indicated antibodies. B, nuclear localization of PARP1 in wild-type DT40 cells and loss of signal in PARP1−/− DT40 cells. Immunofluorescence microscopy was conducted with anti-PARP1 antibody and nuclei were stained with 4′,6-diamidino-2-phenylindole (DAPI). C, PAR level in wild-type and PARP1−/− DT40 cells measured by ELISA. Cells were treated with the indicated olaparib concentrations for 2 hours. D, survival curves after continuous olaparib treatment for 72 hours. Cellular ATP activity was used to measure cell viability. The survival of untreated cells was set as 100%. Error bars represent SD (n = 3). E, cell-cycle analyses after olaparib treatment. The indicated cell lines were treated with 10 μmol/L olaparib for 24 hours; after which, cells were analyzed by flow cytometry analysis. The left and right peaks indicate G1 and late S and G2 populations, respectively. PI, propidium iodide. F, DNA damage measured by immunofluorescence microscopy with anti–γ-H2AX antibody after 2 hours treatments with or without 10 μmol/L olaparib. Left, representative images and right plot, a quantitation of the number of γ-H2AX foci per individual cells. The number of cells analyzed (N) and the average number of γ-H2AX foci (Ave.) are shown above the plot.
Olaparib stabilizes toxic PARP1–DNA complexes
The finding that the cytotoxicity of olaparib was selective for PARP1-proficient cells raised the possibility that olaparib might generate toxic DNA lesions mediated by PARP1 at spontaneously occurring base lesions. To examine such possibility, we fractionated wild-type DT40 cell lysates into nuclear-soluble and chromatin-bound fractions. To increase base damage, we combined olaparib with low dose of the classic alkylating agent MMS. Under drug-free condition, most PARP1 was in the nuclear-soluble fraction (Fig 2A, lanes 1 and 7). Olaparib or MMS alone slightly increased PARP1 chromatin binding (Fig 2A, lanes 8 and 9, and 2B). Notably, in the presence of MMS, increasing the concentration of olaparib greatly induced PARP1 accumulation in the chromatin-bound fraction (Fig. 2A, lanes 10–12, and 2B). Under these conditions, PAR levels were reduced by olaparib even in the presence of 0.01% MMS (Fig. 2A).
Olaparib stabilizes PARP1–DNA complexes that are more toxic than unrepaired SSBs. A, Western blot analysis of nuclear soluble and chromatin-bound fractions prepared from wild-type DT40 cells. Cells were treated for 30 minutes without drug (lanes 1 and 7) or as indicated. Blots were probed with indicated antibodies. Histone H3 and topoisomerase I were used as positive markers for chromatin-bound and nuclear soluble fractions, respectively. The asterisk indicates a nonspecific band. B, quantification of PARP1 levels in chromatin binding fractions 30 minutes after the indicated drug treatments. Signal intensity was quantified using Image J software (NIH) from 4 to 7 independent Western blot analyses (mean values ± SD) and normalized to untreated cells set as one. C, survival curves of wild-type (left) and PARP1−/− (right) cells treated with MMS alone or with the indicated olaparib concentrations (indicated beside each curve in micromolar unit). The survival of treated cells without MMS was set as 100%. Data are mean ± SD (n = 3). D, cell-cycle analysis of wild-type and PARP1−/− DT40 cells after 2- and/or 4-hour continuous treatments. Cells were pulse-labeled with BrdUrd for 10 minutes following the drug treatment to measure DNA synthesis. PI, propidium iodide.
Olaparib stabilizes PARP1–DNA complexes that are more toxic than unrepaired SSBs. A, Western blot analysis of nuclear soluble and chromatin-bound fractions prepared from wild-type DT40 cells. Cells were treated for 30 minutes without drug (lanes 1 and 7) or as indicated. Blots were probed with indicated antibodies. Histone H3 and topoisomerase I were used as positive markers for chromatin-bound and nuclear soluble fractions, respectively. The asterisk indicates a nonspecific band. B, quantification of PARP1 levels in chromatin binding fractions 30 minutes after the indicated drug treatments. Signal intensity was quantified using Image J software (NIH) from 4 to 7 independent Western blot analyses (mean values ± SD) and normalized to untreated cells set as one. C, survival curves of wild-type (left) and PARP1−/− (right) cells treated with MMS alone or with the indicated olaparib concentrations (indicated beside each curve in micromolar unit). The survival of treated cells without MMS was set as 100%. Data are mean ± SD (n = 3). D, cell-cycle analysis of wild-type and PARP1−/− DT40 cells after 2- and/or 4-hour continuous treatments. Cells were pulse-labeled with BrdUrd for 10 minutes following the drug treatment to measure DNA synthesis. PI, propidium iodide.
To examine the cytotoxicity of the olaparib-induced PARP1–DNA complexes, we determined the cytotoxicity of MMS with increasing concentrations of olaparib in wild-type and PARP1−/− cells (Fig. 2C). As expected from the widely accepted concept that PARylation is necessary for BER, PARP1−/− cells were hypersensitive to MMS. Moreover, olaparib had no impact on MMS sensitivity in PARP1−/− cells (Fig. 2C, right), confirming the lack of “off-target effects” of olaparib. In parallel, olaparib synergistically sensitized wild-type cells to MMS in a concentration-dependent manner (Fig 2C, left). Notably, the MMS sensitivity of wild-type cells treated with olaparib was greater than that of PARP1−/− cells (compare left and right in Fig. 2C). These results indicate that PARylation inhibition does not simply explain the hypersensitivity of olaparib-treated cells to MMS. Cell-cycle responses in wild-type cells (Fig. 2D) also showed that the combination of olaparib with MMS suppressed BrdUrd incorporation more than either single treatment, and also induced sub-G1 accumulation after 4 hours, indicative of dying cells. In contrast, olaparib had no additional impact on PARP1−/− cells as compared with MMS alone. Together, these results show that olaparib synergistically increases the cytotoxicity of MMS in a PARP1-dependent manner, as it induces the trapping of PARP1–DNA complexes.
Trapping of human PARP1 and PARP2 by olaparib
To determine whether both human PARP1 and PARP2 act in the cytotoxicity of olaparib, we compared PARP1−/− DT40 cells complemented either with human PARP1 (PARP1−/−; hPARP1) or PARP2 (PARP1−/−; hPARP2; ref. 35). Complementation with either hPARP1 or hPARP2 restored sensitivity to olaparib (Fig. 3A), showing that both hPARP1 and hPARP2 can be poisoned by olaparib.
Stabilization of toxic PARP1- and PARP2–DNA complexes by olaparib. A, survival curves of wild-type, PARP1−/−, PARP1−/− complemented with human PARP1 (PARP1−/−; hPARP1) and PARP2 (PARP1−/−; hPARP2) DT40 cell lines treated with olaparib. Cells were treated with olaparib as in Fig. 1D. Error bars represent SD (n = 3). B–E, Western blotting of nuclear soluble and chromatin-bound fractions prepared from human DU145 cells. Blots were probed with the indicated antibodies. B, cells were treated for 4 hours without drug (lanes 1 and 7) or as indicated. C, rapid reduction of PARP1 and PARP2 binding to chromatin after washing out olaparib. Cells were treated for 4 hours without drug or with 10 μmol/L olaparib in the presence of 0.01% MMS. Drugs were removed by washing cells with cold PBS 3 times, followed by incubation in prewarmed medium for the indicated times. D, control siRNA (siCtrl), PARR1 siRNA (siP1), and PARP2 siRNA (siP2) were transfected into DU145 cells. Three days later, cells were treated for 4 hours with10 μmol/L olaparib in the presence of 0.01% MMS. E, reduced cytotoxicity of olaparib in DU145 cells downregulated for PARP1. Cells were transfected with the indicated siRNAs for 60 hours before adding olaparib for an additional 72 hours. Sensitivity was determined as Fig. 1D. Error bars represent SD (n = 3).
Stabilization of toxic PARP1- and PARP2–DNA complexes by olaparib. A, survival curves of wild-type, PARP1−/−, PARP1−/− complemented with human PARP1 (PARP1−/−; hPARP1) and PARP2 (PARP1−/−; hPARP2) DT40 cell lines treated with olaparib. Cells were treated with olaparib as in Fig. 1D. Error bars represent SD (n = 3). B–E, Western blotting of nuclear soluble and chromatin-bound fractions prepared from human DU145 cells. Blots were probed with the indicated antibodies. B, cells were treated for 4 hours without drug (lanes 1 and 7) or as indicated. C, rapid reduction of PARP1 and PARP2 binding to chromatin after washing out olaparib. Cells were treated for 4 hours without drug or with 10 μmol/L olaparib in the presence of 0.01% MMS. Drugs were removed by washing cells with cold PBS 3 times, followed by incubation in prewarmed medium for the indicated times. D, control siRNA (siCtrl), PARR1 siRNA (siP1), and PARP2 siRNA (siP2) were transfected into DU145 cells. Three days later, cells were treated for 4 hours with10 μmol/L olaparib in the presence of 0.01% MMS. E, reduced cytotoxicity of olaparib in DU145 cells downregulated for PARP1. Cells were transfected with the indicated siRNAs for 60 hours before adding olaparib for an additional 72 hours. Sensitivity was determined as Fig. 1D. Error bars represent SD (n = 3).
In human DU145 prostate cancer cells, olaparib induced chromatin binding of both PARP1 and PARP2 in the presence of 0.01% MMS (Fig. 3B, lanes 10–12). Trapping of PARP1 and PARP2 was rapidly reversible upon drug removal (Fig. 3C) with release of PARP1 and PARP2 from chromatin within 30 minutes. Dual trapping of PARP1 and PARP2 by olaparib was also observed in other human cancer cells: OVCAR4 (ovarian cell line) and SF295 (central nervous system; Supplementary Fig. S3). Suppression of PARP1 expression level by siRNA reduced the trapping of PARP1 but did not affect the trapping of PARP2 (Fig. 3D). Conversely, siRNA against PARP2 abolished the trapping of PARP2 by olaparib without affecting PARP1 trapping (Fig. 3D). Similar results were obtained in SF295 cells (Supplementary Fig. S3B), showing that olaparib traps both PARP1- and PARP2–DNA complexes independently of each other.
Next, we tested whether genetic depletion of PARP1 or PARP2 would make human cells tolerant to olaparib, as observed in DT40 cells (see Fig. 1D). As expected, transfection of DU145 cells with siRNA against PARP1 reduced sensitivity to olaparib (Fig. 3E). However, siRNA against PARP2 had minimal impact despite the efficiency of the knockdown (see Fig. 3D). These results indicate that the presence of PARP1, whose expression is much higher than PARP2 in human cells (5, 37) is the major mediator of olaparib-induced antiproliferative activity in DU145 cells.
Differential potency of PARP inhibitors to poison PARP1 and PARP2 in cells
Having established the dual effect of olaparib on PARP (catalytic inhibition and PARP trapping), we compared olaparib with 2 other clinically relevant PARP inhibitors, veliparib and niraparib (see structures in Supplementary Fig. S1; ref. 27). Catalytic PARP inhibition was examined by 2 methods (Western blotting and ELISA assays; refs. 34, 38). Both assays showed that olaparib was the most potent catalytic PARP inhibitor, followed by veliparib and niraparib (Fig. 4A and 4B). For all 3 drugs, PAR levels were almost undetectable at and above 1 μmol/L drug concentration (Fig. 4A).
Three clinical PARP inhibitors differ in their potency to poison PARP1 and PARP2 irrespective of their potency to inhibit PARP catalytic activity. A, drug-induced PARylation inhibition. Western blotting of PAR levels in whole cell lysates from wild-type DT40 cells treated as indicated for 30 minutes. Blots were probed with anti-PAR antibody. Whole cell lysates from untreated PARP1−/− cells were used as a control. Asterisks indicate nonspecific bands. B, representative PAR ELISA assays in DT40 cells treated with olaparib, veliparib, or niraparib for 120 minutes. PAR level of untreated cells was set as 100%. IC50 (inhibitory concentration 50%) of olaparib, veliparib, and niraparib were 1.2 nmol/L, 10.5 nmol/L, and 50.5 nmol/L, respectively. C–F, differential cytotoxicity of the 3 PARP inhibitors. Viability was measured as in Fig. 1D. Error bars represent SD (n = 3). Survival curves of the wild-type cells (C), PARP1−/− cells (D), wild-type cells treated with a subtoxic MMS concentration (0.00025%; see Fig. 2C; E), and BRCA2-deficient DT40 cells (F). Sensitivity was determined as in Fig. 1D. Error bars represent SD (n = 3). The viability of PARP1−/− cell to 0.00025% MMS (22%) was shown with a horizontal dashed line with annotation (E). G, survival curves of DU145 treated with the indicated PARP inhibitors for 72 hours. H, survival curves of DU145 treated with MMS alone or MMS plus 1 μmol/L olaparib or veliparib or niraparib for 72 hours.
Three clinical PARP inhibitors differ in their potency to poison PARP1 and PARP2 irrespective of their potency to inhibit PARP catalytic activity. A, drug-induced PARylation inhibition. Western blotting of PAR levels in whole cell lysates from wild-type DT40 cells treated as indicated for 30 minutes. Blots were probed with anti-PAR antibody. Whole cell lysates from untreated PARP1−/− cells were used as a control. Asterisks indicate nonspecific bands. B, representative PAR ELISA assays in DT40 cells treated with olaparib, veliparib, or niraparib for 120 minutes. PAR level of untreated cells was set as 100%. IC50 (inhibitory concentration 50%) of olaparib, veliparib, and niraparib were 1.2 nmol/L, 10.5 nmol/L, and 50.5 nmol/L, respectively. C–F, differential cytotoxicity of the 3 PARP inhibitors. Viability was measured as in Fig. 1D. Error bars represent SD (n = 3). Survival curves of the wild-type cells (C), PARP1−/− cells (D), wild-type cells treated with a subtoxic MMS concentration (0.00025%; see Fig. 2C; E), and BRCA2-deficient DT40 cells (F). Sensitivity was determined as in Fig. 1D. Error bars represent SD (n = 3). The viability of PARP1−/− cell to 0.00025% MMS (22%) was shown with a horizontal dashed line with annotation (E). G, survival curves of DU145 treated with the indicated PARP inhibitors for 72 hours. H, survival curves of DU145 treated with MMS alone or MMS plus 1 μmol/L olaparib or veliparib or niraparib for 72 hours.
In terms of cytotoxicity in wild-type DT40 cells, niraparib was the most potent drug, whereas veliparib was the least cytotoxic and olaparib was intermediate (Fig. 4C). The cytotoxicity was mediated by PARP1 as shown by the relative lack of effect of the drugs in PARP1−/− cells (Fig. 4D). Moreover, in the case of veliparib, high drug concentration (50 μmol/L) was required to synergize with MMS beyond the sensitivity of PARP1−/− cells, whereas for olaparib or niraparib, 1 μmol/L (50-fold less) was enough to achieve superior synergism (Fig. 4E). These results suggested that veliparib had a lesser potency to poison PARP1 than olaparib and niraparib. Despite its minimal cytotoxicity in wild-type cells, veliparib sensitized BRCA2tr/− cells, which is consistent with the general concept of synthetic lethality by PARP inhibitors in BRCA-deficient tumors. Yet, veliparib was approximately 10-fold less cytotoxic than olaparib in the BRCA2tr/− cells, and niraparib was even more effective than olaparib as a single agent in the BRCA2tr/− cells (Fig. 4F).
The differential cytotoxicity of the 3 PARP inhibitors was also observed in human DU145 cells (Fig. 4G and H). Single treatment with either olaparib or niraparib killed DU145 in a concentration-dependent manner, whereas veliparib had little cytotoxicity (Fig. 4G). Under conditions in which 1 μmol/L veliparib showed no effect, 1 μmol/L olaparib or niraparib synergized markedly with MMS (Fig. 4H). Similar results were obtained in the human SF295 cell line (Supplementary Fig. S4A and S4B). Therefore, the cytotoxicity of 3 clinically relevant PARP inhibitors shows marked differences across different cell lines with niraparib and olaparib being the most cytotoxic and veliparib the least. Moreover, our results show that the antiproliferative activity of the PARP inhibitors, including in BRCA2tr/− cells, is not directly correlated with their PARP catalytic inhibition potencies.
We next tested whether the differential cytotoxicity of the PARP inhibitors was related to their trapping of PARP1 and PARP2. Even at 50 μmol/L, veliparib induced significantly less trapping of PARP1 and PARP2 on chromatin than olaparib or niraparib at 10 μmol/L across different cell types (DT40, human DU145 and SF295 cells; Fig. 5A and Supplementary Fig. S4C). To examine the tightness of PARP–DNA complexes induced by different PARP inhibitors, further fractionation experiments were carried out (Fig. 5B and C). Most of PARP1 was eluted in fraction B from untreated or veliparib-treated cells, whereas it was mostly eluted in fraction D in olaparib- and niraparib–treated cells. Histone H3 and monoubiquitinated PCNA were also eluted dominantly in fraction D. These results indicate that PARP1–DNA complexes trapped by olaparib and niraparib are tighter than by veliparib. We also compared in our different assays, 4-AN, a commonly used PARP inhibitor (32, 33, 39). Supplementary Fig. S5 shows that its PARP poisoning potency was stronger than veliparib, but was less than olaparib and niraparib, whereas having similar catalytic PARP inhibitory ability as niraparib. Together, these results show that PARP inhibitors are not equivalent with respect to trapping PARP–DNA complexes. Veliparib is the least potent poison of PARP1 and PARP2 at concentrations that are effectively blocking PARP catalytic activity.
Differential cellular trapping of PARP1 and PARP2 by clinical PARP inhibitors. A, Western blotting of nuclear soluble and chromatin-bound fractions prepared from wild-type DT40 cells (top) and DU145 cells (bottom) with or without drug treatment as indicated. Blots were probed with indicated antibodies. B and C, fractionation of PARP–DNA complexes stabilized by PARP inhibitors. DU145 cells were treated for 4 hours without drug (lanes 1–5) or as indicated. B, scheme for the preparation of indicated fractions (P1, A–D) using different stringency buffers (Supplementary Methods). C, the eluted proteins in each fraction were blotted with indicated antibodies. Arrows, monoubiquitinated PCNA.
Differential cellular trapping of PARP1 and PARP2 by clinical PARP inhibitors. A, Western blotting of nuclear soluble and chromatin-bound fractions prepared from wild-type DT40 cells (top) and DU145 cells (bottom) with or without drug treatment as indicated. Blots were probed with indicated antibodies. B and C, fractionation of PARP–DNA complexes stabilized by PARP inhibitors. DU145 cells were treated for 4 hours without drug (lanes 1–5) or as indicated. B, scheme for the preparation of indicated fractions (P1, A–D) using different stringency buffers (Supplementary Methods). C, the eluted proteins in each fraction were blotted with indicated antibodies. Arrows, monoubiquitinated PCNA.
Differential potency of PARP inhibitors to poison PARP1 and PARP2 in biochemical assays
To elucidate whether the differential trapping of PARP–DNA complexes by the 3 PARP inhibitors could be observed at the molecular level, we set up a PARP1–DNA–binding assay based on fluorescence anisotropy. A nicked oligonucleotide duplex DNA with a single 5′-dRP end at the break site was used as fluorescent substrate (Fig. 6A). Its anisotropy was enhanced upon PARP1 binding to the damaged site. PARylation following addition of NAD+ reduced the fluorescence anisotropy signal by freeing the DNA. Figure 6B shows that all 3 PARP inhibitors enhanced the fluorescence anisotropy signal, which reflects the stabilization of PARP1–DNA complexes. Olaparib and niraparib were approximately 10-fold more potent than veliparib. Time-course experiments following NAD+ addition (Fig. 6C) also showed that the PARP inhibitors slowed down the dissociation of PARP1–DNA complexes with olaparib and niraparib being significantly more potent than veliparib. These results are in agreement with the cellular data showing greater stabilization of PARP–DNA complexes by olaparib and niraparib than by veliparib.
Differential biochemical trapping of PARP1 by clinical PARP inhibitors. A, scheme of the fluorescence anisotropy (FA)–binding assay. The star indicates the site labeled on the DNA substrate with Alexa Fluor488. Unbound nicked DNA substrate rotates fast and gives low fluorescence anisotropy value. PARP1 binding to the substrate slows the rotation and gives high fluorescence anisotropy value. Addition of NAD+ leads to PARP1 dissociation from DNA due to autoPARylation. B, concentration-dependent PARP1–DNA association using the 3 different PARP inhibitors. Fluorescence anisotropy value was measured 60 minutes after adding NAD+. C, time-course of PARP1–DNA dissociation in the presence of the 3 different PARP inhibitors (1 μmol/L each). The absence of PARP inhibitors immediately reduces PARP1–DNA complexes (DMSO control). In the absence of NAD+, PARP1–DNA complexes remained stable for at least 120 minutes (no NAD+).
Differential biochemical trapping of PARP1 by clinical PARP inhibitors. A, scheme of the fluorescence anisotropy (FA)–binding assay. The star indicates the site labeled on the DNA substrate with Alexa Fluor488. Unbound nicked DNA substrate rotates fast and gives low fluorescence anisotropy value. PARP1 binding to the substrate slows the rotation and gives high fluorescence anisotropy value. Addition of NAD+ leads to PARP1 dissociation from DNA due to autoPARylation. B, concentration-dependent PARP1–DNA association using the 3 different PARP inhibitors. Fluorescence anisotropy value was measured 60 minutes after adding NAD+. C, time-course of PARP1–DNA dissociation in the presence of the 3 different PARP inhibitors (1 μmol/L each). The absence of PARP inhibitors immediately reduces PARP1–DNA complexes (DMSO control). In the absence of NAD+, PARP1–DNA complexes remained stable for at least 120 minutes (no NAD+).
PARP–DNA complexes are repaired by several repair pathways besides HR
To gain insight into the cellular repair pathways involved in the repair of stabilized PARP1- and PARP2–DNA complexes, we compared the sensitivity of 30 isogenic mutant DT40 cell lines with preestablished genetic DNA repair deficiencies (Fig. 7A and Supplementary Table S1). An obvious difference between the 3 PARP inhibitors was the resistance of PARP1−/− cell lines to olaparib and niraparib (light blue in Fig. 7A), which is consistent with the requirement of PARP1 for the drugs cytotoxic activities (see above). As expected, the HR-deficient cell lines: BRCA2, XRCC3, XRCC2, BRCA1, RAD54, and H2AX (dark blue in Fig. 7A) were hypersensitive to all 3 PARP inhibitors. The greater differential for veliparib (>16-fold) is probably because of veliparib's low cytotoxicity in wild-type cells. Notably, deficiencies in check point including ATM and RAD9 (light pink), postreplication repair pathways including PCNA ubiquitination and its ubiquitin ligases RAD18 and UBC13 (brown), all FA genes tested (magenta), polymerase β (orange), and FEN1 (purple) conferred greater than approximately 2-fold hypersensitivity to all 3 drugs, indicating that DNA lesions induced by PARP inhibitors are not only repaired by classic HR (see Discussion). Involvement of the FA pathway was also inferred from enhanced monoubiquitination of FANCD2, a key regulator of the pathway, in response to olaparib (Supplementary Fig. S6). Interestingly, deletion of PARP1 in RAD18−/− cells (40) restored their resistance to olaparib and niraparib (Fig. 7A, compare RAD18 and PARP1/RAD18), suggesting that toxic PARP–DNA complexes but not catalytic inhibition of PARP kills RAD18−/− cells with defective postreplication repair.
DNA repair pathways for PARP inhibitors. A, sensitivity profiles of olaparib, veliparib, and niraparib in a panel of 30 DT40 isogenic repair mutant DT40 cell lines (see Supplementary Table S1 for additional information on the cell lines). The sensitivity in each mutant cell line was determined relative to wild-type cells (see Fig. 1D). Negative and positive scores indicate that a given cell line is hypersensitive or resistant to the drug, respectively. Each bar length reflects the degree of sensitivity or resistance to the drug (log2 units). Each bar is colored according to the category of DNA repair function (see Supplementary Table S1): black, nonhomologous end joining; pink, checkpoint; dark blue, HR; brown, postreplication repair; light blue, PARP1 and PARP1/RAD18; light green, removal of Top1 or Top2 cleavage complexes; magenta, Fanconi anemia (FA) pathway; orange, DNA polymerases; purple, FEN1; red, nucleotide excision repair; and green, RecQ helicases. The IC90 for olaparib and niraparib and IC50 for veliparib are shown at the bottom (mean ± SD, n = 13). B and C, schematic illustration of the dual inhibitory mechanism of PARP inhibitors (poisoning and catalytic inhibition). B, molecular interaction map showing the regulatory pathways for the formation and dissociation of PARP–DNA complexes. Symbols are: ↔, binding of PARP and DNA with the node corresponding to the PARP–DNA complex; ⇒, activation; ┤, inhibition (50). PARP inhibitors enhance the PARP–DNA complexes by 2 mechanisms: 1, inhibition of NAD+ binding; and 2, binding of PARP inhibitor in the NAD site activates (allosterically) the binding of PARP to DNA. See Discussion for further details. C, dual cytotoxic mechanisms of PARP inhibitors. 1, catalytic inhibition (upper pathway) interferes with the repair of SSBs, leading to replication fork damage that requires HR repair. 2, trapping of PARP–DNA complexes also leads to replication fork damage but uses additional repair pathways including Fanconi pathway (FA), template switching (TS), ATM, FEN1 (replicative flap endonuclease), and polymerase β.
DNA repair pathways for PARP inhibitors. A, sensitivity profiles of olaparib, veliparib, and niraparib in a panel of 30 DT40 isogenic repair mutant DT40 cell lines (see Supplementary Table S1 for additional information on the cell lines). The sensitivity in each mutant cell line was determined relative to wild-type cells (see Fig. 1D). Negative and positive scores indicate that a given cell line is hypersensitive or resistant to the drug, respectively. Each bar length reflects the degree of sensitivity or resistance to the drug (log2 units). Each bar is colored according to the category of DNA repair function (see Supplementary Table S1): black, nonhomologous end joining; pink, checkpoint; dark blue, HR; brown, postreplication repair; light blue, PARP1 and PARP1/RAD18; light green, removal of Top1 or Top2 cleavage complexes; magenta, Fanconi anemia (FA) pathway; orange, DNA polymerases; purple, FEN1; red, nucleotide excision repair; and green, RecQ helicases. The IC90 for olaparib and niraparib and IC50 for veliparib are shown at the bottom (mean ± SD, n = 13). B and C, schematic illustration of the dual inhibitory mechanism of PARP inhibitors (poisoning and catalytic inhibition). B, molecular interaction map showing the regulatory pathways for the formation and dissociation of PARP–DNA complexes. Symbols are: ↔, binding of PARP and DNA with the node corresponding to the PARP–DNA complex; ⇒, activation; ┤, inhibition (50). PARP inhibitors enhance the PARP–DNA complexes by 2 mechanisms: 1, inhibition of NAD+ binding; and 2, binding of PARP inhibitor in the NAD site activates (allosterically) the binding of PARP to DNA. See Discussion for further details. C, dual cytotoxic mechanisms of PARP inhibitors. 1, catalytic inhibition (upper pathway) interferes with the repair of SSBs, leading to replication fork damage that requires HR repair. 2, trapping of PARP–DNA complexes also leads to replication fork damage but uses additional repair pathways including Fanconi pathway (FA), template switching (TS), ATM, FEN1 (replicative flap endonuclease), and polymerase β.
Discussion
To our knowledge, our study provides the first direct evidence that PARP inhibitors trap both PARP1- and PARP2–DNA complexes. It shows the relevance of this mechanism for the antiproliferative activity (and therefore anticancer activity) of PARP inhibitors and provides novel insights into the molecular pathways that repair PARP–DNA complexes. The clinical relevance of our findings stems from the fact that the concentrations (<10 μmol/L) required to readily detect PARP1- and PARP2–DNA complexes are below the peak concentration of olaparib (24 μmol/L) in clinical trials (23).
Automodification is a well-established mechanism by which the highly negatively charged PAR polymers dissociate PARP enzymes from DNA (41; Fig. 7B, step d). In addition, our data show that PARP inhibitors differ in their potency to stabilize PARP–DNA complexes irrespective of their relative potency to inhibit PARP catalytic activity (see Fig 4, Supplementary Figs. S4 and S5). Indeed, PARP–DNA complexes were detectable under conditions in which PAR automodification remained detectable (Fig. 5A, Supplementary Figs. S3 and S4C). Similarly, a recent study showed that nuclear PARylation was observed even after treatment with 10 μmol/L veliparib with MMS (42). Therefore, we propose that PARP inhibitors poison PARP by 2 mechanisms. One is by inhibition of PARylation, which increases the binding of PARP to DNA (Fig. 7B, annotations 1, c, and d; refs. 41, 43); the other is drug binding to the NAD+ site, which allosterically enhances the DNA binding of PARP's N-terminal zinc finger domain (Fig. 7B, annotations 2 and e). The drug allosteric mechanism fits with the greater potency of the bulky inhibitors, niraparib and olaparib to produce PARP–DNA complexes as compared with veliparib and 4-AN (see Supplementary Fig. S1). This hypothetical reverse allosteric mechanism is supported by the fact that tight binding of PARP's N-terminus zinc domain to damaged DNA (Fig. 7B, step a) is known to allosterically activate the C-terminal catalytic pocket of PARP (1–3, 42, 44; Fig. 7B, step b). Hence, the PARP inhibitors can be viewed as blocking with the normal allosteric cycle of PARP and as revealing the normal impact of NAD+ binding to PARP. PARP binding to damaged DNA might be normally primed by NAD+ binding. In turn, PARP binding to damaged DNA induces the known allosteric change that activates PARP's catalytic site (44). The allosteric enhancement of PARP–DNA binding by PARP inhibitors is likely to be more significant for the anticancer activity of PARP inhibitors than catalytic inhibition. Indeed, while PARP1−/− DT40 cells overexpressing catalytic inactive PARP1 are mildly more sensitive to MMS than PARP1−/− cells (45), we find that olaparib, under conditions that induce PARP–DNA complexes, produces a much more drastic hypersensitivity to MMS than PARP−/− cells (Fig. 2C).
Our experiments clearly show that PARP inhibitors that were assumed to be equivalent based on PARP catalytic inhibition, are not equivalent with respect to their potency to poison PARP, suggesting the importance of categorizing PARP inhibitors according to their PARP poisoning potency. Depending on which PARP inhibitor is used, the clinical and experimental results could be very different. Furthermore, this persistent tendency to confuse different PARP inhibitors would make it impossible for the others to recapitulate the results depending on which PARP inhibitor was used. Therefore, our findings are relevant to basic researcher and clinicians who use PARP inhibitors. Specifically, our findings uniquely explain why veliparib has much lower toxicity than olaparib and niraparib although the 3 drugs have similar potency to inhibit PARP catalytic activity.
Our DT40 cell panel data confirm the critical importance of HRs including BRCA1 and BRCA2. They are consistent with the lethality of PARP1 siRNAs or short hairpin RNA (shRNA) in BRCA-deficient cells (21, 22) and with the activity of olaparib in BRCA-deficient tumors (23). BRCA2tr/−PARP1−/− double mutant DT40 cells are also unviable (unpublished data). These lines of evidence strongly support the current anticancer selectivity of PARP inhibitors based on the interference with SSB repair (Fig. 7C, upper pathway). On the other hand, in the case of the RAD18−/− cells, their sensitivity to the PARP inhibitors is unambiguously derived from the trapping of toxic PARP–DNA complexes (Fig. 7C, lower pathway) rather than from catalytic PARP inhibition, as the hypersensitivity of the RAD18−/− cells was rescued in RAD18−/−PARP1−/− double mutant DT40 cells.
Our present data (Fig. 7A) add XRCC2, XRCC3, RAD54, and H2AX to the list of HR genes that confer high sensitivity to PARP inhibitors. The hypersensitivity of check point–deficient cells (ATM−/− and RAD9−/−) and the resistance of nonhomologous end joining–deficient cells (KU70−/− and LIGASE IV−/−) extend recent reports (46, 47). The hypersensitivity of FEN1 and polymerase β–deficient cells may imply involvement of BER pathway or increase of endogenous DNA damage in these cells. RAD18 and UBC13, which are both responsible for PCNA ubiquitination on lysine 164 (48, 49), were identified in our DT40 panel as novel genes for survival to PARP inhibitors. Consistently, the PCNA mutant K164R also showed hypersensitivity to all 3 PARP inhibitors. PCNA ubiquitination has been proposed to play a central role in postreplication repair that uses translesion synthesis (TLS) and template switching to bypass bulky DNA lesions. For the repair of PARP–DNA lesion, template switching seems more important than TLS because the representative TLS polymerases, POLZ and POLH, mutant cells do not show apparent sensitivity to all drugs. The FA pathway seems critical for the repair of PARP–DNA complexes (see Supplementary Fig. S6). These results could rationalize the use of PARP inhibitors in tumors with deficiencies in the genes listed earlier.
In summary, this is the first report showing the toxic PARP–DNA complex by PARP inhibitors, and the greatly different potency to poison PARP among clinical PARP inhibitors irrespective of the potency to inhibit PARP catalytic activity. It also reveals new genetic repair pathways that are critical for the repair of PARP–DNA complexes, and potential synthetic lethality, which can be exploited for additional clinical use of PARP inhibitors.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Authors' Contributions
Conception and design: J. Murai, Y. Pommier
Development of methodology: J. Murai, S.N. Huang, Y. Zhang, J. Ji, Y. Pommier
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): J. Murai, S.N. Huang, B.B. Das, A. Renaud, J. Ji
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): J. Murai, S.N. Huang, B.B. Das, A. Renaud, Y. Zhang, J. Ji, S. Takeda, Y. Pommier
Writing, review, and/or revision of the manuscript: J. Murai, B.B. Das, J.H. Doroshow, S. Takeda, Y. Pommier
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): Y. Zhang, J.H. Doroshow
Study supervision: J. Ji, Y. Pommier
Acknowledgments
The authors thank Drs. Kouji Hirota and Hiroyuki Sasanuma (Kyoto University, Kyoto, Japan) for technical advice, Minoru Takata (Kyoto University, Kyoto, Japan) for sharing FANCC, FANCD2, and FANCG mutant DT40 cell lines, and Keith Caldecott (University of Sussex, UK) for kindly providing TDP2 mutant DT40 cell line. The authors also appreciate Drs. Mitsuko Masutani (National Cancer Center Research Institute, Japan) and Jun Nakamura (University of North Carolina) for a critical discussion.
Grant Support
J. Murai is a recipient of a fellowship from the Japan Society for the Promotion of Science (JSPS). This work was supported by the Intramural Program of the National Cancer Institute, Center for Cancer Research, and by the JSPS Core-to-Core Program.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.