DNA polymerase beta: structure-fidelity relationship from Pre-steady-state kinetic analyses of all possible correct and incorrect base pairs for wild type and R283A mutant

Correct Nucleotide Selection Is Confined at the Binding Site of Polymerase Enzymes

DNA polymerases (Pols) add incoming nucleotides (deoxyribonucleoside triphosphate (dNTPs)) to growing DNA strands, a crucial step for DNA synthesis. The insertion of correct (vs incorrect) nucleotides relates to Pols' fidelity, which defines Pols' ability to faithfully replicate DNA strands in a template-dependent manner. We and others have demonstrated that reactant alignment and correct base pairing at the Pols catalytic site are crucial structural features to fidelity. Here, we first used equilibrium molecular simulations to demonstrate that the local dynamics at the protein-DNA interface in the proximity of the catalytic site is different when correct vs incorrect dNTPs are bound to polymerase β (Pol β). Formation and dynamic stability of specific interatomic interactions around the incoming nucleotide influence the overall binding site architecture. This explains why certain Pols' mutants can affect the local catalytic environment and influence the selection of correct vs incorrect nucleotides. In particular, this is here demonstrated by analyzing the interaction network formed by the residue R283, whose mutant R283A has an experimentally measured lower capacity of differentiating correct (G:dCTP) vs incorrect (G:dATP) base pairing in Pol β. We also used alchemical free-energy calculations to quantify the G:dCTP →G:dATP transformation in Pol β wild-type and mutant R283A. These results correlate well with the experimental trend, thus corroborating our mechanistic insights. Sequence and structural comparisons with other Pols from the same family suggest that these findings may also be valid in similar enzymes.

Dynamic basis for dA•dGTP and dA•d8OGTP misincorporation via Hoogsteen base pairs

Replicative errors contribute to the genetic diversity needed for evolution but in high frequency can lead to genomic instability. Here, we show that DNA dynamics determine the frequency of misincorporating the A•G mismatch, and altered dynamics explain the high frequency of 8-oxoguanine (8OG) A•8OG misincorporation. NMR measurements revealed that Aanti•Ganti (population (pop.) of >91%) transiently forms sparsely populated and short-lived Aanti+•Gsyn (pop. of ~2% and kex = kforward + kreverse of ~137 s-1) and Asyn•Ganti (pop. of ~6% and kex of ~2,200 s-1) Hoogsteen conformations. 8OG redistributed the ensemble, rendering Aanti•8OGsyn the dominant state. A kinetic model in which Aanti+•Gsyn is misincorporated quantitatively predicted the dA•dGTP misincorporation kinetics by human polymerase β, the pH dependence of misincorporation and the impact of the 8OG lesion. Thus, 8OG increases replicative errors relative to G because oxidation of guanine redistributes the ensemble in favor of the mutagenic Aanti•8OGsyn Hoogsteen state, which exists transiently and in low abundance in the A•G mismatch.

Structural Insights into Phosphorylation-Mediated Polymerase Function Loss for DNA Polymerase β Bound to Gapped DNA

DNA polymerase β is a member of the X-family of DNA polymerases, playing a critical role in the base excision repair (BER) pathway in mammalian cells by implementing the nucleotide gap-filling step. In vitro phosphorylation of DNA polymerase β with PKC on S44 causes loss in the enzyme's DNA polymerase activity but not single-strand DNA binding. Although these studies have shown that single-stranded DNA binding is not affected by phosphorylation, the structural basis behind the mechanism underlying phosphorylation-induced activity loss remains poorly understood. Previous modeling studies suggested phosphorylation of S44 was sufficient to induce structural changes that impact the enzyme's polymerase function. However, the S44 phosphorylated-enzyme/DNA complex has not been modeled so far. To address this knowledge gap, we conducted atomistic molecular dynamics simulations of pol β complexed with gapped DNA. Our simulations, which used explicit solvent and lasted for microseconds, revealed that phosphorylation at the S44 site, in the presence of Mg ions, induced significant conformational changes in the enzyme. Specifically, these changes led to the transformation of the enzyme from a closed to an open structure. Additionally, our simulations identified phosphorylation-induced allosteric coupling between the inter-domain region, suggesting the existence of a putative allosteric site. Taken together, our results provide a mechanistic understanding of the conformational transition observed due to phosphorylation in DNA polymerase β interactions with gapped DNA. Our simulations shed light on the mechanisms of phosphorylation-induced activity loss in DNA polymerase β and reveal potential targets for the development of novel therapeutics aimed at mitigating the effects of this post-translational modification.

Phosphorylation Induced Conformational Transitions in DNA Polymerase β

DNA polymerase β (pol β) is a member of the X- family of DNA polymerases that catalyze the distributive addition of nucleoside triphosphates during base excision DNA repair. Previous studies showed that the enzyme was phosphorylated in vitro with PKC at two serines (44 and 55), causing loss of DNA polymerase activity but not DNA binding. In this work, we have investigated the phosphorylation-induced conformational changes in DNA polymerase β in the presence of Mg ions. We report a comprehensive atomic resolution study of wild type and phosphorylated DNA polymerase using molecular dynamics (MD) simulations. The results are examined via novel methods of internal dynamics and energetics analysis to reveal the underlying mechanism of conformational transitions observed in DNA pol β. The results show drastic conformational changes in the structure of DNA polymerase β due to S44 phosphorylation. Phosphorylation-induced conformational changes transform the enzyme from a closed to an open structure. The dynamic cross-correlation shows that phosphorylation enhances the correlated motions between the different domains. Centrality network analysis reveals that the S44 phosphorylation causes structural rearrangements and modulates the information pathway between the Lyase domain and base pair binding domain. Further analysis of our simulations reveals that a critical hydrogen bond (between S44 and E335) disruption and the formation of three additional salt bridges are potential drivers of these conformational changes. In addition, we found that two of these additional salt bridges form in the presence of Mg ions on the active sites of the enzyme. These results agree with our previous study of DNA pol β S44 phosphorylation without Mg ions which predicted the deactivation of DNA pol β. However, the phase space of structural transitions induced by S44 phosphorylation is much richer in the presence of Mg ions.

Two-Metal-Ion Catalysis: Inhibition of DNA Polymerase Activity by a Third Divalent Metal Ion

Almost all DNA polymerases (pols) exhibit bell-shaped activity curves as a function of both pH and Mg²⁺ concentration. The pol activity is reduced when the pH deviates from the optimal value. When the pH is too low the concentration of a deprotonated general base (namely, the attacking 3′-hydroxyl of the 3′ terminal residue of the primer strand) is reduced exponentially. When the pH is too high the concentration of a protonated general acid (i.e., the leaving pyrophosphate group) is reduced. Similarly, the pol activity also decreases when the concentration of the divalent metal ions deviates from its optimal value: when it is too low, the binding of the two catalytic divalent metal ions required for the full activity is incomplete, and when it is too high a third divalent metal ion binds to pyrophosphate, keeping it in the replication complex longer and serving as a substrate for pyrophosphorylysis within the complex. Currently, there is a controversy about the role of the third metal ion which we will address in this review.

The ligation of pol β mismatch insertion products governs the formation of promutagenic base excision DNA repair intermediates

