Modifying the beta,gamma leaving-group bridging oxygen alters nucleotide incorporation efficiency, fidelity, and the catalytic mechanism of DNA polymerase beta

DNA polymerase catalysis and fidelity studies typically compare incorporation of "right" versus "wrong" nucleotide bases where the leaving group is pyrophosphate. Here we use dGTP analogues replacing the beta,gamma-bridging O with CH2, CHF, CF2, or CCl2 to explore leaving-group effects on the nucleotidyl transfer mechanism and fidelity of DNA polymerase (pol) beta. T.G mismatches occur with fidelities similar to dGTP with the exception of the CH2 analogue, which is incorporated with 5-fold higher fidelity. All analogues are observed to bind opposite template C with Kds between 1 and 4 microM, and structural evidence suggests that the analogues bind in essentially the native conformation, making them suitable substrates for probing linear free energy relationships (LFERs) in transient-kinetics experiments. Importantly, Brnsted correlations of log(kpol) versus leaving-group pKa for both right and wrong base incorporation reveal similar sensitivities (betalg approximately -0.8) followed by departures from linearity, suggesting that a chemical step rather than enzyme conformational change is rate-limiting for either process. The location of the breaks relative to pKas of CF2, O, and the sterically bulky CCl2-bridging compounds suggests a modification-induced change in the mechanism by stabilization of leaving-group elimination. The results are addressed theoretically in terms of the energetics of successive primer 3'-O addition (bond forming) and pyrophosphate analogue elimination (bond breaking) reaction energy barriers.     

Arg660Ser mutation in Thermus aquaticus DNA polymerase I suppresses T-->C transitions: implication of wobble base pair formation at the nucleotide incorporation step

We examined the replication fidelity of an Arg660Ser (R660S) mutant of Thermus aquaticus DNA polymerase I (Taq pol I). In a forward mutation assay, R660S showed a marked reduction in T→C transitions, one of the most frequent errors made by the wild-type enzyme. Steady-state kinetics showed that R660S discriminates against dGTP incorporation at a template T 13-fold better than the wild-type. R660S was also 3.2-fold less efficient than the wild-type at extending a T:dG mismatch. These results indicate that R660S has enhanced fidelity during incorporation and extension, which reduces its T→C transition frequency. Interestingly, R660S also discriminated correct from incorrect nucleotides at the incorporation step of C:dATP, A:dATP, G:dATP and C:8-OH-dGTP mispairs 28-, 6.0-, 4.1- and 6.8-fold better, respectively, than the wild-type, although it may not always be as accurate as the wild-type at the extension step. A structural model suggests that Arg660 may participate in two interactions that influence fidelity; the guanidinium group of Arg660 might interact with the incoming guanine base at the major groove and it might compete for forming another interaction with the primer terminus. Substituting Arg with Ser may eliminate or alter these interactions and destabilize the closed complex with incorrect substrates. Our data also suggest that T:dGTP and C:dATP base pairs form ‘wobble’ structures at the incorporation step of Taq pol I.     

Fidelity of the human mitochondrial DNA polymerase

We have quantified the fidelity of polymerization of DNA by human mitochondrial DNA polymerase using synthetic DNA oligonucleotides and recombinant holoenzyme and examining each of the possible 16-base pair combinations. Although the kinetics of incorporation for all correct nucleotides are similar, with an average Kd of 0.8 microM and an average k(pol) of 37 s(-1), the kinetics of misincorporation vary widely. The ground state binding Kd of incorrect bases ranges from a low of 25 microM for a dATP:A mispair to a high of 360 microM for a dCTP:T mispair. Similarly, the rates of incorporation of incorrect bases vary from 0.0031 s(-1) for a dCTP:C mispair to 1.16 s(-1) for a dGTP:T mispair. Due to the variability in the kinetic parameters for misincorporation, the estimates of fidelity range from 1 error in 3563 nucleotides for dGTP:T to 1 error in 2.3 x 10(6) nucleotides for dCTP:C. Interestingly, the discrimination against a dGTP:T mismatch is 16.5 times lower than that of a dTTP:G mismatch due to a tighter Kd for ground state binding and a faster rate of incorporation of the dGTP:T mismatch relative to the dTTP:G mismatch. We calculate an average fidelity of 1 error in 440,000 nucleotides.     

