Slow rate of phosphodiester bond formation accounts for the strong bias that Taq DNA polymerase shows against 2',3'-dideoxynucleotide terminators

Beyond Chelation: EDTA Tightly Binds Taq DNA Polymerase, MutT and dUTPase and Directly Inhibits dNTPase Activity

EDTA is commonly used as an efficient chelator of metal ion enzyme cofactors. It is highly soluble, optically inactive and does not interfere with most chemicals used in standard buffers making EDTA a common choice to generate metal-free conditions for biochemical and biophysical investigations. However, the controversy in the literature on metal-free enzyme activities achieved using EDTA or by other means called our attention to a putative effect of EDTA beyond chelation. Here, we show that EDTA competes for the nucleotide binding site of the nucleotide hydrolase dUTPase by developing an interaction network within the active site similar to that of the substrate. To achieve these findings, we applied kinetics and molecular docking techniques using two different dUTPases. Furthermore, we directly measured the binding of EDTA to dUTPases and to two other dNTPases, the Taq polymerase and MutT using isothermal titration calorimetry. EDTA binding proved to be exothermic and mainly enthalpy driven with a submicromolar dissociation constant considerably lower than that of the enzyme:substrate or the Mg:EDTA complexes. Control proteins, including an ATPase, did not interact with EDTA. Our findings indicate that EDTA may act as a selective inhibitor against dNTP hydrolyzing enzymes and urge the rethinking of the utilization of EDTA in enzymatic experiments.

Sequence motifs capable of forming DNA stem-loop structures act as a replication diode

Calculating peak‐height ratios between single‐nucleotide polymorphisms (SNP) alleles in sequencing chromatograms is a practical method for estimating their copy number proportions (CNPs). However, it is surprising that sequencing DNA from different directions might yield different results. We analyzed three adjacent SNPs within the ovine period circadian‐clock 2 (PER2) gene that displayed such behavior. We compared Sanger and DNA‐seq sequencing for this locus and applied high‐resolution melt and MFOLD analyses to point to the DNA secondary structure that underlined this phenomenon. A synthetic system of oligonucleotides cloned into plasmids was used to further test the effect of such structures on sequencing. Our analyses indicated that a stem–loop structure capable of G–T pairing mediated the orientation bias by stabilizing this structure for specific alleles in heterozygous situations. We propose that this wobble‐like pairing hinders DNA polymerase passage on one strand while, on the complementary strand, the nonpaired A–C nucleotide counterparts allow unobstructed replication. Experimentation with synthetic amplicons that form similar stem–loop structures supported our hypothesis. We coined the term ‘replication diode’ for this effect and demonstrated that we can minimize it by lowering DNA and salt concentration. We also demonstrated that common genomic palindromes might induce the replication diode effect by applying bidirectional sequencing to an amplicon containing the palindrome within the human miRNA 1‐1 gene. Hence, to obtain reliable peak‐height ratios, bidirectional sequencing should be practiced at the lowest possible ionic strength whenever estimating CNPs. Further research is needed to determine whether the observed variable stem–loop structures affect PER2 regulation in vivo.

Pre-steady-state Kinetic Analysis of a Family D DNA Polymerase from Thermococcus sp. 9°N Reveals Mechanisms for Archaeal Genomic Replication and Maintenance

Family D DNA polymerases (polDs) have been implicated as the major replicative polymerase in archaea, excluding the Crenarchaeota branch, and bear little sequence homology to other DNA polymerase families. Here we report a detailed kinetic analysis of nucleotide incorporation and exonuclease activity for a Family D DNA polymerase from Thermococcus sp. 9°N. Pre-steady-state single-turnover nucleotide incorporation assays were performed to obtain the kinetic parameters, kpol and Kd, for correct nucleotide incorporation, incorrect nucleotide incorporation, and ribonucleotide incorporation by exonuclease-deficient polD. Correct nucleotide incorporation kinetics revealed a relatively slow maximal rate of polymerization (kpol ∼ 2.5 s(-1)) and especially tight nucleotide binding (Kd (dNTP) ∼ 1.7 μm), compared with DNA polymerases from Families A, B, C, X, and Y. Furthermore, pre-steady-state nucleotide incorporation assays revealed that polD prevents the incorporation of incorrect nucleotides and ribonucleotides primarily through reduced nucleotide binding affinity. Pre-steady-state single-turnover assays on wild-type 9°N polD were used to examine 3'-5' exonuclease hydrolysis activity in the presence of Mg(2+) and Mn(2+). Interestingly, substituting Mn(2+) for Mg(2+) accelerated hydrolysis rates > 40-fold (kexo ≥ 110 s(-1) versus ≥ 2.5 s(-1)). Preference for Mn(2+) over Mg(2+) in exonuclease hydrolysis activity is a property unique to the polD family. The kinetic assays performed in this work provide critical insight into the mechanisms that polD employs to accurately and efficiently replicate the archaeal genome. Furthermore, despite the unique properties of polD, this work suggests that a conserved polymerase kinetic pathway is present in all known DNA polymerase families.

DNA polymerases drive DNA sequencing-by-synthesis technologies: both past and present

