Perspective: pre-chemistry conformational changes in DNA polymerase mechanisms

Correct Nucleotide Selection Is Confined at the Binding Site of Polymerase Enzymes

DNA polymerases (Pols) add incoming nucleotides (deoxyribonucleoside triphosphate (dNTPs)) to growing DNA strands, a crucial step for DNA synthesis. The insertion of correct (vs incorrect) nucleotides relates to Pols' fidelity, which defines Pols' ability to faithfully replicate DNA strands in a template-dependent manner. We and others have demonstrated that reactant alignment and correct base pairing at the Pols catalytic site are crucial structural features to fidelity. Here, we first used equilibrium molecular simulations to demonstrate that the local dynamics at the protein-DNA interface in the proximity of the catalytic site is different when correct vs incorrect dNTPs are bound to polymerase β (Pol β). Formation and dynamic stability of specific interatomic interactions around the incoming nucleotide influence the overall binding site architecture. This explains why certain Pols' mutants can affect the local catalytic environment and influence the selection of correct vs incorrect nucleotides. In particular, this is here demonstrated by analyzing the interaction network formed by the residue R283, whose mutant R283A has an experimentally measured lower capacity of differentiating correct (G:dCTP) vs incorrect (G:dATP) base pairing in Pol β. We also used alchemical free-energy calculations to quantify the G:dCTP →G:dATP transformation in Pol β wild-type and mutant R283A. These results correlate well with the experimental trend, thus corroborating our mechanistic insights. Sequence and structural comparisons with other Pols from the same family suggest that these findings may also be valid in similar enzymes.

Trapping a non-cognate nucleotide upon initial binding for replication fidelity control in SARS-CoV-2 RNA dependent RNA polymerase

The RNA dependent RNA polymerase (RdRp) in SARS-CoV-2 is a highly conserved enzyme responsible for viral genome replication/transcription. To understand how the viral RdRp achieves fidelity control during such processes, here we computationally investigate the natural non-cognate vs. cognate nucleotide addition and selectivity during viral RdRp elongation. We focus on the nucleotide substrate initial binding (RdRp active site open) to the prechemical insertion (active site closed) of the RdRp. The current studies were first carried out using microsecond ensemble equilibrium all-atom molecular dynamics (MD) simulations. Due to the slow conformational changes (from open to closed) during nucleotide insertion and selection, enhanced or umbrella sampling methods have been further employed to calculate the free energy profiles of the nucleotide insertion. Our studies find notable stability of noncognate dATP and GTP upon initial binding in the active-site open state. The results indicate that while natural cognate ATP and Remdesivir drug analogue (RDV-TP) are biased toward stabilization in the closed state to facilitate insertion, the natural non-cognate dATP and GTP can be well trapped in off-path initial binding configurations and prevented from insertion so that to be further rejected. The current work thus presents the intrinsic nucleotide selectivity of SARS-CoV-2 RdRp for natural substrate fidelity control, which should be considered in antiviral drug design.

DNA polymerase β: Closing the gap between structure and function

DNA polymerase (dpol) β has served as a model for structural, kinetic, and computational characterization of the DNA synthesis reaction. The laboratory directed by Samuel H. Wilson has utilized a multifunctional approach to analyze the function of this enzyme at the biological, chemical, and molecular levels for nearly 50 years. Over this time, it has become evident that correlating static crystallographic structures of dpol β with solution kinetic measurements is a daunting task. However, aided by computational and spectroscopic approaches, novel and unexpected insights have emerged. While dpols generally insert wrong nucleotides with similar poor efficiencies, their capacity to insert the right nucleotide depends on the identity of the dpol. Accordingly, the ability to choose right from wrong depends on the efficiency of right, rather than wrong, nucleotide insertion. Structures of dpol β in various liganded forms published by the Wilson laboratory, and others, have provided molecular insights into the molecular attributes that hasten correct nucleotide insertion and deter incorrect nucleotide insertion. Computational approaches have bridged the gap between structures of intermediate complexes and provided insights into this basic and essential chemical reaction.