A reexamination of the nucleotide incorporation fidelity of DNA polymerases

Intensive study has been devoted to understanding the kinetic and structural bases underlying the exceptionally high fidelity (low error frequencies) of the typical DNA polymerase. Commonly proposed explanations have included (i) the concept of fidelity check points, in which the correctness of a nascent base pair match is tested at multiple points along the reaction pathway, and (ii) an induced-fit fidelity enhancement mechanism based on a rate-limiting, substrate-induced conformational change. In this article, we consider the evidence and theoretical framework for the involvement of such mechanisms in fidelity enhancement. We suggest that a "simplified" model, in which fidelity is derived fundamentally from differential substrate binding at the transition state of a rate-limiting chemical step, is consistent with known data and sufficient to explain the substrate selectivity of these enzymes.     

Altering DNA polymerase incorporation fidelity by distorting the dNTP binding pocket with a bulky carcinogen-damaged template

Fidelity of DNA polymerases is predominantly governed by an induced fit mechanism in which the incoming dNTP in the ternary complex fits tightly into a binding pocket whose geometry is determined by the nature of the templating base. However, modification of the template with a bulky carcinogen may alter the dNTP binding pocket and thereby the polymerase incorporation fidelity. High fidelity DNA polymerases, such as bacteriophage T7 DNA polymerase, are predominantly blocked by bulky chemical lesions on the template strand during DNA replication. However, some mutagenic bypass can occur, which may lead to carcinogenesis. Experimental studies have shown that a DNA covalent adduct derived from (+)-anti-BPDE [(+)-(7R,8S,9S,10R)-7,8-dihydroxy-9,10-epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene], a carcinogenic metabolite of benzo[a]pyrene (BP), primarily blocks Sequenase 2.0, an exo(-) T7 DNA polymerase; however, a mismatched dATP can be preferentially inserted opposite the damaged adenine templating base within the active site of the polymerase [Chary, P., and Lloyd, R. S. (1995) Nucleic Acids Res. 23, 1398-1405]. The goal of this work is to elucidate structural features that contribute to DNA polymerase incorporation fidelity in the presence of this bulky covalent adduct and to interpret the experimental findings on a molecular level. We have carried out molecular modeling and molecular dynamics simulations with AMBER 6.0, investigating a T7 DNA polymerase primer-template closed ternary complex containing this 10S (+)-trans-anti-[BP]-N(6)-dA adduct in the templating position within the polymerase active site. All four incoming dNTPs were studied. The simulations show that the BP ring system fits well into an open pocket on the major groove side of the modified template adenine with anti glycosidic bond conformation, without disturbing critical polymerase-DNA interactions. However, steric hindrance between the BP ring system and the primer-template DNA causes displacement of the modified template adenine, so that the dNTP base binding pocket is enlarged. This alteration can explain the experimentally observed preference for incorporation of dATP opposite this lesion. These studies also rationalize the observed lower probabilities of incorporation of the other three nucleotides. Our results suggest that the differences in incorporation of dGTP, dCTP, and dTTP are due to the effects of imperfect geometric complementarity. Thus, the simulations suggest that altered DNA polymerase incorporation fidelity can result from adduct-induced changes in the dNTP base binding pocket geometry. Furthermore, plausible structural explanations for the observed effects of [BP]-N(6)-dA adduct stereochemistry on the observed stalling patterns are proposed.     

The properties of steric gate mutants reveal different constraints within the active sites of Y-family and A-family DNA polymerases

Y-family (lesion-bypass) DNA polymerases show the same overall structural features seen in other members of the polymerase superfamily, yet their active sites are more open, with fewer contacts to the DNA and nucleotide substrates. This raises the question of whether analogous active-site side chains play equivalent roles in the bypass polymerases and their classical DNA polymerase counterparts. In Klenow fragment, an A-family DNA polymerase, the steric gate side chain (Glu710) not only prevents ribonucleotide incorporation but also plays an important role in discrimination against purine-pyrimidine mispairs. In this work we show that the steric gate (Phe12) of the Y-family polymerase Dbh plays a very minor role in fidelity, despite its analogous role in sugar selection. Using ribonucleotide discrimination to report on the positioning of a mispaired dNTP, we found that the pyrimidine of a Pu-dPyTP nascent mispair occupies a similar position to that of a correctly paired dNTP in the Dbh active site, whereas in Klenow fragment the mispaired dNTP sits higher in the active site pocket. If purine-pyrimidine mispairs adopt the expected wobble geometry, the difference between the two polymerases can be attributed to the binding of the templating base, with the looser binding site of Dbh permitting a variety of template conformations with only minimal adjustment at the incoming dNTP. In Klenow fragment the templating base is more rigidly held, so that changes in base pair geometry would affect the dNTP position, allowing the Glu710 side chain to serve as a sensor of nascent mispairs.     