Next-generation sequencing (NGS) technologies have revolutionized modern biological and biomedical research. The engines responsible for this innovation are DNA polymerases; they catalyze the biochemical reaction for deriving template sequence information. In fact, DNA polymerase has been a cornerstone of DNA sequencing from the very beginning. Escherichia coli DNA polymerase I proteolytic (Klenow) fragment was originally utilized in Sanger’s dideoxy chain-terminating DNA sequencing chemistry. From these humble beginnings followed an explosion of organism-specific, genome sequence information accessible via public database. Family A/B DNA polymerases from mesophilic/thermophilic bacteria/archaea were modified and tested in today’s standard capillary electrophoresis (CE) and NGS sequencing platforms. These enzymes were selected for their efficient incorporation of bulky dye-terminator and reversible dye-terminator nucleotides respectively. Third generation, real-time single molecule sequencing platform requires slightly different enzyme properties. Enterobacterial phage ϕ29 DNA polymerase copies long stretches of DNA and possesses a unique capability to efficiently incorporate terminal phosphate-labeled nucleoside polyphosphates. Furthermore, ϕ29 enzyme has also been utilized in emerging DNA sequencing technologies including nanopore-, and protein-transistor-based sequencing. DNA polymerase is, and will continue to be, a crucial component of sequencing technologies.

Conformational dynamics of Thermus aquaticus DNA polymerase I during catalysis

Despite the fact that DNA polymerases have been investigated for many years and are commonly used as tools in a number of molecular biology assays, many details of the kinetic mechanism they use to catalyze DNA synthesis remain unclear. Structural and kinetic studies have characterized a rapid, pre-catalytic open-to-close conformational change of the Finger domain during nucleotide binding for many DNA polymerases including Thermus aquaticus DNA polymerase I (Taq Pol), a thermostable enzyme commonly used for DNA amplification in PCR. However, little has been performed to characterize the motions of other structural domains of Taq Pol or any other DNA polymerase during catalysis. Here, we used stopped-flow Förster resonance energy transfer to investigate the conformational dynamics of all five structural domains of the full-length Taq Pol relative to the DNA substrate during nucleotide binding and incorporation. Our study provides evidence for a rapid conformational change step induced by dNTP binding and a subsequent global conformational transition involving all domains of Taq Pol during catalysis. Additionally, our study shows that the rate of the global transition was greatly increased with the truncated form of Taq Pol lacking the N-terminal domain. Finally, we utilized a mutant of Taq Pol containing a de novo disulfide bond to demonstrate that limiting protein conformational flexibility greatly reduced the polymerization activity of Taq Pol.

Enthalpic switch-points and temperature dependencies of DNA binding and nucleotide incorporation by Pol I DNA polymerases

This study examines the relationship between the DNA binding thermodynamics and the enzymatic activity of the Klenow and Klentaq Pol I DNA polymerases from Escherichia coli and Thermus aquaticus. Both polymerases bind DNA with nanomolar affinity at temperatures down to at least 5°C, but have lower than 1% enzymatic activity at these lower temperatures. For both polymerases it is found that the temperature of onset of significant enzymatic activity corresponds with the temperature where the enthalpy of binding (ΔHbinding) crosses zero (TH) and becomes favorable (negative). This TH/activity upshift temperature is 15°C for Klenow and 30°C for Klentaq. The results indicate that a negative free energy of DNA binding alone is not sufficient to proceed to catalysis, but that the enthalpic versus entropic balance of binding may be a modulator of the temperature dependence of enzymatic function. Analysis of the temperature dependence of the catalytic activity of Klentaq polymerase using expanded Eyring theory yields thermodynamic patterns for ΔG(‡), ΔH(‡), and TΔS(‡) that are highly analogous to those commonly observed for direct DNA binding. Eyring analysis also finds a significant ΔCp(‡) of formation of the activated complex, which in turn indicates that the temperature of maximal activity, after which incorporation rate slows with increasing temperature, will correspond with the temperature where the activation enthalpy (ΔH(‡)) switches from positive to negative.

Rapid incorporation kinetics and improved fidelity of a novel class of 3'-OH unblocked reversible terminators

Recent developments of unique nucleotide probes have expanded our understanding of DNA polymerase function, providing many benefits to techniques involving next-generation sequencing (NGS) technologies. The cyclic reversible termination (CRT) method depends on efficient base-selective incorporation of reversible terminators by DNA polymerases. Most terminators are designed with 3′-O-blocking groups but are incorporated with low efficiency and fidelity. We have developed a novel class of 3′-OH unblocked nucleotides, called Lightning Terminators™, which have a terminating 2-nitrobenzyl moiety attached to hydroxymethylated nucleobases. A key structural feature of this photocleavable group displays a ‘molecular tuning’ effect with respect to single-base termination and improved nucleotide fidelity. Using Therminator™ DNA polymerase, we demonstrate that these 3′-OH unblocked terminators exhibit superior enzymatic performance compared to two other reversible terminators, 3′-O-amino-TTP and 3′-O-azidomethyl-TTP. Lightning Terminators™ show maximum incorporation rates (k_pol) that range from 35 to 45 nt/s, comparable to the fastest NGS chemistries, yet with catalytic efficiencies (k_pol/K_D) comparable to natural nucleotides. Pre-steady-state kinetic studies of thymidine analogs revealed that the major determinant for improved nucleotide selectivity is a significant reduction in k_pol by >1000-fold over TTP misincorporation. These studies highlight the importance of structure–function relationships of modified nucleotides in dictating polymerase performance.

