Kinetic Mechanism of DNA Polymerases: Contributions of Conformational Dynamics and a Third Divalent Metal Ion

Kinetic Investigation of Resistance to Islatravir Conferred by Mutations in HIV-1 Reverse Transcriptase

Islatravir (EFdA) is a novel nucleoside reverse transcriptase translocation inhibitor (NRTTI) that potently blocks HIV-1 replication in vivo. Its unique structural features in contrast to nucleoside reverse transcriptase inhibitors (NRTIs), particularly the 4'-ethynyl and 3'-hydroxy groups, contribute to its high clinical potency. Once intracellularly activated to EFdA 5'-triphosphate (EFdA-TP), it competes with dATP for incorporation by HIV-1 reverse transcriptase (RT) during HIV-1 genomic replication. The 4'-ethynyl group of incorporated EFdA-MP interacts with a hydrophobic pocket of HIV-1 RT, hindering DNA translocation and terminating DNA synthesis. The M184V mutation, commonly associated with resistance to NRTIs such as lamivudine and emtricitabine, and the M184V/A114S mutations, both located within the hydrophobic pocket, were shown to reduce Islatravir susceptibility in cell-based viral resistance selection assays. To elucidate the mechanisms by which these mutations affect Islatravir inhibition, we employed pre-steady-state kinetics to investigate their impact on EFdA-TP incorporation by HIV-1 RT using both DNA and RNA templates. We found that M184V had a modest effect on EFdA-TP incorporation efficiency, increasing it 2-fold with the DNA template and decreasing it 3-fold with the RNA template. In contrast, M184V/A114S significantly inhibited EFdA-TP incorporation, reducing its incorporation efficiency 5.4-fold with the DNA template and 181-fold with the RNA template. These reductions were primarily attributable to corresponding decreases in EFdA-TP incorporation rate constants of 18-fold and 105-fold, respectively. These results suggest that, unlike FDA-approved NRTIs, the clinical efficacy of Islatravir, may not be substantially compromised by the M184V mutation alone but will be significantly reduced by the M184V/A114S mutations.

An Endo IV-driven signal amplifier for sensitive detection of cancer biomarker TdT

Terminal deoxynucleotidyl transferase (TdT) is a template-independent DNA polymerase and its function is mainly involved in the development of the human immune system. The activity of TdT is normally suppressed in healthy human cells but is overexpressed in certain diseases, such as leukemia and other malignant blood disorders. Therefore, the detection of TdT activity is crucial for clinical diagnosis and disease progression monitoring. In this study, we developed an Endo IV-driven signal amplifier to assess TdT activity. TdT can recognize primers and add dATP to the 3' hydroxyl end of the primers, forming a polyA tail. Fluorescent probes (polyT) containing AP sites can bind to the polyA tail and the AP site will be hydrolyzed by Endo IV. Then the fluorescent probe is cut and emits fluorescence. Subsequently, the intact fluorescent probe rebinds to the polyA tail, creating a cyclic process for signal amplification. The detection limit for TdT was 1 × 10-4 U μL-1, with a linear range of 1 × 10-4 to 1 × 10-1 U μL-1. This biosensor expands the current TdT detection toolkit. Its excellent specificity and sensitivity in various kinds of biological environments provide a powerful detection tool for clinical TdT testing.

Recent Advances in Biocatalytic and Chemoenzymatic Synthesis of Oligonucleotides

Access to synthetic oligonucleotides is crucial for applications in diagnostics, therapeutics, synthetic biology, and nanotechnology. Traditional solid phase synthesis is limited by sequence length and complexities, low yields, high costs and poor sustainability. Similarly, polymerase-based approaches such as in vitro transcription and primer extension reactions do not permit any control on the positioning of modifications and display poor substrate tolerance. In response, biocatalytic and chemoenzymatic strategies have emerged as promising alternatives, offering selective and efficient pathways for oligonucleotide synthesis. These methods leverage the precision and efficiency of enzymes to construct oligonucleotides with high fidelity. Recent advancements have focused on optimized systems and/or engineered enzymes enabling the incorporation of chemically modified nucleotides. Biocatalytic approaches, particularly those using DNA/RNA polymerases provide advantages in milder reaction conditions and enhanced sustainability. Chemoenzymatic methods, combining chemical synthesis and enzymes, have proven to be effective in overcoming limitations of traditional solid phase synthesis. This review summarizes recent developments in biocatalytic and chemoenzymatic strategies to construct oligonucleotides, highlighting innovations in enzyme engineering, substrate and reaction condition optimization for various applications. We address crucial details of the methods, their advantages, and limitations as well as important insights for future research directions in oligonucleotide production.

Pre-Steady-State Kinetic Studies of Nucleotide Incorporation into a Single-Nucleotide Gapped DNA Substrate Catalyzed by Human DNA Polymerase β

DNA polymerase β (Polβ) is a key enzyme in DNA base excision repair (BER). Despite extensive research, several microscopic rate constants within the kinetic mechanism of nucleotide incorporation into single-nucleotide gapped DNA by Polβ have not been determined and the identity of the rate-limiting step remains controversial. Here, we employed pre-steady-state kinetic methods and determined the rate constants for correct dNTP association (k2 = 4.5 × 106 M-1 s-1) and dissociation (k-2 = 118 s-1) as well as DNA product release (k7=0.93 s-1). Previously, uncertainty regarding the transition state of phosphodiester bond formation has led to confusion regarding the interpretation of the sulfur elemental effect between the incorporations of dNTP and its thio analog Sp-dNTPαS. However, recent results from time-resolved X-ray crystallographic studies of three DNA polymerases have allowed us to revise the benchmark of sulfur elemental effect for a rate-limiting chemistry step from 4-11 to 10-160. By using the revised benchmark, we determined the sulfur elemental effects for correct and incorrect nucleotide incorporation to be 3.94 and 64.6, respectively. These suggest the chemistry step limits mismatched, but not matched, nucleotide incorporation. Furthermore, the 2.1-fold difference in the reaction amplitudes of the pulse-quench and pulse-chase assays provides definitive evidence that a protein conformational change step prior to the chemistry step is rate-limiting for matched nucleotide incorporation. These findings unify the kinetic mechanism of correct nucleotide incorporation for Polβ and all other kinetically characterized DNA polymerases and reverse transcriptases, in which the protein conformational change prior to the chemistry step is rate-limiting.

Semi-rational evolution of a recombinant DNA polymerase for modified nucleotide incorporation efficiency