DNA ligase I and DNA ligase III/XRCC1 complex catalyze the ultimate ligation step following DNA polymerase (pol) β nucleotide insertion during base excision repair (BER). Pol β Asn279 and Arg283 are the critical active site residues for the differentiation of an incoming nucleotide and a template base and the N-terminal domain of DNA ligase I mediates its interaction with pol β. Here, we show inefficient ligation of pol β insertion products with mismatched or damaged nucleotides, with the exception of a Watson–Crick-like dGTP insertion opposite T, using BER DNA ligases in vitro. Moreover, pol β N279A and R283A mutants deter the ligation of the promutagenic repair intermediates and the presence of N-terminal domain of DNA ligase I in a coupled reaction governs the channeling of the pol β insertion products. Our results demonstrate that the BER DNA ligases are compromised by subtle changes in all 12 possible noncanonical base pairs at the 3′-end of the nicked repair intermediate. These findings contribute to understanding of how the identity of the mismatch affects the substrate channeling of the repair pathway and the mechanism underlying the coordination between pol β and DNA ligase at the final ligation step to maintain the BER efficiency.

Environmental Effects on Guanine-Thymine Mispair Tautomerization Explored with Quantum Mechanical/Molecular Mechanical Free Energy Simulations

DNA bases can adopt energetically unfavorable tautomeric forms that enable the formation of Watson-Crick-like (WC-like) mispairs, which have been proposed to give rise to spontaneous mutations in DNA and misincorporation errors in DNA replication and translation. Previous NMR and computational studies have indicated that the population of WC-like guanine-thymine (G-T) mispairs depends on the environment, such as the local nucleic acid sequence and solvation. To investigate these environmental effects, herein G-T mispair tautomerization processes are studied computationally in aqueous solution, in A-form and B-form DNA duplexes, and within the active site of a DNA polymerase λ variant. The wobble G-T (wG-T), WC-like G-T*, and WC-like G*-T forms are considered, where * indicates the enol tautomer of the base. The minimum free energy paths for the tautomerization from the wG-T to the WC-like G-T* and from the WC-like G-T* to the WC-like G*-T are computed with mixed quantum mechanical/molecular mechanical (QM/MM) free energy simulations. The reaction free energies and free energy barriers are found to be significantly influenced by the environment. The wG-T→G-T* tautomerization is predicted to be endoergic in aqueous solution and the DNA duplexes but slightly exoergic in the polymerase, with Arg517 and Asn513 providing electrostatic stabilization of G-T*. The G-T*→G*-T tautomerization is also predicted to be slightly more thermodynamically favorable in the polymerase relative to these DNA duplexes. These simulations are consistent with an experimentally driven kinetic misincorporation model suggesting that G-T mispair tautomerization occurs in the ajar polymerase conformation or concertedly with the transition from the ajar to the closed polymerase conformation. Furthermore, the order of the associated two proton transfer reactions is predicted to be different in the polymerase than in aqueous solution and the DNA duplexes. These studies highlight the impact of the environment on the thermodynamics, kinetics, and fundamental mechanisms of G-T mispair tautomerization, which plays a role in a wide range of biochemically important processes.

Transitions in DNA polymerase β μs-ms dynamics related to substrate binding and catalysis

DNA polymerase β (pol β) plays a central role in the DNA base excision repair pathway and also serves as an important model polymerase. Dynamic characterization of pol β from methyl-TROSY 13C-1H multiple quantum CPMG relaxation dispersion experiments of Ile and Met sidechains and previous backbone relaxation dispersion measurements, reveals transitions in μs-ms dynamics in response to highly variable substrates. Recognition of a 1-nt-gapped DNA substrate is accompanied by significant backbone and sidechain motion in the lyase domain and the DNA binding subdomain of the polymerase domain, that may help to facilitate binding of the apoenzyme to the segments of the DNA upstream and downstream from the gap. Backbone μs-ms motion largely disappears after formation of the pol β-DNA complex, giving rise to an increase in uncoupled μs-ms sidechain motion throughout the enzyme. Formation of an abortive ternary complex using a non-hydrolyzable dNTP results in sidechain motions that fit to a single exchange process localized to the catalytic subdomain, suggesting that this motion may play a role in catalysis.

Human DNA Polymerase μ Can Use a Noncanonical Mechanism for Multiple Mn2+-Mediated Functions

Recent research on the structure and mechanism of DNA polymerases has continued to generate fundamentally important features, including a noncanonical pathway involving "prebinding" of metal-bound dNTP (MdNTP) in the absence of DNA. While this noncanonical mechanism was shown to be a possible subset for African swine fever DNA polymerase X (Pol X) and human Pol λ, it remains unknown whether it could be the primary pathway for a DNA polymerase. Pol μ is a unique member of the X-family with multiple functions and with unusual Mn²⁺ preference. Here we report that Pol μ not only prebinds MdNTP in a catalytically active conformation but also exerts a Mn²⁺ over Mg²⁺ preference at this early stage of catalysis, for various functions: incorporation of dNTP into a single nucleotide gapped DNA, incorporation of rNTP in the nonhomologous end joining (NHEJ) repair, incorporation of dNTP to an ssDNA, and incorporation of an 8-oxo-dGTP opposite template dA (mismatched) or dC (matched). The structural basis of this noncanonical mechanism and Mn²⁺ over Mg²⁺ preference in these functions was analyzed by solving 19 structures of prebinding binary complexes, precatalytic ternary complexes, and product complexes. The results suggest that the noncanonical pathway is functionally relevant for the multiple functions of Pol μ. Overall, this work provides the structural and mechanistic basis for the long-standing puzzle in the Mn²⁺ preference of Pol μ and expands the landscape of the possible mechanisms of DNA polymerases to include both mechanistic pathways.

Uniform Free-Energy Profiles of the P-O Bond Formation and Cleavage Reactions Catalyzed by DNA Polymerases β and λ

Human X-family DNA polymerases β (Polβ) and λ (Polλ) catalyze the nucleotidyl-transfer reaction in the base excision repair pathway of the cellular DNA damage response. Using empirical valence bond and free-energy perturbation simulations, we explore the feasibility of various mechanisms for the deprotonation of the 3′-OH group of the primer DNA strand, and the subsequent formation and cleavage of P–O bonds in four Polβ, two truncated Polλ (tPolλ), and two tPolλ Loop1 mutant (tPolλΔL1) systems differing in the initial X-ray crystal structure and nascent base pair. The average calculated activation free energies of 14, 18, and 22 kcal mol^–1 for Polβ, tPolλ, and tPolλΔL1, respectively, reproduce the trend in the observed catalytic rate constants. The most feasible reaction pathway consists of two successive steps: specific base (SB) proton transfer followed by rate-limiting concerted formation and cleavage of the P–O bonds. We identify linear free-energy relationships (LFERs) which show that the differences in the overall activation and reaction free energies among the eight studied systems are determined by the reaction free energy of the SB proton transfer. We discuss the implications of the LFERs and suggest pK_a of the 3′-OH group as a predictor of the catalytic rate of X-family DNA polymerases.

Induced Fit in the Selection of Correct versus Incorrect Nucleotides by DNA Polymerase β

DNA polymerase β (Pol β) repairs single-nucleotide gapped DNA (sngDNA) by enzymatic incorporation of the Watson-Crick partner nucleotide at the gapped position opposite the templating nucleotide. The process by which the matching nucleotide is incorporated into a sngDNA sequence has been relatively well-characterized, but the process of discrimination from nucleotide misincorporation remains unclear. We report here NMR spectroscopic characterization of full-length, uniformly labeled Pol β in apo, sngDNA-bound binary, and ternary complexes containing matching and mismatching nucleotide. Our data indicate that, while binding of the correct nucleotide to the binary complex induces chemical shift changes consistent with the process of enzyme closure, the ternary Pol β complex containing a mismatching nucleotide exhibits no such changes and appears to remain in an open, unstable, binary-like conformation. Our findings support an induced-fit mechanism for polymerases in which a closed ternary complex can only be achieved in the presence of matching nucleotide.

Substrate-induced DNA polymerase β activation

