Local conformations and competitive binding affinities of single- and double-stranded primer-template DNA at the polymerization and editing active sites of DNA polymerases

The thumb subdomain of Pyrococcus furiosus DNA polymerase is responsible for deoxyuracil binding, hydrolysis and polymerization of nucleotides

B-family DNA polymerases, which are found in eukaryotes, archaea, viruses, and some bacteria, participate in DNA replication and repair. Starting from the N-terminus of archaeal and bacterial B-family DNA polymerases, three domains include the N-terminal, exonuclease, and polymerase domains. The N-terminal domain of the archaeal B-family DNA polymerase has a conserved deoxyuracil-binding pocket for specially binding the deoxyuracil base on DNA. The exonuclease domain is responsible for removing the mismatched base pair. The polymerase domain is the core functional domain and takes a highly conserved structure composed of fingers, palm and thumb subdomains. Previous studies have demonstrated that the thumb subdomain mainly functions as a DNA-binding element and has coordination with the exonuclease domain and palm subdomain. To further elucidate the possible functions of the thumb subdomain of archaeal B-family DNA polymerases, the thumb subdomain of Pyrococcus furiosus DNA polymerase was mutated, and the effects on three activities were characterized. Our results demonstrate that the thumb subdomain participates in the three activities of archaeal B-family DNA polymerases as a common structural element. Both the N-terminal deoxyuracil-binding pocket and thumb subdomain are critical for deoxyuracil binding. Moreover, the thumb subdomain assists DNA polymerization and hydrolysis reactions, but it does not contribute to the fidelity of DNA polymerization.

The shaping of a molecular linguist: How a career studying DNA energetics revealed the language of molecular communication

My personal and professional journeys have been far from predictable based on my early childhood. Owing to a range of serendipitous influences, I miraculously transitioned from a rebellious, apathetic teenage street urchin who did poorly in school to a highly motivated, disciplined, and ambitious academic honors student. I was the proverbial “late bloomer.” Ultimately, I earned my PhD in biophysical chemistry at Yale, followed by a postdoc fellowship at Berkeley. These two meccas of thermodynamics, coupled with my deep fascination with biology, instilled in me a passion to pursue an academic career focused on mapping the energy landscapes of biological systems. I viewed differential energetics as the language of molecular communication that would dictate and control biological structures, as well as modulate the modes of action associated with biological functions. I wanted to be a “molecular linguist.” For the next 50 years, my group and I used a combination of spectroscopic and calorimetric techniques to characterize the energy profiles of the polymorphic conformational space of DNA molecules, their differential ligand-binding properties, and the energy landscapes associated with mutagenic DNA damage recognition, repair, and replication. As elaborated below, the resultant energy databases have enabled the development of quantitative molecular biology through the rational design of primers, probes, and arrays for diagnostic, therapeutic, and molecular-profiling protocols, which collectively have contributed to a myriad of biomedical assays. Such profiling is further justified by yielding unique energy-based insights that complement and expand elegant, structure-based understandings of biological processes.

Understanding the Interaction of Polyelectrolyte Architectures with Proteins and Biosystems

The counterions neutralizing the charges on polyelectrolytes such as DNA or heparin may dissociate in water and greatly influence the interaction of such polyelectrolytes with biomolecules, particularly proteins. In this Review we give an overview of studies on the interaction of proteins with polyelectrolytes and how this knowledge can be used for medical applications. Counterion release was identified as the main driving force for the binding of proteins to polyelectrolytes: Patches of positive charge become multivalent counterions of the polyelectrolyte and lead to the release of counterions from the polyelectrolyte and a concomitant increase in entropy. This is shown from investigations on the interaction of proteins with natural and synthetic polyelectrolytes. Special emphasis is paid to sulfated dendritic polyglycerols (dPGS). The Review demonstrates that we are moving to a better understanding of charge-charge interactions in systems of biological relevance. Research along these lines will aid and promote the design of synthetic polyelectrolytes for medical applications.

Single-molecule FRET reveals proofreading complexes in the large fragment of Bacillus stearothermophilus DNA polymerase I