Visualizing Rev1 catalyze protein-template DNA synthesis

During DNA replication, replicative DNA polymerases may encounter DNA lesions, which can stall replication forks. One way to prevent replication fork stalling is through the recruitment of specialized translesion synthesis (TLS) polymerases that have evolved to incorporate nucleotides opposite DNA lesions. Rev1 is a specialized TLS polymerase that bypasses abasic sites, as well as minor-groove and exocyclic guanine adducts. Lesion bypass is accomplished using a unique protein-template mechanism in which the templating base is evicted from the DNA helix and the incoming dCTP hydrogen bonds with an arginine side chain of Rev1. To understand the protein-template mechanism at an atomic level, we employed a combination of time-lapse X-ray crystallography, molecular dynamics simulations, and DNA enzymology on the Saccharomyces cerevisiae Rev1 protein. We find that Rev1 evicts the templating base from the DNA helix prior to binding the incoming nucleotide. Binding the incoming nucleotide changes the conformation of the DNA substrate to orient it for nucleotidyl transfer, although this is not coupled to large structural changes in Rev1 like those observed with other DNA polymerases. Moreover, we found that following nucleotide incorporation, Rev1 converts the pyrophosphate product to two monophosphates, which drives the reaction in the forward direction and prevents pyrophosphorolysis. Following nucleotide incorporation, the hydrogen bonds between the incorporated nucleotide and the arginine side chain are broken, but the templating base remains extrahelical. These postcatalytic changes prevent potentially mutagenic processive synthesis by Rev1 and facilitate dissociation of the DNA product from the enzyme.

Eukaryotic Base Excision Repair: New Approaches Shine Light on Mechanism

Genomic DNA is susceptible to endogenous and environmental stresses that modify DNA structure and its coding potential. Correspondingly, cells have evolved intricate DNA repair systems to deter changes to their genetic material. Base excision DNA repair involves a number of enzymes and protein cofactors that hasten repair of damaged DNA bases. Recent advances have identified macromolecular complexes that assemble at the DNA lesion and mediate repair. The repair of base lesions generally requires five enzymatic activities: glycosylase, endonuclease, lyase, polymerase, and ligase. The protein cofactors and mechanisms for coordinating the sequential enzymatic steps of repair are being revealed through a range of experimental approaches. We discuss the enzymes and protein cofactors involved in eukaryotic base excision repair, emphasizing the challenge of integrating findings from multiple methodologies. The results provide an opportunity to assimilate biochemical findings with cell-based assays to uncover new insights into this deceptively complex repair pathway.

Evidence that a catalytic glutamate and an 'Arginine Toggle' act in concert to mediate ATP hydrolysis and mechanochemical coupling in a viral DNA packaging motor

ASCE ATPases include ring-translocases such as cellular helicases and viral DNA packaging motors (terminases). These motors have conserved Walker A and B motifs that bind Mg²⁺-ATP and a catalytic carboxylate that activates water for hydrolysis. Here we demonstrate that Glu179 serves as the catalytic carboxylate in bacteriophage λ terminase and probe its mechanistic role. All changes of Glu179 are lethal: non-conservative changes abrogate ATP hydrolysis and DNA translocation, while the conservative E179D change attenuates ATP hydrolysis and alters single molecule translocation dynamics, consistent with a slowed chemical hydrolysis step. Molecular dynamics simulations of several homologous terminases suggest a novel mechanism, supported by experiments, wherein the conserved Walker A arginine ‘toggles’ between interacting with a glutamate residue in the ‘lid’ subdomain and the catalytic glutamate upon ATP binding; this switch helps mediate a transition from an ‘open’ state to a ‘closed’ state that tightly binds nucleotide and DNA, and also positions the catalytic glutamate next to the γ-phosphate to align the hydrolysis transition state. Concomitant reorientation of the lid subdomain may mediate mechanochemical coupling of ATP hydrolysis and DNA translocation. Given the strong conservation of these structural elements in terminase enzymes, this mechanism may be universal for viral packaging motors.

Determining selection free energetics from nucleotide pre-insertion to insertion in viral T7 RNA polymerase transcription fidelity control

