DNA polymerase: structural homology, conformational dynamics, and the effects of carcinogenic DNA adducts

Non-canonical amino acids uncover the significant impact of Tyr671 on Taq DNA polymerase catalytic activity

Responsible for synthesizing the complementary strand of the DNA template, DNA polymerase is a crucial enzyme in DNA replication, recombination and repair. A highly conserved tyrosine (Tyr), located at the C-terminus of the O-helix in family A DNA polymerases, plays a critical role in enzyme activity and fidelity. Here, we combined the technology of genetic code extension to incorporate non-canonical amino acids and molecular dynamics (MD) simulations to uncover the mechanisms by which Tyr671 impacts substrate binding and conformation transitions in a DNA polymerase from Thermus aquaticus. Five non-canonical amino acids, namely l-3,4-dihydroxyphenylalanine (l-DOPA), p-aminophenylalanine (pAF), p-acetylphenylalanine (pAcF), p-cyanophenylalanine (pCNF) and p-nitrophenylalanine (pNTF), were individually incorporated at position 671. Strikingly, Y671pAF and Y671DOPA were active, but with lower activity compared to Y671F and wild-type. Y671pAF showed a higher fidelity than the Y671F, despite both possessing lower fidelity than the wild-type. Metadynamics and long-timescale MD simulations were carried out to probe the role of mutations in affecting protein structure, including open conformation, open-to-closed conformation transition, closed conformation, and closed-to-open conformation transition. The MD simulations clearly revealed that the size of the 671 amino acid residue and interactions with substrate or nearby residues were critical for Tyr671 to determine enzyme activity and fidelity.

Effects of the Y432S Cancer-Associated Variant on the Reaction Mechanism of Human DNA Polymerase κ

