Structural and dynamical insights into the PH domain of p62 in human TFIIH

Abstract TFIIH is a crucial transcription and DNA repair factor consisting of the seven-subunit core. The core subunit p62 contains a pleckstrin homology domain (PH-D), which is essential for locating TFIIH at transcription initiation and DNA damage sites, and two BSD (BTF2-like transcription factors, synapse-associated proteins and DOS2-like proteins) domains. A recent cryo-electron microscopy (cryo-EM) structure of human TFIIH visualized most parts of core, except for the PH-D. Here, by nuclear magnetic resonance spectroscopy we have established the solution structure of human p62 PH-D connected to the BSD1 domain by a highly flexible linker, suggesting the flexibility of PH-D in TFIIH. Based on this dynamic character, the PH-D was modeled in the cryo-EM structure to obtain the whole human TFIIH core structure, which indicates that the PH-D moves around the surface of core with a specific but limited spatial distribution; these dynamic structures were refined by molecular dynamics (MD) simulations. Furthermore, we built models, also refined by MD simulations, of TFIIH in complex with five p62-binding partners, including transcription factors TFIIEα, p53 and DP1, and nucleotide excision repair factors XPC and UVSSA. The models explain why the PH-D is crucially targeted by these factors, which use their intrinsically disordered acidic regions for TFIIH recruitment.

p62 is an essential subunit of TFIIH, not only for formation of the core TFIIH complex but also for recruiting TFIIH to appropriate functional sites. Human p62 has four structural domains: an N-terminal pleckstrin homology (PH) domain (PH-D) (residues 1-103); two BSD (BTF2-like transcription factors, synapse-associated proteins and DOS2-like proteins) domains, BSD1 (residues 109-147) and BSD2 (residues 189-225); and a C-terminal three-helix bundle (residues 454-548). The PH-D is responsible for recruitment of TFIIH through its specific interaction with intrinsically disordered acidic regions of various transcription factors and several nucleotide excision repair factors, such as TFIIE␣ (28,29), p53 (30,31), EKLF (32), DP1 (33), VP16 (34,35) and EBNA2 (36) during transcription, and XPC (37,38), UVSSA (39) and XPG (40) during nucleotide excision repair. Recent cryo-electron microscopy (cryo-EM) studies of human apo TFIIH (41) and yeast holo TFIIH in complex with the pre-initiation complex (PIC) and Mediator (42) have provided detailed pictures of the core TFIIH complex and how it is connected to the CAK complex. In the cryo-EM structure of human TFIIH, al-most all parts of the seven subunits in the core can be visualized except for the N-terminal PH-D of p62.
To gain insight into the structure of the p62 PH-D in the core TFIIH complex, here, keeping in mind the fact that the BSD1 domain can be clearly identified in the cryo-EM structure, we have solved the solution structure of the tandem PH and BSD1 domains of human p62 (residues 1-158) by nuclear magnetic resonance (NMR) spectroscopy. These domains are connected by a short linker (residues 104-108) and do not interact with each other. The flexibility of the interdomain region has been evaluated by measuring 15 N relaxation parameters and 1 H exchange rates with solvent. Based on our dynamic NMR structures of the tandem PH and BSD1 domains coupled with the cryo-EM structure (43), we have built reasonable dynamic structural models of human TFIIH including the flexible p62 PH-D by superposing the BSD1 domains from the solution and cryo-EM structures; we have further refined these dynamic structures by using all-atom molecular dynamics (MD) simulations. In addition, we have modeled the structures, also refined by the MD simulations, of TFIIH in complex with five p62binding partners that are responsible for recruiting TFIIH to sites of transcription initiation and DNA damage.