An elongation cycle of a transcribing RNA polymerase (RNAP) usually consists of multiple kinetics steps, so there exist multiple kinetic checkpoints where non-cognate nucleotides can be selected against. We conducted comprehensive free energy calculations on various nucleotide insertions for viral T7 RNAP employing all-atom molecular dynamics simulations. By comparing insertion free energy profiles between the non-cognate nucleotide species (rGTP and dATP) and a cognate one (rATP), we obtained selection free energetics from the nucleotide pre-insertion to the insertion checkpoints, and further inferred the selection energetics down to the catalytic stage. We find that the insertion of base mismatch rGTP proceeds mainly through an off-path along which both pre-insertion screening and insertion inhibition play significant roles. In comparison, the selection against dATP is found to go through an off-path pre-insertion screening along with an on-path insertion inhibition. Interestingly, we notice that two magnesium ions switch roles of leave and stay during the dATP on-path insertion. Finally, we infer that substantial selection energetic is still required to catalytically inhibit the mismatched rGTP to achieve an elongation error rate ∼10-4 or lower; while no catalytic selection seems to be further needed against dATP to obtain an error rate ∼10-2.

C, Da LT, Yu J. A Viral T7 RNA Polymerase Ratcheting Along DNA With Fidelity Control

RNA polymerase (RNAP) from bacteriophage T7 is a representative single-subunit viral RNAP that can transcribe with high promoter activities without assistances from transcription factors. We accordingly studied this small transcription machine computationally as a model system to understand underlying mechanisms of mechano-chemical coupling and fidelity control in the RNAP transcription elongation. Here we summarize our computational work from several recent publications to demonstrate first how T7 RNAP translocates via Brownian alike motions along DNA right after the catalytic product release. Then we show how the backward translocation motions are prevented at post-translocation upon successful nucleotide incorporation, which is also subject to stepwise nucleotide selection and acts as a pawl for "selective ratcheting". The structural dynamics and energetics features revealed from our atomistic molecular dynamics (MD) simulations and related analyses on the single-subunit T7 RNAP thus provided detailed and quantitative characterizations on the Brownian-ratchet working scenario of a prototypical transcription machine with sophisticated nucleotide selectivity for fidelity control. The presented mechanisms can be more or less general for structurally similar viral or mitochondrial RNAPs and some of DNA polymerases, or even for the RNAP engine of the more complicated transcription machinery in higher organisms.

Molecular dynamics simulations suggest changes in electrostatic interactions as a potential mechanism through which serine phosphorylation inhibits DNA polymerase β activity

DNA polymerase β is a 39 kDa enzyme that is a major component of Base Excision Repair in human cells. The enzyme comprises two major domains, a 31 kDa domain responsible for the polymerase activity and an 8 kDa domain, which bind ssDNA and has a deoxyribose phosphate (dRP) lyase activity. DNA polymerase β was shown to be phosphorylated in vitro with protein kinase C (PKC) at serines 44 and 55 (S44 and S55), resulting in loss of its polymerase enzymic activity, but not its ability to bind ssDNA. In this study, we investigate the potential phosphorylation-induced structural changes for DNA polymerase β using molecular dynamics simulations. The simulations show drastic conformational changes of the polymerase structure as a result of S44 phosphorylation. Phosphorylation-induced conformational changes transform the closed (active) enzyme structure into an open one. Further analysis of the results points to a key hydrogen bond and newly formed salt bridges as potential drivers of these structural fluctuations. The changes observed with S55/44 and S55 phosphorylation were less dramatic and the integrity of the H-bond was not compromised. Thus the phosphorylation of S44 is the major contributor to structural fluctuations that lead to loss of enzymatic activity.

Mapping Functional Substrate-Enzyme Interactions in the pol β Active Site through Chemical Biology: Structural Responses to Acidity Modification of Incoming dNTPs