Engineering improved B-family DNA polymerases to catalyze 3'-O-modified nucleotide reversible terminators is limited by an insufficient understanding of the structural determinants that define polymerization efficiency. To explore the key mechanism for unnatural nucleotide incorporation, we engineered a B-family DNA polymerase from Thermococcus Kodakaraenis (KOD pol) by using semi-rational design strategies. We first scanned the active pocket of KOD pol through site-directed saturation mutagenesis and combinatorial mutations and identified a variant Mut_C2 containing five mutation sites (D141A, E143A, L408I, Y409A, A485E) using a high-throughput microwell-based screening method. Mut_C2 demonstrated high catalytic efficiency in incorporating 3'-O-azidomethyl-dATP labeled with a Cy3 dye, whereas the wild-type KOD pol failed to catalyze it. Computational simulations were then conducted of the DNA binding region of KOD pol to predict additional mutations with enhanced catalytic activity, which were subsequently experimentally verified. By a stepwise combinatorial mutagenesis approach, we obtained an eleven-mutation variant, named Mut_E10 by introducing additional mutations to the Mut_C2 variant. Mut_E10, which carried six specific mutations (S383T, Y384F, V389I, V589H, T676K, and V680M) within the DNA-binding region, demonstrated over 20-fold improvement in enzymatic activity as compared to Mut_C2. In addition, Mut_E10 demonstrated satisfactory performance in two different sequencing platforms (BGISEQ-500 and MGISEQ-2000), indicating its potential for commercialization. Our study demonstrates that a significant enhancement in its catalytic efficiency towards modified nucleotides can be achieved efficiently through combinatorial mutagenesis of residues in the active site and DNA binding region of DNA polymerases. These findings contribute to a comprehensive understanding of the mechanisms that underlie the incorporation of modified nucleotides by DNA polymerase. The sites of beneficial mutations, as well as the nucleotide incorporation mechanism identified in this study, can provide valuable guidance for the engineering of other B-family DNA polymerases.

Metals in Motion: Understanding Labile Metal Pools in Bacteria

Metal ions are essential for all life. In microbial cells, potassium (K+) is the most abundant cation and plays a key role in maintaining osmotic balance. Magnesium (Mg2+) is the dominant divalent cation and is required for nucleic acid structure and as an enzyme cofactor. Microbes typically require the transition metals manganese (Mn), iron (Fe), copper (Cu), and zinc (Zn), although the precise set of metal ions needed to sustain life is variable. Intracellular metal pools can be conceptualized as a chemically complex mixture of rapidly exchanging (labile) ions, complemented by those reservoirs that exchange slowly relative to cell metabolism (sequestered). Labile metal pools are buffered by transient interactions with anionic metabolites and macromolecules, with the ribosome playing a major role. Sequestered metal pools include many metalloproteins, cofactors, and storage depots, with some pools redeployed upon metal depletion. Here, I review the size, composition, and dynamics of intracellular metal pools and highlight the major gaps in understanding.

Unveiling the Mechanism of Deprotonation and Proton Transfer of DNA Polymerase Catalysis via Single-Molecule Conductance

DNA polymerases (Pols) play important roles in the transmission of genetic information. Although the function and (de)regulation of Pols are linked to many human diseases, the key mechanism of 3'-OH deprotonation and the PPi formation are not totally clear. In this work, a method is presented to detect the full catalytic cycle of human Pol (hPol β) in graphene-molecule-graphene single-molecule junctions. Real-time in situ monitoring successfully revealed the spatial and temporal properties of the open and closed conformation states of hPol β, distinguishing the reaction states in the Pols catalytic cycle and unveiling 3'-OH deprotonation and pyrophosphate (PPi) formation mechanism of hPol β. Proton inventory experiment demonstrated that the rate-limiting step of PPi formation is deprotonation, which occurs before a reverse conformational change. Additionally, by detecting the acidity (pKa), it is found that MgA-bound OH- acted as a general base and activated the nucleophile of 3'-OH, and that acidic residue D190 or D192 coordinated with MgB as a proton donor to PPi. This work provides useful insights into a fundamental chemical reaction that impacts genome synthesis efficiency and Pol fidelity, which the discovery of Pol-targeting drugs and design of artificial Pols for DNA synthetic applications are expected to accelerated.

Mechanistic Basis for a Single Amino Acid Residue Mutation Causing Human DNA Ligase 1 Deficiency, A Rare Pediatric Disease

In mammalian cells, DNA ligase 1 (LIG1) functions as the primary DNA ligase in both genomic replication and single-strand break repair. Several reported mutations in human LIG1, including R305Q, R641L, and R771W, cause LIG1 syndrome, a primary immunodeficiency. While the R641L and R771W mutations, respectively located in the nucleotidyl transferase and oligonucleotide binding domains, have been biochemically characterized and shown to reduce catalytic efficiency, the recently reported R305Q mutation within the DNA binding domain (DBD) remains mechanistically unexplored. The R641L and R771W mutations are known to decrease the catalytic activity of LIG1 by affecting both interdomain interactions and DNA binding during catalysis, without significantly impacting overall DNA affinity. To elucidate the molecular basis of the LIG1 syndrome-causing R305Q mutation, we purified this single-residue mutant protein and investigated its secondary structure, protein stability, DNA binding affinity, and catalytic efficiency. Our findings reveal that the R305Q mutation significantly impairs the function of LIG1 by disrupting the DBD-DNA interactions, leading to a 7-21-fold lower DNA binding affinity and a 33-300-fold reduced catalytic efficiency of LIG1. Additionally, the R305Q mutation slightly decreases LIG1's protein stability by 2 to 3.6 °C, on par with the effect observed previously with either the R641L or R771W mutant. Collectively, our results uncover a new mechanism whereby the R305Q mutation impairs LIG1-catalyzed nicked DNA ligation, resulting in LIG1 syndrome, and highlight the crucial roles of the DBD-DNA interactions in tight DNA binding and efficient LIG1 catalysis.

Third Metal Ion Dictates the Catalytic Activity of the Two-Metal-Ion Pre-Ribosomal RNA-Processing Machinery

Nucleic acid processing enzymes use a two-Mg2+-ion motif to promote the formation and cleavage of phosphodiester bonds. Yet, recent evidence demonstrates the presence of spatially conserved second-shell cations surrounding the catalytic architecture of proteinaceous and RNA-dependent enzymes. The RNase mitochondrial RNA processing (MRP) complex, which cleaves the ribosomal RNA (rRNA) precursor at the A3 cleavage site to yield mature 5'-end of 5.8S rRNA, hosts in the catalytic core one atypically-located Mg2+ ion, in addition to the ions forming the canonical catalytic motif. Here, we employ biased quantum classical molecular dynamics simulations of RNase MRP to discover that the third Mg2+ ion inhibits the catalytic process. Instead, its displacement in favour of a second-shell monovalent K+ ion propels phosphodiester bond cleavage by enabling the formation of a specific hydrogen bonding network that mediates the essential proton transfer step. This study points to a direct involvement of a transient K+ ion in the catalytic cleavage of the phosphodiester bond and implicates cation trafficking as a general mechanism in nucleic acid processing enzymes and ribozymes.

Mechanism of Nucleic Acid Phosphodiester Bond Cleavage by Human Endonuclease V: MD and QM/MM Calculations Reveal a Versatile Metal Dependence

Human endonuclease V (EndoV) catalytically removes deaminated nucleobases by cleaving the phosphodiester bond as part of RNA metabolism. Despite being implicated in several diseases (cancers, cardiovascular diseases, and neurological disorders) and potentially being a useful tool in biotechnology, details of the human EndoV catalytic pathway remain unclear due to limited experimental information beyond a crystal structure of the apoenzyme and select mutational data. Since a mechanistic understanding is critical for further deciphering the central roles and expanding applications of human EndoV in medicine and biotechnology, molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) calculations were used to unveil the atomistic details of the catalytic pathway. Due to controversies surrounding the number of metals required for nuclease activity, enzyme-substrate models with different numbers of active site metals and various metal-substrate binding configurations were built based on structural data for other nucleases. Subsequent MD simulations revealed the structure and stability of the human EndoV-substrate complex for a range of active site metal binding architectures. Four unique pathways were then characterized using QM/MM that vary in metal number (one versus two) and modes of substrate coordination [direct versus indirect (water-mediated)], with several mechanisms being fully consistent with experimental structural, kinetic, and mutational data for related nucleases, including members of the EndoV family. Beyond uncovering key roles for several active site amino acids (D240 and K155), our calculations highlight that while one metal is essential for human EndoV activity, the enzyme can benefit from using two metals due to the presence of two suitable metal binding sites. By directly comparing one- versus two-metal-mediated P-O bond cleavage reactions within the confines of the same active site, our work brings a fresh perspective to the "number of metals" controversy.

