Studies on the biochemical basis of spontaneous mutation. III. Rate model for DNA polymerase-effected nucleotide misincorporation

We propose a model to investigate the relation between insertion and excision activities of polymerases involved in DNA synthesis, and the frequency of errors resulting from substituting either mismatched bases or base analogues into a DNA molecule. An analytical equation is derived which expresses the error frequency as a function of nucleotide insertion and removal rates. For the general case, given arbitrary rates of insertion and removal, and allowing the enzyme to peel back by excising previously incorporated nucleotides, we have developed a computer simulation for the synthesis of a DNA molecule. In the special case, where insertion and removal frequencies are within the biologically interesting range for spontaneous mutations, the effect of “peelback” on error correction can be obtained analytically. Our results suggest that the magnitude of the removal frequency (3′-exonuclease activity) is the parameter that exerts the greatest influence on error correction capability; the frequency of errors is less sensitive to either the specificity for removal of mismatched relative to correctly matched bases, or to peelback.     

Conformational coupling in DNA polymerase information transfer.

The extraordinary fidelity of DNA replication during forward polymerization and exonuclease error correction is largely a function of a conformational change that occurs in response to a correct dNTP binding to properly base-paired duplex DNA. The conformational change serves as a kinetic barrier to effect the rapid incorporation of correct bases while minimizing the rate of polymerization with incorrect bases and allowing for selective removal of mismatches. However, in spite of the number of attractive features to the conformational change model, the evidence in support of such a rate-limiting step is still subject to significant uncertainty. It is the challenge of further work on DNA polymerases as well as many other enzyme systems to devise new methods to define the transient state of the enzyme during catalysis and to relate the kinetic and thermodynamic parameters to the enzyme structure.     

Substrate specificity and proposed structure of the proofreading complex of T7 DNA polymerase

Faithful replication of genomic DNA by high-fidelity DNA polymerases is crucial for the survival of most living organisms. While high-fidelity DNA polymerases favor canonical base pairs over mismatches by a factor of ∼1 × 10⁵, fidelity is further enhanced several orders of magnitude by a 3′–5′ proofreading exonuclease that selectively removes mispaired bases in the primer strand. Despite the importance of proofreading to maintaining genome stability, it remains much less studied than the fidelity mechanisms employed at the polymerase active site. Here we characterize the substrate specificity for the proofreading exonuclease of a high-fidelity DNA polymerase by investigating the proofreading kinetics on various DNA substrates. The contribution of the exonuclease to net fidelity is a function of the kinetic partitioning between extension and excision. We show that while proofreading of a terminal mismatch is efficient, proofreading a mismatch buried by one or two correct bases is even more efficient. Because the polymerase stalls after incorporation of a mismatch and after incorporation of one or two correct bases on top of a mismatch, the net contribution of the exonuclease is a function of multiple opportunities to correct mistakes. We also characterize the exonuclease stereospecificity using phosphorothioate-modified DNA, provide a homology model for the DNA primer strand in the exonuclease active site, and propose a dynamic structural model for the transfer of DNA from the polymerase to the exonuclease active site based on MD simulations.     

Function and assembly of the bacteriophage T4 DNA replication complex: interactions of the T4 polymerase with various model DNA constructs