We report high-resolution crystal structures of DNA polymerase (pol) β in ternary complex with a panel of incoming dNTPs carrying acidity-modified 5'-triphosphate groups. These novel dNTP analogues have a variety of halomethylene substitutions replacing the bridging oxygen between Pβ and Pγ of the incoming dNTP, whereas other analogues have alkaline substitutions at the bridging oxygen. Use of these analogues allows the first systematic comparison of effects of 5'-triphosphate acidity modification on active site structures and the rate constant of DNA synthesis. These ternary complex structures with incoming dATP, dTTP, and dCTP analogues reveal the enzyme's active site is not grossly altered by the acidity modifications of the triphosphate group, yet with analogues of all three incoming dNTP bases, subtle structural differences are apparent in interactions around the nascent base pair and at the guanidinium groups of active site arginine residues. These results are important for understanding how acidity modification of the incoming dNTP's 5'-triphosphate can influence DNA polymerase activity and the significance of interactions at arginines 183 and 149 in the active site.

Balancing Non-Equilibrium Driving with Nucleotide Selectivity at Kinetic Checkpoints in Polymerase Fidelity Control

High fidelity gene transcription and replication require kinetic discrimination of nucleotide substrate species by RNA and DNA polymerases under chemical non-equilibrium conditions. It is known that sufficiently large free energy driving force is needed for each polymerization or elongation cycle to maintain far-from-equilibrium to achieve low error rates. Considering that each cycle consists of multiple kinetic steps with different transition rates, one expects that the kinetic modulations by polymerases are not evenly conducted at each step. We show that accelerations at different kinetic steps impact quite differently to the overall elongation characteristics. In particular, for forward transitions that discriminate cognate and non-cognate nucleotide species to serve as kinetic selection checkpoints, the transition cannot be accelerated too quickly nor retained too slowly to obtain low error rates, as balancing is needed between the nucleotide selectivity and the non-equilibrium driving. Such a balance is not the same as the speed-accuracy tradeoff in which high accuracy is always obtained at sacrifice of speed. For illustration purposes, we used three-state and five-state models of nucleotide addition in the polymerase elongation and show how the non-equilibrium steady state characteristics change upon variations on stepwise forward or backward kinetics. Notably, by using the multi-step elongation schemes and parameters from T7 RNA polymerase transcription elongation, we demonstrate that individual transitions serving as selection checkpoints need to proceed at moderate rates in order to sustain the necessary non-equilibrium drives as well as to allow nucleotide selections for an optimal error control. We also illustrate why rate-limiting conformational transitions of the enzyme likely play a significant role in the error reduction.

Molecular dynamics simulations suggest changes in electrostatic interactions as a potential mechanism through which serine phosphorylation inhibits DNA Polymerase β's activity

DNA polymerase β is a 39kDa enzyme that is a major component of Base Excision Repair in human cells. The enzyme comprises two major domains, a 31kDa domain responsible for the polymerase activity and an 8kDa domain, which bind ssDNA and has a deoxyribose phosphate (dRP) lyase activity. DNA polymerase β was shown to be phosphorylated in vitro with protein kinase C (PKC) at serines 44 and 55 (S44 and S55), resulting in loss of its polymerase enzymic activity, but not its ability to bind ssDNA. In this study, we investigate the potential phosphorylation-induced structural changes for DNA polymerase β using molecular dynamics. The simulations show drastic conformational changes of the polymerase structure as a result of S44 phosphorylation. Phosphorylation-induced conformational changes transform the closed (active) enzyme structure into an open one. Further analysis of the results points to a key hydrogen bond and newly formed salt bridges as potential drivers of these structural fluctuations. The changes observed with S44/55 and S55 phosphorylation were less dramatic than S44 and the integrity of the H-bond was not compromised. Thus the phosphorylation of S44 is likely the major contributor to structural fluctuations that lead to loss of enzymatic activity.

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.

Molecular Dynamics Simulations and Structural Analysis to Decipher Functional Impact of a Twenty Residue Insert in the Ternary Complex of Mus musculus TdT Isoform