DNA polymerases and substrates undergo conformational changes upon forming protein-ligand complexes. These conformational adjustments can hasten or deter DNA synthesis and influence substrate discrimination. From structural comparison of binary DNA and ternary DNA-dNTP complexes of DNA polymerase β, several side chains have been implicated in facilitating formation of an active ternary complex poised for chemistry. Site-directed mutagenesis of these highly conserved residues (Asp-192, Arg-258, Phe-272, Glu-295, and Tyr-296) and kinetic characterization provides insight into the role these residues play during correct and incorrect insertion as well as their role in conformational activation. The catalytic efficiencies for correct nucleotide insertion for alanine mutants were wild type ∼ R258A > F272A ∼ Y296A > E295A > D192A. Because the efficiencies for incorrect insertion were affected to about the same extent for each mutant, the effects on fidelity were modest (<5-fold). The R258A mutant exhibited an increase in the single-turnover rate of correct nucleotide insertion. This suggests that the wild-type Arg-258 side chain generates a population of non-productive ternary complexes. Structures of binary and ternary substrate complexes of the R258A mutant and a mutant associated with gastric carcinomas, E295K, provide molecular insight into intermediate structural conformations not appreciated previously. Although the R258A mutant crystal structures were similar to wild-type enzyme, the open ternary complex structure of E295K indicates that Arg-258 stabilizes a non-productive conformation of the primer terminus that would decrease catalysis. Significantly, the open E295K ternary complex binds two metal ions indicating that metal binding cannot overcome the modified interactions that have interrupted the closure of the N-subdomain.

Phylogenetic analysis and evolutionary origins of DNA polymerase X-family members

Mammalian DNA polymerase (pol) β is the founding member of a large group of DNA polymerases now termed the X-family. DNA polymerase β has been kinetically, structurally, and biologically well characterized and can serve as a phylogenetic reference. Accordingly, we have performed a phylogenetic analysis to understand the relationship between pol β and other members of the X-family of DNA polymerases. The bacterial X-family DNA polymerases, Saccharomyces cerevisiae pol IV, and four mammalian X-family polymerases appear to be directly related. These enzymes originated from an ancient common ancestor characterized in two Bacillus species. Understanding distinct functions for each of the X-family polymerases, evolving from a common bacterial ancestor is of significant interest in light of the specialized roles of these enzymes in DNA metabolism.

Fluorescence resonance energy transfer studies of DNA polymerase β: the critical role of fingers domain movements and a novel non-covalent step during nucleotide selection

During DNA repair, DNA polymerase β (Pol β) is a highly dynamic enzyme that is able to select the correct nucleotide opposite a templating base from a pool of four different deoxynucleoside triphosphates (dNTPs). To gain insight into nucleotide selection, we use a fluorescence resonance energy transfer (FRET)-based system to monitor movement of the Pol β fingers domain during catalysis in the presence of either correct or incorrect dNTPs. By labeling the fingers domain with ((((2-iodoacetyl)amino)ethyl)amino)naphthalene-1-sulfonic acid (IAEDANS) and the DNA substrate with Dabcyl, we are able to observe rapid fingers closing in the presence of correct dNTPs as the IAEDANS comes into contact with a Dabcyl-labeled, one-base gapped DNA. Our findings show that not only do the fingers close after binding to the correct dNTP, but that there is a second conformational change associated with a non-covalent step not previously reported for Pol β. Further analyses suggest that this conformational change corresponds to the binding of the catalytic metal into the polymerase active site. FRET studies with incorrect dNTP result in no changes in fluorescence, indicating that the fingers do not close in the presence of incorrect dNTP. Together, our results show that nucleotide selection initially occurs in an open fingers conformation and that the catalytic pathways of correct and incorrect dNTPs differ from each other. Overall, this study provides new insight into the mechanism of substrate choice by a polymerase that plays a critical role in maintaining genome stability.

Recent insight into the kinetic mechanisms and conformational dynamics of Y-Family DNA polymerases

The kinetic mechanisms by which DNA polymerases catalyze DNA replication and repair have long been areas of active research. Recently discovered Y-family DNA polymerases catalyze the bypass of damaged DNA bases that would otherwise block replicative DNA polymerases and stall replication forks. Unlike DNA polymerases from the five other families, the Y-family DNA polymerases have flexible, solvent-accessible active sites that are able to tolerate various types of damaged template bases and allow for efficient lesion bypass. Their promiscuous active sites, however, also lead to fidelities that are much lower than those observed for other DNA polymerases and give rise to interesting mechanistic properties. Additionally, the Y-family DNA polymerases have several other unique structural features and undergo a set of conformational changes during substrate binding and catalysis different from those observed for replicative DNA polymerases. In recent years, pre-steady-state kinetic methods have been extensively employed to reveal a wealth of information about the catalytic properties of these fascinating noncanonical DNA polymerases. Here, we review many of the recent findings on the kinetic mechanisms of DNA polymerization with undamaged and damaged DNA substrates by the Y-family DNA polymerases, and the conformational dynamics employed by these error-prone enzymes during catalysis.

C(α) torsion angles as a flexible criterion to extract secrets from a molecular dynamics simulation

Given the increasing complexity of simulated molecular systems, and the fact that simulation times have now reached milliseconds to seconds, immense amounts of data (in the gigabyte to terabyte range) are produced in current molecular dynamics simulations. Manual analysis of these data is a very time-consuming task, and important events that lead from one intermediate structure to another can become occluded in the noise resulting from random thermal fluctuations. To overcome these problems and facilitate a semi-automated data analysis, we introduce in this work a measure based on C(α) torsion angles: torsion angles formed by four consecutive C(α) atoms. This measure describes changes in the backbones of large systems on a residual length scale (i.e., a small number of residues at a time). Cluster analysis of individual C(α) torsion angles and its fuzzification led to continuous time patches representing (meta)stable conformations and to the identification of events acting as transitions between these conformations. The importance of a change in torsion angle to structural integrity is assessed by comparing this change to the average fluctuations in the same torsion angle over the complete simulation. Using this novel measure in combination with other measures such as the root mean square deviation (RMSD) and time series of distance measures, we performed an in-depth analysis of a simulation of the open form of DNA polymerase I. The times at which major conformational changes occur and the most important parts of the molecule and their interrelations were pinpointed in this analysis. The simultaneous determination of the time points and localizations of major events is a significant advantage of the new bottom-up approach presented here, as compared to many other (top-down) approaches in which only the similarity of the complete structure is analyzed.

The E295K cancer variant of human polymerase β favors the mismatch conformational pathway during nucleotide selection

DNA polymerase β (pol β) is responsible for gap filling synthesis during repair of damaged DNA as part of the base excision repair pathway. Human pol β mutations were recently identified in a high percentage (∼30%) of tumors. Characterization of specific cancer variants is particularly useful to further the understanding of the general mechanism of pol β while providing context to disease contribution. We showed that expression of the carcinoma variant E295K induces cellular transformation. The poor polymerase activity exhibited by the variant was hypothesized to be caused by the destabilization of proper active site assembly by the glutamate to lysine mutation. Here, we show that this variant exhibits an unusual preference for binding dCTP opposite a templating adenine over the cognate dTTP. Biochemical studies indicate that the noncognate competes with the cognate nucleotide for binding to the polymerase active site with the noncognate incorporation a function of higher affinity and not increased activity. In the crystal structure of the variant bound to dA:dCTP, the fingers domain closes around the mismatched base pair. Nucleotide incorporation is hindered because key residues in the polymerase active site are not properly positioned for nucleotidyl transfer. In contrast to the noncognate dCTP, neither the cognate dTTP nor its nonhydrolyzable analog induced fingers closure, as isomorphous difference Fourier maps show that the cognate nucleotides are bound to the open state of the polymerase. Comparison with published structures provides insight into the structural rearrangements within pol β that occur during the process of nucleotide discrimination.