Complexes formed between DNA polymerase and genomic DNA at the replication fork are key elements of the replication machinery. We used sedimentation velocity, fluorescence anisotropy, and surface plasmon resonance to measure the binding interactions between bacteriophage T4 DNA polymerase (gp43) and various model DNA constructs. These results provide quantitative insight into how this replication polymerase performs template-directed 5' --> 3' DNA synthesis and how this function is coordinated with the activities of the other proteins of the replication complex. We find that short (single- and double-stranded) DNA molecules bind a single gp43 polymerase in a nonspecific (overlap) binding mode with moderate affinity (Kd approximately 150 nm) and a binding site size of approximately 10 nucleotides for single-stranded DNA and approximately 13 bp for double-stranded DNA. In contrast, gp43 binds in a site-specific (nonoverlap) mode and significantly more tightly (Kd approximately 5 nm) to DNA constructs carrying a primer-template junction, with the polymerase covering approximately 5 nucleotides downstream and approximately 6-7 bp upstream of the 3'-primer terminus. The rate of this specific binding interaction is close to diffusion-controlled. The affinity of gp43 for the primer-template junction is modulated specifically by dNTP substrates, with the next "correct" dNTP strengthening the interaction and an incorrect dNTP weakening the observed binding. These results are discussed in terms of the individual steps of the polymerase-catalyzed single nucleotide addition cycle and the replication complex assembly process. We suggest that changes in the kinetics and thermodynamics of these steps by auxiliary replication proteins constitute a basic mechanism for protein coupling within the replication complex.     

Effect of oxidatively damaged DNA on the active site preorganization during nucleotide incorporation in a high fidelity polymerase from Bacillus stearothermophilus

