Hydrogen-bonding capability of a templating difluorotoluene nucleotide residue in an RB69 DNA polymerase ternary complex

RB69 DNA polymerase structure, kinetics, and fidelity

This review will summarize our structural and kinetic studies of RB69 DNA polymerase (RB69pol) as well as selected variants of the wild-type enzyme that were undertaken to obtain a deeper understanding of the exquisitely high fidelity of B family replicative DNA polymerases. We discuss how the structures of the various RB69pol ternary complexes can be used to rationalize the results obtained from pre-steady-state kinetic assays. Our main findings can be summarized as follows. (i) Interbase hydrogen bond interactions can increase catalytic efficiency by 5000-fold; meanwhile, base selectivity is not solely determined by the number of hydrogen bonds between the incoming dNTP and the templating base. (ii) Minor-groove hydrogen bond interactions at positions n – 1 and n – 2 of the primer strand and position n – 1 of the template strand in RB69pol ternary complexes are essential for efficient primer extension and base selectivity. (iii) Partial charge interactions among the incoming dNTP, the penultimate base pair, and the hydration shell surrounding the incoming dNTP modulate nucleotide insertion efficiency and base selectivity. (iv) Steric clashes between mismatched incoming dNTPs and templating bases with amino acid side chains in the nascent base pair binding pocket (NBP) as well as weak interactions and large gaps between the incoming dNTPs and the templating base are some of the reasons that incorrect dNTPs are incorporated so inefficiently by wild-type RB69pol. In addition, we developed a tC°–tC_nitro Förster resonance energy transfer assay to monitor partitioning of the primer terminus between the polymerase and exonuclease subdomains.

Alteration in the cavity size adjacent to the active site of RB69 DNA polymerase changes its conformational dynamics

Internal cavities are a common feature of many proteins, often having profound effects on the dynamics of their interactions with substrate and binding partners. RB69 DNA polymerase (pol) has a hydrophobic cavity right below the nucleotide binding pocket at the tip of highly conserved L415 side chain. Replacement of this residue with Gly or Met in other B family pols resulted in higher mutation rates. When similar substitutions for L415 were introduced into RB69pol, only L415A and L415G had dramatic effects on pre-steady-state kinetic parameters, reducing base selectivity by several hundred fold. On the other hand, the L415M variant behaved like the wild-type. Using a novel tC^o-tC_nitro Förster Resonance Energy Transfer (FRET) assay, we were able to show that the partition of the primer terminus between pol and exonuclease (exo) domains was compromised with the L415A and L415G mutants, but not with the L415M variant. These results could be rationalized by changes in their structures as determined by high resolution X-ray crystallography.

Partial base flipping is sufficient for strand slippage near DNA duplex termini

Strand slippage is a structural mechanism by which insertion-deletion (indel) mutations are introduced during replication by polymerases. Three-dimensional atomic-resolution structural pathways are still not known for the decades-old template slippage description. The dynamic nature of the process and the higher energy intermediates involved increase the difficulty of studying these processes experimentally. In the present study, restrained and unrestrained molecular dynamics simulations, carried out using multiple nucleic acid force fields, are used to demonstrate that partial base-flipping can be sufficient for strand slippage at DNA duplex termini. Such strand slippage can occur in either strand, i.e. near either the 3' or the 5' terminus of a DNA strand, which suggests that similar structural flipping mechanisms can cause both primer and template slippage. In the repetitive mutation hot-spot sequence studied, non-canonical base-pairing with exposed DNA groove atoms of a neighboring G:C base-pair stabilizes a partially flipped state of the cytosine. For its base-pair partner guanine, a similar partially flipped metastable intermediate was not detected, and the propensity for sustained slippage was also found to be lower. This illustrates that a relatively small metastable DNA structural distortion in polymerase active sites could allow single base insertion or deletion mutations to occur, and stringent DNA groove molecular recognition may be required to maintain intrinsic DNA polymerase fidelity. The implications of a close relationship between base-pair dissociation, base unstacking, and strand slippage are discussed in the context of sequence dependence of indel mutations.

DNA mismatch synthesis complexes provide insights into base selectivity of a B family DNA polymerase

