Solution study of the Escherichia coli DNA polymerase III clamp loader reveals the location of the dynamic ψχ heterodimer

Several X-ray crystal structures of the E. coli core clamp loader containing the five core (δ′, δ, and three truncated γ) subunits have been determined, but they lack the ψ and χ subunits. We report the first solution structure of the complete seven-subunit clamp loader complex using small angle X-ray scattering. This structure not only provides information about the location of the χ and ψ subunits but also provides a model of the dynamic nature of the clamp loader complex.


I. INTRODUCTION
The rate of DNA synthesis in E. coli ($1000 bp/s) is dependent on the function of the b sliding clamp and its associated clamp loader assembly. [1][2][3] The b sliding clamp is a homodimer ring-shaped protein 4 which confers processivity upon the replicative polymerase III and provides a mobile platform for interaction with various proteins during DNA damage repair, in addition to the polymerase. [5][6][7][8] The clamp loader complex is a member of the AAAþ family of ATPases (adenosine 5 0 -triphosphatases). The clamp loader utilizes ATP binding and hydrolysis to load the sliding clamp onto the DNA at the primer template junction [9][10][11][12] and is comprised of seven subunits (d, 3s or c, d 0 , v and w). [13][14][15][16] The d subunit is responsible for clamp binding and opening. The d 0 subunit acts as a stator and stabilizes the interaction of d subunit with the sliding clamp. [17][18][19][20] The s and c subunits are the active ATPases, and both c and s are products of the dnaX gene. [21][22][23] The s and c clamp loading function is interchangeable and the major difference is in the length of the two proteins; c is a shorter version of s subunits that is created by a translational frame shift. 24,25 Each DnaX subunit (either s or c) binds one molecule of ATP, and the clamp loader binds and hydrolyses three ATP molecules for each loading cycle. 26 The E. coli clamp loader complex contains two other subunits w and v. These two subunits form a tight 1:1 elongated heterodimeric complex. 15,27 The w subunit interacts with the C-terminal region of c subunit. 15,28 These two subunits are essential for bridging the interaction between the clamp loader and single strand DNA binding protein (SSB) in E. coli. The w subunit plays a role in stabilizing the conformational changes induced by ATP binding, the v subunit directly interacts with the C terminus of SSB (8 amino acid residues), and the interaction of SSB with v subunit is increased thousand of folds when SSB is bound to DNA. [29][30][31] The vw complex also plays a role in increasing the affinity of s and c for d.d 0 to a physiologically relevant range. 32 Several X-ray crystal structures of the E. coli clamp loader complex have been determined but all of these structures lack the w and v subunits. 13,14,33 The X-ray crystal structure of the clamp loader complex bound to N-terminal 28 amino acids of w subunit 33 showed that this interaction is restricted to the collar region of the three DnaX (c) subunits, and previous crosslinking studies have shown that the w subunit interacts with the c subunit closest to the d 0 . 29,34 Mass spectroscopy data have revealed a geometrical arrangement of the heterodimer where wv extends over the two c subunits and w adjacent to d. 35 In these studies, we have obtained for the first time using small angle X-ray scattering a solution structure of the complete seven subunit E. coli clamp loader (d, 3c, d 0 , w and v) complex in the presence of ATP. The Small Angle X-ray Scattering (SAXS) model reveals the location of the interaction of the v and w subunit with the C-terminal region of the c subunit away from the single stranded DNA exit site on d subunit and closer to the d 0 subunit in agreement with the crosslinking data. This structure also provides insight about the location of interaction of the vw heterodimer with the clamp loader complex in the absence of DNA and provides possible mechanism for the interaction of the clamp loader complex with SSB.