Insertions/deletions are common evolutionary tools employed to alter the structural and functional repertoire of protein domains. An insert situated proximal to the active site or ligand binding site frequently impacts protein function; however, the effect of distal indels on protein activity and/or stability are often not studied. In this paper, we have investigated a distal insert, which influences the function and stability of a unique DNA polymerase, called terminal deoxynucleotidyl transferase (TdT). TdT (EC:2.7.7.31) is a monomeric 58 kDa protein belonging to family X of eukaryotic DNA polymerases and known for its role in V(D)J recombination as well as in non-homologous end-joining (NHEJ) pathways. Two murine isoforms of TdT, with a length difference of twenty residues and having different biochemical properties, have been studied. All-atom molecular dynamics simulations at different temperatures and interaction network analyses were performed on the short and long-length isoforms. We observed conformational changes in the regions distal to the insert position (thumb subdomain) in the longer isoform, which indirectly affects the activity and stability of the enzyme through a mediating loop (Loop1). A structural rationale could be provided to explain the reduced polymerization rate as well as increased thermosensitivity of the longer isoform caused by peripherally located length variations within a DNA polymerase. These observations increase our understanding of the roles of length variants in introducing functional diversity in protein families in general.

Insertion of oxidized nucleotide triggers rapid DNA polymerase opening

A novel mechanism is unveiled to explain why a pro-mutagenic nucleotide lesion (oxidized guanine, 8-oxoG) causes the mammalian DNA repair polymerase-β (pol-β) to rapidly transition to an inactive open conformation. The mechanism involves unexpected features revealed recently in time-lapse crystallography. Specifically, a delicate water network associated with a lesion-stabilizing auxilliary product ion Mg(p) triggers a cascade of events that leads to poor active site geometry and the rupture of crucial molecular interactions between key residues in both the anti(8-oxoG:C) and syn(8-oxoG:A) systems. Once the base pairs in these lesioned systems are broken, dislocation of both Asp192 (a metal coordinating ligand) and the oxoG phosphate group (PO₄) interfere with the hydrogen bonding between Asp192 and Arg258, whose rotation toward Asp192 is crucial to the closed-to-open enzyme transition. Energetically, the lesioned open states are similar in energy to those of the corresponding closed complexes after chemistry, in marked contrast to the unlesioned pol-β anti(G:C) system, whose open state is energetically higher than the closed state. The delicate surveillance system offers a fundamental protective mechanism in the cell that triggers DNA repair events which help deter insertion of oxidized lesions.

A Jump-from-Cavity Pyrophosphate Ion Release Assisted by a Key Lysine Residue in T7 RNA Polymerase Transcription Elongation

Pyrophosphate ion (PPi) release during transcription elongation is a signature step in each nucleotide addition cycle. The kinetics and energetics of the process as well as how it proceeds with substantial conformational changes of the polymerase complex determine the mechano-chemical coupling mechanism of the transcription elongation. Here we investigated detailed dynamics of the PPi release process in a single-subunit RNA polymerase (RNAP) from bacteriophage T7, implementing all-atom molecular dynamics (MD) simulations. We obtained a jump-from-cavity kinetic model of the PPi release utilizing extensive nanosecond MD simulations. We found that the PPi release in T7 RNAP is initiated by the PPi dissociation from two catalytic aspartic acids, followed by a comparatively slow jump-from-cavity activation process. Combining with a number of microsecond long MD simulations, we also found that the activation process is hindered by charged residue associations as well as by local steric and hydrogen bond interactions. On the other hand, the activation is greatly assisted by a highly flexible lysine residue Lys472 that swings its side chain to pull PPi out. The mechanism can apply in general to single subunit RNA and DNA polymerases with similar molecular structures and conserved key residues. Remarkably, the flexible lysine or arginine residue appears to be a universal module that assists the PPi release even in multi-subunit RNAPs with charge facilitated hopping mechanisms. We also noticed that the PPi release is not tightly coupled to opening motions of an O-helix on the fingers domain of T7 RNAP according to the microsecond MD simulations. Our study thus supports the Brownian ratchet scenario of the mechano-chemical coupling in the transcription elongation of the single-subunit polymerase.

The Closing Mechanism of DNA Polymerase I at Atomic Resolution