Exploring the molecular aspect and updating evolutionary approaches to the DNA polymerase enzymes for biotechnological needs: A comprehensive review

DNA polymerases are essential enzymes that play a key role in living organisms, as they participate in the synthesis and maintenance of the DNA molecule. The intrinsic properties of these enzymes have been widely observed and studied to understand their functions, activities, and behavior, which has allowed their natural power in DNA synthesis to be exploited in modern biotechnology, to the point of making them true pillars of the field. In this context, the laboratory evolution of these enzymes, either by directed evolution or rational design, has led to the generation of a wide range of new DNA polymerases with novel properties, suitable for a variety of biotechnological needs. In this review, we examine DNA polymerases at the molecular level, their biotechnological use, and their evolutionary methods in relation to the novel properties sought, providing a chronological selection of evolved DNA polymerases cited in the literature that we consider to be of great interest. To our knowledge, this work is the first to bring together the molecular, functional and evolutionary aspects of the DNA polymerase enzyme. We believe it will be of great interest to researchers whose aim is to produce new lines of evolved DNA polymerases.

Melanoma-Derived DNA Polymerase Theta Variants Exhibit Altered DNA Polymerase Activity

DNA polymerase θ (Pol θ or POLQ) is primarily involved in repairing double-stranded breaks in DNA through an alternative pathway known as microhomology-mediated end joining (MMEJ) or theta-mediated end joining (TMEJ). Unlike other DNA repair polymerases, Pol θ is thought to be highly error-prone yet critical for cell survival. We have identified several POLQ gene variants from human melanoma tumors that experience altered DNA polymerase activity, including a propensity for incorrect nucleotide selection and reduced polymerization rates compared to WT Pol θ. Variants are 30-fold less efficient at incorporating a nucleotide during repair and up to 70-fold less accurate at selecting the correct nucleotide opposite a templating base. This suggests that aberrant Pol θ has reduced DNA repair capabilities and may also contribute to increased mutagenesis. Moreover, the variants were identified in established tumors, suggesting that cancer cells may use mutated polymerases to promote metastasis and drug resistance.

Three-metal ion mechanism of cross-linked and uncross-linked DNA polymerase β: A theoretical study

In our recent publication, we have proposed a revised base excision repair pathway in which DNA polymerase β (Polβ) catalyzes Schiff base formation prior to the gap-filling DNA synthesis followed by β-elimination. In addition, the polymerase activity of Polβ employs the "three-metal ion mechanism" instead of the long-standing "two-metal ion mechanism" to catalyze phosphodiester bond formation based on the fact derived from time-resolved x-ray crystallography that a third Mg2+ was captured in the polymerase active site after the chemical reaction was initiated. In this study, we develop the models of the uncross-linked and cross-linked Polβ complexes and investigate the "three-metal ion mechanism" vs the "two-metal ion mechanism" by using the quantum mechanics/molecular mechanics molecular dynamics simulations. Our results suggest that the presence of the third Mg2+ ion stabilizes the reaction-state structures, strengthens correct nucleotide binding, and accelerates phosphodiester bond formation. The improved understanding of Polβ's catalytic mechanism provides valuable insights into DNA replication and damage repair.

Viral RNA Replication Suppression of SARS-CoV-2: Atomistic Insights into Inhibition Mechanisms of RdRp Machinery by ddhCTP

The nonstructural protein 12, known as RNA-dependent RNA polymerase (RdRp), is essential for both replication and repair of the viral genome. The RdRp of SARS-CoV-2 has been used as a promising candidate for drug development since the inception of the COVID-19 spread. In this work, we performed an in silico investigation on the insertion of the naturally modified pyrimidine nucleobase ddhCTP into the SARS-CoV-2 RdRp active site, in a comparative analysis with the natural one (CTP). The modification in ddhCTP involves the removal of the 3'-hydroxyl group that prevents the addition of subsequent nucleotides into the nascent strand, acting as an RNA chain terminator inhibitor. Quantum mechanical investigations helped to shed light on the mechanistic source of RdRp activity on the selected nucleobases, and comprehensive all-atom simulations provided insights about the structural rearrangements occurring in the active-site region when inorganic pyrophosphate (PPi) is formed. Subsequently, the intricate pathways for the release of PPi, the catalytic product of RdRp, were investigated using Umbrella Sampling simulations. The results are in line with the available experimental data and contribute to a more comprehensive point of view on such an important viral enzyme.

A wide survey of heavy metals-induced in-vitro DNA replication stress characterized by rate-limited replication

Heavy metals (HMs) are environmental pollutants that pose a threat to human health and have been accepted to cause various diseases, including cancer and developmental disorders. DNA replication stress has been identified to be associated with such diseases. However, the effect of HMs exclusively on DNA replication stress is still not well understood. In this study, DNA replication stress induced by thirteen HMs was assessed using a simplified in-vitro DNA replication model. Two parameters, Cte/Ctc reflecting the cycle threshold value alteration and Ke/Kc reflecting the linear phase slope change, were calculated based on the DNA replication amplification curve to evaluate the rate of exponential and linear phases. These parameters were used to detect the replication rate reflecting in-vitro DNA replication stress induced by tested HMs. According to the effective concentrations and rate-limiting degree, HMs were ranked as follows: Hg, Ce > Pb > Zn > Cr > Cd > Co > Fe > Mn, Cu, Bi, Sr, Ni. Additionally, EDTA could relieve the DNA replication stress induced by some HMs. In conclusion, this study highlights the potential danger of HMs themselves on DNA replication and provides new insight into the possible links between HMs and DNA replication-related diseases.

Substrate Specificity Diversity of Human Terminal Deoxynucleotidyltransferase May Be a Naturally Programmed Feature Facilitating Its Biological Function

Terminal 2'-deoxynucleotidyl transferase (TdT) is a unique enzyme capable of catalysing template-independent elongation of DNA 3' ends during V(D)J recombination. The mechanism controlling the enzyme's substrate specificity, which is necessary for its biological function, remains unknown. Accordingly, in this work, kinetic and mutational analyses of human TdT were performed and allowed to determine quantitative characteristics of individual stages of the enzyme-substrate interaction, which overall may ensure the enzyme's operation either in the distributive or processive mode of primer extension. It was found that conformational dynamics of TdT play an important role in the formation of the catalytic complex. Meanwhile, the nature of the nitrogenous base significantly affected both the dNTP-binding and catalytic-reaction efficiency. The results indicated that neutralisation of the charge and an increase in the internal volume of the active site caused a substantial increase in the activity of the enzyme and induced a transition to the processive mode in the presence of Mg2+ ions. Surrogate metal ions Co2+ or Mn2+ also may regulate the switching of the enzymatic process to the processive mode. Thus, the totality of individual factors affecting the activity of TdT ensures effective execution of its biological function.

Structural Insight into Polymerase Mechanism via a Chiral Center Generated with a Single Selenium Atom