A. Protein expression and purification
The E. coli clamp loader complex subunits d, c, d 0 , v, and w were expressed, purified, and reconstituted as described previously, and the clamp loader was reconstituted from purified subunits. 32,[36][37][38][39] The storage buffer for c complex contained: 20 mM Tris-HCl, 0.5 mM EDTA, 2 mM DTT, 10% glycerol, and 50 mM NaCl. The final ATP and MgCl 2 concentrations for clamp loader complex samples with ATP were 500 lM and 8 mM, respectively.
B. Small angle x-ray scattering SAXS data (X-ray 9.881 keV, k ¼ 1.256 Å , (tantalum edge) at 1 Â 10 10 photons/s) were collected on the F2 station at the Cornell High Energy Synchrotron Source (CHESS) in Ithaca, New York. Prior to data collection, all samples were centrifuged at 15 000 Â g for 10 min to remove any aggregates. For data collection, 30 ll of sample, which was in the storage buffer, was loaded manually and robotically oscillated in the beam. The storage buffer was used as buffer background (30 ll). The X-ray beam was collimated to 250 Â 250 lm diameter and centered on a 2 mm diameter horizontal quartz capillary tube with 10 lm thick walls (Hampton Research, Aliso Viejo, CA). The capillary tube and full X-ray flight path, including beamstop, were kept in vacuo to eliminate air scatter. Data were captured on a Pilatus 100 K-S detector (Dectris, Baden, Switzerland) at 1500 mm distance. Scattering profiles were measured over the range of momentum transfer 0.01 < q ¼ 4p sin(h)/k < 0.27 Å À1 . The exposure time was 30 s, and 10 data sets were collected for possible radiation damage detection. Undamaged exposures were averaged and buffer subtracted using the RAW software. 40 Scattering data were also measured at different protein concentrations to check for presence of aggregation as a function of increasing concentration. For the clamp loader complex, the R g values at varying concentrations (2.0, 1.0, 0.5, and 0.25 mg/ml) were 54.6 Å 6 1.4, 53.8 Å 6 1.0, 51.0 Å 6 1.3, and 54.3 Å 6 1.4. Data collection and scattering derived parameters are summarized in Table I.

C. Software and modelling of SAXS data
The ATSAS suite of programs (Hamburg, EMBL (European Molecular Biology Labs)) was used for data processing. The program PRIMUS was used to obtain the Guinier plot and to calculate the radius of gyration (R g ). 41 Linear fitting was performed on the data having R g .q < 1.3 for Guinier analysis. The program GNOM 42 was used to calculate the particle distance distribution function P(r). The program DAMMIF (which is a rapid ab initio shape determination program using a single phase dummy sphere model) 43 was used to calculate the SAXS envelope by simulated annealing, and DAMMAVER was used to average the models. 44 Programs COOT 45 (crystallographic object oriented toolkit) and Chimera 46 were used to manually model the previously determined X-ray crystal structures of the E. coli d, 3c, d 0 (PDB: 3GLF) 33 complex (with the additionally modelled DnaX regions (369-431)) and heterodimeric v and w complex (PDB:1EM8) 16 into the clamp loader SAX envelope. Chimera was then used to rigid body fit the individual protein subunits into the SAX envelope. The program CRYSOL 47 was used to calculate theoretical scattering curves for comparison of the clamp loader complex crystal models and the experimental 1D scattering curve.