"Gate-keeper" residues and active-site rearrangements in DNA polymerase μ help discriminate non-cognate nucleotides

Incorporating the cognate instead of non-cognate substrates is crucial for DNA polymerase function. Here we analyze molecular dynamics simulations of DNA polymerase μ (pol μ) bound to different non-cognate incoming nucleotides including A:dCTP, A:dGTP, A(syn):dGTP, A:dATP, A(syn):dATP, T:dCTP, and T:dGTP to study the structure-function relationships involved with aberrant base pairs in the conformational pathway; while a pol μ complex with the A:dTTP base pair is available, no solved non-cognate structures are available. We observe distinct differences of the non-cognate systems compared to the cognate system. Specifically, the motions of active-site residue His329 and Asp330 distort the active site, and Trp436, Gln440, Glu443 and Arg444 tend to tighten the nucleotide-binding pocket when non-cognate nucleotides are bound; the latter effect may further lead to an altered electrostatic potential within the active site. That most of these “gate-keeper” residues are located farther apart from the upstream primer in pol μ, compared to other X family members, also suggests an interesting relation to pol μ's ability to incorporate nucleotides when the upstream primer is not paired. By examining the correlated motions within pol μ complexes, we also observe different patterns of correlations between non-cognate systems and the cognate system, especially decreased interactions between the incoming nucleotides and the nucleotide-binding pocket. Altered correlated motions in non-cognate systems agree with our recently proposed hybrid conformational selection/induced-fit models. Taken together, our studies propose the following order for difficulty of non-cognate system insertions by pol μ: T:dGTP<A(syn):dATP<T:dCTP<A:dGTP<A(syn):dGTP<A:dCTP<A:dATP. This sequence agrees with available kinetic data for non-cognate nucleotide insertions, with the exception of A:dGTP, which may be more sensitive to the template sequence. The structures and conformational aspects predicted here are experimentally testable.

DNA polymerase μ (pol μ) is an enzyme that participates in DNA repair and thus has a central role in maintaining the integrity of genetic information. To efficiently repair the DNA, discriminating the cognate instead of non-cognate nucleotides (“fidelity-checking”) is required. Here we analyze molecular dynamics simulations of pol μ bound to different non-cognate nucleotides to study the structure-function relationships involved in the fidelity-checking mechanism of pol μ on the atomic level. Our results suggest that His329, Asp330, Trp436, Gln440, Glu443, and Arg444 are of great importance for pol μ's fidelity-checking mechanism. We also observe altered patterns of correlated motions within pol μ complex when non-cognate instead of cognate nucleotides are bound, which agrees with our recently proposed hybrid conformational selection/induced-fit models. Taken together, our studies help interpret the available kinetic data of various non-cognate nucleotide insertions by pol μ. We also suggest experimentally testable predictions; for example, a point mutation like E443M may reduce the ability of pol μ to insert the cognate more than of non-cognate nucleotides. Our studies suggest an interesting relation to pol μ's unique ability to incorporate nucleotides when the upstream primer is not paired.

Structures of dNTP intermediate states during DNA polymerase active site assembly

DNA polymerase and substrate conformational changes are essential for high-fidelity DNA synthesis. Structures of DNA polymerase (pol) β in complex with DNA show the enzyme in an "open" conformation. Subsequent to binding the nucleotide, the polymerase "closes" around the nascent base pair with two metals positioned for chemistry. However, structures of substrate/active site intermediates prior to closure are lacking. By destabilizing the closed complex, we determined unique ternary complex structures of pol β with correct and incorrect incoming nucleotides bound to the open conformation. These structures reveal that Watson-Crick hydrogen bonding is assessed upon initial complex formation. Importantly, nucleotide-bound states representing intermediate metal coordination states occur with active site assembly. The correct, but not incorrect, nucleotide maintains Watson-Crick hydrogen bonds during interconversion of these states. These structures indicate that the triphosphate of the incoming nucleotide undergoes rearrangement prior to closure, providing an opportunity to deter misinsertion and increase fidelity.

Catalytic effects of mutations of distant protein residues in human DNA polymerase β: theory and experiment

We carried out free-energy calculations and transient kinetic experiments for the insertion of the right (dC) and wrong (dA) nucleotides by wild-type (WT) and six mutant variants of human DNA polymerase β (Pol β). Since the mutated residues in the point mutants, I174S, I260Q, M282L, H285D, E288K, and K289M, were not located in the Pol β catalytic site, we assumed that the WT and its point mutants share the same dianionic phosphorane transition-state structure of the triphosphate moiety of deoxyribonucleotide 5'-triphosphate (dNTP) substrate. On the basis of this assumption, we have formulated a thermodynamic cycle for calculating relative dNTP insertion efficiencies, Ω = (k(pol)/K(D))(mut)/(k(pol)/K(D))(WT) using free-energy perturbation (FEP) and linear interaction energy (LIE) methods. Kinetic studies on five of the mutants have been published previously using different experimental conditions, e.g., primer-template sequences. We have performed a presteady kinetic analysis for the six mutants for comparison with wild-type Pol β using the same conditions, including the same primer/template DNA sequence proximal to the dNTP insertion site used for X-ray crystallographic studies. This consistent set of kinetic and structural data allowed us to eliminate the DNA sequence from the list of factors that can adversely affect calculated Ω values. The calculations using the FEP free energies scaled by 0.5 yielded 0.9 and 1.1 standard deviations from the experimental log Ω values for the insertion of the right and wrong dNTP, respectively. We examined a hybrid FEP/LIE method in which the FEP van der Waals term for the interaction of the mutated amino acid residue with its surrounding environment was replaced by the corresponding van der Waals term calculated using the LIE method, resulting in improved 0.4 and 1.0 standard deviations from the experimental log Ω values. These scaled FEP and FEP/LIE methods were also used to predict log Ω for R283A and R283L Pol β mutants.

Residues Arg283, Arg285, and Ile287 in the nucleotide binding pocket of bovine viral diarrhea virus NS5B RNA polymerase affect catalysis and fidelity

Residues Arg283, Arg285, and Ile287 are highly conserved amino acids in bovine viral diarrhea virus RNA polymerase (BVDV RdRp) and RdRps from related positive-strand RNA viruses. This motif is an important part of the binding pocket for the nascent RNA base pair during initiation and elongation. We found that replacement of the arginines with alanines or more conserved lysines or replacement of isoleucine with alanine or valine alters the ability of the mutant RdRps to incorporate ribonucleotides efficiently. The reduced RdRp activity stems from both decreased ribonucleotide binding and decreased catalytic efficiency in both primer-dependent and de novo initiation, as shown by kinetic studies. In line with other studies on flaviviral RdRps, our data suggest that Arg283 and Ile287 may be implicated in ribonucleotide binding and positioning of the template base in the active site. Arg285 appears to be involved directly in the selection of cognate nucleotide. The findings for Arg285 and Ile287 mutants also agree with similar data from picornavirus RdRps.

Unfavorable electrostatic and steric interactions in DNA polymerase β E295K mutant interfere with the enzyme's pathway

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 as 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 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 as 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 as 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.

Thermodynamic basis for the emergence of genomes during prebiotic evolution