DNA polymerases must quickly and accurately distinguish between similar nucleic acids to form Watson-Crick base pairs and avoid DNA replication errors. Deoxynucleoside triphosphate (dNTP) binding to the DNA polymerase active site induces a large conformational change that is difficult to characterize experimentally on an atomic level. Here, we report an X-ray crystal structure of DNA polymerase I bound to DNA in the open conformation with a dNTP present in the active site. We use this structure to computationally simulate the open to closed transition of DNA polymerase in the presence of a Watson-Crick base pair. Our microsecond simulations allowed us to characterize the key steps involved in active site assembly, and propose the sequence of events involved in the prechemistry steps of DNA polymerase catalysis. They also reveal new features of the polymerase mechanism, such as a conserved histidine as a potential proton acceptor from the primer 3'-hydroxyl.

Noncognate DNA damage prevents the formation of the active conformation of the Y-family DNA polymerases DinB and DNA polymerase κ

Y-family DNA polymerases are specialized to copy damaged DNA, and are associated with increased mutagenesis, owing to their low fidelity. It is believed that the mechanism of nucleotide selection by Y-family DNA polymerases involves conformational changes preceding nucleotidyl transfer, but there is limited experimental evidence for such structural changes. In particular, nucleotide-induced conformational changes in bacterial or eukaryotic Y-family DNA polymerases have, to date, not been extensively characterized. Using hydrogen-deuterium exchange mass spectrometry, we demonstrate here that the Escherichia coli Y-family DNA polymerase DinB and its human ortholog DNA polymerase κ undergo a conserved nucleotide-induced conformational change in the presence of undamaged DNA and the correct incoming nucleotide. Notably, this holds true for damaged DNA containing N(2) -furfuryl-deoxyguanosine, which is efficiently copied by these two polymerases, but not for damaged DNA containing the major groove modification O(6) -methyl-deoxyguanosine, which is a poor substrate. Our observations suggest that DinB and DNA polymerase κ utilize a common mechanism for nucleotide selection involving a conserved prechemical conformational transition promoted by the correct nucleotide and only preferred DNA substrates.

DPT tautomerisation of the G·A(syn) and A*·G*(syn) DNA mismatches: a QM/QTAIM combined atomistic investigation

By applying a combined QM and QTAIM atomistic computational approach we have established for the first time that the G·A(syn) and A*·G*(syn) DNA mismatches (rare tautomers are marked with an asterisk), causing spontaneous transversions with substantially various probabilities, radically differ from each other in their ability to tautomerise through the double proton transfer (DPT). The A*·G*(syn) mismatch tautomerises quite easily (ΔΔG(TS) ≈ 4·kT at room temperature) into the A·G*(syn) mismatch through the asynchronous concerted mechanism, whereas the G·A(syn) base mispair does not tautomerise via the DPT at all, since there is no local minimum corresponding to the tautomerised G*·A*(syn) mismatch on the potential energy surface. It was established that the A·G*(syn) base mispair is a dynamically unstable H-bonded complex with an extremely short lifetime of 2.17 × 10(-13) s. Consequently, the obtained results allow us to believe that spontaneous or forced dissociation of both the G·A(syn) and A*·G*(syn) DNA mismatches by the DNA-polymerase occurs with the preservation of the tautomeric status of the bases.

Applications of quantum mechanical/molecular mechanical methods to the chemical insertion step of DNA and RNA polymerization

We review theoretical attempts to model the chemical insertion reactions of nucleoside triphosphates catalyzed by the nucleic acid polymerases using combined quantum mechanical/molecular mechanical methodology. Due to an existing excellent database of high-resolution X-ray crystal structures, the DNA polymerase β system serves as a useful template for discussion and comparison. The convergence of structures of high-quality complexes and continued developments of theoretical techniques suggest a bright future for understanding the global features of nucleic acid polymerization.

How DNA polymerase X preferentially accommodates incoming dATP opposite 8-oxoguanine on the template