Current hypotheses that attempt to rationalize the high degree of base selectivity exhibited by replicative DNA polymerases (pols) concur that ternary complexes formed with incorrect dNTPs are destabilized. Knowing what accounts for this destabilization is likely to be the key to understanding base discrimination. To address this issue, we have determined crystal structures of ternary complexes with all 12 mismatches using an engineered RB69 pol quadruple mutant (qm, L415A/L561A/S565G/Y567A) that enabled us to capture nascent mispaired dNTPs. These structures show that mismatches in the nascent base-pair binding pocket (NBP) of the qm pol differ markedly from mismatches embedded in binary pol-DNA complexes. Surprisingly, only 3 of 12 mismatches clash with the NBP when they are modeled into the wild-type (wt) pol. The remaining can fit into a wt pol ternary complex but deviate from normal Watson-Crick base-pairs. Repositioning of the templating nucleotide residue and the enlarged NBP in qm ternary complex play important roles in accommodating incorrect incoming dNTPs. From these structures, we propose additional reasons as to why incorrect dNTPs are incorporated so inefficiently by wt RB69 pol: (i) steric clashes with side chains in the NBP after Fingers closing; (ii) weak interactions or large gaps between the incoming dNTP and the templating base; and (iii) burying a protonated base in the hydrophobic environment of the NBP. All of these possibilities would be expected to destabilize the closed ternary complex so that incorporation of incorrect dNTP would be a rare event.

Expanding the horizon of the thymine isostere biochemistry: unique cyclobutane dimers formed by photoreaction between a thymine and a toluene residue in the dinucleotide framework

Substituted toluenyl groups are considered as close isosteres of the thymine residue. They can be recognized by DNA polymerases as if they were thymine. These toluene derivatives are generally inert toward radical additions, including the [2+2] photo-cycloadditions, due to the stable structure of the aromatic ring and are usually used as solvents for radical reactions. Surprisingly, after incorporating toluene into the dinucleotide framework, we found that the UV excited thymine residue readily dimerizes with the toluenyl moiety through a [2+2] photo-addition reaction. Furthermore, the reaction site on the toluenyl moiety is not the C5=C6 bond, as commonly observed in cyclobutane pyrimidine dimers, but the C4=C5 or C3=C4 instead. Such a reaction pattern suggests that in the stacked structure, it is one of these bonds, not the C5=C6, that is close to the thymine C5=C6 bond. A similar structural feature is found in DNA duplex with a thymine replaced by a 2,4-difluorotoluene. Our results argue that although the substituted toluenyl moieties closely mimic the size and shape of the thymine residue, their more hydrophobic nature determines that they stack on DNA bases differently from the natural thymine residue and likely cause local conformational changes in duplex DNA.

Synthetic nucleotides as probes of DNA polymerase specificity

The genetic code is continuously expanding with new nucleobases designed to suit specific research needs. These synthetic nucleotides are used to study DNA polymerase dynamics and specificity and may even inhibit DNA polymerase activity. The availability of an increasing chemical diversity of nucleotides allows questions of utilization by different DNA polymerases to be addressed. Much of the work in this area deals with the A family DNA polymerases, for example, Escherichia coli DNA polymerase I, which are DNA polymerases involved in replication and whose fidelity is relatively high, but more recent work includes other families of polymerases, including the Y family, whose members are known to be error prone. This paper focuses on the ability of DNA polymerases to utilize nonnatural nucleotides in DNA templates or as the incoming nucleoside triphosphates. Beyond the utility of nonnatural nucleotides as probes of DNA polymerase specificity, such entities can also provide insight into the functions of DNA polymerases when encountering DNA that is damaged by natural agents. Thus, synthetic nucleotides provide insight into how polymerases deal with nonnatural nucleotides as well as into the mutagenic potential of nonnatural nucleotides.

Contribution of partial charge interactions and base stacking to the efficiency of primer extension at and beyond abasic sites in DNA