DNA synthesis catalyzed by DNA polymerase is essential for all life forms, and phosphodiester bond formation with phosphorus center inversion is a key step in this process. Herein, by using a single-selenium-atom-modified dNTP probe, we report a novel strategy to visualize the reaction stereochemistry and catalysis. We capture the before- and after-reaction states and provide explicit evidence of the center inversion and in-line attacking SN2 mechanism of DNA polymerization, while solving the diastereomer absolute configurations. Further, our kinetic and thermodynamic studies demonstrate that in the presence of Mg2+ ions (or Mn2+), the binding affinity (Km) and reaction selectivity (kcat/Km) of dGTPαSe-Rp were 51.1-fold (or 19.5-fold) stronger and 21.8-fold (or 11.3-fold) higher than those of dGTPαSe-Sp, respectively, indicating that the diastereomeric Se-Sp atom was quite disruptive of the binding and catalysis. Our findings reveal that the third metal ion is much more critical than the other two metal ions in both substrate recognition and bond formation, providing insights into how to better design the polymerase inhibitors and discover the therapeutics.

Hybrid Quantum Mechanical/Molecular Mechanical Methods For Studying Energy Transduction in Biomolecular Machines

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods have become indispensable tools for the study of biomolecules. In this article, we briefly review the basic methodological details of QM/MM approaches and discuss their applications to various energy transduction problems in biomolecular machines, such as long-range proton transports, fast electron transfers, and mechanochemical coupling. We highlight the particular importance for these applications of balancing computational efficiency and accuracy. Using several recent examples, we illustrate the value and limitations of QM/MM methodologies for both ground and excited states, as well as strategies for calibrating them in specific applications. We conclude with brief comments on several areas that can benefit from further efforts to make QM/MM analyses more quantitative and applicable to increasingly complex biological problems.

In crystallo observation of active site dynamics and transient metal ion binding within DNA polymerases

DNA polymerases are the enzymatic catalysts that synthesize DNA during DNA replication and repair. Kinetic studies and x-ray crystallography have uncovered the overall kinetic pathway and led to a two-metal-ion dependent catalytic mechanism. Diffusion-based time-resolved crystallography has permitted the visualization of the catalytic reaction at atomic resolution and made it possible to capture transient events and metal ion binding that have eluded static polymerase structures. This review discusses past static structures and recent time-resolved structures that emphasize the crucial importance of primer alignment and different metal ions binding during catalysis and substrate discrimination.

Bacterial thermophilic DNA polymerases: A focus on prominent biotechnological applications

DNA polymerases are the enzymes able to replicate the genetic information in nucleic acid. As a result, they are necessary to copy the complete genome of every living creature before cell division and sustain the integrity of the genetic information throughout the life of each cell. Any organism that uses DNA as its genetic information, whether unicellular or multicellular, requires one or more thermostable DNA polymerases to thrive. Thermostable DNA polymerase is important in modern biotechnology and molecular biology because it results in methods such as DNA cloning, DNA sequencing, whole genome amplification, molecular diagnostics, polymerase chain reaction, synthetic biology, and single nucleotide polymorphism detection. There are at least 14 DNA-dependent DNA polymerases in the human genome, which is remarkable. These include the widely accepted, high-fidelity enzymes responsible for replicating the vast majority of genomic DNA and eight or more specialized DNA polymerases discovered in the last decade. The newly discovered polymerases' functions are still being elucidated. Still, one of its crucial tasks is to permit synthesis to resume despite the DNA damage that stops the progression of replication-fork. One of the primary areas of interest in the research field has been the quest for novel DNA polymerase since the unique features of each thermostable DNA polymerase may lead to the prospective creation of novel reagents. Furthermore, protein engineering strategies for generating mutant or artificial DNA polymerases have successfully generated potent DNA polymerases for various applications. In molecular biology, thermostable DNA polymerases are extremely useful for PCR-related methods. This article examines the role and importance of DNA polymerase in a variety of techniques.

Frozen in time: analyzing molecular dynamics with time-resolved cryo-EM

Molecular machines, such as polymerases, ribosomes, or proteasomes, fulfill complex tasks requiring the thermal energy of their environment. They achieve this by restricting random motion along a path of possible conformational changes. These changes are often directed through engagement with different cofactors, which can best be compared to a Brownian ratchet. Many molecular machines undergo three major steps throughout their functional cycles, including initialization, repetitive processing, and termination. Several of these major states have been elucidated by cryogenic electron microscopy (cryo-EM). However, the individual steps for these machines are unique and multistep processes themselves, and their coordination in time is still elusive. To measure these short-lived intermediate events by cryo-EM, the total reaction time needs to be shortened to enrich for the respective pre-equilibrium states. This approach is termed time-resolved cryo-EM (trEM). In this review, we sum up the methodological development of trEM and its application to a range of biological questions.

The metal dependence of single-metal mediated phosphodiester bond cleavage: a QM/MM study of a multifaceted human enzyme

Nucleases catalyze the cleavage of phosphodiester bonds in nucleic acids using a range of metal cofactors. Although it is well accepted that many nucleases rely on two metal ions, the one-metal mediated pathway is debated. Furthermore, one-metal mediated nucleases maintain activity in the presence of many different metals, but the underlying reasons for this broad metal specificity are unknown. The human apurinic/apyrimidinic endonuclease (APE1), which plays a key role in DNA repair, transcription regulation, and gene expression, is a prototypical example of a one-metal dependent nuclease. Although Mg²⁺ is the native metal cofactor, APE1 remains catalytically active in the presence of several metals, with the rate decreasing as Mg²⁺ > Mn²⁺ > Ni²⁺ > Zn²⁺, while Ca²⁺ completely abolished the activity. The present work uses quantum mechanics-molecular mechanics techniques to map APE1-facilitated phosphodiester bond hydrolysis in the presence of these metals. The structural differences in stationary points along the reaction pathway shed light on the interplay between several factors that allow APE1 to remain catalytically active for various metals, with the trend in the barrier heights correlating with the experimentally reported APE1 catalytic activity. In contrast, Ca²⁺ significantly changes the metal coordination and active site geometry, and thus completely inhibits catalysis. Our work thereby provides support for the controversial single-metal mediated phosphodiester bond cleavage and clarifies uncertainties regarding the role of the metal and metal identity in this important reaction. This information is key for future medicinal and biotechnological applications including disease diagnosis and treatment, and protein engineering.

New insights into DNA polymerase mechanisms provided by time-lapse crystallography

DNA polymerases play central roles in DNA replication and repair by catalyzing template-directed nucleotide incorporation. Recently time-lapse X-ray crystallography, which allows one to observe reaction intermediates, has revealed numerous and unexpected mechanistic features of DNA polymerases. In this article, we will examine recent new discoveries that have come from time-lapse crystallography that are currently transforming our understanding of the structural mechanisms used by DNA polymerases. Among these new discoveries are the binding of a third metal ion within the polymerase active site, the mechanisms of translocation along the DNA, the presence of new fidelity checkpoints, a novel pyrophosphatase activity within the active site, and the mechanisms of pyrophosphorolysis.

Kinetic and Mechanistic Studies of the Terminal Uridylyltransferase, Zcchc11 (TUT4)