The modified base 8-oxo-7,8-dihydro-2'-deoxyguanosine (oxoG) is a common DNA adduct produced by the oxidation of DNA by reactive oxygen species. Kinetic data reveal that DNA polymerase X (pol X) from the African swine fever virus incorporates adenine (dATP) opposite to oxoG with higher efficiency than the non-damaged G:C basepair. To help interpret the kinetic data, we perform molecular dynamics simulations of pol X/DNA complexes, in which the template base opposite to the incoming dNTP (dCTP, dATP, dGTP) is oxoG. Our results suggest that pol X accommodates the oxoGsyn:A mispair by sampling closed active conformations that mirror those observed in traditional Watson-Crick complexes. Moreover, for both the oxoGsyn:A and oxoG:C ternary complexes, conformational sampling of the polymerase follows previously described large subdomain movements, local residue motions, and active site reorganization. Interestingly, the oxoGsyn:A system exhibits superior active site geometry in comparison to the oxoG:C system. Simulations for the other mismatch basepair complexes reveal large protein subdomain movement for all systems, except for oxoG:G, which samples conformations close to the open state. In addition, active site geometry and basepairing of the template base with the incoming nucleotide, reveal distortions and misalignments that range from moderate (i.e., oxoG:Asyn) to extreme (i.e., oxoGanti/syn:G). These results agree with the available kinetic data for pol X and provide structural insights regarding the mechanism by which this polymerase can accommodate incoming nucleotides opposite oxoG. Our simulations also support the notion that α-helix E is involved both in DNA binding and active site stabilization. Our proposed mechanism by which pol X can preferentially accommodate dATP opposite template oxoG further underscores the role that enzyme dynamics and conformational sampling operate in polymerase fidelity and function.

Conformational dynamics of Thermus aquaticus DNA polymerase I during catalysis

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

Structure and mechanism of DNA polymerase β

DNA polymerase (pol) β is a small eukaryotic DNA polymerase composed of two domains. Each domain contributes an enzymatic activity (DNA synthesis and deoxyribose phosphate lyase) during the repair of simple base lesions. These domains are termed the polymerase and lyase domains, respectively. Pol β has been an excellent model enzyme for studying the nucleotidyl transferase reaction and substrate discrimination at a molecular level. In this review, recent crystallographic studies of pol β in various liganded and conformational states during the insertion of right and wrong nucleotides as well as during the bypass of damaged DNA (apurinic sites and 8-oxoguanine) are described. Structures of these catalytic intermediates provide unexpected insights into mechanisms by which DNA polymerases enhance genome stability. These structures also provide an improved framework that permits computational studies to facilitate the interpretation of detailed kinetic analyses of this model enzyme.

Optimal and variant metal-ion routes in DNA polymerase β's conformational pathways

To interpret recent structures of the R283K mutant of human DNA repair enzyme DNA polymerase β (pol β) differing in the number of Mg(2+) ions, we apply transition path sampling (TPS) to assess the effect of differing ion placement on the transition from the open one-metal to the closed two-metal state. We find that the closing pathway depends on the initial ion position, both in terms of the individual transition states and associated energies. The energy barrier of the conformational pathway varies from 25 to 58 kJ/mol, compared to the conformational energy barrier of 42 kJ/mol for the wild-type pol β reported previously. Moreover, we find a preferred ion route located in the center of the enzyme, parallel to the DNA. Within this route, the conformational pathway is similar to that of the overall open to closed transition of pol β, but outside it, especially when the ion starts near active site residues Arg258 and Asp190, the conformational pathway diverges significantly. Our findings should apply generally to pol β, since R283K is relatively far from the active site; further experimental and computational work are required to confirm this. Our studies also underscore the common feature that less active mutants have less stable closed states than their open states, in marked contrast to the wild-type enzyme, where the closed state is significantly more stable than the open form.

The E295K cancer variant of human polymerase β favors the mismatch conformational pathway during nucleotide selection

