Steric and electrostatic effects at the C2 atom substituent influence replication and miscoding of the DNA deamination product deoxyxanthosine and analogs by DNA polymerases

A human high-fidelity DNA polymerase holoenzyme has a wide range of lesion bypass activities

During replication, lagging strand lesions are initially encountered by high-fidelity DNA polymerase (pol) holoenzymes comprised of pol δ and the PCNA sliding clamp. To proceed unhindered, pol δ holoenzymes must bypass lesions without stalling. This entails dNMP incorporation opposite the lesion (insertion) and the 5' template nucleotide (extension). Historically, it was viewed that high-fidelity pol holoenzymes stall upon encountering lesions, activating DNA damage tolerance pathways that are ultimately responsible for lesion bypass. Our recent study of 4 prominent lesions revealed that human pol δ holoenzymes support insertion and/or bypass for multiple lesions and the extents of these activities depends on the lesion and pol δ proofreading. In the present study, we expand these analyses to other prominent lesions. Collectively, analyses of 10 lesions from both studies reveal that the insertion and bypass efficiencies of pol δ holoenzymes each span a complete range (0 - 100%). Consequently, the fates of pol δ holoenzymes upon encountering lesions are quite diverse. Furthermore, pol δ proofreading promoted holoenzyme progression at 7 of the 10 lesions and did not deter progression at any. Altogether, the results significantly alter our understanding of the replicative capacity of high-fidelity pol holoenzymes and their functional role(s) in lesion bypass.

Identifying Metal Binding Sites in Proteins Using Homologous Structures, the MADE Approach

In order to identify the locations of metal ions in the binding sites of proteins, we have developed a method named the MADE (MAcromolecular DEnsity and Structure Analysis) approach. The MADE approach represents an evolution of our previous toolset, the ProBiS H2O (MD) methodology, for the identification of conserved water molecules. Our method uses experimental structures of proteins homologous to a query, which are subsequently superimposed upon it. Areas with a particular species present in a similar location among many homologous protein structures are identified using a clustering algorithm. Dense clusters likely represent positions containing species important to the query protein structure or function. We analyze well-characterized apo protein structures and show that the MADE approach can identify clusters corresponding to the expected positions of metal ions in their binding sites. The greatest advantage of our method lies in its generality. It can in principle be applied to any species found in protein records; it is not only limited to metal ions. We additionally demonstrate that the MADE approach can be successfully applied to predict the location of cofactors in computer-modeled structures, e.g., via AlphaFold. We also conduct a careful protein superposition method comparison and find our methodology robust and the results largely independent of the selected protein superposition algorithm. We postulate that with increasing structural data availability, additional applications of the MADE approach will be possible such as non-protein systems, water network identification, protein binding site elaboration, and analysis of binding events, all in a dynamic manner. We have implemented the MADE approach as a plugin for the PyMOL molecular visualization tool. The MADE plugin is available free of charge at https://gitlab.com/Jukic/made_software.

Strand displacement DNA synthesis by DNA polymerase gp90 exo- of Pseudomonas aeruginosa phage 1

Strand displacement DNA synthesis is essential for DNA replication. Gp90, the sole DNA polymerase of Pseudomonas aeruginosa phage 1, can bypass multiply DNA lesions. However, whether it can perform strand displacement synthesis is still unknown. In this work, we found that gp90 exo^- could perform strand displacement synthesis, albeit its activity and processivity were lower than those of primer extension. Gp90 exo^- itself could not unwind Y-shaped or fork DNA. Tail and gap at DNA fork were necessary for efficient synthesis. High GC content obviously inhibited strand displacement synthesis. Consecutive GC sequence at the entrance of fork showed more inhibition effect on DNA synthesis than that in the downstream DNA fork. The fraction of productive polymerase and DNA complex (A values) was higher for fork than gap; while their average extension rates (k_p values) were similar. However, both A and k_p values were lower than those for the primer/template (P/T) substrate. The binding of gp90 exo^- to fork was tighter than P/T or gap in the absence of dATP. In the presence of dATP to form ternary complex, the binding affinity of gp90 exo^- to P/T or gap was increased compared with that in the binary complex. Abasic site, 8-oxoG, and O⁶-MeG inhibited and even blocked strand displacement synthesis. This work shows that gp90 exo^- could perform strand displacement DNA synthesis at DNA fork, discovering the presence of new functions of PaP1 DNA polymerase in DNA replication and propagation of PaP1.