The mechanistic architecture of thermostable Pyrococcus furiosus family B DNA polymerase motif A and its interaction with the dNTP substrate

Thermostable DNA polymerases isolated from archaeal organisms have not been completely characterized kinetically and require further study if we are to understand both their dNTP binding mechanism and their role within the organism. Here, we demonstrate that the thermostable family B DNA polymerase from Pyrococcus furiosus (Pfu Pol) contains sensitive determinants of both dNTP binding and replicational fidelity within the highly conserved motif A. Site-directed mutagenesis of the motif A SYLP region revealed that small shifts in side chain volume result in significant changes in the dNTP binding affinity, steady state kinetics, and fidelity of the enzyme. Mutants of Y410 show high fidelity in both misincorporation assays and forward mutation assays, but display a substantially higher K(m) than the wild type. In contrast, mutations of upstream residue L409 result in a drastically reduced affinity for the correct dNTP, a much higher efficiency of both misincorporation and mismatch extension, and substantially lower fidelity as demonstrated by a PCR-based forward mutation assay. The A408S mutant, however, displayed a significant increase in both dNTP binding affinity and fidelity. In summary, these data show that modulation of motif A can greatly shift both the steady and pre-steady state kinetics of the enzyme as well as the fidelity of Pfu Pol.

Gln151 of HIV-1 Reverse Transcriptase Acts as a Steric Gate Towards Clinically Relevant Acyclic Phosphonate Nucleotide Analogues

Background

In the treatment of HIV, the loose active site of the HIV-1 reverse transcriptase (RT) allows numerous nucleotide analogues to act as proviral DNA ‘chain-terminators’. Acyclic nucleotide phosphonate analogues (ANPs) represent a particular class of nucleotide analogue that does not possess a ribose moiety. The structural basis for their substrate efficiency regarding viral DNA polymerases is poorly understood.

Methods

Pre-steady-state kinetics on HIV-1 RT together with molecular modelling, were used to evaluate the relative characteristics of both the initial binding and incorporation into DNA of three different ANP diphosphates with progressively increasing steric demands on the acyclic linker: adefovir-diphosphate (DP), tenofovir-DP, and cidofovir-DP.

Results

The increase of steric demand in ANPs induced a proportional loss of the binding affinity to wild-type HIV-1 RT ( K d cidofovir-D P>> K d tenofovir-D P> K d adefovir-DP∼ K d dNTPs), consistent with the lack of HIV-1 inhibitory activity for cidofovir. We show that, starting from adefovir-DP, the steric constraints mainly map to Gln151, as its mutation to alanine provides cidofovir-DP sensitivity. Interactions between the Gln151 residue and the methyl group of tenofovir-DP further increase with the mutation Gln151Met, resulting in a specific discrimination and low-level resistance to tenofovir-DP. This alteration is the result of a dual decrease in the binding affinity ( K d) and the catalytic rate ( k pol) of incorporation of tenofovir-DP. By contrast, the tenofovir resistance mutation K65R induces a broad ‘ k pol-dependent’ nonspecific discrimination towards the three ANPs.

Conclusions

Overall, our results show that the efficiency of ANPs to compete against natural nucleotides as substrates for RT is determined by their close interaction with specific amino acids such as Gln151 within the RT active site. These results should help us to map and predict ANP sensitivity determinants in cellular and viral DNA polymerase active sites for which the understanding of different ANP sensitivity patterns are of medical importance.

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.

Enzymatic therapeutic index of acyclovir. Viral versus human polymerase gamma specificity

We have examined the kinetics of incorporation of acyclovir triphosphate by the herpes simplex virus-1 DNA polymerase holoenzyme (Pol-UL42) and the human mitochondrial DNA polymerase using transient kinetic methods. For each enzyme, we compared the kinetic parameters for acyclovir to those governing incorporation of dGTP. The favorable ground state dissociation constant (6 microM) and rate of polymerization (10 s(-1)) afford efficient incorporation of acyclovir triphosphate by the Pol-UL42 enzyme. A discrimination factor of approximately 50 favors dGTP over acyclovir triphosphate, mostly due to a faster maximum rate of dGTP incorporation. Once incorporated, acyclovir is removed with a half-life of approximately 1 h in the presence of a normal concentration of deoxynucleoside triphosphates, leading to a high toxicity index (16,000) toward viral replication. To assess the potential for toxicity toward the host we examined the incorporation and removal of acyclovir triphosphate by the human mitochondrial DNA polymerase. These results suggest moderate inhibition of mitochondrial DNA replication defining a toxicity index of 380. This value is much higher than the value of 1.5 determined for tenofovir, another acyclic nucleoside analog. The enzymatic therapeutic index is only 42 in favoring inhibition of the viral polymerase over polymerase gamma, whereas that for tenofovir is greater than 1,200. Mitochondrial toxicity is relatively low because acyclovir is activated only in infected cells by the promiscuous viral thymidine kinase and otherwise, mitochondrial toxicity would accumulate during long term treatment.

Motions of the fingers subdomain of klentaq1 are fast and not rate limiting: implications for the molecular basis of fidelity in DNA polymerases