The RNA world hypothesis views modern organisms as descendants of RNA molecules. The earliest RNA molecules must have been random sequences, from which the first genomes that coded for polymerase ribozymes emerged. The quasispecies theory by Eigen predicts the existence of an error threshold limiting genomic stability during such transitions, but does not address the spontaneity of changes. Following a recent theoretical approach, we applied the quasispecies theory combined with kinetic/thermodynamic descriptions of RNA replication to analyze the collective behavior of RNA replicators based on known experimental kinetics data. We find that, with increasing fidelity (relative rate of base-extension for Watson-Crick versus mismatched base pairs), replications without enzymes, with ribozymes, and with protein-based polymerases are above, near, and below a critical point, respectively. The prebiotic evolution therefore must have crossed this critical region. Over large regions of the phase diagram, fitness increases with increasing fidelity, biasing random drifts in sequence space toward ‘crystallization.’ This region encloses the experimental nonenzymatic fidelity value, favoring evolutions toward polymerase sequences with ever higher fidelity, despite error rates above the error catastrophe threshold. Our work shows that experimentally characterized kinetics and thermodynamics of RNA replication allow us to determine the physicochemical conditions required for the spontaneous crystallization of biological information. Our findings also suggest that among many potential oligomers capable of templated replication, RNAs may have evolved to form prebiotic genomes due to the value of their nonenzymatic fidelity.

A leading hypothesis for the origin of life describes a prebiotic world where RNA molecules started carrying genetic information for catalyzing their own replication. This origin of biological information is akin to the crystallization of ice from water, where ‘order’ emerges from ‘disorder.’ What does the science of such phase transformations tell us about the emergence of genomes? In this paper, we show that such thermodynamic considerations of RNA synthesis, when combined with kinetics and population dynamics, lead to the conclusion that the ‘crystallization’ of genomes from its basic elements would have been spontaneous for RNAs, but not necessarily for other potential building blocks of genomes in the prebiotic soup.

DNA polymerase activity at the single-molecule level

DNA polymerases are essential enzymes responsible for replication and repair of DNA in all organisms. To replicate DNA with high fidelity, DNA polymerases must select the correct incoming nucleotide substrate during each cycle of nucleotide incorporation, in accordance with the templating base. When an incorrect nucleotide is sometimes inserted, the polymerase uses a separate 3'→5' exonuclease to remove the misincorporated base (proofreading). Large conformational rearrangements of the polymerase-DNA complex occur during both the nucleotide incorporation and proofreading steps. Single-molecule fluorescence spectroscopy provides a unique tool for observation of these dynamic conformational changes in real-time, without the need to synchronize a population of DNA-protein complexes.

Human DNA polymerase beta mutations allowing efficient abasic site bypass

The DNA of every cell in the human body gets damaged more than 50,000 times a day. The most frequent damages are abasic sites. This kind of damage blocks proceeding DNA synthesis by several DNA polymerases that are involved in DNA replication and repair. The mechanistic basis for the incapability of these DNA polymerases to bypass abasic sites is not clarified. To gain insights into the mechanistic basis, we intended to identify amino acid residues that govern for the pausing of DNA polymerase β when incorporating a nucleotide opposite to abasic sites. Human DNA polymerase β was chosen because it is a well characterized DNA polymerase and serves as model enzyme for studies of DNA polymerase mechanisms. Moreover, it acts as the main gap-filling enzyme in base excision repair, and human tumor studies suggest a link between DNA polymerase β and cancer. In this study we employed high throughput screening of a library of more than 11,000 human DNA polymerase β variants. We identified two mutants that have increased ability to incorporate a nucleotide opposite to an abasic site. We found that the substitutions E232K and T233I promote incorporation opposite the lesion. In addition to this feature, the variants have an increased activity and a lower fidelity when processing nondamaged DNA. The mutations described in this work are located in well characterized regions but have not been reported before. A crystallographic structure of one of the mutants was obtained, providing structural insights.

Characterization of the elongation complex of dengue virus RNA polymerase: assembly, kinetics of nucleotide incorporation, and fidelity

Dengue virus (DENV) infects 50–100 million people worldwide per year, causing severe public health problems. DENV RNA-dependent RNA polymerase, an attractive target for drug development, catalyzes de novo replication of the viral genome in three phases: initiation, transition, and elongation. The aim of this work was to characterize the mechanism of nucleotide addition catalyzed by the polymerase domain of DENV serotype 2 during elongation using transient kinetic methods. We measured the kinetics of formation of the elongation complex containing the polymerase and a double-stranded RNA by preincubation experiments. The elongation complex assembly is slow, following a one-step binding mechanism with an association rate of 0.0016 ± 0.0001 μm⁻¹s⁻¹ and a dissociation rate of 0.00020 ± 0.00005 s⁻¹ at 37 °C. The elongation complex assembly is 6 times slower at 30 °C and requires Mg²⁺ during preincubation. The assembled elongation complex incorporates a correct nucleotide, GTP, to the primer with a K_d of 275 ± 52 μm and k_pol of 18 ± 1 s⁻¹. The fidelity of the polymerase is 1/34,000, 1/59,000, 1/135,000 for misincorporation of UTP, ATP, and CTP opposite CMP in the template, respectively. The fidelity of DENV polymerase is comparable with HIV reverse transcriptase and the poliovirus polymerase. This work reports the first description of presteady-state kinetics and fidelity for an RNA-dependent RNA polymerase from the Flaviviridae family.

Identification of critical residues for the tight binding of both correct and incorrect nucleotides to human DNA polymerase λ

DNA polymerase λ (Pol λ) is a novel X-family DNA polymerase that shares 34% sequence identity with DNA polymerase β. Pre-steady-state kinetic studies have shown that the Pol λ-DNA complex binds both correct and incorrect nucleotides 130-fold tighter, on average, than the DNA polymerase β-DNA complex, although the base substitution fidelity of both polymerases is 10(-)(4) to 10(-5). To better understand Pol λ's tight nucleotide binding affinity, we created single-substitution and double-substitution mutants of Pol λ to disrupt the interactions between active-site residues and an incoming nucleotide or a template base. Single-turnover kinetic assays showed that Pol λ binds to an incoming nucleotide via cooperative interactions with active-site residues (R386, R420, K422, Y505, F506, A510, and R514). Disrupting protein interactions with an incoming correct or incorrect nucleotide impacted binding to each of the common structural moieties in the following order: triphosphate≫base>ribose. In addition, the loss of Watson-Crick hydrogen bonding between the nucleotide and the template base led to a moderate increase in K(d). The fidelity of Pol λ was maintained predominantly by a single residue, R517, which has minor groove interactions with the DNA template.

Pre-steady state kinetic studies of the fidelity of nucleotide incorporation by yeast DNA polymerase delta

Eukaryotic DNA polymerase delta (pol delta) is a member of the B family of polymerases and synthesizes most of the lagging strand during DNA replication. Yeast pol delta is a heterotrimer comprised of three subunits: the catalytic subunit (Pol3) and two accessory subunits (Pol31 and Pol32). Although pol delta is one of the major eukaryotic replicative polymerase, the mechanism by which it incorporates nucleotides is unknown. Here we report both steady state and pre-steady state kinetic studies of the fidelity of pol delta. We found that pol delta incorporates nucleotides with an error frequency of 10(-4) to 10(-5). Furthermore, we showed that for correct versus incorrect nucleotide incorporation, there are significant differences between both pre-steady state kinetic parameters (apparent K(d)(dNTP) and k(pol)). Somewhat surprisingly, we found that pol delta synthesizes DNA at a slow rate with a k(pol) of approximately 1 s(-1). We suggest that, unlike its prokaryotic counterparts, pol delta requires replication accessory factors like proliferating cell nuclear antigen to achieve rapid rates of nucleotide incorporation.

Hinge residue I174 is critical for proper dNTP selection by DNA polymerase beta