Zcchc11 (TUT4, TENT3A, Z11) is a nucleotidyltransferase that catalyzes the 3'-polyuridylation of RNA. Our interest in this enzyme stems from its role in blocking the biogenesis of let-7, a family of microRNAs whose members act as tumor suppressors. Z11 polyuridylates pre-let-7, the precursor of let-7, when pre-let-7 is complexed with LIN28, an RNA-binding protein. Polyuridylation of pre-let-7 marks it for degradation. In addition to this LIN28-dependent activity, Z11 also has LIN28-independent activities. In this paper, we report the results of experiments that characterize LIN28-independent activities of Z11. Significant observations include the following. (1) Z11 uridylates not only mature let-7 species but also substrates as small as dinucleotides. (2) For both let-7i and the diribonucleotide AG, Z11 follows a steady-state ordered mechanism, with UTP adding before RNA. (3) Uridylation kinetics of let-7i (UGAGGUAGUAGUUUGUGCUGUU) and two truncated derivatives, GCUGUU and UU, indicate that Z11 manifests selectivity in K_m,RNA; k_cat,RNA values for the three substrates are nearly identical. (4) Z11 preferentially uridylates RNA lacking base-pairing near the 3' terminus. (5) Selectivity of Z11 toward ribonucleoside triphosphates is similar for let-7i and AG, with XTP preference: UTP > CTP > ATP ≫ GTP. Selectivity is manifested in K_m,XTP, with k_cat,XTP values being similar for UTP, CTP, and ATP. (6) Kinetic parameters for RNA turnover are dependent on the structure of the nucleoside triphosphate, consistent with recent structural data indicating stacking of the nucleoside triphosphate base with the base of the 3'-nucleotide of the substrate RNA (Faehnle et al., Nat. Struct. Mol. Biol. 2017, 24, 658).

Insight into the mechanism of DNA synthesis by human terminal deoxynucleotidyltransferase

Terminal deoxynucleotidyltransferase (TdT) is a member of the DNA polymerase X family that is responsible for random addition of nucleotides to single-stranded DNA. We present investigation into the role of metal ions and specific interactions of dNTP with active-site amino acid residues in the mechanisms underlying the recognition of nucleoside triphosphates by human TdT under pre-steady-state conditions. In the elongation mode, the ratios of translocation and dissociation rate constants, as well as the catalytic rate constant were dependent on the nature of the nucleobase. Preferences of TdT in dNTP incorporation were researched by molecular dynamics simulations of complexes of TdT with a primer and dNTP or with the elongated primer. Purine nucleotides lost the "summarised" H-bonding network after the attachment of the nucleotide to the primer, whereas pyrimidine nucleotides increased the number and relative lifetime of H-bonds in the post-catalytic complex. The effect of divalent metal ions on the primer elongation revealed that Me²⁺ cofactor can significantly change parameters of the primer elongation by strongly affecting the rate of nucleotide attachment and the polymerisation mode.

Structural and Molecular Kinetic Features of Activities of DNA Polymerases

DNA polymerases catalyze DNA synthesis during the replication, repair, and recombination of DNA. Based on phylogenetic analysis and primary protein sequences, DNA polymerases have been categorized into seven families: A, B, C, D, X, Y, and RT. This review presents generalized data on the catalytic mechanism of action of DNA polymerases. The structural features of different DNA polymerase families are described in detail. The discussion highlights the kinetics and conformational dynamics of DNA polymerases from all known polymerase families during DNA synthesis.

In crystallo observation of three metal ion promoted DNA polymerase misincorporation

Error-free replication of DNA is essential for life. Despite the proofreading capability of several polymerases, intrinsic polymerase fidelity is in general much higher than what base-pairing energies can provide. Although researchers have investigated this long-standing question with kinetics, structural determination, and computational simulations, the structural factors that dictate polymerase fidelity are not fully resolved. Time-resolved crystallography has elucidated correct nucleotide incorporation and established a three-metal-ion-dependent catalytic mechanism for polymerases. Using X-ray time-resolved crystallography, we visualize the complete DNA misincorporation process catalyzed by DNA polymerase η. The resulting molecular snapshots suggest primer 3´-OH alignment mediated by A-site metal ion binding is the key step in substrate discrimination. Moreover, we observe that C-site metal ion binding preceded the nucleotidyl transfer reaction and demonstrate that the C-site metal ion is strictly required for misincorporation. Our results highlight the essential but separate roles of the three metal ions in DNA synthesis.

By observing DNA polymerase misincorporation with time-resolved crystallography, the authors visualize three-metal ion dependent polymerase catalysis and identify A-site metal-mediated primer alignment as a key step in nucleotide discrimination.

Single-molecule Taq DNA polymerase dynamics

Taq DNA polymerase functions at elevated temperatures with fast conformational dynamics-regimes previously inaccessible to mechanistic, single-molecule studies. Here, single-walled carbon nanotube transistors recorded the motions of Taq molecules processing matched or mismatched template-deoxynucleotide triphosphate pairs from 22° to 85°C. By using four enzyme orientations, the whole-enzyme closures of nucleotide incorporations were distinguished from more rapid, 20-μs closures of Taq's fingers domain testing complementarity and orientation. On average, one transient closure was observed for every nucleotide binding event; even complementary substrate pairs averaged five transient closures between each catalytic incorporation at 72°C. The rate and duration of the transient closures and the catalytic events had almost no temperature dependence, leaving all of Taq's temperature sensitivity to its rate-determining open state.

Two-Metal-Ion Catalysis: Inhibition of DNA Polymerase Activity by a Third Divalent Metal Ion

Almost all DNA polymerases (pols) exhibit bell-shaped activity curves as a function of both pH and Mg²⁺ concentration. The pol activity is reduced when the pH deviates from the optimal value. When the pH is too low the concentration of a deprotonated general base (namely, the attacking 3′-hydroxyl of the 3′ terminal residue of the primer strand) is reduced exponentially. When the pH is too high the concentration of a protonated general acid (i.e., the leaving pyrophosphate group) is reduced. Similarly, the pol activity also decreases when the concentration of the divalent metal ions deviates from its optimal value: when it is too low, the binding of the two catalytic divalent metal ions required for the full activity is incomplete, and when it is too high a third divalent metal ion binds to pyrophosphate, keeping it in the replication complex longer and serving as a substrate for pyrophosphorylysis within the complex. Currently, there is a controversy about the role of the third metal ion which we will address in this review.

Investigating the Conformational Dynamics of a Y-Family DNA Polymerase during Its Folding and Binding to DNA and a Nucleotide

During DNA polymerization, the Y-family DNA polymerases are capable of bypassing various DNA damage, which can stall the replication fork progression. It has been well acknowledged that the structures of the Y-family DNA polymerases have been naturally evolved to undertake this vital task. However, the mechanisms of how these proteins utilize their unique structural and conformational dynamical features to perform the translesion DNA synthesis are less understood. Here, we developed structure-based models to study the precatalytic DNA polymerization process, including DNA and nucleotide binding to DPO4, a paradigmatic Y-family polymerase from Sulfolobus solfataricus. We studied the interplay between the folding and the conformational dynamics of DPO4 and found that DPO4 undergoes first unraveling (unfolding) and then folding for accomplishing the functional "open-to-closed" conformational transition. DNA binding dynamically modulates the conformational equilibrium in DPO4 during the stepwise binding through different types of interactions, leading to different conformational distributions of DPO4 at different DNA binding stages. We observed that nucleotide binding induces modulation of a few contacts surrounding the active site of the DPO4-DNA complex associated with a high free energy barrier. Our simulation results resonate with the experimental evidence that the conformational change at the active site led by nucleotide is the rate-limiting step of nucleotide incorporation. In combination with localized frustration analyses, we underlined the importance of DPO4 conformational dynamics and fluctuations in facilitating DNA and nucleotide binding. Our findings offer mechanistic insights into the processes of DPO4 conformational dynamics associated with the substrate binding and contribute to the understanding of the "structure-dynamics-function" relationship in the Y-family DNA polymerases.