Various kinetic studies on nucleotide incorporation by DNA polymerases have established that a rate-limiting step occurs that is crucial in the mechanism of discrimination between correct versus incorrect nucleotide. Crystallographic studies have indicated that this step may be due to a large open-to-closed conformational transition affecting the fingers subdomain. However, there is no direct evidence to support this hypothesis. In order to investigate whether or not the open-to-closed conformational transition affecting the fingers subdomain is rate limiting, we have developed a fluorescence resonance energy transfer (FRET) system, which monitors motions of the fingers subdomain. We establish that the closing of the fingers subdomain is significantly faster than the kinetically determined rate-limiting step. We propose that the rate-limiting step occurs after the closing of the fingers subdomain and is caused by local reorganization events in the active site.

Poliovirus RNA-dependent RNA polymerase (3Dpol): pre-steady-state kinetic analysis of ribonucleotide incorporation in the presence of Mg2+

We have solved the complete kinetic mechanism for correct nucleotide incorporation catalyzed by the RNA-dependent RNA polymerase from poliovirus, 3D(pol). The phosphoryl-transfer step is flanked by two isomerization steps. The first conformational change may be related to reorientation of the triphosphate moiety of the bound nucleotide, and the second conformational change may be translocation of the enzyme into position for the next round of nucleotide incorporation. The observed rate constant for nucleotide incorporation by 3D(pol) (86 s(-1)) is dictated by the rate constants for both the first conformational change (300 s(-1)) and phosphoryl transfer (520 s(-1)). Changes in the stability of the "activated" ternary complex correlate best with changes in the observed rate constant for incorporation resulting from modification of the nucleotide. With the exception of UTP, the K(d) values for nucleotides are at least 10-fold lower than the cellular concentration of the corresponding nucleotide. Our data predict that transition mutations should occur at a frequency of 1/15000, transversion mutations should occur at a frequency of less than 1/150000, and incorporation of a 2'-deoxyribonucleotide with a correct base should occur at a frequency 1/7500. Together, these data support the conclusion that 3D(pol) is actually as faithful as an exonuclease-deficient, replicative DNA polymerase. We discuss the implications of this work on the development of RNA-dependent RNA polymerase inhibitors for use as antiviral agents.

Comparative kinetics of nucleotide analog incorporation by vent DNA polymerase

Comparative kinetic and structural analyses of a variety of polymerases have revealed both common and divergent elements of nucleotide discrimination. Although the parameters for dNTP incorporation by the hyperthermophilic archaeal Family B Vent DNA polymerase are similar to those previously derived for Family A and B DNA polymerases, parameters for analog incorporation reveal alternative strategies for discrimination by this enzyme. Discrimination against ribonucleotides was characterized by a decrease in the affinity of NTP binding and a lower rate of phosphoryl transfer, whereas discrimination against ddNTPs was almost exclusively due to a slower rate of phosphodiester bond formation. Unlike Family A DNA polymerases, incorporation of 9-[(2-hydroxyethoxy)methyl]X triphosphates (where X is adenine, cytosine, guanine, or thymine; acyNTPs) by Vent DNA polymerase was enhanced over ddNTPs via a 50-fold increase in phosphoryl transfer rate. Furthermore, a mutant with increased propensity for nucleotide analog incorporation (Vent(A488L) DNA polymerase) had unaltered dNTP incorporation while displaying enhanced nucleotide analog binding affinity and rates of phosphoryl transfer. Based on kinetic data and available structural information from other DNA polymerases, we propose active site models for dNTP, ddNTP, and acyNTP selection by hyperthermophilic archaeal DNA polymerases to rationalize structural and functional differences between polymerases.

DNA Polymerase Template Interactions Probed by Degenerate Isosteric Nucleobase Analogs

The development of novel artificial nucleobases and detailed X-ray crystal structures for primer/template/DNA polymerase complexes provide opportunities to assess DNA-protein interactions that dictate specificity. Recent results have shown that base pair shape recognition in the context of DNA polymerase must be considered a significant component. The isosteric azole carboxamide nucleobases (compounds 1-5; ) differ only in the number and placement of nitrogen atoms within a common shape and therefore present unique electronic distributions that are shown to dictate the selectivity of template-directed nucleotide incorporation by DNA polymerases. The results demonstrate how nucleoside triphosphate substrate selection by DNA polymerase is a complex phenomenon involving electrostatic interactions in addition to hydrogen bonding and shape recognition. These azole nucleobase analogs offer unique molecular tools for probing nonbonded interactions dictating substrate selection and fidelity of DNA polymerases.

Structural insights into the origins of DNA polymerase fidelity

DNA polymerases discriminate from a pool of structurally similar molecules to insert the correct nucleotide to preserve Watson-Crick base pairing rules. The ability to choose between "right and wrong" is highly dependent on the identity of the polymerase. Because naturally occurring polymerases with divergent fidelities insert incorrect nucleotides with comparable efficiencies, fidelity is primarily governed by the ability to insert the correct nucleotide. DNA polymerases generally bind the correct nucleotide with similar affinities, but low-fidelity polymerases insert correct nucleotides more slowly than higher fidelity enzymes. A comparison of crystallographic ternary substrate complexes of DNA polymerases from five families exhibiting a range of nucleotide insertion rates reveals possible structural features that lead to rapid, efficient, and faithful DNA synthesis.