Ribonucleoside triphosphates promote T7 DNA replication and the lysis of T7-Infected Escherichia coli

rNTPs are structurally similar to dNTPs, but their concentrations are much higher than those of dNTPs in cells. rNTPs in solutions or rNMP at the primer terminus or embedded in template always inhibit or block DNA replication, due to the reduced Mg²⁺ apparent concentration, competition of rNTPs with dNTPs, and the extra repulsive interaction of rNTP or rNMP with polymerase active site. In this work, unexpectedly, we found rNTPs can promote T7 DNA replication with the maximal promotion at rNTPs/dNTPs concentration ratio of 20. This promotion was not due to the optimized Mg²⁺ apparent concentration or the direct incorporation of extra rNMPs into DNA. This promotion was dependent on the concentrations and types of rNTPs. Kinetic analysis showed that this promotion was originated from the increased fraction of polymerase-DNA productive complex and the accelerated DNA polymerization. Further evidence showed that more polymerase-DNA complex was formed and their binding affinity was also enhanced in the presence of extra rNTPs. Moreover, this promotion in T7 DNA replication also accelerated the lysis of T7-infected host Escherichia coli. This work discovered that rNTPs could promote DNA replication, completely different from the traditional concept that rNTPs always inhibit DNA replication.

Epigenetic DNA Modification N6-Methyladenine Inhibits DNA Replication by DNA Polymerase of Pseudomonas aeruginosa Phage PaP1

N⁶-methyladenine (6mA), a newly identified epigenetic modification, plays important roles in regulation of various biological processes. However, the effect of 6mA on DNA replication has been little addressed. In this work, we investigated how 6mA affected DNA replication by DNA polymerase of Pseudomonas aeruginosa Phage PaP1 (gp90 exo^-). The presence of 6mA, as well as its intermediate hypoxanthine (Hyp), inhibited DNA replication by gp90 exo^-. The 6mA reduced dTTP incorporation efficiency by 10-fold and inhibited next-base extension efficiency by 100-fold. Differently, dCTP was preferentially incorporated opposite Hyp among four dNTPs. Gp90 exo^- reduced the extension priority beyond the 6mA:T pair rather than the 6mA:C mispair and preferred to extend beyond Hyp:C rather than the Hyp:T pair. Incorporation of dTTP opposite 6mA and dCTP opposite Hyp showed fast burst phases. The burst rate and burst amplitude were both reduced for 6mA compared with unmodified A. Moreover, the total incorporation efficiency ( k_pol/ K_d,dNTP) was decreased for dTTP incorporation opposite 6mA and dCTP incorporation opposite Hyp compared with dTTP incorporation opposite A. 6mA reduced the incorporation rate ( k_pol), and Hyp increased the dissociation constant ( K_d,dNTP). However, 6mA or Hyp on template did not affect the binding of DNA polymerase to DNA in binary or ternary complexes. This work provides new insight into the inhibited effects of epigenetic modification of 6mA on DNA replication in PaP1.

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

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

Error-prone bypass of O6-methylguanine by DNA polymerase of Pseudomonas aeruginosa phage PaP1

O⁶-Methylguanine (O⁶-MeG) is highly mutagenic and is commonly found in DNA exposed to methylating agents, generally leads to G:C to A:T mutagenesis. To study DNA replication encountering O⁶-MeG by the DNA polymerase (gp90) of P. aeruginosa phage PaP1, we analyzed steady-state and pre-steady-state kinetics of nucleotide incorporation opposite O⁶-MeG by gp90 exo^-. O⁶-MeG partially inhibited full-length extension by gp90 exo^-. O⁶-MeG greatly reduces dNTP incorporation efficiency, resulting in 67-fold preferential error-prone incorporation of dTTP than dCTP. Gp90 exo^- extends beyond T:O⁶-MeG 2-fold more efficiently than C:O⁶-MeG. Incorporation of dCTP opposite G and incorporation of dCTP or dTTP opposite O⁶-MeG show fast burst phases. The pre-steady-state incorporation efficiency (k_pol/K_d,dNTP) is decreased in the order of dCTP:G>dTTP:O⁶-MeG>dCTP:O⁶-MeG. The presence of O⁶-MeG at template does not affect the binding affinity of polymerase to DNA but it weakened their binding in the presence of dCTP and Mg²⁺. Misincorporation of dTTP opposite O⁶-MeG further weakens the binding affinity of polymerase to DNA. The priority of dTTP incorporation opposite O⁶-MeG is originated from the fact that dTTP can induce a faster conformational change step and a faster chemical step than dCTP. This study reveals that gp90 bypasses O⁶-MeG in an error-prone manner and provides further understanding in DNA replication encountering mutagenic alkylation DNA damage for P. aeruginosa phage PaP1.