Human DNA polymerases are vital for genetic information management. Their function involves catalyzing the synthesis of DNA strands with unparalleled accuracy, which ensures the fidelity and stability of the human genomic blueprint. Several disease-associated mutations and their functional impact on DNA polymerases have been reported. One particular polymerase, human DNA polymerase kappa (Pol κ), has been reported to be susceptible to several cancer-associated mutations. The Y432S mutation in Pol κ, associated with various cancers, is of interest due to its impact on polymerization activity and markedly reduced thermal stability. Here, we have used computational simulations to investigate the functional consequences of the Y432S using classical molecular dynamics (MD) and coupled quantum mechanics/molecular mechanics (QM/MM) methods. Our findings suggest that Y432S induces structural alterations in domains responsible for nucleotide addition and ternary complex stabilization while retaining structural features consistent with possible catalysis in the active site. Calculations of the minimum energy path associated with the reaction mechanism of the wild type (WT) and Y432S Pol κ indicate that, while both enzymes are catalytically competent (in terms of energetics and the active site's geometries), the cancer mutation results in an endoergic reaction and an increase in the catalytic barrier. Interactions with a third magnesium ion and environmental effects on nonbonded interactions, particularly involving key residues, contribute to the kinetic and thermodynamic distinctions between the WT and mutant during the catalytic reaction. The energetics and electronic findings suggest that active site residues favor the catalytic reaction with dCTP3- over dCTP4-.

Unlocking The Mysteries of DNA Adducts with Artificial Intelligence

Cellular genome is considered a dynamic blueprint of a cell since it encodes genetic information that gets temporally altered due to various endogenous and exogenous insults. Largely, the extent of genomic dynamicity is controlled by the trade-off between DNA repair processes and the genotoxic potential of the causative agent (genotoxins or potential carcinogens). A subset of genotoxins form DNA adducts by covalently binding to the cellular DNA, triggering structural or functional changes that lead to significant alterations in cellular processes via genetic (e. g., mutations) or non-genetic (e. g., epigenome) routes. Identification, quantification, and characterization of DNA adducts are indispensable for their comprehensive understanding and could expedite the ongoing efforts in predicting carcinogenicity and their mode of action. In this review, we elaborate on using Artificial Intelligence (AI)-based modeling in adducts biology and present multiple computational strategies to gain advancements in decoding DNA adducts. The proposed AI-based strategies encompass predictive modeling for adduct formation via metabolic activation, novel adducts' identification, prediction of biochemical routes for adduct formation, adducts' half-life predictions within biological ecosystems, and, establishing methods to predict the link between adducts chemistry and its location within the genomic DNA. In summary, we discuss some futuristic AI-based approaches in DNA adduct biology.

Understanding the Effect of Multiple Domain Deletion in DNA Polymerase I from Geobacillus Sp. Strain SK72

The molecular structure of DNA polymerase I or family A polymerases is made up of three major domains that consist of a single polymerase domain with two extra exonuclease domains. When the N-terminal was deleted, the enzyme was still able to perform basic polymerase activity with additional traits that used isothermal amplification. However, the 3′-5′ exonuclease domain that carries a proofreading activity was disabled. Yet, the structure remained attached to the 5′-3′ polymerization domain without affecting its ability. The purpose of this non-functional domain still remains scarce. It either gives negative effects or provides structural support to the DNA polymerase. Here, we compared the effect of deleting each domain against the polymerase activity. The recombinant wild type and its variants were successfully purified and characterized. Interestingly, SK72-Exo (a large fragment excluding the 5′-3′ exonuclease domain) exhibited better catalytic activity than the native SK72 (with all three domains) at similar optimum temperature and pH profile, and it showed longer stability at 70 °C. Meanwhile, SK72-Exo2 (polymerization domain without both the 5′-3′ and 3′-5′ exonuclease domain) displayed the lowest activity with an optimum at 40 °C and favored a more neutral environment. It was also the least stable among the variants, with almost no activity at 50 °C for the first 10 min. In conclusion, cutting both exonuclease domains in DNA polymerase I has a detrimental effect on the polymerization activity and structural stability.

Computational Simulations of DNA Polymerases: Detailed Insights on Structure/Function/Mechanism from Native Proteins to Cancer Variants

Genetic information is vital in the cell cycle of DNA-based organisms. DNA polymerases (DNA Pols) are crucial players in transactions dealing with these processes. Therefore, the detailed understanding of the structure, function, and mechanism of these proteins has been the focus of significant effort. Computational simulations have been applied to investigate various facets of DNA polymerase structure and function. These simulations have provided significant insights over the years. This perspective presents the results of various computational studies that have been employed to research different aspects of DNA polymerases including detailed reaction mechanism investigation, mutagenicity of different metal cations, possible factors for fidelity synthesis, and discovery/functional characterization of cancer-related mutations on DNA polymerases.

Investigation of Intradomain Motions of a Y-Family DNA Polymerase during Substrate Binding and Catalysis

DNA polymerases catalyze DNA synthesis through a stepwise kinetic mechanism that begins with binding to DNA, followed by selection, binding, and incorporation of a nucleotide into an elongating primer. It is hypothesized that subtle active site adjustments in a polymerase to align reactive moieties limit the rate of correct nucleotide incorporation. DNA damage can impede this process for many DNA polymerases, causing replication fork stalling, genetic mutations, and potentially cell death. However, specialized Y-family DNA polymerases are structurally evolved to efficiently bypass DNA damage in vivo, albeit at the expense of replication fidelity. Dpo4, a model Y-family polymerase from Sulfolobus solfataricus, has been well-studied kinetically, structurally, and computationally, which yielded a mechanistic understanding of how the Y-family DNA polymerases achieve their unique catalytic properties. We previously employed a real-time Förster resonance energy transfer (FRET) technique to characterize the global conformational motions of Dpo4 during DNA binding as well as nucleotide binding and incorporation by monitoring changes in distance between sites on the polymerase and DNA, and even between domains of Dpo4. Here, we extend the utility of our FRET methodology to observe conformational transitions within individual domains of Dpo4 during DNA binding and nucleotide incorporation. The results of this novel, intradomain FRET approach unify findings from many studies to fully clarify the complex DNA binding mechanism of Dpo4. Furthermore, intradomain motions in the Finger domain during nucleotide binding and incorporation, for the first time, report on the rate-limiting step of a single-nucleotide addition catalyzed by Dpo4.

Steric gate residues of Y-family DNA polymerases DinB and pol kappa are crucial for dNTP-induced conformational change

Discrimination against ribonucleotides by DNA polymerases is critical to preserve DNA integrity. For many DNA polymerases, including those of the Y family, rNTP discrimination has been attributed to steric clashes between a residue near the active site, the steric gate, and the 2'-hydroxyl of the incoming rNTP. Here we used hydrogen/deuterium exchange (HDX) mass spectrometry (MS) to probe the effects of the steric gate in the Y-family DNA polymerases Escherichia coli DinB and human DNA pol κ. Formation of a ternary complex with a G:dCTP base pair in the active site resulted in slower hydrogen exchange relative to a ternary complex with G:rCTP in the active site. The protection from exchange was localized to regions both distal and proximal to the active site, suggesting that DinB and DNA pol κ adopt different conformations depending on the sugar of the incoming nucleotide. In contrast, when the respective steric gate residues were mutated to alanine, the differences in HDX between the dNTP- and rNTP-bound ternary complexes were attenuated such that for DinB(F13A) and pol κ(Y112A), ternary complexes with either G:dCTP or G:rCTP base pairs had similar HDX profiles. Furthermore, the HDX in these ternary complexes resembled that of the rCTP-bound state rather than the dCTP-bound state of the wild-type enzymes. Primer extension assays confirmed that DinB(F13A) and pol κ(Y112A) do not discriminate against rNTPs to the same extent as the wild-type enzymes. Our observations indicate that the steric gate is crucial for rNTP discrimination because of its role in specifically promoting a dNTP-induced conformational change and that rNTP discrimination occurs in a relatively closed state of the polymerases.

Real-time surface plasmon resonance study of biomolecular interactions between polymerase and bulky mutagenic DNA lesions

Surface plasmon resonance (SPR) was used to measure polymerase-binding interactions of the bulky mutagenic DNA lesions N-(2′-deoxyguanosin-8-yl)-4′-fluoro-4-aminobiphenyl (FABP) or N-(2′-deoxyguanosin-8-yl)-7-fluoro-2-acetylaminofluorene (FAAF) in the context of two unique 5′-flanking bases (CG*A and TG*A). The enzymes used were exo-nuclease-deficient Klenow fragment (Kf-exo^–) or polymerase β (pol β). Specific binary and ternary DNA binding affinities of the enzymes were characterized at subnanomolar concentrations. The SPR results showed that Kf-exo^– binds strongly to a double strand/single strand template/primer junction, whereas pol β binds preferentially to double-stranded DNA having a one-nucleotide gap. Both enzymes exhibited tight binding to native DNA, with high nucleotide selectivity, where the K_D values for each base pair increased in the order dCTP ≪ dTTP ∼ dATP ≪ dGTP. In contrast to that for pol β, Kf-exo^– binds tightly to lesion-modified templates; however, both polymerases exhibited minimal nucleotide selectivity toward adducted DNA. Primer steady-state kinetics and ¹⁹F NMR results support the SPR data. The relative insertion efficiency f_ins of dCTP opposite FABP was significantly higher in the TG*A sequence compared to that in CG*A. Although Kf-exo^– was not sensitive to the presence of a DNA lesion, FAAF-induced conformational heterogeneity perturbed the active site of pol β, weakening the enzyme’s ability to bind to FAAF adducts compared to FABP adducts. The present study demonstrates the effectiveness of SPR for elucidating how lesion-induced conformational heterogeneity affects the binding capability of polymerases and ultimately the nucleotide insertion efficiency.

Conformational insights into the lesion and sequence effects for arylamine-induced translesion DNA synthesis: 19F NMR, surface plasmon resonance, and primer kinetic studies

Adduct-induced DNA damage can affect transcription efficiency and DNA replication and repair. We previously investigated the effects of the 3′-next flanking base (G*CT vs G*CA; G*, FABP, N-(2′-deoxyguanosin-8-yl)-4′-fluoro-4-aminobiphenyl; FAF, N-(2′-deoxyguanosin-8-yl)-7-fluoro-2-aminofluorene) on the conformation of arylamine-DNA lesions in relation to E. coli nucleotide excision repair (JainV., HiltonB., LinB., PatnaikS., LiangF., DarianE., ZouY., MackerellA. D.Jr., and ChoB. P. (2013) Nucleic Acids Res., 41, 869−88023180767). Here, we report the differential effects of the same pair of sequences on DNA replication in vitro by the polymerases exofree Klenow fragment (Kf-exo^–) and Dpo4. We obtained dynamic ¹⁹F NMR spectra for two 19-mer modified templates during primer elongation: G*CA [d(5′-CTTACCATCG*CAACCATTC-3′)] and G*CT [d(5′-CTTACCATCG*CTACCATTC-3′)]. We found that lesion stacking is favored in the G*CT sequence compared to the G*CA counterpart. Surface plasmon resonance binding results showed consistently weaker affinities for the modified DNA with the binding strength in the order of FABP > FAF and G*CA > G*CT. Primer extension was stalled at (n) and near (n – 1 and n + 1) the lesion site, and the extent of blockage and the extension rates across the lesion were influenced by not only the DNA sequences but also the nature of the adduct’s chemical structure (FAF vs FABP) and the polymerase employed (Kf-exo^– vs Dpo4). Steady-state kinetics analysis with Kf-exo^– revealed the most dramatic sequence and lesion effects at the lesion (n) and postinsertion (n + 1) sites, respectively. Taken together, these results provide insights into the important role of lesion-induced conformational heterogeneity in modulating translesion DNA synthesis.

C(α) torsion angles as a flexible criterion to extract secrets from a molecular dynamics simulation

Given the increasing complexity of simulated molecular systems, and the fact that simulation times have now reached milliseconds to seconds, immense amounts of data (in the gigabyte to terabyte range) are produced in current molecular dynamics simulations. Manual analysis of these data is a very time-consuming task, and important events that lead from one intermediate structure to another can become occluded in the noise resulting from random thermal fluctuations. To overcome these problems and facilitate a semi-automated data analysis, we introduce in this work a measure based on C(α) torsion angles: torsion angles formed by four consecutive C(α) atoms. This measure describes changes in the backbones of large systems on a residual length scale (i.e., a small number of residues at a time). Cluster analysis of individual C(α) torsion angles and its fuzzification led to continuous time patches representing (meta)stable conformations and to the identification of events acting as transitions between these conformations. The importance of a change in torsion angle to structural integrity is assessed by comparing this change to the average fluctuations in the same torsion angle over the complete simulation. Using this novel measure in combination with other measures such as the root mean square deviation (RMSD) and time series of distance measures, we performed an in-depth analysis of a simulation of the open form of DNA polymerase I. The times at which major conformational changes occur and the most important parts of the molecule and their interrelations were pinpointed in this analysis. The simultaneous determination of the time points and localizations of major events is a significant advantage of the new bottom-up approach presented here, as compared to many other (top-down) approaches in which only the similarity of the complete structure is analyzed.

Carcinogenic adducts induce distinct DNA polymerase binding orientations

DNA polymerases must accurately replicate DNA to maintain genome integrity. Carcinogenic adducts, such as 2-aminofluorene (AF) and N-acetyl-2-aminofluorene (AAF), covalently bind DNA bases and promote mutagenesis near the adduct site. The mechanism by which carcinogenic adducts inhibit DNA synthesis and cause mutagenesis remains unclear. Here, we measure interactions between a DNA polymerase and carcinogenic DNA adducts in real-time by single-molecule fluorescence. We find the degree to which an adduct affects polymerase binding to the DNA depends on the adduct location with respect to the primer terminus, the adduct structure and the nucleotides present in the solution. Not only do the adducts influence the polymerase dwell time on the DNA but also its binding position and orientation. Finally, we have directly observed an adduct- and mismatch-induced intermediate state, which may be an obligatory step in the DNA polymerase proofreading mechanism.

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.