The molecular mechanism of multidrug resistance by the Q151M human immunodeficiency virus type 1 reverse transcriptase and its suppression using alpha-boranophosphate nucleotide analogues

Nucleoside analogues are currently used to treat human immunodeficiency virus infections. The appearance of up to five substitutions (A62V, V75I, F77L, F116Y, and Q151M) in the viral reverse transcriptase promotes resistance to these drugs, and reduces efficiency of the antiretroviral chemotherapy. Using pre-steady state kinetics, we show that Q151M and A62V/V75I/F77L/F116Y/Q151M substitutions confer to reverse transcriptase (RT) the ability to discriminate an analogue relative to its natural counterpart, and have no effect on repair of the analogue-terminated DNA primer. Discrimination results from a selective decrease of the catalytic rate constant k(pol): 18-fold (from 7 to 0.3 s(-1)), 13-fold (from 1.9 to 0.14 s(-1)), and 12-fold (from 13 to 1 s(-1)) in the case of ddATP, ddCTP, and 3'-azido-3'-deoxythymidine 5'-triphosphate (AZTTP), respectively. The binding affinities of the triphosphate analogues for RT remain unchanged. Molecular modeling explains drug resistance by a selective loss of electrostatic interactions between the analogue and RT. Resistance was overcome using alpha-boranophosphate nucleotide analogues. Using A62V/V75I/F77L/F116Y/Q151M RT, k(pol) increases up to 70- and 13-fold using alpha-boranophosphate-ddATP and alpha-boranophosphate AZTTP, respectively. These results highlight the general capacity of such analogues to circumvent multidrug resistance when RT-mediated nucleotide resistance originates from the selective decrease of the catalytic rate constant k(pol).

Acyclic and dideoxy terminator preferences denote divergent sugar recognition by archaeon and Taq DNA polymerases

The nucleotide-binding site in a variety of DNA polymerases was probed by analyzing incorporation of dideoxy and acyclic chain terminators. Family B archaeon DNA polymerases Vent, Deep Vent, 9 degrees N and Pfu incorporated acyclic in preference to dideoxy terminators, while the Family A DNA polymerases Taq and Klenow preferred dideoxy terminators. These divergent biases suggest that significant differences exist in sugar recognition in these two classes of polymerases. Mutants of Vent (A488L) and Taq (F667Y) that increase incorporation of dideoxy terminators maintained the acyclic/dideoxy bias of the parent enzyme, while more efficiently incorporating both dideoxy and acyclic terminators. The preference of archaeon DNA polymerases for acyclic analogs was exploited in chain terminator DNA sequence and genotype analysis. This technology was additionally aided by identification of specific dye-modified bases that improve terminator incorporation over that of the unmodified terminator. These three enhancing effects, (i) acyclic analogs, (ii) archaeon variants and (iii) specific dyes, appear to act additively and independently to increase terminator incorporation efficiency, collectively enhancing incorporation approximately 8000-fold over the wild-type incorporation of dideoxynucleotides. Fluorescent chain terminator DNA sequence traces demonstrate the applicability of these advances in improving terminator incorporation, as required in DNA sequence and genotype determinations.

Mechanism-based suppression of dideoxynucleotide resistance by K65R human immunodeficiency virus reverse transcriptase using an alpha-boranophosphate nucleoside analogue

The amino acid change K65R in human immunodeficiency virus type 1-reverse transcriptase (RT) confers viral resistance to various 2',3'-dideoxynucleoside drugs in vivo. Using pre-steady state kinetic methods, we found that K65R-reverse transcriptase is 3.2-14-fold resistant to 2',3'-dideoxynucleotides in vitro relative to wild-type reverse transcriptase, in agreement with resistance levels observed in vivo. A decreased catalytic rate constant k(pol) mostly accounts for the lower incorporation efficiency observed for 2',3'-dideoxynucleotides. Examination of the crystal structure of the RT.DNA.dNTP complex suggested that both the charge at position 65 and the 3'-OH of the incoming nucleotide act in synergy during the creation of the phosphodiester bond, resulting in a more pronounced decreased catalytic rate constant for 2',3'-dideoxynucleotides than for dNTPs. This type of intramolecular activation of the leaving phosphate by the 3'-OH group appears to be conserved in several nucleotide phosphotransferases. These data were used to design dideoxynucleotide analogues targeting K65R RT specifically. alpha-Boranophosphate ddATP was found to be a 2-fold better substrate than dATP and inhibited DNA synthesis by K65R RT 153-fold better than ddATP. This complete suppression of drug resistance at the nucleotide level could serve for other reverse transcriptases for which drug resistance is achieved at the catalytic step.

Crystal structures of a ddATP-, ddTTP-, ddCTP, and ddGTP- trapped ternary complex of Klentaq1: insights into nucleotide incorporation and selectivity