Error-Free Bypass of 7,8-dihydro-8-oxo-2'-deoxyguanosineby DNA Polymerase of Pseudomonas aeruginosa Phage PaP1

As one of the most common forms of oxidative DNA damage, 7,8-dihydro-8-oxo-2′-deoxyguanosine (8-oxoG) generally leads to G:C to T:A mutagenesis. To study DNA replication encountering 8-oxoG by the sole DNA polymerase (Gp90) of Pseudomonasaeruginosa phage PaP1, we performed steady-state and pre-steady-state kinetic analyses of nucleotide incorporation opposite 8-oxoG by Gp90 D234A that lacks exonuclease activities on ssDNA and dsDNA substrates. Gp90 D234A could bypass 8-oxoG in an error-free manner, preferentially incorporate dCTP opposite 8-oxoG, and yield similar misincorporation frequency to unmodified G. Gp90 D234A could extend beyond C:8-oxoG or A:8-oxoG base pairs with the same efficiency. dCTP incorporation opposite G and dCTP or dATP incorporation opposite 8-oxoG showed fast burst phases. The burst of incorporation efficiency (k_pol/K_d,_dNTP) is decreased as dCTP:G > dCTP:8-oxoG > dATP:8-oxoG. The presence of 8-oxoG in DNA does not affect its binding to Gp90 D234A in a binary complex but it does affect it in a ternary complex with dNTP and Mg²⁺, and dATP misincorporation opposite 8-oxoG further weakens the binding of Gp90 D234A to DNA. This study reveals Gp90 D234A can bypass 8-oxoG in an error-free manner, providing further understanding in DNA replication encountering oxidation lesion for P.aeruginosa phage PaP1.

Role of magnesium ions in the reaction mechanism at the interface between Tm1631 protein and its DNA ligand

A protein, Tm1631 from the hyperthermophilic organism Thermotoga maritima belongs to a domain of unknown function protein family. It was predicted that Tm1631 binds with the DNA and that the Tm1631–DNA complex is an endonuclease repair system with a DNA repair function (Konc et al. PLoS Comput Biol 9(11): e1003341, 2013). We observed that the severely bent, strained DNA binds to the protein for the entire 90 ns of classical molecular dynamics (MD) performed; we could observe no significant changes in the most distorted region of the DNA, where the cleavage of phosphodiester bond occurs. In this article, we modeled the reaction mechanism at the interface between Tm1631 and its proposed ligand, the DNA molecule, focusing on cleavage of the phosphodiester bond. After addition of two Mg²⁺ ions to the reaction center and extension of classical MD by 50 ns (totaling 140 ns), the DNA ligand stayed bolted to the protein. Results from density functional theory quantum mechanics/molecular mechanics (QM/MM) calculations suggest that the reaction is analogous to known endonuclease mechanisms: an enzyme reaction mechanism with two Mg²⁺ ions in the reaction center and a pentacovalent intermediate. The minimum energy pathway profile shows that the phosphodiester bond cleavage step of the reaction is kinetically controlled and not thermodynamically because of a lack of any energy barrier above the accuracy of the energy profile calculation. The role of ions is shown by comparing the results with the reaction mechanisms in the absence of the Mg²⁺ ions where there is a significantly higher reaction barrier than in the presence of the Mg²⁺ ions.Graphical abstract

Mechanisms of mutagenesis: DNA replication in the presence of DNA damage