Structural basis for differential insertion kinetics of dNMPs opposite a difluorotoluene nucleotide residue

We have recently challenged the widely held view that 2,4-difluorotoluene (dF) is a nonpolar isosteric analogue of the nucleotide dT, incapable of forming hydrogen bonds (HBs). To gain a further understanding for the kinetic preference that favors dAMP insertion opposite a templating dF, a result that mirrors the base selectivity that favors dAMP insertion opposite dT by RB69 DNA polymerase (RB69pol), we determined presteady-state kinetic parameters for incorporation of four dNMPs opposite dF by RB69pol and solved the structures of corresponding ternary complexes. We observed that both the F2 and F4 substituent of dF in these structures serve as HB acceptors forming HBs either directly with dTTP and dGTP or indirectly with dATP and dCTP via ordered water molecules. We have defined the shape and chemical features of each dF/dNTP pair in the RB69pol active site without the corresponding phosphodiester-linkage constraints of dF/dNs when they are embedded in isolated DNA duplexes. These features can explain the kinetic preferences exhibited by the templating dF when the nucleotide incorporation is catalyzed by wild type RB69pol or its mutants. We further show that the shapes of the dNTP/dF nascent base pair differ markedly from the corresponding dNTP/dT in the pol active site and that these differences have a profound effect on their incorporation efficiencies.     

DNA structure and polymerase fidelity.

The accuracy of DNA replication results from both the intrinsic DNA polymerase fidelity and the DNA sequence. Although the recent structural studies on polymerases have brought new insights on polymerase fidelity, the role of DNA sequence and structure is less well understood. Here, the analysis of the crystal structures of hotspots for polymerase slippage including (CA)n and (A)n tracts in different intermolecular contexts reveals that, in the B-form, these sequences share common structural alterations which may explain the high rate of replication errors. In particular, a two-faced "Janus-like" structure with shifted base-pairs in the major groove but an apparent normal geometry in the minor groove constitutes a molecular decoy specifically suitable to mislead the polymerases. A model of the rat polymerase beta bound to this structure suggests that an altered conformation of the nascent template-primer duplex can interfere with correct nucleotide incorporation by affecting the geometry of the active site and breaking the rules of base-pairing, while at the same time escaping enzymatic mechanisms of error discrimination which scan for the correct geometry of the minor groove.In contrast, by showing that the A-form greatly attenuates the sequence-dependent structural alterations in hotspots, this study suggests that the A-conformation of the nascent template-primer duplex at the vicinity of the polymerase active site will contribute to fidelity. The A-form may play the role of a structural buffer which preserves the correct geometry of the active site for all sequences. The detailed comparison of the conformation of the nascent template-primer duplex in the available crystal structures of DNA polymerase-DNA complexes shows that polymerase beta, the least accurate enzyme, is unique in binding to a B-DNA duplex even close to its active site. This model leads to several predictions which are discussed in the light of published experimental data.     

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.     

Energetic mapping of transition state analogue interactions with human and Plasmodium falciparum purine nucleoside phosphorylases