The mechanism by which DNA polymerase I enzymes function has been the subject of extensive biochemical and structural studies. We previously determined the structure of a ternary complex of the large fragment of DNA polymerase I from Thermus aquaticus (Klentaq1) bound to a primer/template DNA and a dideoxycytidine 5'-triphosphate (ddCTP). In this report, we present the details of the 2.3-A resolution crystal structures of three additional ternary complexes of Klentaq1 bound to a primer/template DNA and a dideoxyguanosine 5'-triphosphate (ddGTP), a dideoxythymidine 5'-triphosphate (ddTTP), or a dideoxyadenosine 5'-triphosphate (ddATP). Comparison of the active site of the four ternary complexes reveals that the protein residues around the nascent base pair (that formed between the incoming dideoxynucleoside triphosphate [ddNTP] and the template base) form a snug binding pocket into which only a correct Watson-Crick base pair can fit. Except in the ternary complex bound to dideoxyguanosine 5'-triphosphate, there are no sequence specific contacts between the protein side chains and the nascent base pair, suggesting that steric constraints imposed by the protein onto the nascent base pair is the major contributor to nucleotide selectivity at the polymerase active site. The protein around the polymerase active site also shows plasticity, which may be responsible for the substrate diversity of the enzyme. Two conserved side chains, Q754 and R573, form hydrogen bonds with the N3 atom in the purine base and O2 atom in the pyrimidine base at the minor groove side of the base pair formed by the incorporated ddNMP and the corresponding template base in all the four ternary complexes. These hydrogen-bonding interactions may provide a means of detecting misincorporation at this position.

Prokaryotic DNA polymerase I: evolution, structure, and "base flipping" mechanism for nucleotide selection

Accurate transmission of DNA material from one generation to the next is crucial for prolonged cell survival. Following the discovery of DNA polymerse I in Escherichia coli, the DNA polymerase I class of enzymes has served as the prototype for studies on structural and biochemical mechanisms of DNA replication. Recently, a series of genomic, mutagenesis and structural investigations have provided key insights into how Pol I class of enzymes function and evolve. X-ray crystal structures of at least three Pol I class of enzymes have been solved in the presence of DNA and dNTP, thus allowing a detailed description of a productive replication complex. Rapid-quench stop-flow studies have helped define individual steps during nucleotide incorporation and conformational changes that are rate limiting during catalysis. Studies in our laboratory have generated large libraries of active mutant enzymes (8000) containing a variety of substitutions within the active site, some of which exhibit altered biochemical properties. Extensive genomic information of Pol I has recently become available, as over 50 polA genes from different prokaryotic species have been sequenced. In light of these advancements, we review here the structure-function relationships of Pol I, and we highlight those interactions that are responsible for the high fidelity of DNA synthesis. We present a mechanism for "flipping" of the complementary template base to enhance interactions with the incoming nucleotide substrate during DNA synthesis.

Enhanced ribonucleotide incorporation by an O-helix mutant of Thermus aquaticus DNA polymerase I

The O-helix of DNA polymerases has been implicated in substrate discrimination and replication fidelity. In this study, wild-type Thermus aquaticus DNA polymerase I (Taq pol I) and an O-helix mutant A661E was examined for their ability to discriminate between ribonucleotides and deoxyribonucleotides. Steady-state nucleotide extension kinetics were carried out using a template cytidine and each nucleotide dNTP and rGTP. Wild-type Taq pol I and A661E demonstrated similar Vmax and Km values for the correct nucleotide dGTP. However, A661E discriminated between incorrect and correct nucleotide less well than wild-type; discrimination was reduced by factors of 9.5-, 5.6- and 15-fold for dATP, dTTP and rGTP, respectively. These data suggest that A661E is efficient polymerases in the presence of the correct deoxynucleotide, dGTP, but it is impaired in ability to discriminate between correct and incorrect deoxyribonucleotides or between ribo- and deoxyribonucleotides. A structural model of Taq pol I is described in which the mutation A661E alters the interactions between the O-helix and the terminal two phosphate groups in the primer strand.

Role of Q190 of MuLV RT in ddNTP resistance and fidelity of DNA synthesis: a molecular model of interactions with substrates

Gln190 of MuLV reverse transcriptase (RT) plays an important role in the catalytic mechanism of MuLV RT for its conservative and non-conservative mutant derivatives exhibit low catalytic activity. We now report that both Q190N and Q190A MuLV RTs are more efficient in their activity to incorporate ddNTPs and exhibit higher fidelity than the wild-type (WT) enzyme of DNA synthesis in both RNA- and DNA-directed reactions. To obtain some insight into the structural basis for the differential utilization of dNTP and ddNTP by the mutant enzymes, we modeled the binary and the ternary complexes of MuLV RT using corresponding HIV-1 RT structures and available structure of the fragment of MuLV RT. Q190 of MuLV RT appears to be essential for the interaction with 3'OH of dNTP. The lack of a 3'OH moiety in ddNTP does not permit the binding of ddNTPs to WT MuLV RT. However, the shorter side chain of Q190N (or A) mutant MuLV RT and the absence of 3'OH in ddNTP result in the rearrangement of hydrophobic interactions favoring the binding and limited incorporation of ddNTPs. In addition, while modeling the binary and ternary complexes of MuLV RT, we noted that in the formation of the ternary complex, an interaction of Q190 with dNTP substrate requires a shift from its interaction with the template base. This may be achieved by a small conformational change or motion of the loop between beta9 and alphaH containing Q190, which may correspond to the conformational change step requiring participation of Q190 during the catalytic reaction as reported in an earlier biochemical investigation.