There is increasing interest in the use of DNA polymerases (DNA pols) in next-generation sequencing strategies. These methodologies typically rely on members of the A and B family of DNA polymerases that are classified as high-fidelity DNA polymerases. These enzymes possess the ability to selectively incorporate the correct nucleotide opposite a templating base with an error frequency of only 1 in 10⁶ insertion events. How they achieve this remarkable fidelity has been the subject of numerous investigations, yet the mechanism by which these enzymes achieve this level of accuracy remains elusive. Several smFRET assays were designed to monitor the conformational changes associated with the nucleotide selection mechanism(s) employed by DNA pols. smFRET has also been used to monitor the movement of DNA pols along a DNA substrate as well as to observe the formation of proof-reading complexes. One member among this class of enzymes, the large fragment of Bacillus stearothermophilus DNA polymerase I (Bst pol I LF), contains both 5'→3' polymerase and 3'→5' exonuclease domains, but reportedly lacks exonuclease activity. We have designed a smFRET assay showing that Bst pol I LF forms proofreading complexes. The formation of proofreading complexes at the single molecule level is strongly influenced by the presence of the 3' hydroxyl at the primer-terminus of the DNA substrate. Our assays also identify an additional state, observed in the presence of a mismatched primer-template terminus, that may be involved in the transfer of the primer-terminus from the polymerase to the exonuclease active site.

Analysis of amplicon-based NGS data from neurological disease gene panels: a new method for allele drop-out management

BACKGROUND

Amplicon-based targeted resequencing is a commonly adopted solution for next-generation sequencing applications focused on specific genomic regions. The reliability of such approaches rests on the high specificity and deep coverage, although sequencing artifacts attributable to PCR-like amplification can be encountered. Between these artifacts, allele drop-out, which is the preferential amplification of one allele, causes an artificial increase in homozygosity when heterozygous mutations fall on a primer pairing region. Here, a procedure to manage such artifacts, based on a pipeline composed of two steps of alignment and variant calling, is proposed. This methodology has been compared to the Illumina Custom Amplicon workflow, available on Illumina MiSeq, on the analysis of data obtained with four newly designed TruSeq Custom Amplicon gene panels.

RESULTS

Four gene panels, specific for Parkinson disease, for Intracerebral Hemorrhage Diseases (COL4A1 and COL4A2 genes) and for Familial Hemiplegic Migraine (CACNA1A and ATP1A2 genes) were designed. A total of 119 samples were re-sequenced with Illumina MiSeq sequencer and panel characterization in terms of coverage, number of variants found and allele drop-out potential impact has been carried out. Results show that 14 % of identified variants is potentially affected by allele drop-out artifacts and that both the Custom Amplicon workflow and the procedure proposed here could correctly identify them. Furthermore, a more complex configuration in presence of two mutations was simulated in silico. In this configuration, our proposed methodology outperforms Custom Amplicon workflow, being able to correctly identify two mutations in all the studied configurations.

CONCLUSIONS

Allele drop-out plays a crucial role in amplicon-based targeted re-sequencing and specific procedures in data analysis of amplicon data should be adopted. Although a consensus has been established in the elimination of primer sequences from aligned data (e.g., via primer sequence trimming or soft clipping), more complex configurations need to be managed in order to increase the retrieved information from available data. Our method shows how to manage one of these complex configurations, when two mutations occur.

Time Course Analysis of Enzyme-Catalyzed DNA Polymerization