DNA polymerase beta (pol beta) is the key gap-filling polymerase in base excision repair, the DNA repair pathway responsible for repairing up to 20000 endogenous lesions per cell per day. Pol beta is also widely used as a model polymerase for structure and function studies, and several structural regions have been identified as being critical for the fidelity of the enzyme. One of these regions is the hydrophobic hinge, a network of hydrophobic residues located between the palm and fingers subdomains. Previous work by our lab has shown that hinge residues Y265, I260, and F272 are critical for polymerase fidelity by functioning in discrimination of the correct from incorrect dNTP during ground state binding. Our work aimed to elucidate the role of hinge residue I174 in polymerase fidelity. To study this residue, we conducted a genetic screen to identify mutants with a substitution at residue I174 that resulted in a mutator polymerase. We then chose the mutator mutant I174S for further study and found that it follows the same general kinetic pathway as and has an overall protein folding similar to that of wild-type (WT) pol beta. Using single-turnover kinetic analysis, we found that I174S exhibits decreased fidelity when inserting a nucleotide opposite a template base G, and this loss of fidelity is due primarily to a loss of discrimination during ground state dNTP binding. Molecular dynamics simulations show that mutation of residue I174 to serine results in an overall tightening of the hinge region, resulting in aberrant protein dynamics and fidelity. These results point to the hinge region as being critical in the maintenance of the proper geometry of the dNTP binding pocket.

A binding free energy decomposition approach for accurate calculations of the fidelity of DNA polymerases

DNA polymerase beta (pol beta) is a small eukaryotic enzyme with the ability to repair short single-stranded DNA gaps that has found use as a model system for larger replicative DNA polymerases. For all DNA polymerases, the factors determining their catalytic power and fidelity are the interactions between the bases of the base pair, amino acids near the active site, and the two magnesium ions. In this report, we study effects of all three aspects on human pol beta transition state (TS) binding free energies by reproducing a consistent set of experimentally determined data for different structures. Our calculations comprise the combination of four different base pairs (incoming pyrimidine nucleotides incorporated opposite both matched and mismatched purines) with four different pol beta structures (wild type and three mutants). We generate three fragments of the incoming deoxynucleoside 5'-triphosphate-TS and run separate calculations for the neutral base part and the highly charged triphosphate part, using different dielectric constants in order to account for the specific dielectric response. This new approach improves our ability to predict the effect of matched and mismatched base pairing and of mutations in DNA polymerases on fidelity and may be a useful tool in studying the potential of DNA polymerase mutations in the development of cancer. It also supports our point of view with regards to the origin of the structural control of fidelity, allowing for a quantified description of the fidelity of DNA polymerases.

Relationship between conformational changes in pol lambda's active site upon binding incorrect nucleotides and mismatch incorporation rates

The correct replication and repair of DNA is critical for a cell's survival. Here, we investigate the fidelity of mammalian DNA polymerase lambda (pol lambda) utilizing dynamics simulation of the enzyme bound to incorrect incoming nucleotides including A:C, A:G, A(syn):G, A:A, A(syn):A, and T:G, all of which exhibit differing incorporation rates for pol lambda as compared to A:T bound to pol lambda. The wide range of DNA motion and protein residue side-chain motions observed in the mismatched systems demonstrates distinct differences when compared to the reference (correct base pair) system. Notably, Arg517's interactions with the DNA template strand bases in the active site are more limited, and Arg517 displays increased interactions with the incorrect dNTPs. This effect suggests that Arg517 helps provide a base-checking mechanism to discriminate correct from incorrect dNTPs. In addition, we find Tyr505 and Phe506 also play key roles in this base checking. A survey of the electrostatic potential landscape of the active sites and concomitant changes in electrostatic interaction energy between Arg517 and the dNTPs reveals that pol lambda binds incorrect dNTPs less tightly than the correct dNTP. These trends lead us to propose the following order for mismatch insertion by pol lambda: A:C > A:G > A(syn):G > T:G > A(syn):A > A:A. This sequence agrees with available kinetic data for incorrect nucleotide insertion opposite template adenine, with the exception of T:G, which may be more sensitive to the insertion context.

Dissection of the sequence-specific DNA binding and exonuclease activities reveals a superactive yet apoptotically impaired mutant p53 protein

Both sequence-specific DNA binding and exonuclease activities have been mapped to the central conserved core domain of p53. To gain more information about these two activities a series of mutants were generated that changed core domain histidine residues. Of these mutants, only one, H115N p53, showed markedly reduced exonuclease activity (ca. 15% of wild-type). Surprisingly, purified H115N p53 protein was found to be significantly more potent than wild-type p53 in binding to DNA by several criteria including gel mobility shift assay, filter binding and DNase I footprinting. Interestingly as well, non-specific DNA binding by the core domain of H115N p53 is superior to that of wild-type p53. To study H115N p53 in vivo, clones of H1299 cells expressing tetracycline regulated wild-type or H115N p53 were generated. H115N was both more potent than wild-type p53 in inducing p53 target genes such as p21 and PIG3 and was also more effective in arresting cells in G1. Unexpectedly, in contrast to wild-type p53, H115N p53 was markedly impaired in causing apoptosis when cells were subjected to DNA damage. Our results indicate that the exonuclease activity and transcriptional activation functions of p53 can be separated. They also extend previous findings showing that cell cycle arrest and apoptosis are separable functions of p53. Finally, these experiments confirm that DNA binding and xonuclease activities are distinct features of the p53 core domain.

RB69 DNA polymerase mutants with expanded nascent base-pair-binding pockets are highly efficient but have reduced base selectivity

We have investigated the effect of systematically enlarging the nascent base-pair-binding pocket (NBP) of a replicative DNA polymerase from bacteriophage RB69 (RB69 pol) on the incorporation efficiency (k(pol)/K(d,app)) for both correct and incorrect dNMPs. Accordingly, we replaced residues L561, Y567, and S565 in the NBP with Ala, Ala, and Gly, respectively. We combined L561A and Y567A to give a double mutant and then introduced the S565G mutation to give a triple mutant. The efficiency of incorrect dNMP insertion increased markedly relative to the wild type with the single mutants and increased further as the number of substitutions in the NBP increased. The difference in incorporation efficiency for mispairs between the mutants and the wild-type RB69 pol was due mainly to k(pol). Unexpectedly, enlarging the NBP had a minimal effect on the incorporation efficiency of correct dNMPs. Our kinetic data suggest that replicative DNA pols exert base discrimination via "negative selection" against mispairs by using residues in the NBP, particularly the three residues analyzed in this study, to allow rapid incorporation of only correct base pairs. This proposal differs from how geometry and "tightness of fit" of the NBP is often invoked to account for rapid incorporation of correct base pairs, namely, that a tighter fit within the NBP leads to an increase in insertion rates [Kool, E. T. (2002) Annu. Rev. Biochem. 71, 191-219]. We related our findings to that of a model translesion DNA pol, Sulfolobus solfataricus Dpo4. We concur with the main conclusion of a previous study [Mizukami, S., et al. (2006) Biochemistry 45, 2772-2778], namely, that lesion bypass pols exhibit low incorporation efficiencies for correct dNMPs (leading to relative promiscuity) not because of a more open NBP but because of a loose fit of substrates bound in the catalytic centers. This is a property not shared by RB69 pol and its mutants.

DNA polymerase family X: function, structure, and cellular roles