Poliovirus RNA-dependent RNA polymerase (3Dpol): structural, biochemical, and biological analysis of conserved structural motifs A and B

We have constructed a structural model for poliovirus RNA-dependent RNA polymerase (3D(pol)) in complex with a primer-template (sym/sub) and ATP. Residues found in conserved structural motifs A (Asp-238) and B (Asn-297) are involved in nucleotide selection. Asp-238 appears to couple binding of nucleotides with the correct sugar configuration to catalytic efficiency at the active site of the enzyme. Asn-297 is involved in selection of ribonucleoside triphosphates over 2'-dNTPs, a role mediated most likely via a hydrogen bond between the side chain of this residue and the 2'-OH of the ribonucleoside triphosphate. Substitutions at position 238 or 297 of 3D(pol) produced derivatives exhibiting a range of catalytic efficiencies when assayed in vitro for poly(rU) polymerase activity or sym/sub elongation activity. A direct correlation existed between activity on sym/sub and biological phenotypes; a 2.5-fold reduction in polymerase elongation rate produced virus with a temperature-sensitive growth phenotype. These data permit us to propose a detailed, structural model for nucleotide selection by 3D(pol), confirm the biological relevance of the sym/sub system, and provide additional evidence for kinetic coupling between RNA synthesis and subsequent steps in the virus life cycle.

Coupling ribose selection to fidelity of DNA synthesis. The role of Tyr-115 of human immunodeficiency virus type 1 reverse transcriptase

The catalytic efficiency of incorporation of deoxyribonucleotides by wild-type human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) was around 100-fold higher than for dideoxyribonucleotides, in Mg(2+)-catalyzed reactions, and more than 10,000-fold higher than for nucleotides having a 2'-hydroxyl group in Mg(2+)- and Mn(2+)-catalyzed reactions. Mutant RTs with nonconservative substitutions affecting Tyr-115 (Y115V, Y115A, and Y115G) showed a dramatic reduction in their ability to discriminate against ribonucleotides in the presence of both cations. However, selectivity of deoxyribonucleotides versus ribonucleotides was not affected in mutants Y115W and F160W. The substitution of Tyr-115 with Val or Gly had no effect on discrimination against dideoxyribonucleotides, but these mutants were less efficient than the wild-type RT in discriminating against cordycepin 5'-triphosphate. We also show that Tyr-115 is involved in fidelity of DNA synthesis, but substitutions at this position have different effects depending on the metal cofactor used in the assays. In Mg(2+)-catalyzed reactions, removal of the side chain of Tyr-115 reduced misinsertion and mispair extension fidelity, while opposite effects were observed in Mn(2+)-catalyzed reactions. Our results indicate that the aromatic side chain of Tyr-115 plays a role as a "steric gate" preventing the incorporation of nucleotides with a 2'-hydroxyl group in a cation-independent manner, while its influence on fidelity could be modulated by Mg(2+) or Mn(2+).

Structure-based design of Taq DNA polymerases with improved properties of dideoxynucleotide incorporation

The Taq DNA polymerase is the most commonly used enzyme in DNA sequencing. However, all versions of Taq polymerase are deficient in two respects: (i) these enzymes incorporate each of the four dideoxynucleoside 5' triphosphates (ddNTPs) at widely different rates during sequencing (ddGTP, for example, is incorporated 10 times faster than the other three ddNTPs), and (ii) these enzymes show uneven band-intensity or peak-height patterns in radio-labeled or dye-labeled DNA sequence profiles, respectively. We have determined the crystal structures of all four ddNTP-trapped closed ternary complexes of the large fragment of the Taq DNA polymerase (Klentaq1). The ddGTP-trapped complex structure differs from the other three ternary complex structures by a large shift in the position of the side chain of residue 660 in the O helix, resulting in additional hydrogen bonds being formed between the guanidinium group of this residue and the base of ddGTP. When Arg-660 is mutated to Asp, Ser, Phe, Tyr, or Leu, the enzyme has a marked and selective reduction in ddGTP incorporation rate. As a result, the G track generated during DNA sequencing by these Taq polymerase variants does not terminate prematurely, and higher molecular-mass G bands are detected. Another property of these Taq polymerase variants is that the sequencing patterns produced by these enzymes are remarkably even in band-intensity and peak-height distribution, thus resulting in a significant improvement in the accuracy of DNA sequencing.

Improved DNA Sequencing Accuracy and Detection of Heterozygous Alleles Using Manganese Citrate and Different Fluorescent Dye Terminators

The use of dideoxynucleotide triphosphates labeled with different fluorescent dyes (dye terminators) is the most versatile method for automated DNA sequencing. However, variation in peak heights reduces base-calling accuracy and limits heterozygous allele detection, favoring use of dye-labeled primers for this purpose. We have discovered that the addition of a manganese salt to the PE Applied Biosystems dye-terminator sequencing kits overcomes these limitations for the older rhodamine dyes as well as the more recent dichloro-rhodamine dyes (dRhodamine and BigDyes). Addition of manganese to reactions containing dRhodamine-based dye terminators produced the highest base-calling accuracy. This combination resulted in the most uniform electropherogram profiles, superior to those produced by BigDye terminators and published for dye primers, and facilitated detection of heterozygous alleles.