Extracting kinetic parameters from DNA polymerase-catalyzed processive polymerization data using traditional initial-rate analysis has proven to be problematic for multiple reasons. The first substrate, DNA template, is a heterogeneous polymer and binds tightly to DNA polymerase. Further, the affinity and speed of incorporation of the second substrate, deoxynucleoside triphosphate (dNTP), vary greatly depending on the nature of the templating base and surrounding sequence. Here, we present a mathematical model consisting of the DNA template-binding step and a Michaelis-Menten-type nucleotide incorporation step acting on a DNA template with a finite length. The model was numerically integrated and globally fitted to experimental reaction time courses. The time courses were determined by monitoring the processive synthesis of oligonucleotides of lengths between 50 and 120 nucleotides by DNA polymerase I (Klenow fragment exo-) using the fluorophore PicoGreen. For processive polymerization, we were able to estimate an enzyme-template association rate k1 of 7.4 μM-1 s-1, a disassociation rate k-1 of 0.07 s-1, and a Kd of 10 nM, and the steady-state parameters for correct dNTP incorporation give kcat values of 2.5-3.3 s-1 and Km values of 0.51-0.86 μM. From the analysis of time courses measured between 5 and 25 °C, an activation energy for kcat of 82 kJ mol-1 was calculated, and it was found that up to 73% of Klenow fragment becomes inactivated or involved in unproductive binding at lower temperatures. Finally, a solvent deuterium kinetic isotope effect (KIE) of 3.0-3.2 was observed under processive synthesis conditions, which suggests that either the intrinsic KIE is unusually high, at least 30-40, or previous findings, showing that the phosphoryl transfer step occurs rapidly and is flanked by two slow conformational changes, need to be re-evaluated. We suggest that the numerical integration of rate equations provides a high level of flexibility and generally produces superior results compared to those of initial-rate analysis in the study of DNA polymerase kinetics and, by extension, other complex enzyme systems.

Mapping the interactions of the single-stranded DNA binding protein of bacteriophage T4 (gp32) with DNA lattices at single nucleotide resolution: polynucleotide binding and cooperativity

We here use our site-specific base analog mapping approach to study the interactions and binding equilibria of cooperatively-bound clusters of the single-stranded DNA binding protein (gp32) of the T4 DNA replication complex with longer ssDNA (and dsDNA) lattices. We show that in cooperatively bound clusters the binding free energy appears to be equi-partitioned between the gp32 monomers of the cluster, so that all bind to the ssDNA lattice with comparable affinity, but also that the outer domains of the gp32 monomers at the ends of the cluster can fluctuate on and off the lattice and that the clusters of gp32 monomers can slide along the ssDNA. We also show that at very low binding densities gp32 monomers bind to the ssDNA lattice at random, but that cooperatively bound gp32 clusters bind preferentially at the 5′-end of the ssDNA lattice. We use these results and the gp32 monomer-binding results of the companion paper to propose a detailed model for how gp32 might bind to and interact with ssDNA lattices in its various binding modes, and also consider how these clusters might interact with other components of the T4 DNA replication complex.

Characterization of PA-N terminal domain of Influenza A polymerase reveals sequence specific RNA cleavage

Influenza virus uses a unique cap-snatching mechanism characterized by hijacking and cleavage of host capped pre-mRNAs, resulting in short capped RNAs, which are used as primers for viral mRNA synthesis. The PA subunit of influenza polymerase carries the endonuclease activity that catalyzes the host mRNA cleavage reaction. Here, we show that PA is a sequence selective endonuclease with distinct preference to cleave at the 3′ end of a guanine (G) base in RNA. The G specificity is exhibited by the native influenza polymerase complex associated with viral ribonucleoprotein particles and is conferred by an intrinsic G specificity of the isolated PA endonuclease domain PA-Nter. In addition, RNA cleavage site choice by the full polymerase is also guided by cap binding to the PB2 subunit, from which RNA cleavage preferentially occurs at the 12th nt downstream of the cap. However, if a G residue is present in the region of 10–13 nucleotides from the cap, cleavage preferentially occurs at G. This is the first biochemical evidence of influenza polymerase PA showing intrinsic sequence selective endonuclease activity.

Dynamics of site switching in DNA polymerase