The X family of DNA polymerases in eukaryotic cells consists of terminal transferase and DNA polymerases beta, lambda, and mu. These enzymes have similar structural portraits, yet different biochemical properties, especially in their interactions with DNA. None of these enzymes possesses a proofreading subdomain, and their intrinsic fidelity of DNA synthesis is much lower than that of a polymerase that functions in cellular DNA replication. In this review, we discuss the similarities and differences of three members of Family X: polymerases beta, lambda, and mu. We focus on biochemical mechanisms, structural variation, fidelity and lesion bypass mechanisms, and cellular roles. Remarkably, although these enzymes have similar three-dimensional structures, their biochemical properties and cellular functions differ in important ways that impact cellular function.

Multifactor dimensionality reduction analysis identifies specific nucleotide patterns promoting genetic polymorphisms

Background

The fidelity of DNA replication serves as the nidus for both genetic evolution and genomic instability fostering disease. Single nucleotide polymorphisms (SNPs) constitute greater than 80% of the genetic variation between individuals. A new theory regarding DNA replication fidelity has emerged in which selectivity is governed by base-pair geometry through interactions between the selected nucleotide, the complementary strand, and the polymerase active site. We hypothesize that specific nucleotide combinations in the flanking regions of SNP fragments are associated with mutation.

Results

We modeled the relationship between DNA sequence and observed polymorphisms using the novel multifactor dimensionality reduction (MDR) approach. MDR was originally developed to detect synergistic interactions between multiple SNPs that are predictive of disease susceptibility. We initially assembled data from the Broad Institute as a pilot test for the hypothesis that flanking region patterns associate with mutagenesis (n = 2194). We then confirmed and expanded our inquiry with human SNPs within coding regions and their flanking sequences collected from the National Center for Biotechnology Information (NCBI) database (n = 29967) and a control set of sequences (coding region) not associated with SNP sites randomly selected from the NCBI database (n = 29967). We discovered seven flanking region pattern associations in the Broad dataset which reached a minimum significance level of p ≤ 0.05. Significant models (p << 0.001) were detected for each SNP type examined in the larger NCBI dataset. Importantly, the flanking region models were elongated or truncated depending on the nucleotide change. Additionally, nucleotide distributions differed significantly at motif sites relative to the type of variation observed. The MDR approach effectively discerned specific sites within the flanking regions of observed SNPs and their respective identities, supporting the collective contribution of these sites to SNP genesis.

Conclusion

The present study represents the first use of this computational methodology for modeling nonlinear patterns in molecular genetics. MDR was able to identify distinct nucleotide patterning around sites of mutations dependent upon the observed nucleotide change. We discovered one flanking region set that included five nucleotides clustered around a specific type of SNP site. Based on the strongly associated patterns identified in this study, it may become possible to scan genomic databases for such clustering of nucleotides in order to predict likely sites of future SNPs, and even the type of polymorphism most likely to occur.

Mismatched and matched dNTP incorporation by DNA polymerase beta proceed via analogous kinetic pathways

While matched nucleotide incorporation by DNA polymerase beta (Pol beta) has been well-studied, a true understanding of polymerase fidelity requires comparison of both matched and mismatched dNTP incorporation pathways. Here we examine the mechanism of misincorporation for wild-type (WT) Pol beta and an error-prone I260Q variant using stopped-flow fluorescence assays and steady-state fluorescence spectroscopy. In stopped-flow, a biphasic fluorescence trace is observed for both enzymes during mismatched dNTP incorporation. The fluorescence transitions are in the same direction as that observed for matched dNTP, albeit with lower amplitude. Assignments of the fast and slow fluorescence phases are designated to the same mechanistic steps previously determined for matched dNTP incorporation. For both WT and I260Q mismatched dNTP incorporation, the rate of the fast phase, reflecting subdomain closing, is comparable to that induced by correct dNTP. Pre-steady-state kinetic evaluation reveals that both enzymes display similar correct dNTP insertion profiles, and the lower fidelity intrinsic to the I260Q mutant results from enhanced efficiency of mismatched incorporation. Notably, in comparison to WT, I260Q demonstrates enhanced intensity of fluorescence emission upon mismatched ternary complex formation. Both kinetic and steady-state fluorescence data suggest that relaxed discrimination against incorrect dNTP by I260Q is a consequence of a loss in ability to destabilize the mismatched ternary complex. Overall, our results provide first direct evidence that mismatched and matched dNTP incorporations proceed via analogous kinetic pathways, and support our standing hypothesis that the fidelity of Pol beta originates from destabilization of the mismatched closed ternary complex and chemical transition state.

An intramolecular FRET system monitors fingers subdomain opening in Klentaq1

A major goal of polymerase research is to determine the mechanism through which a nucleotide complementary to a templating DNA base is selected and delivered to the polymerase active site. Structural evidence suggests a large open-to-closed conformational change affecting the fingers subdomain as being crucial to the process. We previously designed a FRET system capable of measuring the rate of fingers subdomain closure in the presence of correct nucleotide. However, this FRET system was limited in that it could not directly measure the rate of fingers subdomain opening by FRET after polymerization or in the absence of DNA. Here we report the development of a new system capable of measuring both fingers subdomain closure and reopening by FRET, and show that the rate of fingers subdomain opening is limited only by the rate of polymerization. We anticipate that this system will scale down to the single molecule level, allowing measurement of fingers subdomain movements in the presence of incorrect nucleotide and in the absence of DNA.

Incorrect nucleotide insertion at the active site of a G:A mismatch catalyzed by DNA polymerase beta

Based on a recent ternary complex crystal structure of human DNA polymerase beta with a G:A mismatch in the active site, we carried out a theoretical investigation of the catalytic mechanism of incorrect nucleotide incorporation using molecular dynamics simulation, quantum mechanics, combined quantum mechanics, and molecular mechanics methods. A two-stage mechanism is proposed with a nonreactive active-site structural rearrangement prechemistry step occurring before the nucleotidyl transfer reaction. The free energy required for formation of the prechemistry state is found to be the major factor contributing to the decrease in the rate of incorrect nucleotide incorporation compared with correct insertion and therefore to fidelity enhancement. Hence, the transition state and reaction barrier for phosphodiester bond formation after the prechemistry state are similar to that for correct insertion reaction. Key residues that provide electrostatic stabilization of the transition state are identified.

Mechanistic consequences of temperature on DNA polymerization catalyzed by a Y-family DNA polymerase

Our previous publication shows that Sulfolobus solfataricus Dpo4 utilizes an ‘induced-fit’ mechanism to select correct incoming nucleotides at 37°C. Here, we provide a comprehensive report elucidating the kinetic mechanism of a DNA polymerase at a reaction temperature higher than 37°C in an attempt to determine the effect of temperature on enzyme fidelity and mechanism. The fidelity of Dpo4 did not change considerably with a 30°C increase in reaction temperature, suggesting that the fidelity of Dpo4 at 80°C is similar to that determined here at 56°C. Amazingly, the incorporation rate for correct nucleotides increased by 18 900-fold from 2°C to 56°C, similar in magnitude to that observed for incorrect nucleotides, thus not perturbing fidelity. Three independent lines of kinetic evidence indicate that a protein conformational change limits correct nucleotide incorporations at 56°C. Furthermore, the activation energy for the incorporation of a correct nucleotide was determined to be 32.9 kcal/mol, a value considerably larger than those values estimated for a rate-limiting chemistry step, providing a fourth line of evidence to further substantiate this conclusion. These results herein provide evidence that Dpo4 utilizes the ‘induced-fit’ mechanism to select a correct nucleotide at all temperatures.

A pre-equilibrium before nucleotide binding limits fingers subdomain closure by Klentaq1

