Mutations in DNA polymerase β (pol β) have been associated with approximately 30% of human tumors. The E295K mutation of pol β has been linked to gastric carcinoma via interference with base excision repair. To interpret the different behavior of E295K compared to wild-type pol β in atomic and energetic detail, we resolve a binary crystal complex of E295K at 2.5 Å and apply transition path sampling (TPS) to delineate the closing pathway of the E295K pol β mutant. Conformational changes are important components in the enzymatic pathway that lead to and ready the enzyme for the chemical reaction. Our analyses show that the closing pathway of E295K mutant differs from the wild-type pol β in terms of the individual transition states along the pathway, associated energies, and the active site conformation in the final closed form of the mutant. In particular, the closed state of E295K has a more distorted active site than the active site in the wild-type pol β. In addition, the total energy barrier in the conformational closing pathway is 65 ± 11 kJ/mol, much higher than that estimated for both correct (e.g., G:C) and incorrect (e.g., G:A) wild-type pol β systems (42 ± 8 kJ/mol and 45 ± 7 kJ/mol, respectively). In particular, the rotation of Arg258 is the rate-limiting step in the conformational pathway of E295K due to unfavorable electrostatic and steric interactions. The distorted active site in the closed relative to open state and the high energy barrier in the conformational pathway may explain in part why the E295K mutant is observed to be inactive. Interestingly, however, following the closing of the thumb but prior to the rotation of Arg258, the E295K mutant complex has a similar energy level compared to the wild-type pol β. This suggests that the E295K mutant may associate with DNA with similar affinity, but it may be hampered in continuing the process of chemistry. Supporting experimental data come from the observation that the catalytic activity of wild-type pol β is hampered when E295K is present: this may arise from the competition between E295K and wild-type enzyme for the DNA. These combined results suggest that the low insertion efficiency of E295K mutant compared to wild-type pol β may be related to a closed-form distorted by unfavorable electrostatic and steric interactions between Arg258 and other key residues. The active site is thus less competent for proceeding to the chemical reaction, which may also involve a higher reaction barrier than the wild-type or may not be possible in this mutant. Our analysis also suggests further experiments for other mutants to test the above hypothesis and dissect the roles of steric and electrostatic factors on enzyme behavior.

To investigate whether an open-to-closed transition before the chemical step and induced-fit mechanism exist in DNA polymerase μ (pol μ), we analyze a series of molecular-dynamics simulations with and without the incoming nucleotide in various forms, including mutant systems, based on pol μ's crystal ternary structure. Our simulations capture no significant large-scale motion in either the DNA or the protein domains of pol μ. However, subtle residue motions can be distinguished, specifically of His(329) and Asp(330) to assemble in pol μ's active site, and of Gln(440) and Glu(443) to help accommodate the incoming nucleotide. Mutant simulations capture a DNA frameshift pairing and indicate the importance of Arg(444) and Arg(447) in stacking with the DNA template, and of Arg(448) and Gln(440) in helping to stabilize the position of both the DNA template and the incoming nucleotide. Although limited sampling in the molecular-dynamics simulations cannot be ruled out, our studies suggest an absence of a large-scale motion in pol μ. Together with the known crystallization difficulties of capturing the open form of pol μ, our studies also raise the possibility that a well-defined open form may not exist. Moreover, we suggest that residues Arg(448) and Gln(440) may be crucial for preventing insertion frameshift errors in pol μ.

To interpret recent structures of the R283K mutant of human DNA repair enzyme DNA polymerase β (pol β) differing in the number of Mg(2+) ions, we apply transition path sampling (TPS) to assess the effect of differing ion placement on the transition from the open one-metal to the closed two-metal state. We find that the closing pathway depends on the initial ion position, both in terms of the individual transition states and associated energies. The energy barrier of the conformational pathway varies from 25 to 58 kJ/mol, compared to the conformational energy barrier of 42 kJ/mol for the wild-type pol β reported previously. Moreover, we find a preferred ion route located in the center of the enzyme, parallel to the DNA. Within this route, the conformational pathway is similar to that of the overall open to closed transition of pol β, but outside it, especially when the ion starts near active site residues Arg258 and Asp190, the conformational pathway diverges significantly. Our findings should apply generally to pol β, since R283K is relatively far from the active site; further experimental and computational work are required to confirm this. Our studies also underscore the common feature that less active mutants have less stable closed states than their open states, in marked contrast to the wild-type enzyme, where the closed state is significantly more stable than the open form.