DNA polymerases replicate DNA by catalyzing the template-directed polymerization of deoxynucleoside triphosphate (dNTP) substrates onto the 3' end of a growing DNA primer strand. Many DNA polymerases also possess a separate 3'-5' exonuclease activity that is used to remove misincorporated nucleotides from the nascent DNA (proofreading). The polymerase (pol) and exonuclease (exo) activities are spatially separated in different enzyme domains, indicating that a mechanism must exist to transfer the growing primer terminus from one site to the other. Here we report a single-molecule Förster resonance energy transfer (smFRET) system that directly monitors the movement of a DNA substrate between the pol and exo sites of DNA polymerase I Klenow fragment (KF). FRET trajectories recorded during the encounter between single polymerase and DNA molecules reveal that DNA can channel between the pol and exo sites in both directions while remaining closely associated with the enzyme (intramolecular transfer). In addition, it is evident from the trajectories that DNA can also dissociate from one site and subsequently rebind at the other (intermolecular transfer). Rate constants for each pathway have been determined by dwell-time analysis, revealing that intramolecular transfer is the faster of the two pathways. Unexpectedly, a mispaired primer terminus accesses the exo site more frequently when dNTP substrates are also present in solution, which is expected to enhance proofreading. Together, these results explain how the separate pol and exo activities of KF are physically coordinated to achieve efficient proofreading.

Unwinding of primer-templates by archaeal family-B DNA polymerases in response to template-strand uracil

Archaeal family-B DNA polymerases bind tightly to deaminated bases and stall replication on encountering uracil in template strands, four bases ahead of the primer-template junction. Should the polymerase progress further towards the uracil, for example, to position uracil only two bases in front of the junction, 3′–5′ proof-reading exonuclease activity becomes stimulated, trimming the primer and re-setting uracil to the +4 position. Uracil sensing prevents copying of the deaminated base and permanent mutation in 50% of the progeny. This publication uses both steady-state and time-resolved 2-aminopurine fluorescence to show pronounced unwinding of primer-templates with Pyrococcus furiosus (Pfu) polymerase–DNA complexes containing uracil at +2; much less strand separation is seen with uracil at +4. DNA unwinding has long been recognized as necessary for proof-reading exonuclease activity. The roles of M247 and Y261, amino acids suggested by structural studies to play a role in primer-template unwinding, have been probed. M247 appears to be unimportant, but 2-aminopurine fluorescence measurements show that Y261 plays a role in primer-template strand separation. Y261 is also required for full exonuclease activity and contributes to the fidelity of the polymerase.

[The fidelity mechanism of DNA synthesis]

Accurate DNA synthesis is vital to maintain genome stability and ensure propagation of species. Synthetic errors have far reaching consequences. Therefore, DNA synthesis is remarkably accurate. The high fidelity is mainly achieved through three steps: ① nucleotide selection, which is based on hydrogen, base pair shape, or some other elements; ② 3'→5' exonuclease proofreading, which removes mis-incorporated nucleotides in cis or trans; ③ repair process, which could correct mismatched nucleotides escaping from proofreading, such as mismatch repair, excission repair, homologous recombination repair, and translesion DNA synthesis. Because all polymerases are suitable targets for anticancer or antiviral drugs, their fidelity is involved in drug resistance and side effects. Understanding the molecular basis of synthesis fidelity is of vital importance. In this review, the fidelity mechanisms of DNA synthesis will be discussed in detail. Furthermore, their application perspective was discussed.

Differential temperature-dependent multimeric assemblies of replication and repair polymerases on DNA increase processivity

Differentiation of binding accurate DNA replication polymerases over error prone DNA lesion bypass polymerases is essential for the proper maintenance of the genome. The hyperthermophilic archaeal organism Sulfolobus solfataricus (Sso) contains both a B-family replication (Dpo1) and a Y-family repair (Dpo4) polymerase and serves as a model system for understanding molecular mechanisms and assemblies for DNA replication and repair protein complexes. Protein cross-linking, isothermal titration calorimetry, and analytical ultracentrifugation have confirmed a previously unrecognized dimeric Dpo4 complex bound to DNA. Binding discrimination between these polymerases on model DNA templates is complicated by the fact that multiple oligomeric species are influenced by concentration and temperature. Temperature-dependent fluorescence anisotropy equilibrium binding experiments were used to separate discrete binding events for the formation of trimeric Dpo1 and dimeric Dpo4 complexes on DNA. The associated equilibria are found to be temperature-dependent, generally leading to improved binding at higher temperatures for both polymerases. At high temperatures, DNA binding of Dpo1 monomer is favored over binding of Dpo4 monomer, but binding of Dpo1 trimer is even more strongly favored over binding of Dpo4 dimer, thus providing thermodynamic selection. Greater processivities of nucleotide incorporation for trimeric Dpo1 and dimeric Dpo4 are also observed at higher temperatures, providing biochemical validation for the influence of tightly bound oligomeric polymerases. These results separate, quantify, and confirm individual and sequential processes leading to the formation of oligomeric Dpo1 and Dpo4 assemblies on DNA and provide for a concentration- and temperature-dependent discrimination of binding undamaged DNA templates at physiological temperatures.