Preparation of human p62 1-158
The 13 C/ 15 N-labeled fragment of human p62 1-158 (residues 1-158) was expressed as a hexa-histidine-tagged product in a pET15b vector (Merck Millipore) in Escherichia coli BL21 (DE3) Gold (Agilent Technologies). The cells were grown at 37 • C in M9 minimal medium containing [ 15 N]ammonium chloride and [ 13 C]glucose. The product was expressed by induction with 1 mM isopropyl-␤-D-thiogalactopyranoside. After 21 h of growth at 20 • C, the cells were collected, resuspended in buffer A [20 mM Tris-HCl (pH 8.0), 10% glycerol, 1 M NaCl], lysed by sonication and then centrifuged. The supernatant was loaded onto a His-Tag Purification Resin (Roche) column equilibrated with buffer A, and the column was washed with buffer A containing 10 mM imidazole-HCl. The sample was eluted by 500 mM imidazole-HCl, peak fractions were pooled and the buffer was changed to thrombin cleavage buffer (10 mM Na 2 HPO 4 , 1.8 mM KH 2 PO 4 , 500 mM NaCl, 2.7 mM KCl, pH 7.3) by using an Amicon Ultra device (Merck Millipore). To remove the histidine tag, the sample was digested with thrombin for 63 h at 20 • C. After concentration via an Amicon Ultra, the sample was applied to a Superdex75 (GE Healthcare) column equilibrated with 20 mM potassium phosphate (pH 6.8) and 500 mM NaCl.

NMR structure determination
For structure determination, we used 0.9-1.2 mM 13 C/ 15 Nlabeled p62 1-158 in 20 mM potassium phosphate (pH 6.8) and 5 mM deuterated DTT prepared in either 90% H 2 O/10% D 2 O or 99.9% D 2 O. NMR experiments were performed at 20 • C on Bruker AVANCE III HD 600-MHz and 950-MHz spectrometers, each equipped with a cryo-TCI probe. Backbone and side-chain resonances were assigned by using standard triple-resonance NMR experiments (44).