During DNA synthesis, base stacking and Watson-Crick (WC) hydrogen bonding increase the stability of nascent base pairs when they are in a ternary complex. To evaluate the contribution of base stacking to the incorporation efficiency of dNTPs when a DNA polymerase encounters an abasic site, we varied the penultimate base pairs (PBs) adjacent to the abasic site using all 16 possible combinations. We then determined pre-steady-state kinetic parameters with an RB69 DNA polymerase variant and solved nine structures of the corresponding ternary complexes. The efficiency of incorporation for incoming dNTPs opposite an abasic site varied between 2- and 210-fold depending on the identity of the PB. We propose that the A rule can be extended to encompass the fact that DNA polymerase can bypass dA/abasic sites more efficiently than other dN/abasic sites. Crystal structures of the ternary complexes show that the surface of the incoming base was stacked against the PB's interface and that the kinetic parameters for dNMP incorporation were consistent with specific features of base stacking, such as surface area and partial charge-charge interactions between the incoming base and the PB. Without a templating nucleotide residue, an incoming dNTP has no base with which it can hydrogen bond and cannot be desolvated, so that these surrounding water molecules become ordered and remain on the PB's surface in the ternary complex. When these water molecules are on top of a hydrophobic patch on the PB, they destabilize the ternary complex, and the incorporation efficiency of incoming dNTPs is reduced.

Using a fluorescent cytosine analogue tC(o) to probe the effect of the Y567 to Ala substitution on the preinsertion steps of dNMP incorporation by RB69 DNA polymerase

Residues in the nascent base pair binding pocket (NBP) of bacteriophage RB69 DNA polymerase (RB69pol) are responsible for base discrimination. Replacing Tyr567 with Ala leads to greater flexibility in the NBP, increasing the probability of misincorporation. We used the fluorescent cytosine analogue, 1,3-diaza-2-oxophenoxazine (tC(o)), to identify preinsertion step(s) altered by NBP flexibility. When tC(o) is the templating base in a wild-type (wt) RB69pol ternary complex, its fluorescence is quenched only in the presence of dGTP. However, with the RB69pol Y567A mutant, the fluorescence of tC(o) is also quenched in the presence of dATP. We determined the crystal structure of the dATP/tC(o)-containing ternary complex of the RB69pol Y567A mutant at 1.9 Å resolution and found that the incoming dATP formed two hydrogen bonds with an imino-tautomerized form of tC(o). Stabilization of the dATP/tC(o) base pair involved movement of the tC(o) backbone sugar into the DNA minor groove and required tilting of the tC(o) tricyclic ring to prevent a steric clash with L561. This structure, together with the pre-steady-state kinetic parameters and dNTP binding affinity, estimated from equilibrium fluorescence titrations, suggested that the flexibility of the NBP, provided by the Y567 to Ala substitution, led to a more favorable forward isomerization step resulting in an increase in dNTP binding affinity.

Probing minor groove hydrogen bonding interactions between RB69 DNA polymerase and DNA

Minor groove hydrogen bonding (HB) interactions between DNA polymerases (pols) and N3 of purines or O2 of pyrimidines have been proposed to be essential for DNA synthesis from results obtained using various nucleoside analogues lacking the N3 or O2 contacts that interfered with primer extension. Because there has been no direct structural evidence to support this proposal, we decided to evaluate the contribution of minor groove HB interactions with family B pols. We have used RB69 DNA pol and 3-deaza-2'-deoxyadenosine (3DA), an analogue of 2-deoxyadenosine, which has the same HB pattern opposite T but with N3 replaced with a carbon atom. We then determined pre-steady-state kinetic parameters for the insertion of dAMP opposite dT using primer/templates (P/T)-containing 3DA. We also determined three structures of ternary complexes with 3DA at various positions in the duplex DNA substrate. We found that the incorporation efficiency of dAMP opposite dT decreased 10(2)-10(3)-fold even when only one minor groove HB interaction was missing. Our structures show that the HB pattern and base pair geometry of 3DA/dT is exactly the same as those of dA/dT, which makes 3DA an optimal analogue for probing minor groove HB interactions between a DNA polymerase and a nucleobase. In addition, our structures provide a rationale for the observed 10(2)-10(3)-fold decrease in the rate of nucleotide incorporation. The minor groove HB interactions between position n - 2 of the primer strand and RB69pol fix the rotomer conformations of the K706 and D621 side chains, as well as the position of metal ion A and its coordinating ligands, so that they are in the optinal orientation for DNA synthesis.

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.