Human purine nucleoside phosphorylase (huPNP) is essential for human T-cell division by removing deoxyguanosine and preventing dGTP imbalance. Plasmodium falciparum expresses a distinct PNP (PfPNP) with a unique substrate specificity that includes 5'-methylthioinosine. The PfPNP functions both in purine salvage and in recycling purine groups from the polyamine synthetic pathway. Immucillin-H is an inhibitor of both huPNP and PfPNPs. It kills activated human T-cells and induces purine-less death in P. falciparum. Immucillin-H is a transition state analogue designed to mimic the early transition state of bovine PNP. The DADMe-Immucillins are second generation transition state analogues designed to match the fully dissociated transition states of huPNP and PfPNP. Immucillins, DADMe-Immucillins and related analogues are compared for their energetic interactions with human and P. falciparum PNPs. Immucillin-H and DADMe-Immucillin-H are 860 and 500 pM inhibitors against P. falciparum PNP but bind human PNP 15-35 times more tightly. This common pattern is a result of kcat for huPNP being 18-fold greater than kcat for PfPNP. This energetic binding difference between huPNP and PfPNP supports the k(chem)/kcat binding argument for transition state analogues. Preferential PfPNP inhibition is gained in the Immucillins by 5'-methylthio substitution which exploits the unique substrate specificity of PfPNP. Human PNP achieves part of its catalytic potential from 5'-OH neighboring group participation. When PfPNP acts on 5'-methylthioinosine, this interaction is not possible. Compensation for the 5'-OH effect in the P. falciparum enzyme is provided by improved leaving group interactions with Asp206 as a general acid compared with Asn at this position in huPNP. Specific atomic modifications in the transition state analogues cause disproportionate binding differences between huPNP and PfPNPs and pinpoint energetic binding differences despite similar transition states.     

Conformational dynamics of DNA polymerase probed with a novel fluorescent DNA base analogue

DNA polymerases discriminate between correct and incorrect nucleotide substrates during a "nonchemical" step that precedes phosphodiester bond formation in the enzymatic cycle of nucleotide incorporation. Despite the importance of this process in polymerase fidelity, the precise nature of the molecular events involved remains unknown. Here we report a fluorescence resonance energy transfer (FRET) system that monitors conformational changes of a polymerase-DNA complex during selection and binding of nucleotide substrates. This system utilizes the fluorescent base analogue 1,3-diaza-2-oxophenothiazine (tC) as the FRET donor and Alexa-555 (A555) as the acceptor. The tC donor was incorporated within a model DNA primer/template in place of a normal base, adjacent to the primer 3' terminus, while the A555 acceptor was attached to an engineered cysteine residue (C751) located in the fingers subdomain of the Klenow fragment (KF) polymerase. The FRET efficiency increased significantly following binding of a correct nucleotide substrate to the KF-DNA complex, showing that the fingers had closed over the active site. Fluorescence anisotropy titrations utilizing tC as a reporter indicated that the DNA was more tightly bound by the polymerase under these conditions, consistent with the formation of a closed ternary complex. The rate of the nucleotide-induced conformational transition, measured in stopped-flow FRET experiments, closely matched the rate of correct nucleotide incorporation, measured in rapid quench-flow experiments, indicating that the conformational change was the rate-limiting step in the overall cycle of nucleotide incorporation for the labeled KF-DNA system. Taken together, these results indicate that the FRET system can be used to probe enzyme conformational changes that are linked to the biochemical function of DNA polymerase.     

DNA polymerase catalysis in the absence of Watson-Crick hydrogen bonds: analysis by single-turnover kinetics

We report the first pre-steady-state kinetic studies of DNA replication in the absence of hydrogen bonds. We have used nonpolar nucleotide analogues that mimic the shape of a Watson-Crick base pair to investigate the kinetic consequences of a lack of hydrogen bonds in the polymerase reaction catalyzed by the Klenow fragment of DNA polymerase I from Escherichia coli. With a thymine isostere lacking hydrogen-bonding ability in the nascent pair, the efficiency (k(pol)/Kd) of the polymerase reaction is decreased by 30-fold, affecting the ground state (Kd) and transition state (k(pol)) approximately equally. When both thymine and adenine analogues in the nascent pair lack hydrogen-bonding ability, the efficiency of the polymerase reaction is decreased by about 1000-fold, with most of the decrease attributable to the transition state. Reactions using nonpolar analogues at the primer-terminal base pair demonstrated the requirement for a hydrogen bond between the polymerase and the minor groove of the primer-terminal base. The R668A mutation of Klenow fragment abolished this requirement, identifying R668 as the probable hydrogen-bond donor. Detailed examination of the kinetic data suggested that Klenow fragment has an extremely low tolerance of even minor deviations of the analogue base pairs from ideal Watson-Crick geometry. Consistent with this idea, some analogue pairings were better tolerated by Klenow fragment mutants having more spacious active sites. In contrast, the Y-family polymerase Dbh was much less sensitive to changes in base pair dimensions and more dependent upon hydrogen bonding between base-paired partners.     

