HPF1 binds to activated ART domain of PARP1
HPF1 has been shown to bind to the CAT domain (HD-ART) of PARP1 in cells in response to DNA damage 23. Since HD is an autoinhibitory domain that blocks productive binding to NAD+ in the resting state, and it undergoes local unfolding to enable NAD+ binding and ADP-ribose transferase activity when PARP1 binds to DNA breaks 19, we asked if unfolding or removing the HD subdomain is the prerequisite for HPF1 binding to PARP1. To answer this question, we tested HPF1 binding to PARP1 CAT domain in the full-length autoinhibited form or in the constitutively activated form through removing the majority (residues 679-786) of the autoinhibitory domain HD (CAT ΔHD) using size-exclusion chromatography (SEC). Interestingly, it is observed that CAT ΔHD, but not full-length CAT, forms a complex with HPF1 (Figure 1a). The interaction was further quantitatively characterized by isothermal titration calorimetry (ITC), which showed that CAT ΔHD binds to HPF1 with a dissociation constant (Kd) of about 1.5μM and 1:1 stoichiometry (Figure 1b, Table 1, Supplemental Figure S2,S3), while full-length CAT showed no binding to HPF1 (Figure 1c). Since these observations confirmed that HPF1 binds to activated CAT, we also tested the binding of HPF1 to full-length PARP1 activated by DNA (Table 1, Supplemental Figure S2). HPF1 binds to full-length PARP1 with a Kd of about 2.8μM, slightly weaker than its binding affinity for CAT ΔHD, probably related to transient folding of the HD domain in solution.
The human HPF1/PARP1-CAT ΔHD Complex crystal structure
In order to understand the assembly and function of the HPF1/PARP1 complex, the crystal structure of human HPF1/PARP1-CAT ΔHD complex bound by benzamide (introduced when expressing the protein) was determined and refined to 1.98 Å (Figure 2a, Supplemental Figure S4 and Table S1). HPF1 and PARP1 were confirmed to form a 1:1 complex since there were one HPF1 and one PARP1-CAT ΔHD molecule in the asymmetric unit of the complex crystal structure.
The high resolution structure of the HPF1/PARP1-CAT ΔHD complex reveals an extensive interface between the two proteins that envelopes the region of the PARP1 active site (see interface I in Figure 2a). The hetero-dimer binding interface was confirmed through extensive mutagenesis combined with ADP-ribosylation assay and ITC assay. It has been well established that HPF1 binding to PARP1 restricts hyper-automodification of PARP1 23, 26. Indeed, hyper-automodification of PARP1 showed up clearly on SDS-PAGE as smeared bands above the un-modified PARP1 band, while wild-type HPF1 binding abolished this effect (Figure 2b, see lanes 3 and 4). On the designated hetero-dimer interface, HPF1 Phe268, Phe280, Asp283, Cys285 and Lys307 are found to directly interact with PARP1 residues (Figure 2a, top-right and bottom-right insets), and these residues are conserved in HPF1 from different species (Supplemental Figure S5). Indeed, mutating most of these HPF1 residues (F268S, F280A, D283H, C285H and K307S) significantly restored hyper-automodification of PARP1 (Figure 2b, see lanes 9, 11, 12, 14 and 15), indicating loss of binding between the HPF1 mutants and PARP1. ITC assays further confirmed that these mutations dramatically reduced the binding affinity between HPF1 and PARP1-CAT ΔHD (Table 1, Supplemental Figure S2, S3). Importantly, HPF1 F268S and D283H mutations completely abolished the interaction between these two proteins, which may be used in following functional studies on the HPF1/PARP1 system. Interestingly, although PARP1, PARP2 and PARP3 share pretty high sequence homology in the CAT domain (25.41% of the sequences are identical and the other 18.85% are strongly similar), PARP3 lacks the key residues corresponding to His826, Leu985, Ser1012, Leu1013 and Trp1014 in PARP1, which according to the HPF1/PARP1-CAT ΔHD complex structure and the mutagenesis data are important for HPF1-PARP interaction, while these key residues are largely conserved between PARP1 and PARP2 (Figure 2c, Supplemental Figure S6). This explains why HPF1 can bind to PARP1 and PARP2, but not to PARP3 23.
In examining the structure, we noted another interface between a different HPF1 and PARP1 molecule in the crystal lattice (see interface II in Figure 2a, top-left inset). However, mutating key residues on this interface (Glu138, Phe139 and Lys216) did not interfere with the binding (Figure 2b, see lanes 5, 6, 7; Table 1, Supplemental Figure S2, S3), and these residues are not conserved in HPF1 (Supplemental Figure S5). Thus this contact does not appear to be relevant to the catalytic function of the HPF1/PARP1 complex.
We also determined structures of human and mouse HPF1 alone, and that the conformation does not change upon formation of the complex with PARP1 (see Supplemental Figure S7). Notably, analysis of the surface electrostatic potential showed that positively and negatively charged areas are evenly distributed on the surface of this protein, except for part of helix α9/α10 and the linker loops connecting helices α6 and α7, where the protein surface is dominantly negatively charged (Supplemental Figure S7B).
HPF1 binding remodels the PARP1 active site for histone binding and serine ADPr
Remarkably, our functional analysis of the structure shows that HPF1 and PARP1 form a composite active site, with HPF1 contributing key active site residues.
We first analyzed the surface electrostatic potential of HPF1, PARP1-CAT ΔHD and the complex (Figure 3a). Though PARP1 has been shown to play a key role in DDR involving the modification of histones 16, 28, the surface of ART subdomain, especially the active site, is positively charged, which is clearly incompatible with the positively charged nucleosome (Supplemental Figure S8). Interestingly, HPF1 binding to PARP1-CAT dramatically changes the situation. As is mentioned above, HPF1 surface covering helix α9, α10 and the loops connecting helices α6 and α7 is negatively charged (Supplemental Figure S7). Many solvent-exposed acidic amino-acid residues are located in this region, including Asp235, Glu240, Glu243, Asp245, Glu273, Asp283, Glu284, Asp286 and Glu292, and Arg239 turns out to be the only basic residue in this region (Figure 3a). Notably, Arg239 and most of these acidic residues are conserved in HPF1 from different species (Supplemental Figure S5). When binding to PARP1, this negatively charged region of HPF1 merges to the active site of PARP1-ART and changes the surface potential of the latter, creating an overall negatively charged joint active site that would accommodate the access of the positively charged histones (Figure 3a). The surface of the HD subdomain facing the active site is also negatively charged due to three acidic residues Glu756, Asp766 and Asp770 and lack of any basic residue here, hence seems to form a negatively charged narrow tunnel with the joint active site (Figure 3b). However, the HD subdomain is believed to undergo partial unfolding upon PARP1 activation 19, and therefore it is unknown whether the tunnel still maintains. We speculate that if the HD subdomain does not completely unfold, the negatively charged tunnel may still maintain (but likely undergoes some conformational changes) that is presumably most accessible for positively charged extended peptide substrate (say histone H3 N-terminal tail), while the folded and/or negatively charged proteins may not readily access this joint active site to get ADP-ribosylated.
A benzamide molecule was found to bind in the active center of the HPF1/PARP1-CAT ΔHD complex structure (Supplemental Figure S4). When superimposed with the previously reported PARP1-CAT ΔHD/BAD complex crystal structure (PDB 6BHV) 22, the benzamide molecule in our HPF1/PARP1-CAT ΔHD complex structure well overlaps the benzamide moiety of BAD, the NAD+ analog (Figure 3c). Importantly, we noted that the HPF1 Glu284 side-chain carboxyl, being positioned/locked by the Arg239-Glu284 salt-bridge, is located beside the ribose of the ADP moiety of carba-NAD+ (representing the acceptor for ADP-ribosylation 29, see Figure 3d), and is about 1.9Å farther away than the latter from the C1” atom of the nicotinamide ribose moiety of BAD/NAD+ (representing the donor for ADP-ribosylation, see Figure 3c), which from the structural biological view would strongly suggest an important functional role of this residue rather as a catalyzer for ADP-ribose transfer to a third party (say a serine residue) than as a direct acceptor itself for ADP-ribose. In NAD+, the positively charged nicotinamide nitrogen trends to deprive an electron from the connecting ribose C1” carbon, rendering the latter an oxacarbenium ion that would then be attacked by a nucleophile such as a negatively charged/deprotonated amino-acid side-chain. Since acidic amino-acids glutamate and aspartate may spontaneously lose a proton to become negatively charged at physiological pH, these residues can be readily ADP-ribosylated once they access the active center of ART (for example, in the absence of HPF1, see Figure 3e). However, in the HPF1/PARP1-CAT ΔHD complex structure the Glu284 side-chain carboxyl is too far away (~4.6Å) from the NAD+ to attack the ribose C1” carbon (Figure 3c). Therefore it seems unlikely that Glu284 itself can get ADP-ribosylated. On the other hand, since the non-acidic serine residue does not spontaneously lose a proton at physiological pH, it must be activated (deprotonated) by a third party prior to ADP-ribosylation. The Glu284 side-chain carboxyl may fulfill this task when it becomes negatively charged (Figure 3f). Indeed, mutating HPF1 Glu 284 to Alanine (E284A) completely abolished ADP-ribosylation of histones (Figure 2b, lane 13). The same observation and proposal of function have been made by Suskiewicz et al. in their structural studies of HPF1/PARP2-CAT ΔHD structure and functional studies including PARP1 and PARP2 27.
Why does HPF1 binding hinder PARP1 automodification and poly-ADP-ribosylation?
Based on the verified interaction mode between HPF1 and PARP1-CAT ΔHD (corresponding to interface I shown in Figure 2a), an overall model of HPF1 binding to full-length activated PARP1 except the Zn2 and AD domains is proposed through overlapping the PARP1-ART domain with the corresponding domain in the PARP1/DNA complex structure (PDB 4DQY) (Figure 4a) 18. HPF1 binds to ART at the active site of the enzyme, while in the meanwhile it is well accommodated with other domains of PARP1 except for HD that undergoes local unfolding upon PARP1 activation. This model immediately hints why HPF1 binding restricts automodification of PARP1. According to the previous model proposed by Langelier et al., the AD/BRCT domain carrying most of the automodification sites is located close to the CAT domain presumably for being ADP-ribosylated 18 (Supplemental Figure S1). This location is roughly the same place where HPF1 binds (compare Figure 4a to Supplemental Figure S1). We therefore speculate that HPF1 binding would dislodge AD from this location, hence limits the automodification of PARP1, at least in the folded BRCT domain.
The complex structure hints that HPF1 not only restricts hyper-automodification of PARP1, but may limit the modification from poly-ADP-ribosylation to mono-ADP-ribosylation. An early study by Ruf et al. observed the ADP moiety of carba-NAD+ binding in the active site, and this ADP moiety was believed to represent the acceptor for the next ADP-ribose during poly-ADP-ribose chain elongation and branching 29. Comparing to the HPF1/PARP1-CAT ΔHD complex structure reported here, we noted that the ADP (i.e. the acceptor for next ADP-ribose in poly-ADP-ribosylation) binding site is partly occupied by HPF1 residues Asp283 and Glu284 (Figure 3d), which would implicate that HPF1 binding is likely to prevent poly-ADP-ribose chain elongation and branching. However, we speculate that poly-ADP-ribosylation may still take place, likely at lower probability, due to the relatively weak and thus transient binding between HPF1 and PARP1. This explains the long smeared bands observed in our ADP-ribosylation assay (Figure 2b, see histone H3/H4ADPr).
The unique role of HPF1 Arg239
Arg239 is highly conserved in HPF1 (Supplemental Figure S5), and has been thought to be important for the interaction between HPF1 and PARPs. Mutating this residue to alanine (R239A) did restore automodification of PARP1 (Figure 2b, lane 8), which was previously interpreted to be loss of binding between HPF1 and PARP1 23.
However, in the human HPF1/PARP1-CAT ΔHD and HPF1/PARP2-CAT ΔHD complex structures, HPF1 Arg239 was found to only interact with Glu284 and Asp286, but not with any residue from PARP1/2 (Figure 4b). Furthermore, ITC assays confirmed that this mutation only mildly weakens the binding affinity between HPF1 and PARP1, indicating that without Arg239 HPF1 still binds pretty well to PARP1 (Table 1). Curiously, in the ADP-ribosylation assays, the band of HPF1 R239A, but not of any other HPF1 mutants, clearly shifted on the SDS-PAGE (Figure 2b, lane 8, black arrow), which looks very likely that HPF1 R239A itself was ADP-ribosylated. We then analyzed the shifted HPF1 R239A band by Mass Spectrometry, which confirmed that HPF1 R239A did get ADP-ribosylated on Asp235 and Glu240, two acidic residues located on both sides of residue 239 in the long loop connecting helices α6/α7 (Figure 4b, 4c). The long loop regions are often difficult to resolve in crystal structure due to their intrinsic flexibility, but in HPF1 the long loop connecting helices α6 and α7 (residues 217-244) is ordered and well resolved. As is shown in Figure 4b, the C-terminal half of this long loop is stabilized by two groups of interactions: one is the salt-bridge network including Arg239-Asp286 and Arg239-Glu284, which also acts to position the catalytically important Glu284 side-chain; the other is the hydrophobic interaction between Tyr238 phenyl and Val218 side-chain, and hydrogen bond between Tyr238 hydroxyl and Glu292 carboxyl. It is inferred that when Arg239 was mutated to alanine that lost the interlocking with Glu284 and Asp286, the loop became so flexible that Asp235 and Glu240 gained the freedom to attack NAD+ in the active center and became ADP-ribosylated. Importantly, although HPF1 R239A still bound PARP1 pretty well, it also restored the automodification of PARP1 (Figure 2b, lane 8). The PARP1 automodification domain (AD) contains a folded subdomain BRCT (residues 380-480) and an unstructured automodification peptide (residues 481-530). In a previous in vitro ADP-ribosylation assay, the serine residues in the Lys-Ser motif within the automodification peptide (Ser499, Ser507and Ser519), but not that in the BRCT subdomain (Lys467-Ser468), were ADP-ribosylated 24 in the presence of HPF1, indicating that the Lys-Ser motif in a flexible peptide, but not in a folded domain may access the remodeled HPF1/PARP1 active center under certain circumstances. In our assays, HPF1 WT binding robustly limited automodification of PARP1, while R239A, although still bound to PARP1, restored the automodification. Taken together, our data indicate that the rigidity of the loop region strengthened by Arg239-Glu284/Arg239-Asp286 interactions is essential for limiting PARP1 automodification in the unstructured automodification peptide region.
In the ADP ribosylation assay we also noted that the HPF1 E284A band mysteriously became weaker (Figure 2b, lane 13, gray arrow), indicating losing a small amount of the un-modified protein (or ADP-ribosylation of a small portion of the protein). However, the modification of this mutant HPF1 protein must be very minor since no shifted band was obviously seen on the gel. We speculate that in this mutant, after losing the Arg239-Glu284 salt-bridge but still maintaining the Arg239-Asp286 salt-bridge, the loop region gained some flexibility but the acidic residues (say Asp239 and Glu240) did not gain enough freedom to attack NAD+ to be ADP-ribosylated. Accordingly, HPF1 E284A only mildly restored PARP1auto-modification, which again supported the notion that the rigidity of this loop region strengthened by Arg239-Glu284/Arg239-Asp286 interactions is essential for limiting PARP1 automodification.
Why is HPF1 Arg239 important for HPF1/PARP1 binding?
The above observations do not explain why mutating Arg239 to alanine weakens the interaction between HPF1 andPARP1, since no residue in the loop region (including Arg239 itself) directly interacts with PARP1. We noted that the HPF1 E284A mutation showed a very unique behavior in our assays. It is the only mutation that enhanced the binding between HPF1 and PARP1 (Table 1). Taken into account the negative electrostatic property of the joint active site discussed above, we propose that the existence of the negatively charged Glu284 in the active center is not favorable for HPF1/PARP1 binding, while the salt-bridge with the positively charged Arg239 side-chain may act to neutralize the negative charge on Glu284 carboxyl, thus facilitate the binding. Mutating Arg239 to directly break Arg239-Glu284 interaction, or mutating Tyr238 to destroy the local conformation supporting Arg239-Glu284 interaction, would restore the negative charge of Glu284, thus weaken HPF1/PARP1 binding. This explains why HPF1 R239A (and/or Y238A, as was found in the previous study 23) mutation weakens the binding between the two proteins, and indicates that although mutating Arg239 (and/or Tyr238) may apparently weaken HPF1/PARP1 binding, they may not be the right choices to act as the non-binding control in functional studies because the mutation(s) lead to much more complicated consequences rather than simply abolishing the binding – and, they do not really abolish the binding.
Taken together, our study indicates that Arg239 is a key residue acting to 1) position Glu284 at the right location for catalyzing serine ADP-ribosylation; 2) stabilize the conformation of the long loop region across Arg239 to limit auto-modification of PARP1; 3) facilitate HPF1/PARP1 binding through neutralizing the negative charge of Glu284 side-chain.