The Role of Natural Polymorphic Variants of DNA Polymerase β in DNA Repair

DNA polymerase β (Polβ) is considered the main repair DNA polymerase involved in the base excision repair (BER) pathway, which plays an important part in the repair of damaged DNA bases usually resulting from alkylation or oxidation. In general, BER involves consecutive actions of DNA glycosylases, AP endonucleases, DNA polymerases, and DNA ligases. It is known that protein-protein interactions of Polβ with enzymes from the BER pathway increase the efficiency of damaged base repair in DNA. However natural single-nucleotide polymorphisms can lead to a substitution of functionally significant amino acid residues and therefore affect the catalytic activity of the enzyme and the accuracy of Polβ action. Up-to-date databases contain information about more than 8000 SNPs in the gene of Polβ. This review summarizes data on the in silico prediction of the effects of Polβ SNPs on DNA repair efficacy; available data on cancers associated with SNPs of Polβ; and experimentally tested variants of Polβ. Analysis of the literature indicates that amino acid substitutions could be important for the maintenance of the native structure of Polβ and contacts with DNA; others affect the catalytic activity of the enzyme or play a part in the precise and correct attachment of the required nucleotide triphosphate. Moreover, the amino acid substitutions in Polβ can disturb interactions with enzymes involved in BER, while the enzymatic activity of the polymorphic variant may not differ significantly from that of the wild-type enzyme. Therefore, investigation regarding the effect of Polβ natural variants occurring in the human population on enzymatic activity and protein-protein interactions is an urgent scientific task.

Beyond the Lesion: Back to High Fidelity DNA Synthesis

High fidelity (HiFi) DNA polymerases (Pols) perform the bulk of DNA synthesis required to duplicate genomes in all forms of life. Their structural features, enzymatic mechanisms, and inherent properties are well-described over several decades of research. HiFi Pols are so accurate that they become stalled at sites of DNA damage or lesions that are not one of the four canonical DNA bases. Once stalled, the replisome becomes compromised and vulnerable to further DNA damage. One mechanism to relieve stalling is to recruit a translesion synthesis (TLS) Pol to rapidly synthesize over and past the damage. These TLS Pols have good specificities for the lesion but are less accurate when synthesizing opposite undamaged DNA, and so, mechanisms are needed to limit TLS Pol synthesis and recruit back a HiFi Pol to reestablish the replisome. The overall TLS process can be complicated with several cellular Pols, multifaceted protein contacts, and variable nucleotide incorporation kinetics all contributing to several discrete substitution (or template hand-off) steps. In this review, we highlight the mechanistic differences between distributive equilibrium exchange events and concerted contact-dependent switching by DNA Pols for insertion, extension, and resumption of high-fidelity synthesis beyond the lesion.

Direct Measurement of Initial Rate of Enzyme Reaction by Electrokinetic Filtration Using a Hydrogel-plugged Capillary Device

A novel electrokinetic filtration device using a plugged hydrogel was developed to directly measure the initial rate of enzyme reactions. In the proposed method, the enzyme reaction proceeded only for a short time when the substrate was passed through a thin layer of enzyme trapped by the hydrogel without any lag times for mixing and detection. In experimental conditions, alkaline phosphatase (enzyme) was filtrated at a cathodic-side interface of the plugged hydrogel by molecular sieving effect, providing the thin enzyme zone whose thickness was approximately 100 μm. When 4-methylumberiferyl phosphate (substrate) was electrokinetically introduced into the device after trapping the enzyme, 4-methylumberiferone (product) was generated by the enzyme reaction for only 1.26 s as the substrate passed through the trapped enzyme zone. As a result, the initial rate of the enzyme reaction could be directly calculated to 31.0 μM/s by simply dividing the concentration of the product by the tunable reaction time. Compared to the initial rate obtained by mixing the enzyme and substrate solutions, the value of the maximum velocity of the enzyme reaction was 30-fold larger than that in the mixing method due to the preconcentration of the enzyme by trapping. The Michaelis-Menten constant in the proposed method was 2.7-fold larger than that in the mixing method, suggesting the variation of changes in the equilibrium of complex formation under the experimental conditions.

Rational Design of Nonbonded Point Charge Models for Divalent Metal Cations with Lennard-Jones 12-6 Potential

Exploring a metal-involved biochemical process at a molecular level often requires a reliable description of metal properties in aqueous solution by classical nonbonded models. An additional C₄ term for considering ion-induced dipole interactions was previously proposed to supplement the widely used Lennard-Jones 12-6 potential (known as the 12-6-4 LJ-type model) with good accuracy. Here, we demonstrate an alternative to modeling divalent metal cations (M²⁺) with the traditional 12-6 LJ potential by developing nonbonded point charge models for use with 11 water models: TIP3P, SPC/E, SPC/E_b, TIP4P-Ew, TIP4P-D, and TIP4P/2005 and the more recent OPC3, TIP3P-FB, OPC, TIP4P-FB, and a99SB-disp. Our designed models simultaneously reproduce the experimental hydration free energy, ion-oxygen distance, and coordination number in the first hydration shell accurately for most of the metal cations, an accuracy equivalent to that of the complex 12-6-4 LJ-type and double exponential potential models. A systematic comparison with the existing M²⁺ models is presented as well in terms of effective ion radii, diffusion constants, water exchange rates, and ion-water interactions. Molecular dynamics simulations of metal substitution in Escherichia coli glyoxalase I variants show the great potential of our new models for metalloproteins.

Selenium atom on phosphate enhances specificity and sensitivity of DNA polymerization and detection

DNA polymerization is of high specificity in vivo. However, its specificity is much lower in vitro, which limits advanced applications of DNA polymerization in ultrasensitive nucleic acid detection. Herein, we report a unique mechanism of single selenium-atom modified dNTP (dNTPαSe) to enhance polymerization specificity. We have found that both dNTPαSe (approximately 660 fold) and Se-DNA (approximately 2.8 fold) have lower binding affinity to DNA polymerase than canonical ones, and the Se-DNA duplex has much lower melting-temperature (Tm) than the corresponding canonical DNA duplex. The reduced affinity and Tm can destabilize the substrate-primer-template-enzyme assembly, thereby largely slowing down the mismatch of DNA polymerization and enhancing the amplification specificity and in turn detection sensitivity. Furthermore, the Se-strategy enables us to develop the selenium enhanced specific isothermal amplification (SEA) for nucleic acid detection with high specificity and sensitivity (up to detection of single-digit copies), allowing convenient detection of clinical HPV and COVID-19 viruses in the low-copy number. Clearly, we have discovered the exciting mechanism for enhancing DNA polymerization accuracy, amplification specificity and detection sensitivity by SEA, up to two orders of magnitude higher.

Pyrophosphate release acts as a kinetic checkpoint during high-fidelity DNA replication by the Staphylococcus aureus replicative polymerase PolC