Kinetic analysis of the unique error signature of human DNA polymerase ν

The fidelity of DNA synthesis by A-family DNA polymerases ranges from very accurate for bacterial, bacteriophage, and mitochondrial family members to very low for certain eukaryotic homologues. The latter include DNA polymerase ν (Pol ν) which, among all A-family polymerases, is uniquely prone to misincorporating dTTP opposite template G in a highly sequence-dependent manner. Here we present a kinetic analysis of this unusual error specificity, in four different sequence contexts and in comparison to Pol ν's more accurate A-family homologue, the Klenow fragment of Escherichia coli DNA polymerase I. The kinetic data strongly correlate with rates of stable misincorporation during gap-filling DNA synthesis. The lower fidelity of Pol ν compared to that of Klenow fragment can be attributed primarily to a much lower catalytic efficiency for correct dNTP incorporation, whereas both enzymes have similar kinetic parameters for G-dTTP misinsertion. The major contributor to sequence-dependent differences in Pol ν error rates is the reaction rate, k(pol). In the sequence context where fidelity is highest, k(pol) for correct G-dCTP incorporation by Pol ν is ~15-fold faster than k(pol) for G-dTTP misinsertion. However, in sequence contexts where the error rate is higher, k(pol) is the same for both correct and mismatched dNTPs, implying that the transition state does not provide additional discrimination against misinsertion. The results suggest that Pol ν may be fine-tuned to function when high enzyme activity is not a priority and may even be disadvantageous and that the relaxed active-site specificity toward the G-dTTP mispair may be associated with its cellular function(s).     

Loop II of DNA polymerase beta is important for discrimination during substrate binding

Loop II of DNA polymerase beta (pol β) consists of 14 amino acid residues and is highly flexible and solvent exposed. Previous research from our laboratory has shown that this loop is important for polymerase activity and fidelity. In the study presented here, we demonstrate that a shortened five amino acid residue loop compromises the fidelity of pol β. This five-residue loop, termed ENEYP, induces one base frameshift errors and A–C transversions within a specific sequence context. We demonstrate that ENEYP misincorporates dGTP opposite template A at higher efficiencies than wild-type pol β. The kinetic basis for misincorporation is a defect in discrimination of the correct from incorrect dNTP substrate at the level of ground-state binding. Our results are consistent with the idea that loop II of pol β functions to maintain accurate DNA synthesis by a direct or indirect influence on the nucleotide binding pocket.     

Effect of the Y955C mutation on mitochondrial DNA polymerase nucleotide incorporation efficiency and fidelity

The human mitochondrial DNA polymerase (pol γ) is responsible for the replication of the mitochondrial genome. Mutation Y955C in the active site of pol γ results in early onset progressive external ophthalmoplegia, premature ovarian failure, and Parkinson's disease. In single-turnover kinetic studies, we show that the Y955C mutation results in a decrease in the maximal rate of polymerization and an increase in the K(m) for correct incorporation. The mutation decreased the specificity constant for correct incorporation of dGTP, TTP, and ATP to values of 1.5, 0.35, and 0.044 μM(-1) s(-1), respectively, representing reductions of 30-, 110-, and 1300-fold, respectively, relative to the value for the wild-type enzyme. The fidelity of incorporation was reduced 6-130-fold, largely because of the significant decrease in the specificity constant for correct dATP:T incorporation. For example, k(cat)/K(m) for forming a TTP:T mismatch was decreased 10-fold from 0.0002 to 0.00002 μM(-1) s(-1) by the Y955C mutant, but the 1300-fold slower incorporation of the correct dATP:T relative to that of the wild type led to a 130-fold lower fidelity. While correct incorporation of 8-oxo-dGTP was largely unchanged, the level of incorporation of 8-oxo-dG with dA in the template strand was reduced 500-fold. These results support a role for Y955 in stabilizing A:T base pairs at the active site of pol γ and suggest that the severe clinical symptoms of patients with this mutation may be due, in part, to the reduced efficiency of incorporation of dATP opposite T, and that the autosomal dominant phenotype may arise from the resulting higher mutation frequency.