Environmental mutagens cause DNA damage that disturbs replication and produces mutations, leading to cancer and other diseases. We discuss mechanisms of mutagenesis resulting from DNA damage, from the level of DNA replication by a single polymerase to the complex DNA replisome of some typical model organisms (including bacteriophage T7, T4, Sulfolobus solfataricus, Escherichia coli, yeast and human). For a single DNA polymerase, DNA damage can affect replication in three major ways: reducing replication fidelity, causing frameshift mutations, and blocking replication. For the DNA replisome, protein interactions and the functions of accessory proteins can yield rather different results even with a single DNA polymerase. The mechanism of mutation during replication performed by the DNA replisome is a long-standing question. Using new methods and techniques, the replisomes of certain organisms and human cell extracts can now be investigated with regard to the bypass of DNA damage. In this review, we consider the molecular mechanism of mutagenesis resulting from DNA damage in replication at the levels of single DNA polymerases and complex DNA replisomes, including translesion DNA synthesis.

Kinetic analysis of bypass of O(6)- methylguanine by the catalytic core of yeast DNA polymerase eta

Alkylating agents can form O(6)-methylguansine (O(6)-MeG). To study the intrinsic kinetic behaviors of bypassing O(6)-MeG, we used the catalytic core of yeast DNA polymerase η (Pol ηcore, residues 1-513), instead of the full-length Pol η, to study their elementary steps, eliminating the effects of the C-terminal C2H2 motif on dNTP incorporation. The misincorporation frequencies were 10(-4) for G and 0.055-0.446 for O(6)-MeG. O(6)-MeG does not affect the extension efficiency. Pol ηcore showed no fast burst phase for any incorporation opposite G or O(6)-MeG. Primer extension was greatly blocked by O(6)-MeG and about 67% dTTP, 31% dCTP and 2% dATP were incorporated opposite O(6)-MeG. This study provides further understanding of the mutation mechanism of alkylated lesion for yeast DNA polymerase η.

Synthetic nucleotides as probes of DNA polymerase specificity

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

Study of Tamiflu sensitivity to variants of A/H5N1 virus using different force fields

An accurate estimation of binding free energy of a ligand to receptor ΔG(bind) is one of the most important problems in drug design. The success of solution of this problem is expected to depend on force fields used for modeling a ligand-receptor complex. In this paper, we consider the impact of four main force fields, AMBER99SB, CHARMM27, GROMOS96 43a1, and OPLS-AA/L, on the binding affinity of Oseltamivir carboxylate to the wild-type and Y252H, N294S, and H274Y mutants of glycoprotein neuraminidase from the pandemic A/H5N1 virus. Having used the molecular mechanic-Poisson-Boltzmann surface area method, we have shown that ΔG(bind), obtained by AMBER99SB, OPLS-AA/L, and CHARMM27, shows the high correlation with the available experimental data. They correctly capture the binding ranking Y252H → WT → N294S → H274Y observed in experiments (Collins, P. J. et al. Nature 2008, 453, 1258). In terms of absolute values of binding scores, results obtained by AMBER99SB are in the nearest range with experiments, while OPLS-AA/L, which is applied to study binding of Oseltamivir to the influenza virus for the first time, gives rather big negative values for ΔG(bind). GROMOS96 43a1 provides a lower correlation as it supports Oseltamivir to be more resistant to N294S than H274Y. Our study suggests that force fields have pronounced influence on theoretical estimations of binding free energy of a ligand to receptor. The effect of all-atom models on dynamics of the binding pocket as well as on the hydrogen-bond network between Oseltamivir and receptors is studied in detail. The hydrogen network, obtained by GROMOS, is weakest among four studied force fields.

Microwave catalysis revisited: an analytical solution

In our previous work [Bren, U., et al. J. Phys. Chem. A 2008, 112, 166] we proposed a novel physical mechanism for microwave catalysis based on rotationally hot reactive species and verified its validity through a Monte Carlo simulation of a realistic chemical reaction: neutral ester hydrolysis. This article represents a continuation of our ongoing effort toward quantitative understanding of the microwave catalytic effect. It provides a derivation of an analytical solution for the microwave catalysis. The obtained expression is compared with the results of the Monte Carlo simulation and is applied to reproduce the microwave catalytic effect experimentally observed in the polyethylene terephthalate solvolysis. Implications for the interactions of microwaves with living organisms in the context of widespread mobile telephony are also discussed.