Bacterial replication is a fast and accurate process, with the bulk of genome duplication being catalyzed by the α subunit of DNA polymerase III within the bacterial replisome. Structural and biochemical studies have elucidated the overall properties of these polymerases, including how they interact with other components of the replisome, but have only begun to define the enzymatic mechanism of nucleotide incorporation. Using transient-state methods, we have determined the kinetic mechanism of accurate replication by PolC, the replicative polymerase from the Gram-positive pathogen Staphylococcus aureus. Remarkably, PolC can recognize the presence of the next correct nucleotide prior to completing the addition of the current nucleotide. By modulating the rate of pyrophosphate byproduct release, PolC can tune the speed of DNA synthesis in response to the concentration of the next incoming nucleotide. The kinetic mechanism described here would allow PolC to perform high fidelity replication in response to diverse cellular environments.

Kinetic investigation of the polymerase and exonuclease activities of human DNA polymerase ε holoenzyme

In eukaryotic DNA replication, DNA polymerase ε (Polε) is responsible for leading strand synthesis, whereas DNA polymerases α and δ synthesize the lagging strand. The human Polε (hPolε) holoenzyme is comprised of the catalytic p261 subunit and the noncatalytic p59, p17, and p12 small subunits. So far, the contribution of the noncatalytic subunits to hPolε function is not well understood. Using pre-steady-state kinetic methods, we established a minimal kinetic mechanism for DNA polymerization and editing catalyzed by the hPolε holoenzyme. Compared with the 140-kDa N-terminal catalytic fragment of p261 (p261N), which we kinetically characterized in our earlier studies, the presence of the p261 C-terminal domain (p261C) and the three small subunits increased the DNA binding affinity and the base substitution fidelity. Although the small subunits enhanced correct nucleotide incorporation efficiency, there was a wide range of rate constants when incorporating a correct nucleotide over a single-base mismatch. Surprisingly, the 3'→5' exonuclease activity of the hPolε holoenzyme was significantly slower than that of p261N when editing both matched and mismatched DNA substrates. This suggests that the presence of p261C and the three small subunits regulates the 3'→5' exonuclease activity of the hPolε holoenzyme. Together, the 3'→5' exonuclease activity and the variable mismatch extension activity modulate the overall fidelity of the hPolε holoenzyme by up to 3 orders of magnitude. Thus, the presence of p261C and the three noncatalytic subunits optimizes the dual enzymatic activities of the catalytic p261 subunit and makes the hPolε holoenzyme an efficient and faithful replicative DNA polymerase.

Building better polymerases: Engineering the replication of expanded genetic alphabets

DNA polymerases are today used throughout scientific research, biotechnology, and medicine, in part for their ability to interact with unnatural forms of DNA created by synthetic biologists. Here especially, natural DNA polymerases often do not have the "performance specifications" needed for transformative technologies. This creates a need for science-guided rational (or semi-rational) engineering to identify variants that replicate unnatural base pairs (UBPs), unnatural backbones, tags, or other evolutionarily novel features of unnatural DNA. In this review, we provide a brief overview of the chemistry and properties of replicative DNA polymerases and their evolved variants, focusing on the Klenow fragment of Taq DNA polymerase (Klentaq). We describe comparative structural, enzymatic, and molecular dynamics studies of WT and Klentaq variants, complexed with natural or noncanonical substrates. Combining these methods provides insight into how specific amino acid substitutions distant from the active site in a Klentaq DNA polymerase variant (ZP Klentaq) contribute to its ability to replicate UBPs with improved efficiency compared with Klentaq. This approach can therefore serve to guide any future rational engineering of replicative DNA polymerases.

Investigating the trade-off between folding and function in a multidomain Y-family DNA polymerase

The way in which multidomain proteins fold has been a puzzling question for decades. Until now, the mechanisms and functions of domain interactions involved in multidomain protein folding have been obscure. Here, we develop structure-based models to investigate the folding and DNA-binding processes of the multidomain Y-family DNA polymerase IV (DPO4). We uncover shifts in the folding mechanism among ordered domain-wise folding, backtracking folding, and cooperative folding, modulated by interdomain interactions. These lead to ‘U-shaped’ DPO4 folding kinetics. We characterize the effects of interdomain flexibility on the promotion of DPO4–DNA (un)binding, which probably contributes to the ability of DPO4 to bypass DNA lesions, which is a known biological role of Y-family polymerases. We suggest that the native topology of DPO4 leads to a trade-off between fast, stable folding and tight functional DNA binding. Our approach provides an effective way to quantitatively correlate the roles of protein interactions in conformational dynamics at the multidomain level.

Environmental Effects on Guanine-Thymine Mispair Tautomerization Explored with Quantum Mechanical/Molecular Mechanical Free Energy Simulations

DNA bases can adopt energetically unfavorable tautomeric forms that enable the formation of Watson-Crick-like (WC-like) mispairs, which have been proposed to give rise to spontaneous mutations in DNA and misincorporation errors in DNA replication and translation. Previous NMR and computational studies have indicated that the population of WC-like guanine-thymine (G-T) mispairs depends on the environment, such as the local nucleic acid sequence and solvation. To investigate these environmental effects, herein G-T mispair tautomerization processes are studied computationally in aqueous solution, in A-form and B-form DNA duplexes, and within the active site of a DNA polymerase λ variant. The wobble G-T (wG-T), WC-like G-T*, and WC-like G*-T forms are considered, where * indicates the enol tautomer of the base. The minimum free energy paths for the tautomerization from the wG-T to the WC-like G-T* and from the WC-like G-T* to the WC-like G*-T are computed with mixed quantum mechanical/molecular mechanical (QM/MM) free energy simulations. The reaction free energies and free energy barriers are found to be significantly influenced by the environment. The wG-T→G-T* tautomerization is predicted to be endoergic in aqueous solution and the DNA duplexes but slightly exoergic in the polymerase, with Arg517 and Asn513 providing electrostatic stabilization of G-T*. The G-T*→G*-T tautomerization is also predicted to be slightly more thermodynamically favorable in the polymerase relative to these DNA duplexes. These simulations are consistent with an experimentally driven kinetic misincorporation model suggesting that G-T mispair tautomerization occurs in the ajar polymerase conformation or concertedly with the transition from the ajar to the closed polymerase conformation. Furthermore, the order of the associated two proton transfer reactions is predicted to be different in the polymerase than in aqueous solution and the DNA duplexes. These studies highlight the impact of the environment on the thermodynamics, kinetics, and fundamental mechanisms of G-T mispair tautomerization, which plays a role in a wide range of biochemically important processes.

Catalytically inactive T7 DNA polymerase imposes a lethal replication roadblock

Bacteriophage T7 encodes its own DNA polymerase, the product of gene 5 (gp5). In isolation, gp5 is a DNA polymerase of low processivity. However, gp5 becomes highly processive upon formation of a complex with Escherichia coli thioredoxin, the product of the trxA gene. Expression of a gp5 variant in which aspartate residues in the metal-binding site of the polymerase domain were replaced by alanine is highly toxic to E. coli cells. This toxicity depends on the presence of a functional E. coli trxA allele and T7 RNA polymerase-driven expression but is independent of the exonuclease activity of gp5. In vitro, the purified gp5 variant is devoid of any detectable polymerase activity and inhibited DNA synthesis by the replisomes of E. coli and T7 in the presence of thioredoxin by forming a stable complex with DNA that prevents replication. On the other hand, the highly homologous Klenow fragment of DNA polymerase I containing an engineered gp5 thioredoxin-binding domain did not exhibit toxicity. We conclude that gp5 alleles encoding inactive polymerases, in combination with thioredoxin, could be useful as a shutoff mechanism in the design of a bacterial cell-growth system.