DNA polymerase β (pol β) is responsible for gap filling synthesis during repair of damaged DNA as part of the base excision repair pathway. Human pol β mutations were recently identified in a high percentage (∼30%) of tumors. Characterization of specific cancer variants is particularly useful to further the understanding of the general mechanism of pol β while providing context to disease contribution. We showed that expression of the carcinoma variant E295K induces cellular transformation. The poor polymerase activity exhibited by the variant was hypothesized to be caused by the destabilization of proper active site assembly by the glutamate to lysine mutation. Here, we show that this variant exhibits an unusual preference for binding dCTP opposite a templating adenine over the cognate dTTP. Biochemical studies indicate that the noncognate competes with the cognate nucleotide for binding to the polymerase active site with the noncognate incorporation a function of higher affinity and not increased activity. In the crystal structure of the variant bound to dA:dCTP, the fingers domain closes around the mismatched base pair. Nucleotide incorporation is hindered because key residues in the polymerase active site are not properly positioned for nucleotidyl transfer. In contrast to the noncognate dCTP, neither the cognate dTTP nor its nonhydrolyzable analog induced fingers closure, as isomorphous difference Fourier maps show that the cognate nucleotides are bound to the open state of the polymerase. Comparison with published structures provides insight into the structural rearrangements within pol β that occur during the process of nucleotide discrimination.

"Gate-keeper" residues and active-site rearrangements in DNA polymerase μ help discriminate non-cognate nucleotides

Incorporating the cognate instead of non-cognate substrates is crucial for DNA polymerase function. Here we analyze molecular dynamics simulations of DNA polymerase μ (pol μ) bound to different non-cognate incoming nucleotides including A:dCTP, A:dGTP, A(syn):dGTP, A:dATP, A(syn):dATP, T:dCTP, and T:dGTP to study the structure-function relationships involved with aberrant base pairs in the conformational pathway; while a pol μ complex with the A:dTTP base pair is available, no solved non-cognate structures are available. We observe distinct differences of the non-cognate systems compared to the cognate system. Specifically, the motions of active-site residue His329 and Asp330 distort the active site, and Trp436, Gln440, Glu443 and Arg444 tend to tighten the nucleotide-binding pocket when non-cognate nucleotides are bound; the latter effect may further lead to an altered electrostatic potential within the active site. That most of these “gate-keeper” residues are located farther apart from the upstream primer in pol μ, compared to other X family members, also suggests an interesting relation to pol μ's ability to incorporate nucleotides when the upstream primer is not paired. By examining the correlated motions within pol μ complexes, we also observe different patterns of correlations between non-cognate systems and the cognate system, especially decreased interactions between the incoming nucleotides and the nucleotide-binding pocket. Altered correlated motions in non-cognate systems agree with our recently proposed hybrid conformational selection/induced-fit models. Taken together, our studies propose the following order for difficulty of non-cognate system insertions by pol μ: T:dGTP<A(syn):dATP<T:dCTP<A:dGTP<A(syn):dGTP<A:dCTP<A:dATP. This sequence agrees with available kinetic data for non-cognate nucleotide insertions, with the exception of A:dGTP, which may be more sensitive to the template sequence. The structures and conformational aspects predicted here are experimentally testable.

DNA polymerase μ (pol μ) is an enzyme that participates in DNA repair and thus has a central role in maintaining the integrity of genetic information. To efficiently repair the DNA, discriminating the cognate instead of non-cognate nucleotides (“fidelity-checking”) is required. Here we analyze molecular dynamics simulations of pol μ bound to different non-cognate nucleotides to study the structure-function relationships involved in the fidelity-checking mechanism of pol μ on the atomic level. Our results suggest that His329, Asp330, Trp436, Gln440, Glu443, and Arg444 are of great importance for pol μ's fidelity-checking mechanism. We also observe altered patterns of correlated motions within pol μ complex when non-cognate instead of cognate nucleotides are bound, which agrees with our recently proposed hybrid conformational selection/induced-fit models. Taken together, our studies help interpret the available kinetic data of various non-cognate nucleotide insertions by pol μ. We also suggest experimentally testable predictions; for example, a point mutation like E443M may reduce the ability of pol μ to insert the cognate more than of non-cognate nucleotides. Our studies suggest an interesting relation to pol μ's unique ability to incorporate nucleotides when the upstream primer is not paired.