Single-molecule microscopy reveals new insights into nucleotide selection by DNA polymerase I

The mechanism by which DNA polymerases achieve their extraordinary accuracy has been intensely studied because of the linkage between this process and mutagenesis and carcinogenesis. Here, we have used single-molecule fluorescence microscopy to study the process of nucleotide selection and exonuclease action. Our results show that the binding of Escherichia coli DNA polymerase I (Klenow fragment) to a primer-template is stabilized by the presence of the next correct dNTP, even in the presence of a large excess of the other dNTPs and rNTPs. These results are consistent with a model where nucleotide selection occurs in the open complex prior to the formation of a closed ternary complex. Our assay can also distinguish between primer binding to the polymerase or exonuclease domain and, contrary to ensemble-averaged studies, we find that stable exonuclease binding only occurs with a mismatched primer terminus.

Characterization of the 6-methyl isoxanthopterin (6-MI) base analog dimer, a spectroscopic probe for monitoring guanine base conformations at specific sites in nucleic acids

We here characterize local conformations of site-specifically placed pairs of guanine (G) residues in RNA and DNA, using 6-methyl isoxanthopterin (6-MI) as a conformational probe. 6-MI is a base analog of G and spectroscopic signals obtained from pairs of adjacent 6-MI residues reflect base–base interactions that are sensitive to the sequence context, local DNA conformation and solvent environment of the probe bases. CD signals show strong exciton coupling between stacked 6-MI bases in double-stranded (ds) DNA; this coupling is reduced in single-stranded (ss) DNA sequences. Solvent interactions reduce the fluorescence of the dimer probe more efficiently in ssDNA than dsDNA, while self-quenching between 6-MI bases is enhanced in dsDNA. 6-MI dimer probes closely resemble adjacent GG residues, in that these probes have minimal effects on the stability of dsDNA and on interactions with solvent additive betaine. They also serve as effective template bases, although further polymerase-dependent extension of DNA primers past 6-MI template bases is significantly inhibited. These probes are also used to monitor DNA ‘breathing’ at model replication forks. The 6-MI dimer probe can serve in many contexts as a useful tool to investigate GG conformations at specific sites within the nucleic acid frameworks of functioning macromolecular machines in solution.

Interactions of replication versus repair DNA substrates with the Pol I DNA polymerases from Escherichia coli and Thermus aquaticus

Different DNA polymerases partition differently between replication and repair pathways. In this study we examine if two Pol I family polymerases from evolutionarily distant organisms also differ in their preferences for replication versus repair substrates. The DNA binding preferences of Klenow and Klentaq DNA polymerases, from Escherichia coli and Thermus aquaticus respectively, have been studied using a fluorescence competition binding assay. Klenow polymerase binds primed-template DNA (the replication substrate) with up to 50× higher affinity than it binds to nicked DNA, DNA with a 2 base single-stranded gap, blunt-ended DNA, or to a DNA end with a 3' overhang. In contrast, Klentaq binds all of these DNAs almost identically, indicating that Klenow has a stronger ability to discriminate between replication and repair substrates than Klentaq. In contrast, both polymerases bind mismatched primed-template and blunt-ended DNA tighter than they bind matched primed-template DNA, suggesting that these two proteins may share a similar mechanism to identify mismatched DNA, despite the fact that Klentaq has no proofreading ability. In addition, the presence or absence of 5'- or 3'-phosphates has slightly different effects on DNA binding by the two polymerases, but again reinforce Klenow's more effective substrate discrimination capability.