III. RESULTS
A. 1D scattering curve, shape, and quality characterization We obtained SAXS profiles for the complete clamp loader complex. The 1D scattering curve for the clamp loader complex in the presence of 500 lM ATP is shown ( Fig. 1(a)). The Guinier analysis of the scattering curve (which is analysis of the scattering curve at very low angles, or small q values) is shown ( Fig. 1(b)). The Guinier plot shows no deviation from linearity, and the linearity is consistent with the sample being monodisperse. The R g .q < 1.3 requirement for the globular proteins was satisfied. The radius of gyration R g , which provides information about the rotational volume of the macromolecule, was calculated to be 55.4 6 0.3 Å (sample at 3.7 mg/ml). The consistency of R g values at different concentrations, list in material and methods section, indicates that there were no concentration effects.
The pair distance distribution function, denoted as P(r), was calculated using the inverse Fourier method as implemented in the GNOM 42 program. In SAXS, the P(r) function is used to describe the paired-set of distances between all of the points within the macromolecular structure and is a useful tool for predicting the overall shape of the molecule as well as visibly detecting conformational changes within a macromolecule. The P(r) function curve for the clamp loader complex in Fig. 1(c) shows that the clamp loader complex is generally spherical in shape (the bell shape curve is centered around 65 Å with additional extended domains as depicted by the presence of the tail at larger distances). The D max value was determined to be 200 Å . The D max calculated for the crystal structure of seven subunit clamp loader core complex later shown in Fig. 2(d) was roughly 190 Å . The slightly larger D max is expected due to the hydration shell that is absent from the crystal structures of the clamp loader complex and flexibility of the clamp loader in solution. The fit of the P(r) curve to the 1D scattering curve is also shown. 50 To ensure that the protein was folded properly, we generated the Kratky plot using the SASPLOT software (feature of PRIMUS). The Krakty plot for the clamp loader complex in the presence of ATP is shown (Fig. 1(d)). The plot shows the presence of a bell shape peak and convergence to x-axis at higher q values expected from a folded globular protein.