High performance DNA sequencing, and the detection of mutations and polymorphisms, on the Clipper sequencer

The Visible Genetics Clipper sequencer is a new platform for automated DNA sequencing which employs disposable MicroCel cassettes and 50 microm thick polyacrylamide gels. Two DNA ladders can be analyzed simultaneously in each of 16 lanes on a gel, after labeling with far-red absorbing dyes such as Cy5 and Cy5.5. This allows a simultaneous bidirectional sequencing of four templates. We have evaluated the Clipper sequencer, by cycle-sequencing of an M13 single-stranded DNA standard, and by coupled amplification and sequencing (CLIP) of reverse-transcribed human immunodeficiency virus (HIV-1) RNA standards and clinical patient samples. (i) Limitations of instrument. We have examined basic instrument parameters such as detector stability, background, digital sampling rate, and gain. With proper usage, the optical and electronic subsystems of the Clipper sequencer do not limit the data collection or sequence-determination processes. (ii) Limitations of gel performance. We have also examined the physics of DNA band separation on 50 microm thick MicroCel gels. We routinely obtain well-resolved sequence which can be base-called with 98.5% accuracy to position approximately 450 on an 11 cm gel, and to position approximately 900 on a 25 cm gel. Resolution on 5 and 11 cm gels ultimately is limited by a sharp decrease in spacing between adjacent bands, in the biased reptation separation regime. Fick's (thermal) diffusion appears to be of minor importance on 6 cm or 11 cm gels, but becomes an additional resolution-limiting factor on 25 cm gels. (iii) Limitations of enzymology. Template quality, primer nesting, choice of DNA polymerase, and choice between dye primers and dye terminators are key determinants of the ability to detect mutations and polymorphisms on the Clipper sequencer, as on other DNA sequencers. When CLIP is used with dye-labeled primers and a DNA polymerase of the F667Y, delta(5'--> 3' exo) class, we can routinely detect single-nucleotide mutations and polymorphisms over the 0.35-0.65 heterozygosity range. We present an example of detecting therapeutically relevant mutations in a clinical HIV-1 RNA isolate.

Dye structure affects Taq DNA polymerase terminator selectivity

All DNA sequencing methods have benefited from the development of new F667Y versions of Taq DNA polymerase. However, terminator chemistry methods show less uniform peak height patterns when compared to primer chemistry profiles suggesting that the dyes and/or their linker arms affect enzyme selectivity. We have measured elementary nucleotide rate and binding constants for representative rhodamine- and fluorescein-labeled terminators to determine how they interact with F667 versions of Taq Pol I. We have also developed a rapid gel-based selectivity assay that can be used to screen and to quantify dye-enzyme interactions with F667Y versions of the enzyme. Our results show that 6-TAMRA-ddTTP behaves like unlabeled ddTTP, while 6-FAM-ddTTP shows a 40-fold reduction in the rate constant for polymerization without affecting ground-state nucleotide binding. Detailed mechanism studies indicate that both isomers of different fluorescein dyes interfere with a conformational change step which the polymerase undergoes following nucleotide binding but only when these dyes are attached to pyrimidines. When these same dyes are attached to purines by the same propargylamino linker arm, they show no effect on enzyme selectivity. These studies suggest that it may be possible to develop fluorescein terminators for thermocycle DNA sequencing methods for polymerases that do not discriminate between deoxy- and dideoxynucleotides.

How E. coli DNA polymerase I (Klenow fragment) distinguishes between deoxy- and dideoxynucleotides

Deoxy- and dideoxynucleotides differ only in whether they have a hydroxyl substituent at C-3' of the ribose moiety, and yet the Klenow fragment DNA polymerase prefers the natural (dNTP) substrate by several thousandfold. We have used this preference in order to investigate how Klenow fragment interacts with the sugar portion of an incoming dNTP. We screened mutant derivatives of Klenow fragment so as to identify those amino acid residues that play important roles in distinguishing between dNTPs and ddNTPs. Substitution of Phe762 with Ala or Tyr caused a dramatic decrease in the discrimination against ddNTPs, while mutations in Tyr766 and Glu710 had a smaller effect, suggesting that these two side-chains play secondary roles in the selection of dNTPs over ddNTPs. In order to understand the interactions in the enzyme-DNA-dNTP ternary complex, pre-steady-state kinetic parameters for the incorporation of dNTPs and ddNTPs were determined for wild-type Klenow fragment and for mutant derivatives that showed changes in dNTP/ddNTP discrimination. From elemental effect measurements we infer that selection against dideoxynucleotides takes place in the transition state for the conformational change that precedes phosphoryl transfer. The crucial role of the Phe762 side-chain appears to be to constrain the dNTP molecule so that the 3'-OH can make an interaction with another group within the ternary complex. When Tyr is substituted at position 762, the same interactions can take place to position the dNTP, but specificity against the ddNTP is lost because the phenolic OH can compensate for the missing 3'-OH of the nucleotide. Substitution of the smaller Ala side-chain results in a loss in specificity because the dNTP is no longer appropriately constrained. Measurement of reaction rates as a function of magnesium ion concentration suggests that the interaction made with the dNTP 3'-OH may involve a metal ion and the Glu710 side-chain, the simplest scenario being that both the 3'-OH and the carboxylate of Glu710 are ligands to the same metal ion.