Measurement and theory of hydrogen bonding contribution to isosteric DNA base pairs

We address the recent debate surrounding the ability of 2,4-difluorotoluene (F), a low-polarity mimic of thymine (T), to form a hydrogen-bonded complex with adenine in DNA. The hydrogen bonding ability of F has been characterized as small to zero in various experimental studies, and moderate to small in computational studies. However, recent X-ray crystallographic studies of difluorotoluene in DNA/RNA have indicated, based on interatomic distances, possible hydrogen bonding interactions between F and natural bases in nucleic acid duplexes and in a DNA polymerase active site. Since F is widely used to measure electrostatic contributions to pairing and replication, it is important to quantify the impact of this isostere on DNA stability. Here, we studied the pairing stability and selectivity of this compound and a closely related variant, dichlorotoluene deoxyriboside (L), in DNA, using both experimental and computational approaches. We measured the thermodynamics of duplex formation in three sequence contexts and with all possible pairing partners by thermal melting studies using the van't Hoff approach, and for selected cases by isothermal titration calorimetry (ITC). Experimental results showed that internal F-A pairing in DNA is destabilizing by 3.8 kcal/mol (van't Hoff, 37 °C) as compared with T-A pairing. At the end of a duplex, base-base interactions are considerably smaller; however, the net F-A interaction remains repulsive while T-A pairing is attractive. As for selectivity, F is found to be slightly selective for adenine over C, G, T by 0.5 kcal mol, as compared with thymine's selectivity of 2.4 kcal/mol. Interestingly, dichlorotoluene in DNA is slightly less destabilizing and slightly more selective than F, despite the lack of strongly electronegative fluorine atoms. Experimental data were complemented by computational results, evaluated at the M06-2X/6-31+G(d) and MP2/cc-pVTZ levels of theory. These computations suggest that the pairing energy of F to A is ~28% of that of T-A, and most of this interaction does not arise from the F···HN interaction, but rather from the CH···N interaction. The nucleobase analogue shows no inherent selectivity for adenine over other bases, and L-A pairing energies are slightly weaker than for F-A. Overall, the results are consistent with a small favorable noncovalent interaction of F with A offset by a large desolvation cost for the polar partner. We discuss the findings in light of recent structural studies and of DNA replication experiments involving these analogues.

Structural insights into complete metal ion coordination from ternary complexes of B family RB69 DNA polymerase

We have captured a preinsertion ternary complex of RB69 DNA polymerase (RB69pol) containing the 3' hydroxyl group at the terminus of an extendable primer (ptO3') and a nonhydrolyzable 2'-deoxyuridine 5'-α,β-substituted triphosphate, dUpXpp, where X is either NH or CH(2), opposite a complementary templating dA nucleotide residue. Here we report four structures of these complexes formed by three different RB69pol variants with catalytically inert Ca(2+) and four other structures with catalytically competent Mn(2+) or Mg(2+). These structures provide new insights into why the complete divalent metal-ion coordination complexes at the A and B sites are required for nucleotidyl transfer. They show that the metal ion in the A site brings ptO3' close to the α-phosphorus atom (Pα) of the incoming dNTP to enable phosphodiester bond formation through simultaneous coordination of both ptO3' and the nonbridging Sp oxygen of the dNTP's α-phosphate. The coordination bond length of metal ion A as well as its ionic radius determines how close ptO3' can approach Pα. These variables are expected to affect the rate of bond formation. The metal ion in the B site brings the pyrophosphate product close enough to Pα to enable pyrophosphorolysis and assist in the departure of the pyrophosphate. In these dUpXpp-containing complexes, ptO3' occupies the vertex of a distorted metal ion A coordination octahedron. When ptO3' is placed at the vertex of an undistorted, idealized metal ion A octahedron, it is within bond formation distance to Pα. This geometric relationship appears to be conserved among DNA polymerases of known structure.