Abstract
ADP-ribosyltransferases PARP1 and PARP2 play a major role in DNA repair mechanism by detecting the DNA damage and inducing poly-ADP-ribosylation dependent chromatin relaxation and recruitment of repair proteins. Catalytic PARP inhibitors are used as anticancer drugs especially in the case of tumors arising from sensitizing mutations. Recently, a study showed that Histone PARylation Factor (HPF1) forms a joint active site with PARP1/2. The interaction of HPF1 with PARP1/2 alters the modification site from Aspartate/Glutamate to Serine, which has been shown to be a key ADP-ribosylation event in the context of DNA damage. Therefore, disruption of PARP1/2-HPF1 interaction could be an alternative strategy for drug development to block the PARP1/2 activity. In this study, we describe a FRET based high-throughput screening assay to screen inhibitor libraries against PARP-HPF1 interaction. We optimized the conditions for FRET signal and verified the interaction by competing the FRET pair in multiple ways. The assay is robust and easy to automate. Validatory screening showed the robust performance of the assay, and we discovered two compounds Dimethylacrylshikonin and Alkannin, with µM inhibition potency against PARP1/2-HPF1 interaction. The assay will facilitate the discovery of inhibitors against HPF1-PARP1/2 complex and to develop potentially new effective anticancer agents.
Similar content being viewed by others
Introduction
ADP-ribosylation of proteins is a post-translational modification which is involved in regulation of many biological processes including DNA damage repair, transcription, and cell cycle regulation1. There are 17 ADP-ribosyltransferases (ARTs) in human with homology to Diphtheria toxin (ARTDs)2. These PARP enzymes and tankyrases regulate a wide range of signaling events and three of them, PARP1-3, are known to recognize DNA lesion sites and activate DNA repair cascades by transferring ADP-ribosyl moiety from NAD+ to DNA damage site, histones and to enzyme itself2,3,4,5. The ADP-ribose moiety is transferred in the form of single unit, mono-ADP ribosyl (MAR), and can be subsequently elongated to form poly-ADP-ribose chains (PAR). The PARylation of histone tail leads to nucleosome remodeling and PARP automodification is a regulatory mechanism for their own activity6. The modification occurs on aspartate, glutamate and serine residues but in presence of Histone PARylation Factor (HPF1), PARP1/2 are more specific and predominantly modify serine residues7,8,9,10,11. The preference to serine is enabled by HPF1 as it forms a joint active site with PARP1 and PARP2 as demonstrated by a crystal structures of the PARP1/2–HPF1 complexes12,13.
The role of PARP enzymes in DNA repair and development of small molecule inhibitors gained attention and interest when it was discovered that such catalytic PARP inhibitors have a synthetic lethality with BRCA mutations frequent in e.g. breast cancer14,15. Subsequently, multiple PARP inhibitors have been clinically approved for different indications16. The approved drugs compete with NAD+ substrate and thereby block the catalytic site of the enzyme. The small molecules are not often very selective17 and attempts have been made to generate PARP1 selective inhibitors18 and on the other hand compounds with dual activities19.
PARP1 contains seven well-structured domains including three zinc finger domains followed by BRCT domain, WGR domain and C-terminal catalytic fragment consisting of a regulatory helical domain and the ART domain. The zinc finger and WGR domains bind to the damaged DNA20,21. In contrast, PARP2 contains disordered N-terminal followed by structured WGR and catalytic part but lacks the N-terminal domains present in PARP1. It has been shown that the disordered N-terminal helps PARP2 in DNA binding5,22. In both PARP1 and PARP2 the helical domain is autoinhibitory while the ART domain interacts with HPF13,23. In presence of HPF1, PARP1/2 predominantly modify serine residues8,9. As an alternative strategy to already developed PARP inhibitors, it could be possible to prevent binding of HPF1 to PARP1/2 to inhibit serine ADP-ribosylation in DNA repair. HPF1 initiates PARylation on specific serine residues of histone tails subsequently leading to remodeling of the nucleosome with help of other factors such as ALC1 recognizing the modification24. As HPF1 shares a joint active site with PARP1/2 and plays a key role in increasing enzymatic activities, the inhibition of the interaction between HPF1-PARP would reduce histone PARylation upon DNA damage as single agent or possibly in combination with catalytic PARP inhibitors or DNA damaging agents. In the crystal structure, HPF1 and PARP2 interaction is mediated by a 15 amino acid long helix (261–284) of HPF1 packing against the PARP2 catalytic domain and NAD+ binding cleft12 (Fig. 1a). Through this, HPF1 and PARP2 collectively form a joint active site. The last two residues of PARP2 C-terminus interact at the back of HPF1 C-terminal domain and the overall interface buried is 1328.6 Å225.
In order to generate a signal that can be used in high-throughput mode to discover protein–protein interaction inhibitors as potential starting points for drug development, we decided to label the proteins with fluorophores to develop a FRET assay. The helical autoinhibitory domain would interfere in PARP-HPF1 interaction and hence we expressed fusion proteins where mCerulean (CFP) is attached to the N-terminus of the PARP ART domains. The mCitrine (YFP) is attached to the N-terminus of HPF1 as then the tags would not interfere with the complex formation based on the structure. When there is an interaction between labeled HPF1 and PARP1/2 then YFP and CFP come in close proximity and upon excitation of CFP there is an energy transfer from CFP to YFP resulting in FRET signal (Fig. 1b). Protein–protein interaction inhibitor would result in a loss of the FRET signal. We optimized the assay for screening of inhibitors of PARP2-HPF1 interaction and showed that the homogeneous assay is robust and can be run on a 384-well microplates in high-throughput screening campaigns. During the assay validation, we discovered two compounds Dimethylacrylshikonin and Alkannin as inhibitors of the protein–protein interaction.
Results
Assay setup and FRET measurement
To measure the interaction between PARP1/2 and HPF1 and to develop a high throughput method for inhibitor screening we used a FRET-based approach. We cloned constructs of PARP1 and PARP2 ART domains fused with CFP at the N-terminus and HPF1 fused with YFP at the N-terminus. Purified CFP-PARP1/2 and YFP-HPF1 were mixed, and FRET was measured by excitation of CFP at 410 nm wavelength and measuring YFP fluorescence intensity (527 nm) and CFP fluorescence intensity (477 nm). When CFP and YFP interact, the CFP emission at 477 nm decreases while YFP emission at 527 nm increases. The ratio of emission intensities of 527 nm and 477 nm, rFRET, was used throughout the work as a measure of binding interaction when comparing it with controls. As the yields for produced CFP-PARP2 were higher than those for CFP-PARP1 (25 mg/l and 10 mg/l, respectively), we decided to use CFP-PARP2 for optimization and for screening. For the interaction between HPF1 and PARP2, buffer conditions i.e., pH, buffering reagent, buffer and salt concentration, additives were optimized to get a robust and reproducible signal (Fig. 2).
First the effect of pH on rFRET (5 µM YFP-HPF1 and 1 µM CFP-PARP2) was tested using 10 mM Bis–Tris-Propane (BTP), which has pH ranges from 6.0 to 9.5 (Fig. 2a). At lower pH at 6.0, YFP-HPF1 (5 µM) and CFP-PARP2 (1 µM) FRET pairs showed maximum signal but it was evident that at pH at 6.0 and 6.5 proteins were aggregating in the well and this was interpreted as the source of the higher rFRET signal. Increasing pH resulted in diminishing signal, while at neutral pH 7.0 the rFRET signal was present and there were no signs of protein aggregation in the well. Therefore, we decided to further optimize the buffer conditions at pH 7.0, which is also close to physiological pH and therefore the ionization state of the potential inhibitors will be reasonable. After optimizing pH, 10 mM buffering reagents at pH 7.0 were tested for the FRET pair and since the highest signal was observed with BTP we decided to continue with it (Fig. 2b). Increasing concentrations of the buffering agent, however, decreased the signal (Fig. 2c). This could be an indication of sensitivity of the system to increasing salt concentrations and indeed, when we tested NaCl concentrations from no salt to 1 M we observed that even low salt concentration of 50 mM drastically reduced the rFRET signal (Fig. 2d). We therefore decided to use no salt in the buffer and tested the effect of Polyethylene Glycol 20,000 (PEG 20k) as a crowding agent. PEG has shown to enhance the rFRET signal and help in blocking adhesion the biomolecules to plastic walls26. Increasing PEG concentrations resulted in an increase in the rFRET signal (Fig. 2e). At higher PEG concentration rFRET signal increased very rapidly to unusually high rFRET values based on our previous experience27,28 and we interpreted this to be due to unspecific interactions. We decided therefore to add only a small amount of 3% (w/v) PEG 20k to the FRET buffer. Triton X-100 was found to slightly increase the signal and as it could help in the assay setup including the meniscus formation in the experimental plate, we also added 0.01% (w/v) to the final FRET buffer (Fig. 2f). Based on the optimization, 10 mM BTP, 0.01% tritonX-100, 0.5 mM TCEP, 3% PEG 20k, pH 7.0, was selected as an optimal FRET buffer. We next wanted to test the effect of HPF1-PARP2 stoichiometry and whether we would be able to reduce the proteins concentration while maintaining a good enough signal for screening. Increase in the relative HPF1 concentration slightly increased the rFRET signal systematically as more complex would be present in the mixture (Fig. 2g). The difference was however small and to save protein and to keep the sensitivity of the assay reasonable we made a compromise and decided to use 400 nM CFP-PARP2 and 800 nM YFP-HPF1 for further studies.
Signal and protein stability
In the initial part of optimization procedure, we studied the stability of the rFRET signal over time (Fig. 3a). We observed that the rFRET signal was increasing with time indicating non-steady state conditions that could be a result of instability of the proteins. This increase in signal over time would be a problem also in the context of a larger screening campaign as the dispensing and readouts need to be precisely timed. As the preparation of PARP1 and PARP2 without the autoinhibitory helical domain has been challenging we suspected that the protein constructs would not be stable in the used low salt buffer and measured the thermal stability using differential scanning fluorimetry (nanoDSF). PARP2 started unfolding at very low temperatures (Tm 25.6 ± 0.15 °C) with a not very well defined transition (Fig. 3b).
The increasing rFRET signal and nanoDSF data, collectively suggested that the protein was not stable at room temperature and was aggregating with time potentially bringing the fluorophores to proximity. Buffer additives such as glycine, guanidium hydrochloride, polyoxyethylene lauryl ether (Brij-35) were tested to stabilize the CFP-PARP2 protein, but additives did not help in stabilizing the protein in FRET buffer (data not shown)29,30. To stabilize the PARP2 fragment, fivefold excess of PARP catalytic inhibitor Olaparib was added in the FRET buffer. Addition of Olaparib resulted in a sharp unfolding transition for PARP2 with higher Tm value (44.6 ± 0.07 °C) by making the protein stable at room temperature (Fig. 3b). For the fusion protein 2 transitions were expected as it contains 2 independent domains with different Tms. The thermal melting was run up to 60 °C and Tm for CFP is 80 °C, so a single transition was observed for PARP2 ART domain only. In the presence of Olaparib, CFP-PARP2 and YFP-HPF1 interaction was stable in the FRET assay and simultaneously confirmed that Olaparib was not interfering in the PARP2-HPF1 interaction and could be used in the interaction assay (Fig. 3c).
In the presence of Olaparib, 400 nM CFP-PARP2 and 800 nM YFP-HPF1 were mixed resulting in a stable signal within an experiment, but the rFRET values were not reproducible in different experiments. This could indicate variation in complex formation between experiments depending on the order of addition and salt concentration that comes from the independently purified proteins. After trial and error, FRET pair proteins were mixed first in higher concentration (25 µM CFP-PARP2, 50 µM YFP-HPF1) in the buffer containing excess Olaparib (100 µM) in a so-called intermediate mixture. Final FRET mixture with 400 nM CFP-PARP2 and 800 nM YFP-HPF1 was then prepared by diluting the intermediate mixture in FRET buffer. Mixing of the proteins in high concentration allows formation of complex in saturating amounts before the dilution to the lower assay conditions. This procedure allowed us to get a stable, reproducible, and higher rFRET signal in comparison to a direct mixing to the final FRET conditions (Fig. 3d).
Signal verification
In the optimization phase we used rFRET as a direct readout for the interaction and secondary information like protein aggregation to judge whether the signal resulted from the interaction between the proteins. To test that the interaction was not with a fluorescent tag we carried out a control titration experiment (Fig. S3). To further confirm the specific interaction between PARP2 and HPF1, we next tested if unlabeled HPF1 would compete with the FRET pair and expected that the competition would result in the loss of rFRET signal. Increasing concentration of unlabeled HPF1 (35 nM–500 mM) in the reaction indeed decreased the rFRET signal in a concentration dependent manner, but it should be noted here that there is salt added also to the mixture coming from the unlabeled HPF1 protein stock (Fig. 4a).To rule out the effect of added salt, we repeated the experiment in the absence of HPF1 but using equivalent salt concentrations (25 nM–350 mM) which resulted in higher IC50 (pIC50: 4.2 ± 0.01, IC50: 64 µM) as compared to the competition by HPF1 (pIC50: 5.8 ± 0.04, IC50: 1.5 µM). The 1.6-log unit difference in the apparent IC50 values clearly indicates that unlabeled HPF1 competes with the FRET pair (Fig. 4a). To further strengthen the analysis and to confirm the competition we also exchanged the HPF1 buffer. The desalted unlabeled HPF1 also showed concentration dependent lowering (pIC50: 4.62, IC50: 24.0 µM) of the rFRET signal (Fig. 4b). We however observed aggregation of HPF1 in no salt condition at high concentrations and therefore the active unlabeled HPF1 protein concentration in the reaction was likely lower than the calculated total concentration while the competition was still evident. These results confirmed that unlabeled HPF1 (without salt) competed with the CFP-PARP2 and YFP-HPF1 FRET pair and validated the interaction between HPF1 and PARP2 in the FRET reaction.
In addition to the FRET pair competition with unlabeled HPF1, we also measured dynamic light scattering (DLS) for the CFP-PARP2 (0.5 mg/ml), YFP-HPF1 (1.0 mg/ml), and CPF-PARP2/YFP-HPF1 complex (Fig. S4a,c,d respectively). In DLS, even at the high protein concentration, scattering starts above 30 °C temperature confirming that there is no aggregation at room temperature (25 °C) and the FRET signal is the result of PARP2-HPF1 interaction.
Dissociation constant for PARP2-HPF1 interaction
The binding affinity of PARP2-HPF1 was measured using FRET according to the method described by Song et al.31. An increasing concentration of YFP-HPF1 (250 nM to 10.5 µM) was mixed with a constant concentration of CFP-PARP2 (400 nM) to saturate the CFP with YFP signal. Here we did not consider the effect of salt as ionic strength does not impact the interaction at high protein concentrations. Based on the raw fluorescence data, the FRET emission (EmFRET) was calculated, and the saturation curve showed Kd of 2.0 ± 0.18 µM (Fig. 4c). This is in line with the values and estimates from the literature 1.5 µM for HPF1-PARP1(ART only)13, while the Kd of 3 µM was reported for full length wild type PARP232. The Kd value for HPF1-PARP2 (ART) indicates that in the preparation of the FRET mixture the complex is formed in the intermediate mixing (25 µM CFP-PARP2, 50 µM YFP-HPF1) and this is then partially dissociated when diluted to the final FRET buffer (400 nM CFP-PARP2 and 800 nM YFP-HPF1).
Signal validation for screening
Once the CFP-PARP2 and YFP-HPF1 interaction was verified, the rFRET signal performance was validated by repeating the experiments on different days using different plates. The experiments were performed in 10 µl and 20 µl reaction volumes in 384-well plates, which contained 176 replicates of FRET pair (positive control) and 176 replicates of FRET pair with 1 M NaCl (negative control) (Fig. 5). The Zʹ-factors for 10 µl and 20 µl FRET reaction volumes were determined in a range of 0.85 to 0.92 and 0.75 to 0.88, respectively (Table S1 and S2). The Z’ factor value above 0.5 indicates the separation between negative and positive control is large, and assay with Z’ value above 0.5 is considered as excellent33. The statistical analysis suggests that even in a low assay volume of only 10 µl the quality of the signal was consistent, and the assay can be used for the screening of inhibitors.
Validatory screening
To validate the assay for inhibitor screening, Targetmol inhibitor compound library (1832 compounds), which contains biologically active compounds, was screened (Fig. 5c). For the screening, 10 nl compounds (10 µM concentration in reaction) and 10 µl FRET pair protein mixture (400 nM CFP-PARP2 and 800 nM YFP-HPF1) was dispensed and measured. In the presence of inhibitors, an additional excitation wavelength 430 nm was also used for the FRET measurement to rule out spurious effects due to compound fluorescence. This is based on the idea that for true protein–protein interaction inhibitors the inhibition at different excitation wavelengths is expected to be similar. On the other hand, artifacts due to small molecules fluorescence are expected to have different effects at different excitation wavelengths since they usually have sharper emission and excitation spectra. To analyze inhibition data, multiple filters were applied to rule out non-specific FRET signal coming from intrinsic fluorescence of individual compounds. Raw fluorescence signal with more than 20% fluorescence value (> 1.2* fluorescence at 527 nm; < 1.2*fluorescence at 477 nm) from the positive control FRET pair fluorescence signal, were considered as outlier and were not included for further analysis (Table S3). The % rFRET activity was calculated for each compound based on the controls. The hit limit for library compounds was set to 70.25% activity (5* standard deviations from mean) (Table S3). In the screening, a total of 3 compounds Dimethylacrylshikonin (CAS no. 24502-79-2), Alkannin (CAS no. 517-88-4), Crystal violet (CAS no. 548-62-9) were identified as hits for PARP2-HPF1 interaction (Fig S2 d-f).
Hit validation
To validate the hit compounds from CFP-PARP2 and YFP-HPF1 FRET pair screening, the binding of the compounds (Dimethylacrylshikonin, Alkannin, Crystal violet) with PARP2 and HPF1 were measured using nanoDSF (Fig. S2). When compound binds to the protein it results in the increase in the melting temperature of the protein. The Tm values for PARP2 and HPF1 in absence of hit compound were 44.2 ± 0.10 °C and 46.8 ± 0.02 °C, respectively. In the presence of hit compounds Dimethylacrylshikonin, Alkannin, and Crystal violet, the Tm value of PARP2 were measured as 46.0 ± 0.4 °C, 46.0 ± 0.2 °C, and 40.0 ± 0.5 °C, respectively (Fig. S2a–c). The increase in Tm in presence of Dimethylacrylshikonin and Alkannin, suggested their binding to PARP2 ART domain, but Crystal violet destabilizes the protein. The compounds did not cause a thermal shift for HPF1.
HPF1 has been crystallized with ART domain of PARP1 (PDB id: 6M3I)13 and PARP2 (PDB id: 6TX3)12. The ART domains of both PARP1 and PARP2 share a similar tertiary structure and interact with HPF1 approximately with the same interface and form a joint active site (Fig. S1). To study also the interaction between PARP1 ART domain and HPF1, CFP-PARP1 was mixed with YFP-HPF1 in 1:2 concentration ratio in optimized buffer for PARP2-HPF1 FRET pair (Fig. 6a). The FRET was measured using the same parameters as used for CFP-PARP2 and YFP-HPF1 FRET pair and a stable FRET signal was observed (Fig. 6b), which is not due to aggregation of CPF-PARP1/YFP-HPF1 complex (Fig. S4e). The assay thus allowed us to measure potencies of inhibitors also for the PARP1-HPF1 interaction.
To confirm the inhibition and to identify the potency of Dimethylacrylshikonin and Alkannin against PARP1/2 and HPF1 interaction, their IC50 values were measured. The IC50 values for Dimethylacrylshikonin, Alkannin against PARP1-HPF1 complex were 29.30 µM (pIC50: 4.53 ± 0.40) and 16.23 µM (pIC50: 4.80 ± 0.20), respectively (Fig. 7a,b) while IC50 values against PARP2-HPF1 complex were 15.73 µM (pIC50: 4.80 ± 0.17) and 21.20 µM (pIC50: 4.67 ± 0.06), respectively (Fig. 7c,d). Based on the binding of compounds to PARP ART domain and IC50 measurements, Dimethylacrylshikonin and Alkannin were identified as first described inhibitors of the PARP1/2-HPF1 interaction.
Inhibitory effect of compounds on PARP1/2 activity
Activities of full-length PARP1 or PARP2 were measured in presence of Dimethylacrylshikonin and Alkannin to study their inhibitory effect on PARP1/2 activity (Fig. 8). In presence of HPF1, PARP1/2 activity increases and treatment with hydroxylamine selectively removes PAR chains linked to Asp and Glu but do not affect Ser ADP-ribosylation34,35 and we confirmed this as a control (Fig. S5). The addition of 100 µM inhibitor compounds Dimethylacrylshikonin and Alkannin to the PARP activity reaction resulted in a decrease in PAR smear and increase of unmodified PARP1 (Fig. 8a, lane 6, 7 and 12, 13) and PARP2 (Fig. 8b, lane 6, 7, and 13, 14). The higher concentration of inhibitor compounds (1 mM) drastically reduced the amount of PAR smear and most PARP1 (Fig. 8a, lane 8, 14) and PARP2 (Fig. 8b, lane 8, 15) remained unmodified. The reduction in PARP1/2 modification confirmed that compounds interfere with PARP-HPF1 interaction in the context of full-length proteins and decrease the PARP activity resulting in the decrease in PAR formation.
Discussion
Current PARP inhibitors target the catalytic site and compete with NAD+. In the context of DNA damage the main form of ADP-ribosylation is on serine residues and dependent on HPF1. Interfering with the PARP-HPF1 interaction provides an alternative and unique way to inhibit PARP activity8,12,13. Protein–protein interactions are known to be challenging to target and so far there are no inhibitors reported for the PARP-HPF1 interaction. In this study, we have developed a FRET-based high-throughput screening assay that allows direct monitoring of HPF1-PARP interaction. The assay can be used for screening of small inhibitor molecules against PARP-HPF1 interaction. The proteins are expressed as fusion proteins with fluorescent proteins and when PARP interacts with HPF1 the FRET signal can be monitored using a common plate reader fluorimeter. We validated that the signal results from the specific interaction by competing it out with unlabeled HPF1. The signal and interaction in nanomolar concentrations is sensitive to salt making it convenient to use just high NaCl concentration for the control in lack of known inhibitors. Extensive efforts in optimization resulted in an assay which showed consistently a stable and robust signal (Z’ > 0.7) and can be run with automation in low 10 µl volumes in 384-well plates. The assay requires > 100 nM protein concentrations, but, as specialized reagents are not needed and the recombinant fusion proteins are produced in E. coli, the assay is cost effective and can be applied to large screening campaigns.
Recently, there was a study to find inhibitors against HPF1-PARP1 complex, where the interference of compounds between HPF1 and PARP1 was indirectly monitored by measuring the retention time of PARP1 on DNA36. This fluorescence polarization assay also identified catalytic PARP inhibitors and therefore the mechanism of the identified compounds may not be evident. In contrast, in the FRET-assay we describe here we measure directly the interaction between HPF1 and PARP and the hits, if not interfering with the signal or destabilizing proteins like crystal violet, are expected to act as protein–protein interaction inhibitors. The assay works even in the presence of a catalytic clinically used inhibitor and therefore selects the inhibitors mechanistically and avoids compounds only binding to the catalytic site unless the hit compounds extend from the catalytic site to the HPF1 binding interface. This feature extends the utility of the assay to test the mechanism of inhibitors identified by other means36,37.
During the validatory screening we discovered two hit compounds Dimethylacrylshikonin and Alkannin. These compounds stabilized PARP2 indicating direct binding and had µM potencies for both PARP1-HPF1 and PARP2-HPF1 interaction. The inhibitor compounds were also tested against full length PARP1/2 activity in presence of HPF1 and they were showing inhibitory effect on PARP activity. The screening method will hopefully facilitate the discovery of potent inhibitors against HPF1-PARP1/2 complexes and help to develop effective therapeutic molecules for cancer with a unique mechanism.
Materials and methods
Cloning
Inserts of PARP1/2 ART domain and full length HPF1 were prepared using gene specific primer in PCR. To make expression construct, PARP1 ART fragment [T661-F676 + GSGSGSGG + (D788-W1014)]38 or PARP2 ART domain were cloned into pNIC-CFP plasmid (addgene # 173083) and HPF1 was cloned into pNIC-YFP plasmid (addgene # 173080) using SLIC cloning method39. The pNIC-CFP or pNIC-YFP plasmids were linearized in PCR using site specific primers. 100 ng linearized plasmid was mixed in 1:3 molar ratio with gene PCR products and incubated with T4 DNA polymerase for 2–3 min at room temperature. The mixture was transformed into E. coli strain NEB5α cells and colonies were grown at 37 °C overnight on LB agar media plates containing, Kanamycin as antibiotic, 10 mM Benzamide, and 5% sucrose for SacB-based negative selection40,41. Plasmid of CFP tagged PARP1 contains “His6 tag-CFP-TEV protease site- ART domain [T661-F676 + GSGSGSGG + (D788-W1014)]”, PARP2 construct contains “His6 tag-CFP-TEV protease site- ART domain” and YFP tagged HPF1 contains His6 tag- CFP-TEV protease site followed by full length HPF1.
Protein expression
Plasmid containing CFP-PARP2 and YFP-HPF1 were transformed into E. coli BL21(DE3) and CFP-PARP1 was transformed to E. coli Rosetta2 (DE3) cells for the protein expression and incubated at 37 °C overnight. Transformed colonies were further inoculated in 500 ml autoinduction media which contains trace elements and supplied with 0.8% (w/v) glycerol, 50 μg/ml kanamycin antibiotic. PARP inhibitors 10 mM Benzamide and 2 mM 3-Amino Benzamide were added to the culture medium of CFP-PARP2 and CFP-PARP1, respectively42. Cultures were incubated at 200 rpm at 37 °C until an OD600 of 1.2, and then incubated for 18 h at 16 °C. Next day, cells were harvested by centrifugation at 4200 × g for 45 min at 4 °C. Pellets were resuspended in lysis buffer (50 mM HEPES, 350 mM NaCl, 10 mM imidazole, 0.5 mM TCEP, 10% (w/v) glycerol, pH 7.5) and 100 µM pefablock [4-(2-Aminoethyl) benzenesulfonyl fluoride hydrochloride] protease inhibitor was added in it. Resuspended cells were flash frozen in liquid nitrogen and stored at − 20 °C.
Protein purification
All the constructs (CFP-PARP1, CFP-PARP2, YFP-HPF1) were purified using two step purification i.e., immobilized metal affinity chromatography (IMAC) and size exclusion chromatography (SEC). Harvested cells were thawed, DNase (2 µg/ml) was added, and cells were lysed by sonication. Lysate was centrifuged at 39,000 × g RCF, for 1 h at 4 °C, supernatant was filtered and loaded on lysis buffer equilibrated, Ni–NTA column. The column was washed with 3 column volumes wash buffer 1 (30 mM HEPES, 1 M NaCl, 10% glycerol, 0.5 mM TCEP, 10 mM imidazole, pH 7.5) and wash buffer 2 (30 mM HEPES, 500 mM NaCl, 10% glycerol, 0.5 mM TCEP, 50 mM imidazole, pH 7.5). Protein was eluted in elution buffer (30 mM HEPES, 500 mM NaCl, 10% glycerol, 0.5 mM TCEP, 250 mM Imidazole, pH 7.5). The eluted protein was loaded on pre-equilibrated (30 mM HEPES, 500 mM NaCl, 10% glycerol, 0.5 mM TCEP, pH 7.5) 16/600 superdex-75 size-exclusion chromatography column. Purified proteins were concentrated using 30 kDa membrane filter, aliquoted in small volumes, flash-frozen in liquid nitrogen and stored at − 70 °C. The yield of proteins for YFP-HPF1, CFP-PARP2 and CFP-PARP1 were approximately 70 mg/liter, 25 mg/liter and 10 mg/liter, respectively. Full-length PARP1 expressed and purified as described previously by Venkannagari et al.43. The codon optimized (E. coli) cDNA of the human full length PARP2 isoform 1 (NP_005475) was expressed and purified as described previously by Obaji et al.44. For YFP-HPF1, CFP-PARP1/2, and full-length PARP1/2 all steps of protein purification were performed at 4 °C in a cold room.
FRET measurement
FRET signal measurements were performed using Tecan Spark multimode plate reader. Measurement was done using 384-well plate (proxiplate-384 F plus from Perkin Elmer). CFP was excited using monochromator at 410 nm (20 nm bandwidth) excitation wavelength, emission was recorded at emission wavelengths of 477 nm (10 nm bandwidth) and 527 nm (10 nm bandwidth). Other parameters such as number of flashes 50, Z position 20,500, signal integration time 40 µs, settle time 10 ms, manual gain 90 were optimized and used for the FRET measurement. To calculate the ratiometric FRET (rFRET) value, fluorescence intensity at 527 nm was divided by fluorescence intensity at 477 nm. All the rFRET measurements were done in optimized buffer 10 mM Bis–Tris-Propane, 0.5 mM TCEP, 3% (w/v) PEG 20k, and 0.01% (v/v) TritonX-100, pH 7.5, unless otherwise stated in negative control samples, 1 M NaCl was added.
Protein stability
Protein stability was tested using nanoDSF (differential scanning fluorimetry), where protein unfolding is monitored by tryptophan fluorescence. Different protein stabilizing regents, additives and PARP inhibitors were tested with CFP-PARP2 and YFP-HPF1 to see the effect on their melting temperatures (Tm values) and stability of CFP-PARP2 protein at room temperature. In addition to thermal stability, scattering was also monitored using the differential light scattering (DLS) data from nanoDSF to rule out any possibility of aggregation in PARP and HPF1 (Fig. S4). Any aggregation of protein leads to increase in particle size which results in increase in dynamic light scattering.
Competition assay
Competition assay was performed with an increasing concentration of unlabeled HPF1 and NaCl. In this experiment, the increasing log concentration of unlabeled HPF1 or NaCl was mixed with CFP-PARP2 (400 nM) and YFP-HPF1 (800 nM) FRET pair. FRET pair with 1 M NaCl and without unlabeled HPF1 were used as lowest and highest FRET signal controls, respectively. Data were measured in 4 replicates with 20 µl reaction volume in 384-well microplate. Statistics of the data and graph were analyzed using GraphPad Prism 9.0 using nonlinear regression analysis and data was plotted as inhibitor concentration v/s rFRET signal.
Determination of dissociation constant
To determine the dissociation constants of the FRET pair CFP-PARP2 and YFP-HPF1, constant CFP-PARP2 (400 nM) concentrations were mixed with increasing concentration of YFP-HPF1. First CFP-PARP2 and HPF1 were mixed in high concentration with increasing concentration of YFP-HPF1 in presence of Olaparib. The high concentration of YFP-HPF1 adds salt to the reaction from stock. This additional salt was counter balanced in all the reactions by backfilling of salt in buffer. The final reaction concentrations were achieved by mixing the dilutions to the buffer. The emission fluorescence was measured at 477 nm (10 nm bandwidth), 527 nm (10 nm bandwidth) by excitation at 430 nm (20 nm bandwidth), and emission at 527 nm (10 nm bandwidth) by excitation at 477 nm (20 nm bandwidth) wavelength. All other parameters for fluorescence measurement were kept the same as described in FRET measurements. Reaction was monitored in 10 µl volume with 4 replicates in 384-well plates. For the confidence in the Kd value, the experiments were performed 3 times, independently. Based on the raw fluorescence data, the FRET emission (EmFRET) was calculated using a formula “EmFRET = F527 at 430 excitation – α (F477 at 430 excitation) – β (F527 at 477 excitation)”. The fluorescence correction factor for CFP (α) and for YFP (β) were calculated as 0.477 and 0.058, respectively, using pure CFP and YFP constructs. The calculated EmFRET intensities were fitted in GraphPad Prism 9.0 and Kd value was calculated using nonlinear regression method with a the equation “Y = EmFRETmax-((2*EmFRETmax*Kd)/(X − A + Kd + sqrt(sqr(X − A − Kd) + (4*Kd*X))))” as described by Song et al.31.
Assay signal validation
The assay signal verification was done with 10 µl and 20 µl reaction volume. To verify the signal, FRET pair protein mixture (400 nM CFP-PARP2, 800 nM YFP-HPF1) in optimized buffer (10 mM Bis–Tris-Propane, 0.5 mM TCEP, 3% (w/v) PEG 20 k, 0.01% (v/v) TritonX-100, pH 7.5 and 1 µM Olaparib) with and without 1 M NaCl (192 replicates each condition) was dispensed to the 384-well plate using liquid dispenser robot (Formulatrix Mantis). This experiment was repeated for different days and with different plates on the same day. Data was analyzed in GraphPad Prism 9.0 using nonlinear regression curve fit and Z’-factor was calculated for the confidence in the assay33.
Validatory screening
For the validatory screening, 10 nl of the compounds from Targetmol inhibitor compound library (10 mM stock) were transferred to the 384-well plates using small volume dispensing robot (Echo 650, Beckmann Coulter). 10 µl reaction volume of FRET pair protein mixture (400 nM CFP-PARP2, 800 nM YFP-HPF1) in optimized buffer (10 mM Bis–Tris-Propane, 0.5 mM TCEP, 3% (w/v) PEG 20k, 0.01%(v/v) TritonX-100, pH 7.0 and 1 µM olaparib) was dispensed in wells containing compounds using liquid dispenser robot Formulatrix Mantis. Fluorescence intensities at 477 nm and 527 nm were measured by exciting the reaction sequentially at 410 nm and 430 nm wavelengths. To analyze inhibition data from libraries, filters were applied to the remove non-specific FRET signal coming from intrinsic fluorescence of individual compounds. The %rFRET activity was calculated from rFRET signal (Emission527 / Emission477) using formula “%rFRET = [{(rFRETreaction – rFRET negative control) / (rFRETpositive control – rFRETnegative control)}*100]”. The statistical analysis of inhibition was performed in GraphPad Prism 9.0 using nonlinear regression curve fit.
Hit validation
To study the binding study of hit compounds from screening was done using nanoDSF. For the experiment, unlabelled full length HPF1 and PARP2 ART domain with 1.0 mg/ml concentration were mixed in assay buffer (10 mM Bis–Tris-Propane, 0.5 mM TCEP, 3% (w/v) PEG 20k, and 0.01% (v/v) Triton X-100, pH 7.5). Additionally, PARP2 buffer was supplied with excess of PARP catalytic inhibitor Olaparib (100 µM PARP conc). Screened hit compounds (100 µM) from the FRET assay were added to the reaction mixture. The thermal melting scans were recorded in 3 replicates in three independent experiments from 20 °C to 60 °C with the temperature increase of 1 °C/min using Prometheus nanoDSF. The data was analyzed, and graphs were plotted using GraphPad Prism 9.0.
To validate the inhibition of FRET pair by hit compounds, the 50% inhibition concentration (IC50) values inhibitors against PARP1/2-HPF1 interaction, were measured using dilution series (100 µM to 3.25 nM) of compounds. Data was plotted against log inhibitor concentrations vs rFRET signal, analyzed in GraphPad Prism 9.0 using nonlinear regression analysis and IC50 values for compounds were calculated.
PARP1/2 activity assay on SDS-PAGE
PARP1/2 activities were performed on SDS-PAGE gel as described by Langelier et al.34 using the reaction buffer containing 50 mM Tris, 5 mM MgCl2 (10 × concentration), pH 8.0. For the inhibitory effect on the PARP1/2 activity, 1 µM PARP1 or 1 µM PARP was mixed with 1 µM central nicked double stranded dumbbell DNA (5'p GGG TCT TTT GAC CCT CGA GCT TTT GCT CGA 3') and incubated for 10 min at room temperature. 1 µM HPF1 followed by 500 µM NAD+ were added to the DNA protein mixture were incubated for 5 min. Reaction was stopped by adding 500 µM Olaparib to the reaction. Selected reactions were treated with 1 M hydroxylamine (NH2OH) for 1 h at room temperature. Inhibitory compounds were added to the buffer before starting the reaction. SDS loading buffer was added to the reactions, samples were heated at 95 °C and run on the 4–20% SDS-PAGE gel stained subsequently with page blue staining solution.
Data availability
The datasets generated during the current study are available from the corresponding author on reasonable request.
References
Gibson, B. A. & Kraus, W. L. New insights into the molecular and cellular functions of poly(ADP-ribose) and PARPs. Nat. Rev. Mol. Cell Biol. 13, 411–424 (2012).
Lüscher, B. et al. ADP-ribosyltransferases, an update on function and nomenclature. FEBS J. 289, 7399–7410 (2022).
Schreiber, V., Dantzer, F., Ame, J.-C. & de Murcia, G. Poly(ADP-ribose): Novel functions for an old molecule. Nat. Rev. Mol. Cell Biol. 7, 517–528 (2006).
Javle, M. & Curtin, N. J. The role of PARP in DNA repair and its therapeutic exploitation. Br. J. Cancer 105, 1114–1122 (2011).
Amé, J. C. et al. PARP-2, A novel mammalian DNA damage-dependent poly(ADP-ribose) polymerase. J. Biol. Chem. 274, 17860–17868 (1999).
De Vos, M., Schreiber, V. & Dantzer, F. The diverse roles and clinical relevance of PARPs in DNA damage repair: Current state of the art. Biochem. Pharmacol. 84, 137–146 (2012).
Leidecker, O. et al. Serine is a new target residue for endogenous ADP-ribosylation on histones. Nat. Chem. Biol. 12, 998–1000 (2016).
Bonfiglio, J. J. et al. Serine ADP-ribosylation depends on HPF1. Mol. Cell 65, 932-940.e6 (2017).
Palazzo, L. et al. Serine is the major residue for ADP-ribosylation upon DNA damage. eLife 7, e34334 (2018).
Crawford, K., Bonfiglio, J. J., Mikoč, A., Matic, I. & Ahel, I. Specificity of reversible ADP-ribosylation and regulation of cellular processes. Crit. Rev. Biochem. Mol. Biol. 53, 64–82 (2018).
Larsen, S. C., Hendriks, I. A., Lyon, D., Jensen, L. J. & Nielsen, M. L. Systems-wide analysis of serine ADP-ribosylation reveals widespread occurrence and site-specific overlap with phosphorylation. Cell Rep. 24, 2493-2505.e4 (2018).
Suskiewicz, M. J. et al. HPF1 completes the PARP active site for DNA damage-induced ADP-ribosylation. Nature 579, 598–602 (2020).
Sun, F.-H. et al. HPF1 remodels the active site of PARP1 to enable the serine ADP-ribosylation of histones. Nat. Commun. 12, 1028 (2021).
Farmer, H. et al. Targeting the DNA repair defect in BRCA mutant cells as a therapeutic strategy. Nature 434, 917–921 (2005).
Bryant, H. E. et al. Specific killing of BRCA2-deficient tumours with inhibitors of poly(ADP-ribose) polymerase. Nature 434, 913–917 (2005).
Li, Q. et al. A new wave of innovations within the DNA damage response. Sig. Transduct. Target Ther. 8, 338 (2023).
Wahlberg, E. et al. Family-wide chemical profiling and structural analysis of PARP and tankyrase inhibitors. Nat. Biotechnol. 30, 283–288 (2012).
Jamal, K. et al. Abstract 2609: AZD9574 is a novel, brain penetrant PARP-1 selective inhibitor with activity in an orthotopic, intracranial xenograft model with aberrant DNA repair. Cancer Res. 82, 2609–2609 (2022).
Plummer, R. et al. First-in-human study of the PARP/tankyrase inhibitor E7449 in patients with advanced solid tumours and evaluation of a novel drug-response predictor. Br. J. Cancer 123, 525–533 (2020).
Langelier, M.-F. & Pascal, J. M. PARP-1 mechanism for coupling DNA damage detection to poly(ADP-ribose) synthesis. Curr. Opin. Struct. Biol. 23, 134–143 (2013).
Langelier, M.-F., Planck, J. L., Roy, S. & Pascal, J. M. Structural basis for DNA damage–dependent Poly(ADP-ribosyl)ation by human PARP-1. Science 336, 728–732 (2012).
Meder, V. S., Boeglin, M., De Murcia, G. & Schreiber, V. PARP-1 and PARP-2 interact with nucleophosmin/B23 and accumulate in transcriptionally active nucleoli. J. Cell Sci. 118, 211–222 (2005).
Hakmé, A., Wong, H., Dantzer, F. & Schreiber, V. The expanding field of poly(ADP-ribosyl)ation reactions. EMBO Rep. 9, 1094–1100 (2008).
Bacic, L. et al. Structure and dynamics of the chromatin remodeler ALC1 bound to a PARylated nucleosome. eLife 10, e71420 (2021).
Krissinel, E. & Henrick, K. Inference of macromolecular assemblies from crystalline state. J. Mol. Biol. 372, 774–797 (2007).
Liu, B. et al. Parts-per-million of polyethylene glycol as a non-interfering blocking agent for homogeneous biosensor development. Anal. Chem. 85, 10045–10050 (2013).
Sowa, S. T. et al. A FRET-based high-throughput screening platform for the discovery of chemical probes targeting the scaffolding functions of human tankyrases. Sci. Rep. 10, 12357 (2020).
Sowa, S. T. et al. A molecular toolbox for ADP-ribosyl binding proteins. Cell Rep. Methods 1, 100121 (2021).
Bye, J. W., Platts, L. & Falconer, R. J. Biopharmaceutical liquid formulation: A review of the science of protein stability and solubility in aqueous environments. Biotechnol. Lett. 36, 869–875 (2014).
Schwarz, D. et al. Preparative scale cell-free expression systems: New tools for the large scale preparation of integral membrane proteins for functional and structural studies. Methods 41, 355–369 (2007).
Song, Y., Rodgers, V. G. J., Schultz, J. S. & Liao, J. Protein interaction affinity determination by quantitative FRET technology. Biotechnol. Bioeng. 109, 2875–2883 (2012).
Gaullier, G. et al. Bridging of nucleosome-proximal DNA double-strand breaks by PARP2 enhances its interaction with HPF1. PLoS One 15, e0240932 (2020).
Zhang, J.-H., Chung, T. D. Y. & Oldenburg, K. R. A simple statistical parameter for use in evaluation and validation of high throughput screening assays. SLAS Discov. 4, 67–73 (1999).
Langelier, M.-F., Billur, R., Sverzhinsky, A., Black, B. E. & Pascal, J. M. HPF1 dynamically controls the PARP1/2 balance between initiating and elongating ADP-ribose modifications. Nat. Commun. 12, 6675 (2021).
Kurgina, T. A. et al. Dual function of HPF1 in the modulation of PARP1 and PARP2 activities. Commun. Biol. 4, 1259 (2021).
Kellett, T. et al. HTS discovery of PARP1-HPF1 complex inhibitors in cancer. SLAS Discov. https://doi.org/10.1016/j.slasd.2023.10.003 (2023).
Choo, N. et al. High-throughput imaging assay for drug screening of 3D prostate cancer organoids. SLAS Discov. 26, 1107–1124 (2021).
Dawicki-McKenna, J. M. et al. PARP-1 activation requires local unfolding of an autoinhibitory domain. Mol. Cell 60, 755–768 (2015).
Jeong, J.-Y. et al. One-step sequence- and ligation-independent cloning as a rapid and versatile cloning method for functional genomics studies. Appl. Environ. Microbiol. 78, 5440–5443 (2012).
Hynes, M. F., Quandt, J., O’Connell, M. P. & Pühler, A. Direct selection for curing and deletion of Rhizobium plasmids using transposons carrying the Bacillus subtilis sacB gene. Gene 78, 111–120 (1989).
Reyrat, J. M., Pelicic, V., Gicquel, B. & Rappuoli, R. Counterselectable markers: Untapped tools for bacterial genetics and pathogenesis. Infect. Immun. 66, 4011–4017 (1998).
Banasik, M., Komura, H., Shimoyama, M. & Ueda, K. Specific inhibitors of poly(ADP-ribose) synthetase and mono(ADP-ribosyl)transferase. J. Biol. Chem. 267, 1569–1575 (1992).
Venkannagari, H. et al. Small-molecule chemical probe rescues cells from mono-ADP-ribosyltransferase ARTD10/PARP10-induced apoptosis and sensitizes cancer cells to DNA damage. Cell Chem. Biol. 23, 1251–1260 (2016).
Obaji, E., Haikarainen, T. & Lehtiö, L. Characterization of the DNA dependent activation of human ARTD2/PARP2. Sci. Rep. 6, 34487 (2016).
Acknowledgements
The use of the facilities of the Biocenter Oulu Structural Biology core facility, member of Biocenter Finland, Instruct-ERIC Centre Finland and FINStruct, is gratefully acknowledged. The research was supported by the Research Council of Finland (Grant No. 347026) and by the Biocenter Oulu spearhead project funding.
Author information
Authors and Affiliations
Contributions
S.S.D. carried out most of the experiments and compiled the manuscript. A.G.-P. tested the initial assay setup. L.L. conceived the study and L.L. and A.G.-P. supervised the work. All authors contributed to the final version of the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Dhakar, S.S., Galera-Prat, A. & Lehtiö, L. High-throughput screening assay for PARP-HPF1 interaction inhibitors to affect DNA damage repair. Sci Rep 14, 3875 (2024). https://doi.org/10.1038/s41598-024-54123-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-024-54123-8
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.