Salt effect on thermodynamics and kinetics of a single RNA base pair

Due to the polyanionic nature of RNAs, the structural folding of RNAs are sensitive to solution salt conditions, while there is still lack of a deep understanding of the salt effect on the thermodynamics and kinetics of RNAs at a single base-pair level. In this work, the thermodynamic and the kinetic parameters for the base-pair AU closing/opening at different salt concentrations were calculated by 3-µsec all-atom molecular dynamics (MD) simulations at different temperatures. It was found that for the base-pair formation, the enthalpy change [Formula: see text] is nearly independent of salt concentration, while the entropy change [Formula: see text] exhibits a linear dependence on the logarithm of salt concentration, verifying the empirical assumption based on thermodynamic experiments. Our analyses revealed that such salt concentration dependence of the entropy change mainly results from the dependence of ion translational entropy change for the base pair closing/opening on salt concentration. Furthermore, the closing rate increases with the increasing of salt concentration, while the opening rate is nearly independent of salt concentration. Additionally, our analyses revealed that the free energy surface for describing the base-pair opening and closing dynamics becomes more rugged with the decrease of salt concentration.

Cations in motion: QM/MM studies of the dynamic and electrostatic roles of H+ and Mg2+ ions in enzyme reactions

Here we discuss current trends in the simulations of enzymatic reactions focusing on phosphate catalysis. The mechanistic details of the proton transfers coupled to the phosphate cleavage is one of the key challenges in QM/MM calculations of these and other enzyme catalyzed reactions. The lack of experimental information offers both an opportunity for computations as well as often unresolved controversies. We discuss the example of small GTPases including the important human Ras protein. The high dimensionality and chemical complexity of these reactions demand carefully chosen computational techniques both in terms of the underlying quantum chemical theory and the sampling of the conformational ensemble. We also point out the important role of Mg²⁺ ions, and recent advances in their transient involvement in the catalytic mechanisms.

Confinement and Crowding Effects on Folding of a Multidomain Y-Family DNA Polymerase

Proteins in vivo endure highly various interactions from the luxuriant surrounding macromolecular cosolutes. Confinement and macromolecular crowding are the two major effects that should be considered while comparing the results of protein dynamics from in vitro to in vivo. However, efforts have been largely focused on single domain protein folding up to now, and the quantifications of the in vivo effects in terms of confinements and crowders on modulating the structure and dynamics as well as the physical understanding of the underlying mechanisms on multidomain protein folding are still challenging. Here we developed a topology-based model to investigate folding of a multidomain Y-family DNA polymerase (DPO4) within spherical confined space and in the presence of repulsive and attractive crowders. We uncovered that the entropic component of the thermodynamic driving force led by confinements and repulsive crowders increases the stability of folded states relative to the folding intermediates and unfolded states, while the enthalpic component of the thermodynamic driving force led by attractive crowders gives rise to the opposite effects with less stability. We found that the shapes of DPO4 conformations influenced by the confinements and the crowders are quite different even when only the entropic component of the thermodynamic driving force is considered. We uncovered that under all in vivo conditions, the folding cooperativity of DPO4 decreases compared to that in bulk. We showed that the loss of folding cooperativity can promote the sequential domain-wise folding, which was widely found in cotranslational multidomain protein folding, and effectively prohibit the backtracking led by topological frustrations during multidomain protein folding processes.

Recruiting Mechanism and Functional Role of a Third Metal Ion in the Enzymatic Activity of 5' Structure-Specific Nucleases

Enzymes of the 5′ structure-specific nuclease family are crucial for DNA repair, replication, and recombination. One such enzyme is the human exonuclease 1 (hExo1) metalloenzyme, which cleaves DNA strands, acting primarily as a processive 5′-3′ exonuclease and secondarily as a 5′-flap endonuclease. Recently, in crystallo reaction intermediates have elucidated how hExo1 exerts hydrolysis of DNA phosphodiester bonds. These hExo1 structures show a third metal ion intermittently bound close to the two-metal-ion active site, to which recessed ends or 5′-flap substrates bind. Evidence of this third ion has been observed in several nucleic-acid-processing metalloenzymes. However, there is still debate over what triggers the (un)binding of this transient third ion during catalysis and whether this ion has a catalytic function. Using extended molecular dynamics and enhanced sampling free-energy simulations, we observed that the carboxyl side chain of Glu89 (located along the arch motif in hExo1) flips frequently from the reactant state to the product state. The conformational flipping of Glu89 allows one metal ion to be recruited from the bulk and promptly positioned near the catalytic center. This is in line with the structural evidence. Additionally, our simulations show that the third metal ion assists the departure, through the mobile arch, of the nucleotide monophosphate product from the catalytic site. Structural comparisons of nuclease enzymes suggest that this Glu(Asp)-mediated mechanism for third ion recruitment and nucleic acid hydrolysis may be shared by other 5′ structure-specific nucleases.

DNA polymerase activity on synthetic N3'→P5' phosphoramidate DNA templates

Genetic polymers that could plausibly govern life in the universe might inhabit a broad swath of chemical space. A subset of these genetic systems can exchange information with RNA and DNA and could therefore form the basis for model protocells in the laboratory. N3'→P5' phosphoramidate (NP) DNA is defined by a conservative linkage substitution and has shown promise as a protocellular genetic material, but much remains unknown about its functionality and fidelity due to limited enzymatic tools. Conveniently, we find widespread NP-DNA-dependent DNA polymerase activity among reverse transcriptases, an observation consistent with structural studies of the RNA-like conformation of NP-DNA duplexes. Here, we analyze the consequences of this unnatural template linkage on the kinetics and fidelity of DNA polymerization activity catalyzed by wild-type and variant reverse transcriptases. Template-associated deficits in kinetics and fidelity suggest that even highly conservative template modifications give rise to error-prone DNA polymerase activity. Enzymatic copying of NP-DNA sequences is nevertheless an important step toward the future study and engineering of this synthetic genetic polymer.

Extensive free-energy simulations identify water as the base in nucleotide addition by DNA polymerase

Transphosphorylation of nucleotide triphosphates is the central reaction in DNA replication by DNA polymerase as well as many other biological processes. Despite its importance, the microscopic chemical mechanism of transphosphorylation of nucleotide triphosphates is, in most cases, unknown. Here we use extensive simulations of DNA polymerase η to test mechanistic hypotheses. We systematically survey the reactive space by calculating 2D free-energy surfaces for 10 different plausible mechanisms that have been proposed. We supplement these free-energy surfaces with calculations of pKa for a number of potentially acidic protons in different states relevant to the catalytic cycle. We find that among all of the conditions that we test, the smallest activation barrier occurs for a reaction where a Mg2+-coordinated water deprotonates the nucleophilic 3'-OH, and this deprotonation is concerted with the phosphoryl transfer. The presence of a third Mg2+ in the active site lowers the activation barrier for the water-as-base mechanism, as does protonation of the pyrophosphate leaving group, which is consistent with general acid catalysis. The results demonstrate the value of simulations, when used in conjunction with experimental data, to help establish a microscopic chemical mechanism in a complex environment.