We study the effect of the oxidative lesion 8-oxoguanine (8oxoG) on the preorganization of the active site for DNA replication in the closed (active) state of the Bacillus fragment (BF), a Klenow analog from Bacillus stearothermophilus. Our molecular dynamics and free energy simulations of explicitly solvated model ternary complexes of BF bound to correct dCTP/incorrect dATP opposite guanine (G) and 8oxoG bases in DNA suggest that the lesion introduces structural and energetic changes at the catalytic site to favor dATP insertion. Despite the formation of a stable Watson-Crick pairing in the 8oxoG:dCTP system, the catalytic geometry is severely distorted to possibly slow down catalysis. Indeed, our calculated free energy landscapes associated with active site preorganization suggest additional barriers to assemble an efficient catalytic site, which need to be overcome during dCTP incorporation opposite 8oxoG relative to that opposite undamaged G. In contrast, the catalytic geometry for the Hoogsteen pairing in the 8oxoG:dATP system is highly organized and poised for efficient nucleotide incorporation via the "two-metal-ion" catalyzed phosphoryl transfer mechanism. However, the free energy calculations suggest that the catalytic geometry during dATP incorporation opposite 8oxoG is considerably less plastic than that during dCTP incorporation opposite G despite a very similar, well organized catalytic site for both systems. A correlation analysis of the dynamics trajectories suggests the presence of significant coupling between motions of the polymerase fingers and the primary distance for nucleophilic attack (i.e., between the terminal primer O3' and the dNTP P(alpha.) atoms) during correct dCTP incorporation opposite undamaged G. This coupling is shown to be disrupted during nucleotide incorporation by the polymerase with oxidatively damaged DNA/dNTP substrates. We also suggest that the lesion affects DNA interactions with key polymerase residues, thereby affecting the enzymes ability to discriminate against non-complementary DNA/dNTP substrates. Taken together, our results provide a unified structural, energetic, and dynamic platform to rationalize experimentally observed relative nucleotide incorporation rates for correct dCTP/incorrect dATP insertion opposite an undamaged/oxidatively damaged template G by BF.     

Implications for damage recognition during Dpo4-mediated mutagenic bypass of m1G and m3C lesions

DNA is susceptible to alkylation damage by a number of environmental agents that modify the Watson-Crick edge of the bases. Such lesions, if not repaired, may be bypassed by Y-family DNA polymerases. The bypass polymerase Dpo4 is strongly inhibited by 1-methylguanine (m1G) and 3-methylcytosine (m3C), with nucleotide incorporation opposite these lesions being predominantly mutagenic. Further, extension after insertion of both correct and incorrect bases, introduces additional base substitution and deletion errors. Crystal structures of the Dpo4 ternary extension complexes with correct and mismatched 3'-terminal primer bases opposite the lesions reveal that both m1G and m3C remain positioned within the DNA template/primer helix. However, both correct and incorrect pairing partners exhibit pronounced primer terminal nucleotide distortion, being primarily evicted from the DNA helix when opposite m1G or misaligned when pairing with m3C. Our studies provide insights into mechanisms related to hindered and mutagenic bypass of methylated lesions and models associated with damage recognition by repair demethylases.     

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.     

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.     

Epigenetically modified N6-methyladenine inhibits DNA replication by human DNA polymerase η

N⁶-methyladenine (6mA), as a newly reported epigenetic marker, plays significant roles in regulation of various biological processes in eukaryotes. However, the effect of 6mA on human DNA replication remain elusive. In this work, we used Y-family human DNA polymerase η as a model to investigate the kinetics of bypass of 6mA by hPol η. We found 6mA and its intermediate hypoxanthine (I) on template partially inhibited DNA replication by hPol η. dTMP incorporation opposite 6mA and dCMP incorporation opposite I can be considered as correct incorporation. However, both 6mA and I reduced correct incorporation efficiency, next-base extension efficiency, and the priority in extension beyond correct base pair. Both dTMP incorporation opposite 6mA and dCTP opposite I showed fast burst phases. However, 6mA and I reduced the burst incorporation rates (k_pol) and increased the dissociation constant (K_d,dNTP), compared with that of dTMP incorporation opposite unmodified A. Biophysical binding assays revealed that both 6mA and I on template reduced the binding affinity of hPol η to DNA in binary or ternary complex compared with unmodified A. All the results explain the inhibition effects of 6mA and I on DNA replication by hPol η, providing new insight in the effects of epigenetically modified 6mA on human DNA replication.     

The determination of equilibrium constants for heterogeneous macromolecular interactions. Systems forming 2:1 complexes

In developing a method for analyzing the heterogeneous association nA + mB in equilibrium AnBm, we have specifically investigated the case of n = 2, m = 1 for both the specific case of no appreciable intermediates and the more general case allowing intermediates. Computer-simulated three-dimensional surfaces of the 2:1 model generated from total concentrations of species A and B and the resulting weight-average molecular weights were analyzed with a Gauss-Newton nonlinear least-squares minimization routine. The surfaces generated included normalized random error of varying standard deviations imposed upon both the concentrations and weight-average molecular weights. For comparison purposes, these surfaces were analyzed not only by using the correct 2:1 model, but also by an incorrect (1:1) model and by the other (incorrect) 2:1 model. Except for those situations where the 'experimental' noise was consistently higher than the concentration of one of the species, correct K values were obtained and the correct model was easily distinguished from the incorrect model. The computer routine similarly distinguished between data correctly described as 1:1 and the same data incorrectly analyzed as either 2:1 model. For those cases in which a microscopic Ki value predicts an association such that all species involved for that particular Ki are in appreciable amounts, the Ki value is returned correctly. Correct overall equilibrium constants are also converged upon as long as adequate amounts of A2B, B and A are present.     

An induced-fit kinetic mechanism for DNA replication fidelity: direct measurement by single-turnover kinetics.

An exonuclease-deficient mutant of T7 DNA polymerase was constructed and utilized in a series of kinetic studies on misincorporation and next correct dNTP incorporation. By using a synthetic oligonucleotide template-primer system for which the kinetic pathway for correct incorporation has been solved [Patel, S.S., Wong, I., & Johnson, K. A. (1991) Biochemistry (first of three papers in this issue)], the kinetic parameters for the incorporation of the incorrect triphosphates dATP, dCTP, and dGTP were determined, giving, respectively, kcat/Km values of 91, 23, and 4.3 M-1 s-1 and a discrimination in the polymerization step of 10(5)-10(6). The rates of misincorporation in all cases were linearly dependent on substrate concentration up to 4 mM, beyond which severe inhibition was observed. Competition of correct incorporation versus dCTP revealed an estimated Ki of approximately 6-8 mM, suggesting a corresponding kcat of 0.14s-1. Moderate elemental effects of 19-, 17-, and 34-fold reduction in rates were measured by substituting the alpha-thiotriphosphate analogues for dATP, dCTP, and dGTP, respectively, indicating that the chemistry step is partially rate-limiting. The absence of a burst of incorporation during the first turnover places the rate-limiting step at a triphosphate binding induced conformational change before chemistry. In contrast, the incorporation of the next correct triphosphate, dCTP, on a mismatched DNA substrate was saturable with a Km of 87 microM for dCTP, 4-fold higher than the Kd for the correct incorporation on duplex DNA, and a kcat of 0.025 s-1. A larger elemental effect of 60, however, suggests a rate-limiting chemistry step. The rate of pyrophosphorolysis on a mismatched 3'-end is undetectable, indicating that pyrophosphorolysis does not play a proofreading role in replication. These results show convincingly that the T7 DNA polymerase discriminates against the incorrect triphosphate by an induced-fit conformational change and that, following misincorporation, the enzyme then selects against the resultant mismatched DNA by a slow, rate-limiting chemistry step, thereby allowing sufficient time for the release of the mismatched DNA from the polymerase active site to be followed by exonucleolytic error correction.     

A mutant of DNA polymerase I (Klenow fragment) with reduced fidelity

The kinetic parameters governing incorporation of correct and incorrect bases into synthetic DNA duplexes have been investigated for Escherichia coli DNA polymerase I [Klenow fragment (KF)] and for two mutants, Tyr766Ser and Tyr766Phe. Tyr766 is located at the C-terminus of helix O in the DNA-binding cleft of KF. The catalytic efficiency for correct incorporation of dNTP is reduced 5-fold for Tyr766Ser. The catalytic efficiencies of all 12 possible misincorporations have been determined for both KF and Tyr766Ser by using single-turnover kinetic conditions and a form of the enzyme that is devoid of the 3'-5' exonuclease activity because of other single amino acid replacements. Tyr766Ser displays an increased efficiency of misincorporation (a reduction in fidelity) for several of the 12 mismatches. The largest increase in efficiency of misincorporation for Tyr766Ser occurs for the misincorporation of TMP opposite template guanosine, a 44-fold increase. In contrast, the efficiencies of misincorporation of dAMP opposite template A, G, or C are little affected by the mutation. A determination of the kinetic parameters associated with a complete kinetic scheme has been made for Tyr766Ser. The rate of addition of the next correct nucleotide onto a preexisting mismatch is decreased for Tyr766Ser. The fidelity of Tyr766Phe was not substantially different from that of KF for the misincorporations examined, indicating that it is the loss of the phenolic ring of the side chain of Tyr766 that leads to the significant decrease in fidelity. The results indicate that KF actively participates in the reduction of misincorporations during the polymerization event and that Tyr766 plays an important role in maintaining the high fidelity of replication by KF.     

Human DNA polymerase iota utilizes different nucleotide incorporation mechanisms dependent upon the template base

Human DNA polymerase ι (Polι) is a member of the Y family of DNA polymerases involved in translesion DNA synthesis. Polι is highly unusual in that it possesses a high fidelity on template A, but has an unprecedented low fidelity on template T, preferring to misincorporate a G instead of an A. To understand the mechanisms of nucleotide incorporation opposite different template bases by Polι, we have carried out pre-steady-state kinetic analyses of nucleotide incorporation opposite templates A and T. These analyses have revealed that opposite template A, the correct nucleotide is preferred because it is bound tighter and is incorporated faster than the incorrect nucleotides. Opposite template T, however, the correct and incorrect nucleotides are incorporated at very similar rates, and interestingly, the greater efficiency of G misincorporation relative to A incorporation opposite T arises predominantly from the tighter binding of G. Based on these results, we propose that the incipient base pair is accommodated differently in the active site of Polι dependent upon the template base and that when T is the templating base, Polι accommodates the wobble base pair better than the Watson-Crick base pair.     

Conformational dynamics during misincorporation and mismatch extension defined using a DNA polymerase with a fluorescent artificial amino acid

High-fidelity DNA polymerases select the correct nucleotide over the structurally similar incorrect nucleotides with extremely high specificity while maintaining fast rates of incorporation. Previous analysis revealed the conformational dynamics and complete kinetic pathway governing correct nucleotide incorporation using a high-fidelity DNA polymerase variant containing a fluorescent unnatural amino acid. Here we extend this analysis to investigate the kinetics of nucleotide misincorporation and mismatch extension. We report the specificity constants for all possible misincorporations and characterize the conformational dynamics of the enzyme during misincorporation and mismatch extension. We present free energy profiles based on the kinetic measurements and discuss the effect of different steps on specificity. During mismatch incorporation and subsequent extension with the correct nucleotide, the rates of the conformational change and chemistry are both greatly reduced. The nucleotide dissociation rate, however, increases to exceed the rate of chemistry. To investigate the structural basis for discrimination against mismatched nucleotides, we performed all atom molecular dynamics simulations on complexes with either the correct or mismatched nucleotide bound at the polymerase active site. The simulations suggest that the closed form of the enzyme with a mismatch bound is greatly destabilized due to weaker interactions with active site residues, nonideal base pairing, and a large increase in the distance from the 3ʹ-OH group of the primer strand to the α-phosphate of the incoming nucleotide, explaining the reduced rates of misincorporation. The observed kinetic and structural mechanisms governing nucleotide misincorporation reveal the general principles likely applicable to other high-fidelity DNA polymerases.     

Mechanism of efficient and accurate nucleotide incorporation opposite 7,8-dihydro-8-oxoguanine by Saccharomyces cerevisiae DNA polymerase eta

Most DNA polymerases incorporate nucleotides opposite template 7,8-dihydro-8-oxoguanine (8-oxoG) lesions with reduced efficiency and accuracy. DNA polymerase (Pol) eta, which catalyzes the error-free replication of template thymine-thymine (TT) dimers, has the unique ability to accurately and efficiently incorporate nucleotides opposite 8-oxoG templates. Here we have used pre-steady-state kinetics to examine the mechanisms of correct and incorrect nucleotide incorporation opposite G and 8-oxoG by Saccharomyces cerevisiae Pol eta. We found that Pol eta binds the incoming correct dCTP opposite both G and 8-oxoG with similar affinities, and it incorporates the correct nucleotide bound opposite both G and 8-oxoG with similar rates. While Pol eta incorporates an incorrect A opposite 8-oxoG with lower efficiency than it incorporates a correct C, it does incorporate A more efficiently opposite 8-oxoG than opposite G. This is mainly due to greater binding affinity for the incorrect incoming dATP opposite 8-oxoG. Overall, these results show that Pol eta replicates through 8-oxoG without any barriers introduced by the presence of the lesion.     

Fidelity of DNA ligation: a novel experimental approach based on the polymerisation of libraries of oligonucleotides.

Complete libraries of oligonucleotides were used as substrates for Thermus thermophilus DNA ligase, on a M13mp18 ssDNA template. A 17mer primer was used to start a polymerisation process. Ladders of ligation products were analysed by gel electrophoresis. Octa-, nona- and decanucleotide libraries were compared. Nonanucleotides were optimum for polymerisation and up to 15 monomers were ligated. The fidelity of incorporation was studied by sequencing 28 clones (2268 bases) of nonanucleotide polymers, 12 monomers in length. Of the ligated monomers, 79% were the correct complementary sequence. In a total of 57 (2.5%) mispaired bases, there was a strong bias to G.T, G.A, G.G and A.G mismatches. Of the mismatches, 86% were found to be purines on the incoming oligonucleotide, of which 71% were G. There is evidence for clustering of mismatches within specific 9mers and at specific positions within these 9mers. The most frequent mismatches were at the 5'-terminus of the oligonucleotide, followed by the central position. We suggest that sequence selection was imposed by the ligase and not just by base pairing interactions. The ligase directs polymerisation in the 3' to 5' direction which we propose is linked to its role in lagging strand DNA replication.