Structure calculation
Interproton distance restraints derived from NOE (nuclear Overhauser effect) intensities were grouped into four distance ranges, 1.8-2.7Å (1.8-2.9Å for NOEs involving HN protons), 1.8-3.3Å (1.8-3.5Å for NOEs involving HN protons), 1.8-5.0 and 1.8-6.0Å, corresponding to strong, medium, weak and very weak NOEs, respectively. The upper limit was corrected for constraints involving methyl groups, aromatic ring protons and non-stereospecifically assigned methylene protons. Dihedral angle restraints for ϕ and ψ were obtained from analysis of the backbone chemical shifts with TALOS+ (47). χ 1 and χ 2 angles were restrained ±30 • for three side-chain rotamers. Structure calculations were performed by distance geometry and simulated annealing using the program Xplor-NIH (48,49). In total, 100 structures were calculated. All structures were then subjected to water refinement (50), in which the structures were immersed in a 7.0-Å layer of water molecules. After minimization with 120 steps, a heating stage from 100 to 500 K with 200 steps of MD for every 100-K increment, a refinement stage with 2500 steps at 500 K and a cooling stage from 500 to 25 K with 200 steps for every 25-K decrement were carried out. The refinement protocol was finished with 200 steps of minimization. Structural statistics for the 20 best structures are summarized in Table 1. Structures were analyzed and displayed by using PROCHECK-NMR (51), MOLMOL (52) and PyMol (http://www.pymol.org).

NMR relaxation analysis
The backbone 15 N relaxation parameters were measured by using a 13  the presence and absence of proton saturation, which was achieved with a 5.0-s duration consisting of 120 • 1 H pulses applied every 5.0 ms. R 1 and R 2 were obtained by fitting peak intensities at a series of relaxation durations to an exponential decay curve using CurveFit (53). The uncertainties of the peak intensities were estimated by using duplicated data from the shortest relaxation delay. The uncertainties of the relaxation rates were determined by CurveFit using a Monte Carlo simulation. The steady-state 15 N-{ 1 H} NOE values were determined from peak intensity ratios obtained from spectra acquired with and without proton saturation. The uncertainties were determined from the standard deviation in background noise levels by using NMRView (46).

HET ex -BEST-TROSY experiment
HET ex -BEST-TROSY (54) experiments were performed at 20 • C on a Bruker AVANCE III HD 950-MHz spectrometer equipped with a cryo-TCI probe. A series of spectra was collected with relaxation delays of 339, 709, 1080, 1450 and 2006 ms with and without water saturation. The water signal intensity was measured at each relaxation delay to determine the polarization of water using the 1 H-15 N BEST-TROSY pulse sequence with an additional small flip angle reading pulse (0.1 s). The hydrogen exchange rate k ex was determined by simple fitting using Equations (1) and (2) for the amide proton polarizations: where R A 1 is the longitudinal relaxation rate of amide protons, A eq z is the equilibrium amide proton polarization and W ref z is the water polarization. All NMR spectra were processed by using the program NMRPipe (45) and analyzed by using the program NMRView (46). Least-squares fitting of the data was performed by using Gnuplot software (http://www.gnuplot.info).

CLEANEX-PM experiment
CLEANEX-PM (phase-modulated CLEAN chemical exchange spectroscopy) (55) experiments were performed at 20 • C on a Bruker AVANCE III HD 950-MHz spectrometer equipped with a cryo-TCI probe using a standard Bruker pulse program. A series of spectra was collected with mixing times for the CLEANEX spin lock of 7.5, 15, 30, 50 and 80 ms. The hydrogen exchange rate k ex was determined by using established methods. All NMR spectra were processed by using the program NMRPipe (45) and analyzed by using the program NMRView (46). Least-squares fitting of the data was performed by using Gnuplot software (http://www.gnuplot.info).

Structural modeling of human TFIIH
Initially, 1000 structures of p62 1-158 were calculated as described in the 'Structure calculation' section, and then 399 structures were selected on the basis of violation criteria (distance violations <0.5Å and dihedral angle violations <5 • ). These 399 structures were aligned to the BSD1 domain region (residues 107-147) of p62 in the cryo-EM structure of human apo TFIIH (PDB code 6NMI) by using the program PyMol (http://www.pymol.org). The BSD1 domain region (residues 107-147) of p62 in the cryo-EM structure was then replaced with the PH-D and BSD1 domain (residues 1-147) of p62 in the 399 structures of p62  . Finally, undesirable models in which the PH-D of p62 clashed against other parts of TFIIH were removed, resulting in 275 model structures.

MD simulations of human TFIIH and its complexes with p62binding partners
The initial structure used for MD simulations was prepared from the whole apo human TFIIH core complex, in which the NMR structure was used for the PH-D and BSD1 domain (residues 1-147) of p62 and the cryo-EM structure (PDB code 6NMI) was used for the other parts of TFIIH. The 275 model structures described above were classified in terms of the position of the PH-D, and 6 representative model structures ( Figure 5B, D, F, H, J and L) were subjected to MD simulations. First, in a typical structure (Figure 5L), the missing regions and side chains in the cryo-EM structure, as well as the region joining the NMR and cryo-EM structures of p62 (residues 145-147), were built by MODELLER (56). The long loops at the termini (residues 1-31 and 731-782 of XPB, residues 1-14 of p44, residues 294-308 of p34, residues 69-71 of p8 and residues 213-309 of MAT1) were not modeled. Next, using the model structure in Figure 5L, the other five whole TFIIH structures ( Figure 5B, D, F, H and J) and five TFIIH complexes (Figure 6) were built: the whole TFIIH structures were aligned to the model structure by using a helix of the three-helix bundle of p62 (residues 487-511), and the PH-D and BSD1 domain (residues 1-143) of each structure were replaced with the same domains of the model structure.
The systems for MD simulations of the model TFIIH structures were prepared by using the psfgen, solvate and autoionize plugins implemented in VMD (57). The protonation states of His, Asp, Glu and Lys were determined by using pK a values calculated with PROPKA (58,59) implemented in PDB2PQR (60) under the pH 7 condition. Missing hydrogen atoms were inserted in accordance with topology in the CHARMM36m force field (61-63). The N-and C-termini were set to NH 3 + and COO − , respectively. The coordinate bonds among anionic cysteine or histidine residues and zinc ions were incorporated, as defined in the CHARMM36m force field, and the oxidized iron-sulfur cluster and the coordinate bonds of anionic cysteine residues with the iron-sulfur cluster were also incorporated by using the topologies and parameters determined by Chang and Kim (64). Here, both phosphorylated serine and threonine in the p62-binding partner p53 were set to the dianionic state. The MD unit cell was set to a rectangular cell, and the cell size was determined such that the distance between the edge of the cell and the centered model structure was at least 15Å. Using the TIP3P water model (65), the cell was filled with water molecules, and counterions (Na + ) and 150 mM NaCl were added. The final systems consisted of ∼504 000 atoms.
The all-atom MD simulations were performed by using the MD program package GROMACS version 2019.6 (66-68) with the CHARMM36m force field under periodic boundary conditions. The topology and parameter files were converted into the GROMACS format by using ParmEd (https://github.com/ParmEd/ParmEd). Electrostatic interactions were handled by the smooth particle mesh Ewald method (69), and the van der Waals interactions were truncated by the switching function with a range of 10-12Å. Bond lengths involving hydrogen atoms were constrained by the P-LINKS algorithm (70). In all simulations, position constraints at the C␣ atoms with a 60 kJ mol −2 nm −2 force constant were applied to the TFIIH structure, except for the PH-D, the interdomain linker and the Nterminal part of the BSD1 domain (residues 1-114) of p62, the modeled loops that were originally missing in the cryo-EM structure, and the p62-binding partners. The temperature and pressure were 298 K and 1 atm, respectively. Before the production run, energy minimization by the steepest descent method and three equilibration runs were performed sequentially: the first 125-ps run with NVT ensemble and a 1-fs timestep, and the second and third 100-ps runs with NPT ensemble and a 2-fs timestep. The thermostat in the first and third equilibration runs and the production run was the Nosé-Hoover scheme (71,72), and that in the second equilibration run was the Berendsen algorithm (73). The barostat in the second equilibration run was the Berendsen algorithm (73), and that in the third equilibration run and the production run was the Parrinello-Rahman approach (74,75). After the equilibration runs, a 100-ns production run was carried out with NPT ensemble and a 2-fs timestep. These MD simulations were performed for each initial model structure of the six whole TFIIH models and the five complex models with p62-binding partners (1.1 s in aggregation of the production runs).
The root-mean-square fluctuation (RMSF) of C␣ atoms, the solvent-accessible surface area (SASA) of the amide hydrogen and the protein secondary structures were calculated by using each 100-ns trajectory. Snapshots were saved every 1 ns (100 snapshots in a trajectory). As an estimate of the RMSF, an averaged structure was obtained by structural alignment of the snapshots using C␣ atoms in the PH-D except for the N-terminus (residues 9-100). The RMSF values of the six apo TFIIH simulations or the five complex simulations were then averaged. The SASA of the amide hydro-gen was calculated by the GROMACS module. The average SASA of each residue was obtained from each trajectory, and then the averaged SASA was calculated from the simulations for each apo TFIIH or each complex. Protein secondary structures were determined by using the DSSP program (76) implemented in CPPTRJ (77). The population of snapshots forming the helix or 3 10 helix structure was obtained for each trajectory, and its averaged value was calculated.

Invisible PH-D of p62 in the cryo-EM structure of human TFIIH
Understanding the structure and dynamics of the PH-D of p62 in TFIIH is indispensable for elucidating the mechanism of recruitment of human TFIIH. At present, highresolution cryo-EM structures of human and yeast TFIIH are available. The structure of human TFIIH has been resolved for an apo form (41,43), while that of yeast TFIIH has been resolved for a holo form in the transcription PIC (42). The corresponding structures of each subunit and the overall architecture of the core TFIIH are relatively similar, except for the corresponding PH-Ds of the yeast homolog Tfb1 (Saccharomyces cerevisiae) and human p62: the Tfb1 PH-D is clearly visible, whereas the human p62 PH-D is invisible (Supplementary Figure S1). Recently, an interesting structural study of human PIC has been reported in which all available cryo-EM data have been integrated to build a structural model of PIC in which the PH-D of p62 is visible (78).
Here, we initially compared the amino acid sequence of human p62 with that of its well-studied homolog Tfb1 (Figure 1A and Supplementary Figure S2). The sequence alignment identified a marked difference in the region between the PH and BSD1 domains. In human p62, the PH-D is immediately followed by the BSD1 domain; in yeast Tfb1, in contrast, a linker of ∼50 residues is present between these domains. This substantial sequence dissimilarity is reflected in their cryo-EM structures ( Figure 1B). In both p62 and Tfb1, the BSD1 and BSD2 domains are aligned at nearly the same position with a beads-on-a-string-like topology, and contribute to stabilizing the core architecture by interacting with other four subunits: XPD, p52, p44 and p34 for p62; and Rad3, Tfb2, Ssl1 and Tfb4 for Tfb1. The PH-D of Tfb1 is also visualized in the yeast holo TFIIH structure ( Figure  1B, right). Its PH-D lies apart from the BSD1 domain owing to its long linker, and interacts with another general transcription factor, Tfa1/TFIIE␣ (not shown in Figure 1B). In contrast, the PH-D of p62 is disordered and invisible in the human apo TFIIH structure ( Figure 1B, left), suggesting that, despite the short interdomain linker, this domain is surface exposed and behaves flexibly in the absence of proteins other than TFIIH. The possible state of p62 in human holo TFIIH is discussed in detail below.
To explore the predicted flexible behavior of the PH-D of p62 in human apo TFIIH, we applied solution NMR because it is a powerful tool for investigating both the dynamics and structure of proteins in solution. Although NMR is not advantageous for investigating large proteins like TFIIH, the PH-D of p62 is likely to be free in apo TFIIH; we considered therefore that we would not need to target the whole TFIIH protein for our purpose. In the cryo-EM structure, the p62 subunit is visualized from Ala107 onward, just after the PH-D and before the BSD1 domain ( Figure 1B, left). Considering this, we focused on the region including both the PH and BSD1 domains (residues 1-158, hereafter designated p62 1-158 ) (Figure 2A), and analyzed the secondary structure of p62 1-158 by NMR.
The p62 1-158 construct exhibited well-dispersed signals on a 1 H-15 N-HSQC spectrum ( Figure 2B). Almost all backbone and side-chain signals were assigned without using any deuterated samples. The chemical shift data (CSI) (79) of the backbone atoms indicated that there were seven ␤strands and one ␣-helix in the PH-D and three ␣-helices in the BSD1 domain, consistent with the structures of their isolated forms (40). No secondary structure elements were found in the interdomain linker region between the ␣1 and ␣2 helices ( Figure 2C). Based on the chemical shift data, it is predicted (RCI) (80) that the interdomain linker region is flexible ( Figure 2D).

Structure of the tandem PH and BSD1 domains of human p62
Next, we determined the tertiary structure of p62 1-158 by using NMR (Figure 3). In total, 3338 NOE-derived distance restraints, 96 hydrogen bond restraints and 381 dihedral angle restraints were collected and used to define the structure ( Figure 3A-D and Table 1). Although the RMSD of the backbone atoms for residues 7-147 in the 20 structures was as high as 8.81 ± 2.55Å, each domain structure was well determined with RMSD values of 0.46 ± 0.10Å for residues 7-103 (PH-D) and 0.76 ± 0.34Å for residues 109-147 (BSD1), suggesting that the interdomain linker (104-108) has high mobility ( Table 1). In the PH-D, the seven ␤strands (␤1-␤7) form a ␤-sandwich, in which strands ␤1-␤4 and ␤5-␤7 form the first and second antiparallel ␤-sheets, respectively. These two ␤-sheets are linked by a contact between strands ␤1 and ␤7, and by a 3 10 helix connecting strands ␤4 and ␤5 ( Figure 3E). The following ␣-helix (␣1) is inserted between the two ␤-sheets and extensively interacts with residues in both, thereby stabilizing the ␤-sandwich. Pro101 is likely to participate in the ␣-helix because Phe103 stabilizes the helix by inserting into the hydrophobic core. The structure of the PH-D in p62 1-158 is essentially identical to that in the isolated form (40). Likewise, the structure of the BSD1 domain in p62 1-158 is essentially the same as that in the isolated form, which has been determined by another research group (PDB code 2DII) although the study has not been published.
Here, we describe the structure in a little more detail. The first helix (␣2) of the BSD1 domain bends nearly at a right angle to the second helix (␣3) at Pro122 ( Figure 3E). After the turn (Val133-Ser138), the third helix (␣4) runs antiparallel to helix ␣3 and perpendicularly to helix ␣2. These three helices interact with one another, rendering the BSD1 domain stable and compact. Lys114, Met117 and Leu118 in ␣2, Val123, Leu124, Phe125, Leu127, Tyr128, Leu131 and Val132 in ␣3, and Ala139, Phe142, Trp143, Arg146 and Leu147 in ␣4 form the hydrophobic core. The turn between helices ␣3 and ␣4, where Val136 and Ile137 closely interact with Leu131, looks rigid. Ile137 also contacts Val132 and Phe142, participating in the formation of the hydrophobic core.

Structural model of human TFIIH
Our results enabled us to complete a structural model of the whole human apo TFIIH core complex by first combining the cryo-EM structure with our NMR structure of p62  and then considering the dynamics of p62. To build the structural model, initially 399 NMR structures of p62  were chosen from the total 1000 calculated structures on the basis of violation criteria. The structure of the BSD1 domain in p62 1-158 is essentially identical to that of the isolated form (PDB code 2DII), which was used to build p62 in the cryo-EM structure of human apo TFIIH, enabling the BSD1 domain of p62  to be placed at the same position in the cryo-EM structure. Hence, the BSD1 domain in the cryo-EM structure could be replaced with the PH and BSD1 domains of our 399 NMR structures of p62  . Ultimately, we selected 275 models in which the PH-D of p62 did not clash against other parts of TFIIH. Of note, nearly 70% (275/399) of the models had no steric clashes.
The structural models allow us to visualize the dynamic behavior of the PH-D of p62 in apo TFIIH ( Figure 5A). The PH-D can move close to both the BSD2 domain ( Figure 5B) and the three-helix bundle ( Figure 5D), as well as the BSD1 domain ( Figure 5F). In addition to its own p62 subunit, it can move close to the C-terminal RecA-like (RecA2) domain of XPB ( Figure 5H) and the von Willebrand factor A domain of p44 ( Figure 5J). This spatially proximity is in accordance with the results of a chemical cross-linking/mass spectrometry study that detected cross-links between Lys38 of p62 and Lys549 of XPB, and between Lys38 of p62 and Lys126 of p44 (81). Notably, in the majority of our models, the PH-D occupies free space ( Figure 5L).
To refine the docking structural models, we performed all-atom MD simulations of the apo whole TFIIH structure. The MD simulations for six representative structural models ( Figure 5B, D, F, H, J and L) were carried out in solution, and the snapshot structures obtained at 100 ns are illustrated in Figure 5C, E, G, I, K and M. In the MD simulations, the PH-D, the interdomain linker and the N-terminal part of BSD1 domain of p62 (residues 1-114) were freely moving, whereas the C␣ atoms in the other part were restrained to the cryo-EM structure (see the 'Materials and Methods' section). The structures of the PH-D itself were retained (Supplementary Figure S4A), but their positions and orientations changed from the initial structure. In MD simulations of the model structures in Figure 5B, D and F, the PH-D fluctuated around its initial position because it interacted with the BSD2 domain ( Figure 5C), the three-helix bundle ( Figure 5E) or both the BSD1 domain and the threehelix bundle ( Figure 5G) of p62. In the MD simulations of the model structures in Figure 5H, J and L, the PH-D moved from its initial position largely because of the weak interactions that held it in its initial place ( Figure 5I, K and M).
This dynamic behavior of the PH-D originates from the high mobility of the interdomain linker (residues 104-108) between the PH-D and BSD1 domain. To estimate the interdomain linker mobility, the RMSF of C␣ atoms was calculated from the MD trajectories ( Supplementary Figure S5A). The fluctuation in the PH-D except for the Nterminus was up to ∼2Å. In contrast, the RMSF of the interdomain linker region was markedly higher, indicating that the conformation of the linker is flexible. We found that the flexibility of the interdomain linker region is related to the SASA of the amide hydrogen in backbones (Supplementary Figure S6A). The amide hydrogens in the linker region are exposed to solvent. The dynamic behavior of the PH-D and the mobility of the interdomain linker region observed in MD simulations are in good agreement with our NMR measurements.
Thus, the PH-D of p62 is able to move dynamically around the surface of the core TFIIH complex, transiently interacting with some regions; however, the spatial distribution is limited to some extent owing to the short interdomain linker. As mentioned above, the structural model of human PIC has been reported (78). In that model, the PH-D of p62 is unambiguously positioned between XPB and XPD (Supplementary Figure S7A Figure S7B, right). During formation of the PIC at a transcription start site, the PH and BSD1 domains extensively interact with another transcription factor, TFIIE␣ (28,29,78). It is not clear whether the partial helix unfolding of the PH-D and BSD1 domain is a result of binding to TFIIE␣. Further studies are needed to explore the possibility.

Complex structural models of human TFIIH and p62-binding partners
We have previously determined the solution structures of the PH-D of p62 in complex with five binding partners: the acidic domain of TFIIE␣ (29); the transactivation domains of the tumor suppressor p53 (31) and the cell cycle regulatory transcription factor DP1 (33); and the acidic regions of the nucleotide excision repair factors XPC (38) and UVSSA (39). All of these proteins bind specifically to the PH-D of p62 via an intrinsically disordered acidic region. Upon binding, this region forms an extended string-like conformation and interacts with a basic surface of the PH-D via extensive electrostatic contacts. Simultaneously, it inserts a phenylalanine or tryptophan residue surrounded by acidic residues into a pocket in the PH-D, enhancing bind-ing specificity. In those studies, we used the isolated form of the PH-D, consisting of residues 1-108 (p62 1-108 ). The current study now demonstrates that this construct design is appropriate because Asn108 is located at the C-terminus of the interdomain linker. Nevertheless, it has been suggested that some binding partners interact with several residues in the C-terminal region of the p62 1-108 construct (29,38,39). Here, based on the complete structural model of the human apo TFIIH core complex ( Figure 5), we have explored whether the p62-binding partners would be able to bind to the PH-D of p62 in TFIIH in the same manner as they bind to p62 1-108 .
To build structural models of TFIIH docked to its p62binding partners, we selected five structures with sufficient variety from the 275 structural models of TFIIH ( Figure  5A). In almost all calculated models, TFIIE␣, p53, DP1, XPC and UVSSA could bind to the PH-D of p62 in TFIIH without steric clashes except for the TFIIE␣-TFIIH complex, in which the C-terminus of the acidic domain of TFIIE␣ collided with the N-terminus of the BSD1 domain of p62 (data not shown).
To refine the docking models of the TFIIH with its p62binding partners, MD simulations were performed for the complex structures. The interactions between the PH-D and the p62-binding partner were maintained for all complexes, and none of the partners dissociated during the simulations ( Figure 6). The relative positions and orientations of the PH-D and the binding partner changed with respect to the whole TFIIH, owing to the flexibility of the linker region of p62. As observed for apo TFIIH, the RMSF (Supplementary Figure S5B) and SASA of the amide hydrogens (Supplementary Figure S6B) indicated that the linker region is flexible, even in complex with binding partners.
It should be noted that TFIIE␣ differs from the other p62-binding partners in possessing some structure in the binding site, although the N-terminal intrinsically disordered acidic region is mainly responsible for binding to the PH-D (29); however, our MD simulation indicates that a stable complex forms between the TFIIE␣ acidic domain and the p62 PH-D in human TFIIH without any helix unfolding of the PH-D or BSD1 domain ( Supplementary Figure S4B). The present complex model of the TFIIE␣ acidic domain bound to p62 PH-D is likely to correspond to the initiation step in formation of the final PIC. As mentioned above, yeast Tfb1 PH-D is followed by a long (∼50-residue) linker, with the result that the PH-D sits at a different location in PIC as compared with the position of the human p62 PH-D. Even if the p62 PH-D could be located at the same position as the Tfb1 PH-D, it seems unlikely that TFIIE␣ could bind to the p62 PH-D due to steric hindrance from the structured part of the acidic domain; both the structure in the binding site of TFIIE␣ and the length of the linker of p62 are likely to determine their positions in PIC.
These results verify the importance of an intrinsically disordered region in recognition of the PH-D of p62 in TFIIH. Namely, the PH-D in TFIIH is sufficiently free to allow access to intrinsically disordered partners, but not completely free to allow access by a structured protein because of the short interdomain linker of p62, which restricts the range of motion of the PH-D in TFIIH as mentioned above. Although the dynamic behavior of the PH-D is essential for TFIIH recruitment, a completely free state would allow approach by not only true partners but also many other irrelevant proteins, resulting in inefficient recruitment. It is likely that this proper spatial restriction serves to enhance target selectivity.
It has recently been reported that monoubiquitination of Lys414 of UVSSA facilitates transfer of TFIIH to RNA polymerase II stalled at DNA lesions during transcriptioncoupled nucleotide excision repair (82). Our MD simulation shows that Lys414, located at the C-terminus of the p62-binding site, remains accessible even though UVSSA is bound to the PH-D of p62 in TFIIH.
The refined complex models also suggest some electrostatic interactions between acidic residues around the Cterminus of the binding sites of partners and lysine residues in the interdomain linker of p62: namely, Asp436 of TFIIE␣ and Lys104 of p62; Glu435 of TFIIE␣ and Lys106 of p62 (Supplementary Figure S8A); Glu62 of p53 and Lys106 of p62 (Supplementary Figure S8B); Glu138 of XPC and Lys104 of p62 (Supplementary Figure S8D); and Glu413 of UVSSA and Lys104 of p62 (Supplementary Figure S8E). However, there were no electrostatic interactions in the DP1-TFIIH complex, even though DP1 has four consecutive acidic residues at the C-terminus ( Supplementary Figure S8C).
Thus, we performed NMR chemical shift perturbation experiments for 15 N-labeled p62 1-158 with the five binding partners and identified significant chemical shift changes of these lysine residues in the interdomain linker of p62 when adding all binding partners except for DP1 (Supplementary Figure S9). The NMR titration experiments were well correlated with the complex models observed in MD simulations. In contrast, Arg105 of p62 displayed a large chemical shift change, but the MD models indicated that this basic residue does not participate in direct binding to the partners except for TFIIE␣, where Glu435 or Asp436 might interact with Arg105 of p62; the large chemical shift change of Arg105 in other cases is likely to be caused by interactions of the neighboring lysine residues, Lys104 and/or Lys106, with binding partners (Supplementary Figure S9).

CONCLUSION
Among the 10 subunits of TFIIH, p62 is relatively unique. In addition to its involvement in both formation of the core architecture and recruitment of TFIIH, the cryo-EM structure revealed that an unstructured and extended segment of p62 (residues 266-287) is inserted into the DNA-binding cavity of XPD, suggesting that p62 may function in regulating the helicase activity of XPD (42,43,78). Again, the dynamical behavior of p62 is crucial for TFIIH function. In the present study, we have revealed the dynamic structural model of human apo TFIIH core by using NMR and MD simulations based the cryo-EM structure, and provided the first structural view of the dynamic behavior of the PH-D of p62 in human apo TFIIH, revealing a different range of motion from the PH-D in the previously modeled human holo TFIIH (78) and yeast holo TFIIH (42) structures. We have also presented models of TFIIH complexed with five p62-binding partners refined by MD simulations. These will be useful for discussing interactions between the fulllength proteins. For example, a previous study presented a model of a cryo-EM structure of the full-length p53-RNA polymerase II complex docked in the PIC; however, the Nterminal transactivation domain of p53 was invisible owing to its structural flexibility (83). We believe that the complementary approach (NMR, MD and cryo-EM) demonstrated here will help us to better understand the dynamics of such key proteins in large, multi-subunit protein complexes.

DATA AVAILABILITY
Coordinates for the TFIIH p62  have been deposited in the RCSB Protein Data Bank (PDB) under ID code 7BUL, and chemical shifts, 15 N relaxation parameters (T 1 = 1/R 1 , T 2 = 1/R 2 and heteronuclear NOE) and hydrogen exchange rates (k ex ) have been deposited in the BioMagRes-Bank (BMRB) under accession number 36340.