DNA conformational changes at the primer-template junction regulate the fidelity of replication by DNA polymerase

Local conformational changes in primer-template (P/T) DNA are involved in the selective incorporation of dNTP by DNA polymerases (DNAP). Here we use near UV CD and fluorescence spectra of pairs of base analogue probes, substituted either at the primer terminus or in the coding region of the template strand, to monitor and interpret conformational changes at and near the coding base of the template in P/T DNA complexes with Klenow fragment (KF) DNAP as the polymerase moves through the nucleotide addition cycle. Incoming dNTPs and rNTPs encounter binary complexes in which the 3'-end of the primer shuttles between the polymerization (pol) and exonuclease (exo) sites of DNAPs, even for perfectly complementary P/T DNA sequences. We have used spectral changes of probes inserted in both strands to monitor this two-state distribution and determine how it depends on the formation of ternary complexes with both complementary ("correct") and noncomplementary ("incorrect") NTPs and on the local sequence of the P/T DNA. The results show that the relative occupancy of the exo and pol sites is coupled to conformational changes in the P/T DNA of the complex that are partially regulated by the incoming NTP. We find that the coding base on the template strand is unperturbed by the binding of incorrect dNTPs, while binding of complementary rNTPs induces a novel template conformation. We conclude that, in addition to its editing function, primer strand occupancy of the 3'-exo site may also serve as a regulatory checkpoint for accurate dNTP selection in DNA synthesis.

Thermodynamics of the DNA structural selectivity of the Pol I DNA polymerases from Escherichia coli and Thermus aquaticus

Understanding the thermodynamics of substrate selection by DNA polymerase I is important for characterizing the balance between replication and repair for this enzyme in vivo. Due to their sequence and structural similarities, Klenow and Klentaq, the large fragments of the Pol I DNA polymerases from Escherichia coli and Thermus aquaticus, are considered functional homologs. Klentaq, however, does not have a functional proofreading site. Examination of the DNA binding thermodynamics of Klenow and Klentaq to different DNA structures: single-stranded DNA (ss-DNA), primer-template DNA (pt-DNA), and blunt-end double-stranded DNA (ds-DNA) show that the binding selectivity pattern is similar when examined across a wide range of salt concentration, but can significantly differ at any individual salt concentration. For both proteins, binding of single-stranded DNA shifts from weakest to tightest binding of the three structures as the salt concentration increases. Both Klenow and Klentaq release two to three more ions when binding to pt-DNA and ds-DNA than when binding to ss-DNA. Klenow exhibits significant differences in the Delta C(p) of binding to pt-DNA versus ds-DNA, and a difference in pI for these two complexes, whereas Klentaq does not, suggesting that Klenow and Klentaq discriminate between these two structures differently. Taken together, the data suggest that the two polymerases bind ds-DNA very differently, but that both bind pt-DNA and ss-DNA similarly, despite the absence of a proofreading site in Klentaq.

Single-molecule measurements of synthesis by DNA polymerase with base-pair resolution

The catalytic mechanism of DNA polymerases involves multiple steps that precede and follow the transfer of a nucleotide to the 3'-hydroxyl of the growing DNA chain. Here we report a single-molecule approach to monitor the movement of E. coli DNA polymerase I (Klenow fragment) on a DNA template during DNA synthesis with single base-pair resolution. As each nucleotide is incorporated, the single-molecule Förster resonance energy transfer intensity drops in discrete steps to values consistent with single-nucleotide incorporations. Purines and pyrimidines are incorporated with comparable rates. A mismatched primer/template junction exhibits dynamics consistent with the primer moving into the exonuclease domain, which was used to determine the fraction of primer-termini bound to the exonuclease and polymerase sites. Most interestingly, we observe a structural change after the incorporation of a correctly paired nucleotide, consistent with transient movement of the polymerase past the preinsertion site or a conformational change in the polymerase. This may represent a previously unobserved step in the mechanism of DNA synthesis that could be part of the proofreading process.