Numerous studies have been undertaken to establish the mechanism of dNTP binding and template-directed incorporation by DNA polymerases. It has been established by kinetic experiments that a rate-limiting step, crucial for dNTP selection, occurs before chemical bond formation. Crystallographic studies indicated that this step may be due to a large open-to-closed conformational transition affecting the fingers subdomain. In previous studies, we established a fluorescence resonance energy transfer system to monitor the open-to-closed transition in the fingers subdomain of Klentaq1. By comparing the rates of the fingers subdomain closure with that of the rate-limiting step for Klentaq1, we showed that fingers subdomain motion was significantly faster than the rate-limiting step. We have now used this system to characterize DNA binding as well as to complete a more extensive characterization of incorporation of all four dNTPs. The data indicate that DNA binding occurs by a two-step association and that dissociation of the DNA is significantly slower in the case of the closed ternary complex. The data for nucleotide incorporation indicate a step occurring before dNTP binding, which differs for all four nucleotides. As the only difference between the (E x p/t) complexes is the templating base, it would suggest an important role for the templating base in initial ground state selection.

A unified kinetic mechanism applicable to multiple DNA polymerases

After extensive studies spanning over half a century, there is little consensus on the kinetic mechanism of DNA polymerases. Using stopped-flow fluorescence assays for mammalian DNA polymerase beta (Pol beta), we have previously identified a fast fluorescence transition corresponding to conformational closing, and a slow fluorescence transition matching the rate of single-nucleotide incorporation. Here, by varying pH and buffer viscosity, we have decoupled the rate of single-nucleotide incorporation from the rate of the slow fluorescence transition, thus confirming our previous hypothesis that this transition represents a conformational event after chemistry, likely subdomain reopening. Analysis of an R258A mutant indicates that rotation of the Arg258 side chain is not rate-limiting in the overall kinetic pathway of Pol beta, yet is kinetically significant in subdomain reopening. We have extended our kinetic analyses to a high-fidelity polymerase, Klenow fragment (KF), and a low-fidelity polymerase, African swine fever virus DNA polymerase X (Pol X), and showed that they follow the same kinetic mechanism as Pol beta, while differing in relative rates of single-nucleotide incorporation and the putative conformational reopening. Our data suggest that the kinetic mechanism of Pol beta is not an exception among polymerases, and furthermore, its delineated kinetic mechanism lends itself as a platform for comparison of the kinetic properties of different DNA polymerases and their mutants.

Single-turnover kinetic analysis of the mutagenic potential of 8-oxo-7,8-dihydro-2'-deoxyguanosine during gap-filling synthesis catalyzed by human DNA polymerases lambda and beta

In the presence of 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxodG) damage, many DNA polymerases exhibit a dual coding potential which facilitates efficient incorporation of matched dCTP or mismatched dATP. This also holds true for the insertion of 8-oxodGTP opposite template bases dC and dA. Employing single-turnover kinetic methods, we examined human DNA polymerase beta and its novel X-family homolog, human DNA polymerase lambda, to determine which nucleotide and template base was preferred when encountering 8-oxodG and 8-oxodGTP, respectively. While DNA polymerase beta preferentially incorporated dCTP over dATP, DNA polymerase lambda did not modulate a preference for either dCTP or dATP when opposite 8-oxodG in single-nucleotide gapped DNA, as incorporation proceeded with essentially equal efficiency and probability. Moreover, DNA polymerase lambda is more efficient than DNA polymerase beta to fill this oxidized single-nucleotide gap. Insertion of 8-oxodGTP by both DNA polymerases lambda and beta occurred predominantly against template dA, thereby reiterating how the asymmetrical design of the polymerase active site differentially accommodated the anti and syn conformations of 8-oxodG and 8-oxodGTP. Although the electronegative oxygen at the C8 position of 8-oxodG may induce DNA structural perturbations, human DNA ligase I was found to effectively ligate the incorporated 8-oxodGMP to a downstream strand, which sealed the nicked DNA. Consequently, the erroneous nucleotide incorporations catalyzed by DNA polymerases lambda and beta as well as the subsequent ligation catalyzed by a DNA ligase during base excision repair are a threat to genomic integrity.

DNA polymerase beta catalytic efficiency mirrors the Asn279-dCTP H-bonding strength

Ternary complexes of wild type or mutant form of human DNA polymerase beta (pol beta) bound to DNA and dCTP substrates were studied by molecular dynamics (MD) simulations. The occurrences of contact configurations (CC) of structurally important atom pairs were sampled along the MD trajectories, and converted into free-energy differences, DeltaG(CC). DeltaG(CC) values were correlated with the experimental binding and catalytic free energies for the wild type pol beta and its Arg183Ala, Tyr271Ala, Asp276Val, Lys280Gly, Arg283Ala, and Glu295Ala mutants. The correlation coefficients show that the strength of the H-bond between dCTP and Asn279 is a strong predictor of the mutation-induced changes in the catalytic efficiency of pol beta. This finding is consistent with the view that enzyme preorganization plays a major role in controlling DNA polymerase specific activity.

DNA polymerase beta catalytic efficiency mirrors the Asn279-dCTP H-bonding strength

Ternary complexes of wild type or mutant form of human DNA polymerase beta (pol beta) bound to DNA and dCTP substrates were studied by molecular dynamics (MD) simulations. The occurrences of contact configurations (CC) of structurally important atom pairs were sampled along the MD trajectories, and converted into free-energy differences, DeltaG(CC). DeltaG(CC) values were correlated with the experimental binding and catalytic free energies for the wild type pol beta and its Arg183Ala, Tyr271Ala, Asp276Val, Lys280Gly, Arg283Ala, and Glu295Ala mutants. The correlation coefficients show that the strength of the H-bond between dCTP and Asn279 is a strong predictor of the mutation-induced changes in the catalytic efficiency of pol beta. This finding is consistent with the view that enzyme preorganization plays a major role in controlling DNA polymerase specific activity.

Regulation of DNA repair fidelity by molecular checkpoints: "gates" in DNA polymerase beta's substrate selection

With an increasing number of structural, kinetic, and modeling studies of diverse DNA polymerases in various contexts, a complex dynamical view of how atomic motions might define molecular "gates" or checkpoints that contribute to polymerase specificity and efficiency is emerging. Such atomic-level information can offer insights into rate-limiting conformational and chemical steps to help piece together mechanistic views of polymerases in action. With recent advances, modeling and dynamics simulations, subject to the well-appreciated limitations, can access transition states and transient intermediates along a reaction pathway, both conformational and chemical, and such information can help bridge the gap between experimentally determined equilibrium structures and mechanistic enzymology data. Focusing on DNA polymerase beta (pol beta), we present an emerging view of the geometric, energetic, and dynamic selection criteria governing insertion rate and fidelity mechanisms of DNA polymerases, as gleaned from various computational studies and based on the large body of existing kinetic and structural data. The landscape of nucleotide insertion for pol beta includes conformational changes, prechemistry, and chemistry "avenues", each with a unique deterministic or stochastic pathway that includes checkpoints for selective control of nucleotide insertion efficiency. For both correct and incorrect incoming nucleotides, pol beta's conformational rearrangements before chemistry include a cascade of slow and subtle side chain rearrangements, followed by active site adjustments to overcome higher chemical barriers, which include critical ion-polymerase geometries; this latter notion of a prechemistry avenue fits well with recent structural and NMR data. The chemical step involves an associative mechanism with several possibilities for the initial proton transfer and for the interaction among the active site residues and bridging water molecules. The conformational and chemical events and associated barriers define checkpoints that control enzymatic efficiency and fidelity. Understanding the nature of such active site rearrangements can facilitate interpretation of existing data and stimulate new experiments that aim to probe enzyme features that contribute to fidelity discrimination across various polymerases via such geometric, dynamic, and energetic selection criteria.