B. SAXS ab initio and rigid body modelling
To obtain three-dimensional molecular envelopes of the clamp loader complex, the program DAMMIF 43 was used. P1 symmetry was imposed for model prediction, ten models were generated for the clamp loader (at 3.7 mg/ml), 50 and the program DAMAVER 44 was used to find the reference model, align all the models with the reference model, and average the dummy sphere models. The mean normalized special discrepancy (NSD) value, which provides agreement between the SAXS envelopes, for individual clamp loader complex bound to ATP was 0.49 6 0.016. The NSD values below 1 are acceptable, values between 0.5 and 0.7 reflect a good model, and a value below 0.5 is an excellent model. The predicted solution structure (the average of these ten models) for c complex clamp loader bound in presence of ATP is shown ( Fig. 2(a)). The central primer template-binding cavity as well as the region in the d subunit where the single strand DNA branches away is highlighted in Fig. 2(a). When compared to the previously determined X-ray crystal structure of the clamp loader complex comprised of the five core subunits (PDB: 3GLF) 33 (Fig. 2(b)), there is extra density forming an "arm" on the surface of the five-subunit core that is not present in the surface representation of the crystal structure of d, 3c, and d 0 complex bound to DNA. This density most likely corresponds to the v and w heterodimeric complex that is absent in all the crystal structures of the E. coli clamp loader complex previously determined.
The C-terminal regions of the c subunits (residues 374-431) are truncated in the crystal structures because of the flexibility of this region. In the X-ray crystal structures, residues 374-431 of the 3c subunits are cleaved by protease treatment or use of truncated protein (residues 1-373). 14,33,48 The extra density is large enough in volume to accommodate the missing flexible regions of the three c subunits (marked 1 in Fig. 1(c)) but too small for accommodating the w and v subunits, and therefore, we modelled as the truncated C-terminal regions of the three c subunits and the extra density opposite the d subunit to v and w. Once the sevensubunit clamp loader model was docked into the SAXS envelope, it was discovered that there is a second region of extra density in the envelope as shown (marked 2 in Fig. 1(c)) which we believe indicates movement of the d subunit possibly in readiness for binding the b sliding clamp and is not captured in crystal structures possibly because of crystal packing effects. A crystal-packing diagram for the c-clamp loader is shown. 50 To check the validity of this SAXS envelope, the seven-subunit clamp loader complex (as shown in Fig. 2(c)) was docked into the SAXS 3D envelope and energy minimized using CHIMERA. 46 In addition, we generated a 3D model for the region of the c subunit that was truncated for crystallography (374-431) using the Robetta server from Baker's lab, 49 and this structure was merged onto the three c subunits (Fig. 2(d)). The correlation coefficient value for this docking was 0.9, which indicates a good fit of the crystal structure into the SAXS envelope. The w subunit was modelled to interact with the C-terminal collar region of the c subunit next to d 0 , based on the previous knowledge of w interacting with the clamp loader and v interacting with SSB.
For comparing the SAXS experimental scattering data to the atomic crystal structure models of the clamp loader complex, CRYSOL 47 was used. Three theoretical scattering curves were generated as shown (Fig. 3). The blue curve is the predicted curve for the five-subunit clamp loader complex using the PDB: 3GLF. 33 The v value was 5.6, and theoretical R g value was 45.5 Å . The green curve is the predicted scattering curve for the seven subunit clamp loader complex as shown in Fig. 1(d) with wv modeled close to the d 0 subunit based on ab inito SAXS structure with v value of 2.2. The red curve is the predicted scattering curve for the seven-subunit clamp loader complex with wv subunits positioned close to the ssDNA exit channel on the d subunit with v value of 4.2 and theoretical R g value of 46.8 Å . It is evident from CRYSOL analysis that the predicted curve of the five subunit clamp loader complex in the absence of wv heterodimer does not fit scattering curve at low and high q values and having the v value of 5.6 and R g value of 45.5 Å . The CYROSL predicted curve for the seven-subunit clamp loader complex with wv positioned close to d 0 subunit has the best fit to experimental data as evident by the lowest v value and better fit to the experimental data as well as better agreement of theoretical R g value of 51.1 Å with the experimental R g value of 53.4 Å . The discrepancy between the theoretical scattering curves for the seven subunit clamp loaders (red and green curve) and the experimental data at q ¼ 0.125 Å À1 could possibly be the result of the intramolecular movement or "breathing" of the clamp loader subunits in solution that are otherwise absent from static state X-ray structures. the interaction of the d subunit with the clamp, while the v and w subunits bridge the clamp loader complex with the SSB. The w N-terminal 28 amino acids are responsible for anchoring the wv to the clamp loader complex. 16,33 The crystal structure of the clamp loader complex with the five core subunits in the presence of the w peptide (residues 2-28) has been determined. 33 The X-ray crystal structure and biochemical studies of the clamp loader complex bound to the w peptide showed that this interaction is restricted to the collar region of the three DnaX subunits, 29 but, to this date, there has been no cryo-EM, X-ray, or SAXS structure of the E. coli seven subunit clamp loader complex. In these studies, we have obtained for the first time a low resolution SAXS model of the complete (all seven subunits) gamma clamp loader complex in the presence of ATP. The SAXS envelope suggests that the w and v subunits mainly interact in an asymmetrical fashion with the collar C-terminal region of one of the c subunits (the c subunit next to d 0 ) ( Fig. 2(a)), which is consistent with the previous crosslinking studies showing that the w subunit interacts with the c subunit closest to the d 0 . 29,34 This finding has implication that there is possible wrapping of ssDNA around the collar region for interaction with SSB at the replication fork (Fig. 4). The second possibility is that upon interaction of the clamp loader complex with the DNA, there is movement of wv possibly through the flexible N-terminal region of w, so that the SSB is located closer to the ssDNA exit channel on the d subunit and this second model is in agreement with ion mobility mass spectroscopy data. 35 It should be noted that the solution structure of the clamp loader complex discussed in this paper could possibly describe one of the several conformations of an ensemble of conformations expected of a highly dynamic macromolecule such as the clamp loader complex and that there might be a switching mechanism between the active and inactive state of the clamp loader upon binding of the DNA. The clamp loader complex SAXS model also reveals that the d subunit possibly undergoes some movements to be able to make contact with the b sliding clamp to open or to stabilize an open clamp and that the additional envelope volume highlighted in Fig. 2(c) is indicative of the location of the d subunit in readiness to bind the clamp. These dynamics could not have been captured in previously determined solid-state X-ray crystal structures.