Efficient recognition of an unpaired lesion by a DNA repair glycosylase

The duplex structure of DNA, with its internal base pairing, protects the nucleobases from chemical damage, but it also poses a barrier to DNA-modifying enzymes, including the enzymes that recognize and repair DNA damage. It is known that unpaired (or bulged) nucleotides are significantly more accessible, but it is not known whether they might be recognized by nucleotide-flipping enzymes. We have investigated this question with human alkyladenine DNA glycosylase (AAG). AAG recognizes a wide variety of structurally disparate lesions, including deoxyinosine (I), which results from the spontaneous oxidative deamination of adenosine, and catalyzes the hydrolysis of the N-glycosidic bond to release the lesion base and initiate the base excision repair pathway. We used single-turnover kinetics to characterize the reactions of AAG with synthetic 25-mer oligonucleotides containing a single I lesion in single-stranded, mismatched, or single-nucleotide bulge contexts. We found that AAG has the highest catalytic efficiency toward a lesion that is presented in a single-nucleotide bulge. In contrast, AAG has more than 2000-fold reduced catalytic efficiency toward a single-stranded I-containing oligonucleotide relative to the duplexes. We have observed 20-fold differences in catalytic efficiency for the excision of the presumed biological target (paired with T) relative to alternative pairings such as C that might be formed by the replication of an unrepaired I. Furthermore, a linear free-energy relationship shows a strong inverse correlation between duplex stability and catalytic efficiency (slope = -0.6 to -1.0), indicating that gaining access to the base lesion provides a substantial barrier to AAG-catalyzed initiation of DNA repair. The observation that AAG recognizes a single-nucleotide bulge as efficiently as a mismatch implies that the recognition of DNA damage is remarkably plastic.

Docking of hydroxamic acids into HDAC1 and HDAC8: a rationalization of activity trends and selectivities

A docking protocol using Gold software was developed to predict the binding disposition of histone deacetylase (HDAC) inhibitors, starting from the X-ray structures of HDAC8. The optimized procedure was subsequently utilized to dock into HDAC8 and into a homology model of HDAC1 nearly 40 compounds that had been tested for their inhibitory activity against the two HDAC isozymes. Evaluation of the best binding poses allowed us to identify the ligand properties and the protein residues important for activity and selectivity. HDACs are important anticancer drug targets, and their study is currently being actively pursued. As such, our results could help design new isozyme-selective HDAC inhibitors. Furthermore, this strategy may also be used for the investigation of other HDACs.

Frameshift deletion by Sulfolobus solfataricus P2 DNA polymerase Dpo4 T239W is selective for purines and involves normal conformational change followed by slow phosphodiester bond formation

The human DNA polymerase kappa homolog Sulfolobus solfataricus DNA polymerase IV (Dpo4) produces "-1" frameshift deletions while copying unmodified DNA and, more frequently, when bypassing DNA adducts. As judged by steady-state kinetics and mass spectrometry, bypass of purine template bases to produce these deletions occurred rarely but with 10-fold higher frequency than with pyrimidines. The DNA adduct 1,N(2)-etheno-2'-deoxyguanosine, with a larger stacking surface than canonical purines, showed the highest frequency of formation of -1 frameshift deletions. Dpo4 T239W, a mutant we had previously shown to produce fluorescence changes attributed to conformational change following dNTP binding opposite cognate bases (Beckman, J. W., Wang, Q., and Guengerich, F. P. (2008) J. Biol. Chem. 283, 36711-36723), reported similar conformational changes when the incoming dNTP complemented the base following a templating purine base or bulky adduct (i.e. the "+1" base). However, in all mispairing cases, phosphodiester bond formation was inefficient. The frequency of -1 frameshift events and the associated conformational changes were not dependent on the context of the remainder of the sequence. Collectively, our results support a mechanism for -1 frameshift deletions by Dpo4 that involves formation of active complexes via a favorable conformational change that skips the templating base, without causing slippage or flipping out of the base, to incorporate a complementary residue opposite the +1 base, in a mechanism previously termed "dNTP-stabilized incorporation." The driving force is attributed to be the stacking potential between the templating base and the incoming dNTP base.