Polymerase beta simulations suggest that Arg258 rotation is a slow step rather than large subdomain motions per se

Machines on Genes through the Computational Microscope

Macromolecular machines acting on genes are at the core of life's fundamental processes, including DNA replication and repair, gene transcription and regulation, chromatin packaging, RNA splicing, and genome editing. Here, we report the increasing role of computational biophysics in characterizing the mechanisms of "machines on genes", focusing on innovative applications of computational methods and their integration with structural and biophysical experiments. We showcase how state-of-the-art computational methods, including classical and ab initio molecular dynamics to enhanced sampling techniques, and coarse-grained approaches are used for understanding and exploring gene machines for real-world applications. As this review unfolds, advanced computational methods describe the biophysical function that is unseen through experimental techniques, accomplishing the power of the "computational microscope", an expression coined by Klaus Schulten to highlight the extraordinary capability of computer simulations. Pushing the frontiers of computational biophysics toward a pragmatic representation of large multimegadalton biomolecular complexes is instrumental in bridging the gap between experimentally obtained macroscopic observables and the molecular principles playing at the microscopic level. This understanding will help harness molecular machines for medical, pharmaceutical, and biotechnological purposes.

DNA Polymerase β in the Context of Cancer

DNA polymerase beta (Pol β) is a 39 kD vertebrate polymerase that lacks proofreading ability, yet still maintains a moderate fidelity of DNA synthesis. Pol β is a key enzyme that functions in the base excision repair and non-homologous end joining pathways of DNA repair. Mechanisms of fidelity for Pol β are still being elucidated but are likely to involve dynamic conformational motions of the enzyme upon its binding to DNA and deoxynucleoside triphosphates. Recent studies have linked germline and somatic variants of Pol β with cancer and autoimmunity. These variants induce genomic instability by a number of mechanisms, including error-prone DNA synthesis and accumulation of single nucleotide gaps that lead to replication stress. Here, we review the structure and function of Pol β, and we provide insights into how structural changes in Pol β variants may contribute to genomic instability, mutagenesis, disease, cancer development, and impacts on treatment outcomes.

Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.

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.

Crystal Structure of the Apicoplast DNA Polymerase from Plasmodium falciparum: The First Look at a Plastidic A-Family DNA Polymerase

Plasmodium falciparum, the primary cause of malaria, contains a non-photosynthetic plastid called the apicoplast. The apicoplast exists in most members of the phylum Apicomplexa and has its own genome along with organelle-specific enzymes for its replication. The only DNA polymerase found in the apicoplast (apPOL) was putatively acquired through horizontal gene transfer from a bacteriophage and is classified as an atypical A-family polymerase. Here, we present its crystal structure at a resolution of 2.9Å. P. falciparum apPOL, the first structural representative of a plastidic A-family polymerase, diverges from typical A-family members in two of three previously identified signature motifs and in a region not implicated by sequence. Moreover, apPOL has an additional N-terminal subdomain, the absence of which severely diminishes its 3' to 5' exonuclease activity. A compound known to be toxic to Plasmodium is a potent inhibitor of apPOL, suggesting that apPOL is a viable drug target. The structure provides new insights into the structural diversity of A-family polymerases and may facilitate structurally guided antimalarial drug design.

Fidelity of Nucleotide Incorporation by the RNA-Dependent RNA Polymerase from Poliovirus

Using poliovirus (PV) and its RNA-dependent RNA polymerase (RdRp) as our primary model system, we have advanced knowledge fundamental to the chemistry and fidelity of nucleotide addition by nucleic acid polymerase. Two fidelity checkpoints exist prior to nucleotide addition. The first toggles the enzyme between a nucleotide binding-occluded state and a nucleotide binding-competent state. The second represents an ensemble of conformational states of conserved structural motifs that permits retention of the incoming nucleotide in a state competent for phosphoryl transfer long enough for chemistry to occur. Nucleophilic attack of the alpha-phosphorous atom of the incoming nucleotide produces a pentavalent transition state, collapse of which is facilitated by protonation of the pyrophosphate leaving group by a general acid. All of the relevant conformational states of the enzyme are controlled by a network of interacting residues that permits remote-site residues to control active-site function. The current state of the art for PV RdRp enzymology is such that mechanisms governing fidelity of this enzyme can now be targeted genetically and chemically for development of attenuated viruses and antiviral agents, respectively. Application of the knowledge obtained with the PV RdRp to the development of vaccines and antivirals for emerging RNA viruses represents an important goal for the future.

A human transcription factor in search mode

Transcription factors (TF) can change shape to bind and recognize DNA, shifting the energy landscape from a weak binding, rapid search mode to a higher affinity recognition mode. However, the mechanism(s) driving this conformational change remains unresolved and in most cases high-resolution structures of the non-specific complexes are unavailable. Here, we investigate the conformational switch of the human mitochondrial transcription termination factor MTERF1, which has a modular, superhelical topology complementary to DNA. Our goal was to characterize the details of the non-specific search mode to complement the crystal structure of the specific binding complex, providing a basis for understanding the recognition mechanism. In the specific complex, MTERF1 binds a significantly distorted and unwound DNA structure, exhibiting a protein conformation incompatible with binding to B-form DNA. In contrast, our simulations of apo MTERF1 revealed significant flexibility, sampling structures with superhelical pitch and radius complementary to the major groove of B-DNA. Docking these structures to B-DNA followed by unrestrained MD simulations led to a stable complex in which MTERF1 was observed to undergo spontaneous diffusion on the DNA. Overall, the data support an MTERF1-DNA binding and recognition mechanism driven by intrinsic dynamics of the MTERF1 superhelical topology.

Substrate-induced DNA polymerase β activation

DNA polymerases and substrates undergo conformational changes upon forming protein-ligand complexes. These conformational adjustments can hasten or deter DNA synthesis and influence substrate discrimination. From structural comparison of binary DNA and ternary DNA-dNTP complexes of DNA polymerase β, several side chains have been implicated in facilitating formation of an active ternary complex poised for chemistry. Site-directed mutagenesis of these highly conserved residues (Asp-192, Arg-258, Phe-272, Glu-295, and Tyr-296) and kinetic characterization provides insight into the role these residues play during correct and incorrect insertion as well as their role in conformational activation. The catalytic efficiencies for correct nucleotide insertion for alanine mutants were wild type ∼ R258A > F272A ∼ Y296A > E295A > D192A. Because the efficiencies for incorrect insertion were affected to about the same extent for each mutant, the effects on fidelity were modest (<5-fold). The R258A mutant exhibited an increase in the single-turnover rate of correct nucleotide insertion. This suggests that the wild-type Arg-258 side chain generates a population of non-productive ternary complexes. Structures of binary and ternary substrate complexes of the R258A mutant and a mutant associated with gastric carcinomas, E295K, provide molecular insight into intermediate structural conformations not appreciated previously. Although the R258A mutant crystal structures were similar to wild-type enzyme, the open ternary complex structure of E295K indicates that Arg-258 stabilizes a non-productive conformation of the primer terminus that would decrease catalysis. Significantly, the open E295K ternary complex binds two metal ions indicating that metal binding cannot overcome the modified interactions that have interrupted the closure of the N-subdomain.

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.

Sustained active site rigidity during synthesis by human DNA polymerase μ

DNA polymerase μ (Pol μ) is the only template-dependent human DNA polymerase capable of repairing double-strand DNA breaks (DSBs) with unpaired 3' ends in nonhomologous end joining (NHEJ). To probe this function, we structurally characterized Pol μ's catalytic cycle for single-nucleotide incorporation. These structures indicate that, unlike other template-dependent DNA polymerases, Pol μ shows no large-scale conformational changes in protein subdomains, amino acid side chains or DNA upon dNTP binding or catalysis. Instead, the only major conformational change is seen earlier in the catalytic cycle, when the flexible loop 1 region repositions upon DNA binding. Pol μ variants with changes in loop 1 have altered catalytic properties and are partially defective in NHEJ. The results indicate that specific loop 1 residues contribute to Pol μ's unique ability to catalyze template-dependent NHEJ of DSBs with unpaired 3' ends.

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.

Magnetic isotope and magnetic field effects on the DNA synthesis

Magnetic isotope and magnetic field effects on the rate of DNA synthesis catalysed by polymerases β with isotopic ions (24)Mg(2+), (25)Mg(2+) and (26)Mg(2+) in the catalytic sites were detected. No difference in enzymatic activity was found between polymerases β carrying (24)Mg(2+) and (26)Mg(2+) ions with spinless, non-magnetic nuclei (24)Mg and (26)Mg. However, (25)Mg(2+) ions with magnetic nucleus (25)Mg were shown to suppress enzymatic activity by two to three times with respect to the enzymatic activity of polymerases β with (24)Mg(2+) and (26)Mg(2+) ions. Such an isotopic dependence directly indicates that in the DNA synthesis magnetic mass-independent isotope effect functions. Similar effect is exhibited by polymerases β with Zn(2+) ions carrying magnetic (67)Zn and non-magnetic (64)Zn nuclei, respectively. A new, ion-radical mechanism of the DNA synthesis is suggested to explain these effects. Magnetic field dependence of the magnesium-catalysed DNA synthesis is in a perfect agreement with the proposed ion-radical mechanism. It is pointed out that the magnetic isotope and magnetic field effects may be used for medicinal purposes (trans-cranial magnetic treatment of cognitive deceases, cell proliferation, control of the cancer cells, etc).

Perspective: pre-chemistry conformational changes in DNA polymerase mechanisms

In recent papers, there has been a lively exchange concerning theories for enzyme catalysis, especially the role of protein dynamics/pre-chemistry conformational changes in the catalytic cycle of enzymes. Of particular interest is the notion that substrate-induced conformational changes that assemble the polymerase active site prior to chemistry are required for DNA synthesis and impact fidelity (i.e., substrate specificity). High-resolution crystal structures of DNA polymerase β representing intermediates of substrate complexes prior to the chemical step are available. These structures indicate that conformational adjustments in both the protein and substrates must occur to achieve the requisite geometry of the reactive participants for catalysis. We discuss computational and kinetic methods to examine possible conformational change pathways that lead from the observed crystal structure intermediates to the final structures poised for chemistry. The results, as well as kinetic data from site-directed mutagenesis studies, are consistent with models requiring pre-chemistry conformational adjustments in order to achieve high fidelity DNA synthesis. Thus, substrate-induced conformational changes that assemble the polymerase active site prior to chemistry contribute to DNA synthesis even when they do not represent actual rate-determining steps for chemistry.

Catalytic effects of mutations of distant protein residues in human DNA polymerase β: theory and experiment

We carried out free-energy calculations and transient kinetic experiments for the insertion of the right (dC) and wrong (dA) nucleotides by wild-type (WT) and six mutant variants of human DNA polymerase β (Pol β). Since the mutated residues in the point mutants, I174S, I260Q, M282L, H285D, E288K, and K289M, were not located in the Pol β catalytic site, we assumed that the WT and its point mutants share the same dianionic phosphorane transition-state structure of the triphosphate moiety of deoxyribonucleotide 5'-triphosphate (dNTP) substrate. On the basis of this assumption, we have formulated a thermodynamic cycle for calculating relative dNTP insertion efficiencies, Ω = (k(pol)/K(D))(mut)/(k(pol)/K(D))(WT) using free-energy perturbation (FEP) and linear interaction energy (LIE) methods. Kinetic studies on five of the mutants have been published previously using different experimental conditions, e.g., primer-template sequences. We have performed a presteady kinetic analysis for the six mutants for comparison with wild-type Pol β using the same conditions, including the same primer/template DNA sequence proximal to the dNTP insertion site used for X-ray crystallographic studies. This consistent set of kinetic and structural data allowed us to eliminate the DNA sequence from the list of factors that can adversely affect calculated Ω values. The calculations using the FEP free energies scaled by 0.5 yielded 0.9 and 1.1 standard deviations from the experimental log Ω values for the insertion of the right and wrong dNTP, respectively. We examined a hybrid FEP/LIE method in which the FEP van der Waals term for the interaction of the mutated amino acid residue with its surrounding environment was replaced by the corresponding van der Waals term calculated using the LIE method, resulting in improved 0.4 and 1.0 standard deviations from the experimental log Ω values. These scaled FEP and FEP/LIE methods were also used to predict log Ω for R283A and R283L Pol β mutants.

Unfavorable electrostatic and steric interactions in DNA polymerase β E295K mutant interfere with the enzyme's pathway

Mutations in DNA polymerase β (pol β) have been associated with approximately 30% of human tumors. The E295K mutation of pol β has been linked to gastric carcinoma via interference with base excision repair. To interpret the different behavior of E295K as compared to wild-type pol β in atomic and energetic detail, we resolve a binary crystal complex of E295K at 2.5 Å and apply transition path sampling (TPS) to delineate the closing pathway of the E295K pol β mutant. Conformational changes are important components in the enzymatic pathway that lead to and ready the enzyme for the chemical reaction. Our analyses show that the closing pathway of E295K mutant differs from the wild-type pol β in terms of the individual transition states along the pathway, associated energies, and the active site conformation in the final closed form of the mutant. In particular, the closed state of E295K has a more distorted active site than the active site in the wild-type pol β. In addition, the total energy barrier in the conformational closing pathway is 65 ± 11 kJ/mol, much higher than that estimated for both correct (e.g., G:C) and incorrect (e.g., G:A) wild-type pol β systems (42 ± 8 and 45 ± 7 kJ/mol, respectively). In particular, the rotation of Arg258 is the rate-limiting step in the conformational pathway of E295K due to unfavorable electrostatic and steric interactions. The distorted active site in the closed relative to open state and the high energy barrier in the conformational pathway may explain in part why the E295K mutant is observed to be inactive. Interestingly, however, following the closing of the thumb but prior to the rotation of Arg258, the E295K mutant complex has a similar energy level as compared to the wild-type pol β. This suggests that the E295K mutant may associate with DNA with similar affinity, but it may be hampered in continuing the process of chemistry. Supporting experimental data come from the observation that the catalytic activity of wild-type pol β is hampered when E295K is present: this may arise from the competition between E295K and wild-type enzyme for the DNA. These combined results suggest that the low insertion efficiency of E295K mutant as compared to wild-type pol β may be related to a closed form distorted by unfavorable electrostatic and steric interactions between Arg258 and other key residues. The active site is thus less competent for proceeding to the chemical reaction, which may also involve a higher reaction barrier than the wild-type or may not be possible in this mutant. Our analysis also suggests further experiments for other mutants to test the above hypothesis and dissect the roles of steric and electrostatic factors on enzyme behavior.

Prechemistry versus preorganization in DNA replication fidelity

The molecular origin of nucleotide insertion catalysis and fidelity of DNA polymerases is explored by means of computational simulations. Special attention is paid to the examination of the validity of proposals that invoke prechemistry effects, checkpoints concepts, and dynamical effects. The simulations reproduce the observed fidelity in Pol β, starting with the relevant observed X-ray structures of the complex with the right (R) and wrong (W) nucleotides. The generation of free energy surfaces for the R and W systems also allowed us to analyze different proposals about the origin of the fidelity and to reach several important conclusions. It is found that the potential of mean force (PMF) obtained by proper sampling does not support QM/MM-based proposals of a large barrier before the prechemistry state. Furthermore, examination of dynamical proposals by the renormalization approach indicates that the motions from open to close configurations do not contribute to catalysis or fidelity. Finally we discuss and analyze the induced fit concept and show that, despite its importance, it does not explain fidelity. That is, the fidelity is apparently due to the change in the preorganization of the chemical site, as a result of the relaxation of the binding site upon binding of the incorrect nucleotide. Finally and importantly, since the issue is the barrier associated with the enzyme-substrate (ES)/DNA complex at the chemical transition state and not the path to this complex formation (unless this path involves rate determining steps), it is also not useful to invoke checkpoints while discussing fidelity.

DNA polymerase activity at the single-molecule level

DNA polymerases are essential enzymes responsible for replication and repair of DNA in all organisms. To replicate DNA with high fidelity, DNA polymerases must select the correct incoming nucleotide substrate during each cycle of nucleotide incorporation, in accordance with the templating base. When an incorrect nucleotide is sometimes inserted, the polymerase uses a separate 3'→5' exonuclease to remove the misincorporated base (proofreading). Large conformational rearrangements of the polymerase-DNA complex occur during both the nucleotide incorporation and proofreading steps. Single-molecule fluorescence spectroscopy provides a unique tool for observation of these dynamic conformational changes in real-time, without the need to synchronize a population of DNA-protein complexes.

A model for the dynamics of mammalian family X DNA polymerases

Based on available structural studies, a model is presented for polymerization dynamics of mammalian family X DNA polymerases, including polymerases β, λ, μ, and terminal deoxynucleotidyl transferase (TdT). Using the model, distinct polymerization activities and processivities of the four polymerases acting on different forms of DNA substrate are analyzed and studied theoretically. A "gradient" of template dependence of polymerases β, λ, μ, and TdT is well explained. The much higher occurrence frequencies of the -1 frameshift DNA synthesis by pols λ and μ than that by pol β are well explained. The theoretical results on the polymerization processivities are also in agreement with the available experimental data.

DNA pol λ's extraordinary ability to stabilize misaligned DNA

DNA polymerases have the venerable task of maintaining genome stability during DNA replication and repair. Errors, nonetheless, occur with error propensities that are polymerase specific. For example, DNA polymerase λ (pol λ) generates single-base deletions through template-strand slippage within short repetitive DNA regions much more readily than does the closely related polymerase β (pol β). Here we present in silico evidence to help interpret pol λ's greater tendency for deletion errors than pol β by its more favorable protein/DNA electrostatic interactions immediately around the extrahelical nucleotide on the template strand. Our molecular dynamics and free energy analyses suggest that pol λ provides greater stabilization to misaligned DNA than aligned DNA. Our study of several pol λ mutants of Lys544 (Ala, Phe, Glu) probes the interactions between the extrahelical nucleotide and the adjacent Lys544 to show that the charge of the 544 residue controls stabilization of the DNA misalignment. In addition, we identify other thumb residues (Arg538, Lys521, Arg517, and Arg514) that play coordinating roles in stabilizing pol λ's interactions with misaligned DNA. Interestingly, their aggregate stabilization effect is more important than that of any one component residue, in contrast to aligned DNA systems, as we determined from mutations of these key residues and energetic analyses. No such comparable network of stabilizing misaligned DNA exists in pol β. Evolutionary needs for DNA repair on substrates with minimal base-pairing, such as those encountered by pol λ in the non-homologous end-joining pathway, may have been solved by a greater tolerance to deletion errors. Other base-flipping proteins share similar binding properties and motions for extrahelical nucleotides.

Conformational coupling, bridge helix dynamics and active site dehydration in catalysis by RNA polymerase

Molecular dynamics simulation of Thermus thermophilus (Tt) RNA polymerase (RNAP) in a catalytic conformation demonstrates that the active site dNMP-NTP base pair must be substantially dehydrated to support full active site closing and optimum conditions for phosphodiester bond synthesis. In silico mutant beta R428A RNAP, which was designed based on substitutions at the homologous position (Rpb2 R512) of Saccharomyces cerevisiae (Sc) RNAP II, was used as a reference structure to compare to Tt RNAP in simulations. Long range conformational coupling linking a dynamic segment of the bridge alpha-helix, the extended fork loop, the active site, and the trigger loop-trigger helix is apparent and adversely affected in beta R428A RNAP. Furthermore, bridge helix bending is detected in the catalytic structure, indicating that bridge helix dynamics may regulate phosphodiester bond synthesis as well as translocation. An active site "latch" assembly that includes a key trigger helix residue Tt beta' H1242 and highly conserved active site residues beta E445 and R557 appears to help regulate active site hydration/dehydration. The potential relevance of these observations in understanding RNAP and DNAP induced fit and fidelity is discussed.

The mechanism of the translocation step in DNA replication by DNA polymerase I: a computer simulation analysis

High-fidelity DNA polymerases copy DNA rapidly and accurately by adding correct deoxynucleotide triphosphates to a growing primer strand of DNA. Following nucleotide incorporation, a series of conformational changes translocate the DNA substrate by one base pair step, readying the polymerase for the next round of incorporation. Molecular dynamics simulations indicate that the translocation consists globally of a polymerase fingers-opening transition, followed by the DNA displacement and the insertion of the template base into the preinsertion site. They also show that the pyrophosphate release facilitates the opening transition and that the universally conserved Y714 plays a key role in coupling polymerase opening to DNA translocation. The transition involves several metastable intermediates in one of which the O helix is bent in the vicinity of G711. Completion of the translocation appears to require a gating motion of the O1 helix, perhaps facilitated by the presence of G715. These roles are consistent with the high level of conservation of Y714 and the two glycine residues at these positions. It is likely that a corresponding mechanism is applicable to other polymerases.

Contribution of the reverse rate of the conformational step to polymerase beta fidelity

A complete understanding of the kinetic mechanism of fidelity requires comparison of correct and incorrect dNTP incorporation pathways in both the forward and reverse directions. The studies presented here focus on the dNTP-induced conformational step, which has historically been proposed by many to be the major determinant of fidelity. As it was recently highlighted [Tsai, Y. C., and Johnson, K. A. (2006) Biochemistry 45, 9675-9687], chemistry can be the slowest step in the forward direction of the correct dNTP incorporation pathway, yet the corresponding microscopic rate constant would not contribute toward fidelity in the case when the reverse rate of the conformational step is slower than chemistry. Here we use a stopped-flow technique to directly measure the reverse rate of the conformational step in the DNA polymerase beta (Pol beta) kinetic pathway. Extensive pre-steady-state kinetic studies presented include the utilization of 2-aminopurine-labeled DNA substrates, 2-aminopurine nucleotide triphosphate, a nonhydrolyzable nucleotide analogue dAMPCPP, and a rapid sequential mixing reaction scheme. Additionally, the effect of mismatched dNTPs, various metal ions, and the presence of the 3'-terminal hydroxyl group of the primer on the rate of the reverse "opening" conformational step were analyzed. Our analyses indicate that reverse "opening" is drastically facilitated in the presence of mismatched ternary complexes, which is in agreement with the hypothesis that the ternary complex is destabilized by the presence of incorrect dNTP. By analysis of the relative magnitudes of chemistry and reverse "opening" in the presence of both matched and mismatched matched ternary complexes, this work further validates that, for Pol beta, fidelity is dictated by the differences in free energy required to reach the highest energy transition state of the chemical step.

Insights into the high fidelity of a DNA polymerase I mutant

Mutants of DNA polymerase I from Thermus aquaticus (Taq) with higher fidelity compared to the wild type enzyme were identified in an earlier study by Summerer et al. (Angew Chem Int Ed 44:4712-4715, 2005). Here, one of these mutants, PLQ (consensus residues 879-881), was analysed using molecular dynamics simulations. This was done by calculating the structures of the ternary complex comprising the enzyme, the DNA primer and template as well as the incoming nucleotide before the chemical reaction for the Watson-Crick and different mismatched base pairings. The results show that the high fidelity of the mutant can be explained partly by different specific interactions between the amino acids of the enzyme and the DNA primer end as well as, in some mismatches, a displacement of the primer relative to the incoming deoxyribonucleoside triphosphate and the catalytic magnesium ion. This displacement is facilitated by reduced steric interactions between the enzyme and the DNA.

Quantum mechanics/molecular mechanics investigation of the chemical reaction in Dpo4 reveals water-dependent pathways and requirements for active site reorganization

The nucleotidyl-transfer reaction coupled with the conformational transitions in DNA polymerases is critical for maintaining the fidelity and efficiency of DNA synthesis. We examine here the possible reaction pathways of a Y-family DNA polymerase, Sulfolobus solfataricus DNA polymerase IV (Dpo4), for the correct insertion of dCTP opposite 8-oxoguanine using the quantum mechanics/molecular mechanics (QM/MM) approach, both from a chemistry-competent state and a crystal closed state. The latter examination is important for understanding pre-chemistry barriers to interpret the entire enzyme mechanism, since the crystal closed state is not an ideal state for initiating the chemical reaction. The most favorable reaction path involves initial deprotonation of O3'H via two bridging water molecules to O1A, overcoming an overall potential energy barrier of approximately 20.0 kcal/mol. The proton on O1A-P(alpha) then migrates to the gamma-phosphate oxygen of the incoming nucleotide as O3' attacks P(alpha), and the P(alpha)-O3A bond breaks. The other possible pathway in which the O3'H proton is transferred directly to O1A on P(alpha) has an overall energy barrier of 25.0 kcal/mol. In both reaction paths, the rate-limiting step is the initial deprotonation, and the trigonal-bipyramidal configuration for P(alpha) occurs during the concerted bond formation (O3'-P(alpha)) and breaking (P(alpha)-O3A), indicating the associative nature of the chemical reaction. In contrast, the Dpo4/DNA complex with an imperfect active-site geometry corresponding to the crystal state must overcome a much higher activation energy barrier (29.0 kcal/mol) to achieve a tightly organized site due to hindered O3'H deprotonation stemming from larger distances and distorted conformation of the proton acceptors. This significant difference demonstrates that the pre-chemistry reorganization in Dpo4 costs approximately 4.0 to 9.0 kcal/mol depending on the primer terminus environment. Compared to the higher fidelity DNA polymerase beta from the X-family, Dpo4 has a higher chemical reaction barrier (20.0 vs 15.0 kcal/mol) due to the more solvent-exposed active site.

Exploring the role of large conformational changes in the fidelity of DNA polymerase beta

The relationships between the conformational landscape, nucleotide insertion catalysis and fidelity of DNA polymerase beta are explored by means of computational simulations. The simulations indicate that the transition states for incorporation of right (R) and wrong (W) nucleotides reside in substantially different protein conformations. The protein conformational changes that reproduce the experimentally observed fidelity are significantly larger than the small rearrangements that usually accompany motions from the reactant state to the transition state in common enzymatic reactions. Once substrate binding has occurred, different constraints imposed on the transition states for insertion of R and W nucleotides render it highly unlikely that both transition states can occur in the same closed structure, because the predicted fidelity would then be many orders of magnitude too large. Since the conformational changes reduce the transition state energy of W incorporation drastically they decrease fidelity rather than increase it. Overall, a better agreement with experimental data is attained when the R is incorporated through a transition state in a closed conformation and W is incorporated through a transition state in one or perhaps several partially open conformations. The generation of free energy surfaces for R and W also allow us to analyze proposals about the relationship between induced fit and fidelity.

DNA polymerase beta fidelity: halomethylene-modified leaving groups in pre-steady-state kinetic analysis reveal differences at the chemical transition state

The mechanism of DNA polymerase beta-catalyzed nucleotidyl transfer consists of chemical steps involving primer 3' OH deprotonation, nucleophilic attack, and pyrophosphate leaving-group elimination, preceded by dNTP binding which induces a large-amplitude conformational change for Watson-Crick nascent base pairs. Ambiguity in the nature of the rate-limiting step and active-site structural differences between correct and incorrect base-paired transition states remain obstacles to understanding DNA replication fidelity. Analogues of dGTP where the beta-gamma bridging oxygen is replaced with fluorine-substituted methylene groups have been shown to probe the contribution of leaving-group elimination to the overall catalytic rate (Biochemistry 46, 461-471). Here, the analysis is expanded substantially to include a broad range of halogen substituents with disparate steric and electronic properties. Evaluation of linear free energy relationships for incorporation of dGTP analogues opposite either template base C or T reveals a strong correlation of log(kpol) to leaving group pKa. Significantly different kpol behavior is observed with a subset of the analogues, with magnitude dependent on the identity of the nascent base pair. This observation, and the absence of an analogous effect on ground state analogue binding (Kd values), points to active-site structural differences at the chemical transition state. Reduced catalysis with bulky halo-containing substrates is manifested in the fidelity of T-G incorporation, where the CCl2-bridging analogue shows a 27-fold increase in fidelity over the natural dGTP. Solvent pH and deuterium isotope-effect data are also used to evaluate mechanistic differences between correct and mispaired incorporation.

Nucleotide-induced DNA polymerase active site motions accommodating a mutagenic DNA intermediate

DNA polymerases occasionally insert the wrong nucleotide. For this error to become a mutation, the mispair must be extended. We report a structure of DNA polymerase beta (pol beta) with a DNA mismatch at the boundary of the polymerase active site. The structure of this complex indicates that the templating adenine of the mispair stacks with the primer terminus adenine while the templating (coding) cytosine is flipped out of the DNA helix. Soaking the crystals of the binary complex with dGTP resulted in crystals of a ternary substrate complex. In this case, the templating cytosine is observed within the DNA helix and forms Watson-Crick hydrogen bonds with the incoming dGTP. The adenine at the primer terminus has rotated into a syn-conformation to interact with the opposite adenine in a planar configuration. Yet, the 3'-hydroxyl on the primer terminus is out of position for efficient nucleotide insertion.

Single-molecule and ensemble fluorescence assays for a functionally important conformational change in T7 DNA polymerase

We report fluorescence assays for a functionally important conformational change in bacteriophage T7 DNA polymerase (T7 pol) that use the environmental sensitivity of a Cy3 dye attached to a DNA substrate. An increase in fluorescence intensity of Cy3 is observed at the single-molecule level, reflecting a conformational change within the T7 pol ternary complex upon binding of a dNTP substrate. This fluorescence change is believed to reflect the closing of the T7 pol fingers domain, which is crucial for polymerase function. The rate of the conformational change induced by a complementary dNTP substrate was determined by both conventional stopped-flow and high-time-resolution continuous-flow fluorescence measurements at the ensemble-averaged level. The rate of this conformational change is much faster than that of DNA synthesis but is significantly reduced for noncomplementary dNTPs, as revealed by single-molecule measurements. The high level of selectivity of incoming dNTPs pertinent to this conformational change is a major contributor to replicative fidelity.

The X family portrait: structural insights into biological functions of X family polymerases

The mammalian family X DNA polymerases (DNA polymerases beta, lambda, mu, and TdT) contribute to base excision repair and double-strand break repair by virtue of their ability to fill short gaps in DNA. Structural information now exists for all four of these enzymes, making this the first mammalian polymerase family whose structural portrait is complete. Here we consider how distinctive structural features of these enzymes contribute to their biological functions in vivo.

A unified kinetic mechanism applicable to multiple DNA polymerases

After extensive studies spanning over half a century, there is little consensus on the kinetic mechanism of DNA polymerases. Using stopped-flow fluorescence assays for mammalian DNA polymerase beta (Pol beta), we have previously identified a fast fluorescence transition corresponding to conformational closing, and a slow fluorescence transition matching the rate of single-nucleotide incorporation. Here, by varying pH and buffer viscosity, we have decoupled the rate of single-nucleotide incorporation from the rate of the slow fluorescence transition, thus confirming our previous hypothesis that this transition represents a conformational event after chemistry, likely subdomain reopening. Analysis of an R258A mutant indicates that rotation of the Arg258 side chain is not rate-limiting in the overall kinetic pathway of Pol beta, yet is kinetically significant in subdomain reopening. We have extended our kinetic analyses to a high-fidelity polymerase, Klenow fragment (KF), and a low-fidelity polymerase, African swine fever virus DNA polymerase X (Pol X), and showed that they follow the same kinetic mechanism as Pol beta, while differing in relative rates of single-nucleotide incorporation and the putative conformational reopening. Our data suggest that the kinetic mechanism of Pol beta is not an exception among polymerases, and furthermore, its delineated kinetic mechanism lends itself as a platform for comparison of the kinetic properties of different DNA polymerases and their mutants.

Distinct energetics and closing pathways for DNA polymerase beta with 8-oxoG template and different incoming nucleotides

Background

8-Oxoguanine (8-oxoG) is a common oxidative lesion frequently encountered by DNA polymerases such as the repair enzyme DNA polymerase β (pol β). To interpret in atomic and energetic detail how pol β processes 8-oxoG, we apply transition path sampling to delineate closing pathways of pol β 8-oxoG complexes with dCTP and dATP incoming nucleotides and compare the results to those of the nonlesioned G:dCTP and G:dATPanalogues.

Results

Our analyses show that the closing pathways of the 8-oxoG complexes are different from one another and from the nonlesioned analogues in terms of the individual transition states along each pathway, associated energies, and the stability of each pathway's closed state relative to the corresponding open state. In particular, the closed-to-open state stability difference in each system establishes a hierarchy of stability (from high to low) as G:C > 8-oxoG:C > 8-oxoG:A > G:A, corresponding to -3, -2, 2, 9 k_BT, respectively. This hierarchy of closed state stability parallels the experimentally observed processing efficiencies for the four pairs. Network models based on the calculated rate constants in each pathway indicate that the closed species are more populated than the open species for 8-oxoG:dCTP, whereas the opposite is true for 8-oxoG:dATP.

Conclusion

These results suggest that the lower insertion efficiency (larger K_m) for dATP compared to dCTP opposite 8-oxoG is caused by a less stable closed-form of pol β, destabilized by unfavorable interactions between Tyr271 and the mispair. This stability of the closed vs. open form can also explain the higher insertion efficiency for 8-oxoG:dATP compared to the nonlesioned G:dATP pair, which also has a higher overall conformational barrier. Our study offers atomic details of the complexes at different states, in addition to helping interpret the different insertion efficiencies of dATP and dCTP opposite 8-oxoG and G.

Differing conformational pathways before and after chemistry for insertion of dATP versus dCTP opposite 8-oxoG in DNA polymerase beta

To elucidate how human DNA polymerase beta (pol beta) discriminates dATP from dCTP when processing 8-oxoguanine (8-oxoG), we analyze a series of dynamics simulations before and after the chemical step with dATP and dCTP opposite an 8-oxoG template started from partially open complexes of pol beta. Analyses reveal that the thumb closing of pol beta before chemistry is hampered when the incorrect nucleotide dATP is bound opposite 8-oxoG; the unfavorable interaction between active-site residue Tyr(271) and dATP that causes an anti to syn change in the 8-oxoG (syn):dATP complex explains this slow motion, in contrast to the 8-oxoG (anti):dCTP system. Such differences in conformational pathways before chemistry for mismatched versus matched complexes help explain the preference for correct insertion across 8-oxoG by pol beta. Together with reference studies with a nonlesioned G template, we propose that 8-oxoG leads to lower efficiency in pol beta's incorporation of dCTP compared with G by affecting the requisite active-site geometry for the chemical reaction before chemistry. Furthermore, because the active site is far from ready for the chemical reaction after partial closing or even full thumb closing, we suggest that pol beta is tightly controlled not only by the chemical step but also by a closely related requirement for subtle active-site rearrangements after thumb movement but before chemistry.

Regulation of DNA repair fidelity by molecular checkpoints: "gates" in DNA polymerase beta's substrate selection

With an increasing number of structural, kinetic, and modeling studies of diverse DNA polymerases in various contexts, a complex dynamical view of how atomic motions might define molecular "gates" or checkpoints that contribute to polymerase specificity and efficiency is emerging. Such atomic-level information can offer insights into rate-limiting conformational and chemical steps to help piece together mechanistic views of polymerases in action. With recent advances, modeling and dynamics simulations, subject to the well-appreciated limitations, can access transition states and transient intermediates along a reaction pathway, both conformational and chemical, and such information can help bridge the gap between experimentally determined equilibrium structures and mechanistic enzymology data. Focusing on DNA polymerase beta (pol beta), we present an emerging view of the geometric, energetic, and dynamic selection criteria governing insertion rate and fidelity mechanisms of DNA polymerases, as gleaned from various computational studies and based on the large body of existing kinetic and structural data. The landscape of nucleotide insertion for pol beta includes conformational changes, prechemistry, and chemistry "avenues", each with a unique deterministic or stochastic pathway that includes checkpoints for selective control of nucleotide insertion efficiency. For both correct and incorrect incoming nucleotides, pol beta's conformational rearrangements before chemistry include a cascade of slow and subtle side chain rearrangements, followed by active site adjustments to overcome higher chemical barriers, which include critical ion-polymerase geometries; this latter notion of a prechemistry avenue fits well with recent structural and NMR data. The chemical step involves an associative mechanism with several possibilities for the initial proton transfer and for the interaction among the active site residues and bridging water molecules. The conformational and chemical events and associated barriers define checkpoints that control enzymatic efficiency and fidelity. Understanding the nature of such active site rearrangements can facilitate interpretation of existing data and stimulate new experiments that aim to probe enzyme features that contribute to fidelity discrimination across various polymerases via such geometric, dynamic, and energetic selection criteria.

Correct and incorrect nucleotide incorporation pathways in DNA polymerase beta

Tracking the structural and energetic changes in the pathways of DNA replication and repair is central to the understanding of these important processes. Here we report favorable mechanisms of the polymerase-catalyzed phosphoryl transfer reactions corresponding to correct and incorrect nucleotide incorporations in the DNA by using a novel protocol involving energy minimizations, dynamics simulations, quasi-harmonic free energy calculations, and mixed quantum mechanics/molecular mechanics dynamics simulations. Though the pathway proposed may not be unique and invites variations, geometric and energetic arguments support the series of transient intermediates in the phosphoryl transfer pathways uncovered here for both the G:C and G:A systems involving a Grotthuss hopping mechanism of proton transfer between water molecules and the three conserved aspartate residues in pol beta's active-site. In the G:C system, the rate-limiting step is the initial proton hop with a free energy of activation of at least 17 kcal/mol, which corresponds closely to measured k(pol) values. Fidelity discrimination in pol beta can be explained by a significant loss of stability of the closed ternary complex of the enzyme in the G:A system and much higher activation energy of the initial step of nucleophilic attack, namely deprotonation of terminal DNA primer O3'H group. Thus, subtle differences in the enzyme active-site between matched and mismatched base pairs generate significant differences in catalytic performance.

Energy analysis of chemistry for correct insertion by DNA polymerase beta

X-ray crystallographic structures of human DNA polymerase beta with nonhydrolyzable analogs containing all atoms in the active site required for catalysis provide a secure starting point for a theoretical analysis (quantum mechanics/molecular mechanics) of the mechanism of chemistry without biasing of modeling assumptions as required in previous studies. These structures provide the basis for a detailed quantum mechanics/molecular mechanics study of the path for the complete transfer of a monophosphate nucleoside donor to the sugar acceptor in the active site. The reaction is largely associative with the main energetic step preceded by proton transfer from the terminal primer deoxyribose O3' to Asp-256. The key residues that provide electrostatic stabilization of the transition state are identified and compared with those identified by mutational studies.

Sequential side-chain residue motions transform the binary into the ternary state of DNA polymerase lambda

The nature of conformational transitions in DNA polymerase lambda (pol lambda), a low-fidelity DNA repair enzyme in the X-family that fills short nucleotide gaps, is investigated. Specifically, to determine whether pol lambda has an induced-fit mechanism and open-to-closed transition before chemistry, we analyze a series of molecular dynamics simulations from both the binary and ternary states before chemistry, with and without the incoming nucleotide, with and without the catalytic Mg(2+) ion in the active site, and with alterations in active site residues Ile(492) and Arg(517). Though flips occurred for several side-chain residues (Ile(492), Tyr(505), Phe(506)) in the active site toward the binary (inactive) conformation and partial DNA motion toward the binary position occurred without the incoming nucleotide, large-scale subdomain motions were not observed in any trajectory from the ternary complex regardless of the presence of the catalytic ion. Simulations from the binary state with incoming nucleotide exhibit more thumb subdomain motion, particularly in the loop containing beta-strand 8 in the thumb, but closing occurred only in the Ile(492)Ala mutant trajectory started from the binary state with incoming nucleotide and both ions. Further connections between active site residues and the DNA position are also revealed through our Ile(492)Ala and Arg(517)Ala mutant studies. Our combined studies suggest that while pol lambda does not demonstrate large-scale subdomain movements as DNA polymerase beta (pol beta), significant DNA motion exists, and there are sequential subtle side chain and other motions-associated with Arg(514), Arg(517), Ile(492), Phe(506), Tyr(505), the DNA, and again Arg(514) and Arg(517)-all coupled to active site divalent ions and the DNA motion. Collectively, these motions transform pol lambda to the chemistry-competent state. Significantly, analogs of these residues in pol beta (Lys(280), Arg(283), Arg(258), Phe(272), and Tyr(271), respectively) have demonstrated roles in determining enzyme efficiency and fidelity. As proposed for pol beta, motions of these residues may serve as gate-keepers by controlling the evolution of the reaction pathway before the chemical reaction.

Magnesium-induced assembly of a complete DNA polymerase catalytic complex

The molecular details of the nucleotidyl transferase reaction have remained speculative, as strategies to trap catalytic intermediates for structure determination utilize substrates lacking the primer terminus 3'-OH and catalytic Mg2+, resulting in an incomplete and distorted active site geometry. Since the geometric arrangement of these essential atoms will impact chemistry, structural insight into fidelity strategies has been hampered. Here, we present a crystal structure of a precatalytic complex of a DNA polymerase with bound substrates that include the primer 3'-OH and catalytic Mg2+. This catalytic intermediate was trapped with a nonhydrolyzable deoxynucleotide analog. Comparison with two new structures of DNA polymerase beta lacking the 3'-OH or catalytic Mg2+ is described. These structures provide direct evidence that both atoms are required to achieve a proper geometry necessary for an in-line nucleophilic attack of O3' on the alphaP of the incoming nucleotide.

Prostate-cancer-associated I260M variant of DNA polymerase beta is a sequence-specific mutator

Studies show that 30% of 189 tumors sequenced to date express variants of the polymerase beta (pol beta) protein that are not present in normal tissue. This raises the possibility that variants of pol beta might be linked to the etiology of cancer. Here, we characterize the I260M prostate-cancer-associated variant of pol beta. Ile260 is a key residue of the hydrophobic hinge that is important for the closing of the polymerase. In this study, we demonstrate that the I260M variant is a sequence context-dependent mutator polymerase. Specifically, I260M is a mutator for misalignment-mediated errors in dipyrimidine sequences. I260M is also a low-fidelity polymerase with regard to the induction of transversions within specific sequence contexts. Our results suggest that the hinge influences the geometry of the DNA within the polymerase active site that is important for accurate DNA synthesis. Importantly, characterization of the I260M variant shows that it has a functional phenotype that could be linked to the etiology or malignant progression of human cancer.

Mismatch-induced conformational distortions in polymerase beta support an induced-fit mechanism for fidelity

Molecular dynamics simulations of DNA polymerase (pol) beta complexed with different incorrect incoming nucleotides (G x G, G x T, and T x T template base x incoming nucleotide combinations) at the template-primer terminus are analyzed to delineate structure-function relationships for aberrant base pairs in a polymerase active site. Comparisons, made to pol beta structure and motions in the presence of a correct base pair, are designed to gain atomically detailed insights into the process of nucleotide selection and discrimination. In the presence of an incorrect incoming nucleotide, alpha-helix N of the thumb subdomain believed to be required for pol beta's catalytic cycling moves toward the open conformation rather than the closed conformation as observed for the correct base pair (G x C) before the chemical reaction. Correspondingly, active-site residues in the microenvironment of the incoming base are in intermediate conformations for non-Watson-Crick pairs. The incorrect incoming nucleotide and the corresponding template residue assume distorted conformations and do not form Watson-Crick bonds. Furthermore, the coordination number and the arrangement of ligands observed around the catalytic and nucleotide binding magnesium ions are mismatch specific. Significantly, the crucial nucleotidyl transferase reaction distance (P(alpha)-O3') for the mismatches between the incoming nucleotide and the primer terminus is not ideally compatible with the chemical reaction of primer extension that follows these conformational changes. Moreover, the extent of active-site distortion can be related to experimentally determined rates of nucleotide misincorporation and to the overall energy barrier associated with polymerase activity. Together, our studies provide structure-function insights into the DNA polymerase-induced constraints (i.e., alpha-helix N conformation, DNA base pair bonding, conformation of protein residues in the vicinity of dNTP, and magnesium ions coordination) during nucleotide discrimination and pol beta-nucleotide interactions specific to each mispair and how they may regulate fidelity. They also lend further support to our recent hypothesis that additional conformational energy barriers are involved following nucleotide binding but prior to the chemical reaction.

Activities and mechanism of DNA polymerase beta

DNA polymerase beta plays an essential role in the base excision repair pathway necessary to cleanse the genome of simple base lesions and abasic sites. Abasic sites arise in DNA from spontaneous base loss (depurination) and DNA-damage specific glycosylases that hydrolyze the N-glycosidic bond between the deoxyribose and the damaged base. DNA polymerase beta contributes two enzymatic activities: DNA synthesis and deoxyribose-phosphate removal through nucleotidyl transferase and lyase mechanisms, respectively. The active site for each of these activities resides on a distinct domain of the protein: 31-kDa polymerase domain and amino-terminal 8-kDa lyase domain. The simple organization of each domain and the ability to assay each activity have hastened our understanding of the faithful replication of DNA during repair synthesis and the flux of intermediates through single nucleotide base excision repair and its alternate pathways.

Environmental Effect on the Relative Contribution of the Charge-Transfer Mechanisms within a Short DNA Sequence

Time evolution of the charge-transfer site population is studied in a short DNA sequence to determine the type of governing charge-transfer mechanism. The system consists of a 5'-GAGGG-3' nucleobase sequence coupled with a dissipative bath that represents the DNA phosphate backbone and solvents. Relative contribution of transfer mechanisms to the whole charge-transfer process has been obtained using the on-the-fly filtered propagator functional path integral method with the density matrix decomposition. Partial density matrixes of the incoherent hopping and coherent superexchange pathways as well as the full reduced density matrix have been evaluated and discussed for both debye and ohmic baths. It was found that the relative contribution of the transfer mechanisms is rather sensitive to the frequency-dependent environmental description.

Subtle but variable conformational rearrangements in the replication cycle of Sulfolobus solfataricus P2 DNA polymerase IV (Dpo4) may accommodate lesion bypass

The possible conformational changes of DNA polymerase IV (Dpo4) before and after the nucleotidyl-transfer reaction are investigated at the atomic level by dynamics simulations to gain insight into the mechanism of low-fidelity polymerases and identify slow and possibly critical steps. The absence of significant conformational changes in Dpo4 before chemistry when the incoming nucleotide is removed supports the notion that the "induced-fit" mechanism employed to interpret fidelity in some replicative and repair DNA polymerases does not exist in Dpo4. However, significant correlated movements in the little finger and finger domains, as well as DNA sliding and subtle catalytic-residue rearrangements, occur after the chemical reaction when both active-site metal ions are released. Subsequently, Dpo4's little finger grips the DNA through two arginine residues and pushes it forward. These metal ion correlated movements may define subtle, and possibly characteristic, conformational adjustments that operate in some Y-family polymerase members in lieu of the prominent subdomain motions required for catalytic cycling in other DNA polymerases like polymerase beta. Such subtle changes do not easily provide a tight fit for correct incoming substrates as in higher-fidelity polymerases, but introduce in low-fidelity polymerases different fidelity checks as well as the variable conformational-mobility potential required to bypass different lesions.

Motions of the fingers subdomain of klentaq1 are fast and not rate limiting: implications for the molecular basis of fidelity in DNA polymerases

Various kinetic studies on nucleotide incorporation by DNA polymerases have established that a rate-limiting step occurs that is crucial in the mechanism of discrimination between correct versus incorrect nucleotide. Crystallographic studies have indicated that this step may be due to a large open-to-closed conformational transition affecting the fingers subdomain. However, there is no direct evidence to support this hypothesis. In order to investigate whether or not the open-to-closed conformational transition affecting the fingers subdomain is rate limiting, we have developed a fluorescence resonance energy transfer (FRET) system, which monitors motions of the fingers subdomain. We establish that the closing of the fingers subdomain is significantly faster than the kinetically determined rate-limiting step. We propose that the rate-limiting step occurs after the closing of the fingers subdomain and is caused by local reorganization events in the active site.

In silico studies of the African swine fever virus DNA polymerase X support an induced-fit mechanism

The African swine fever virus DNA polymerase X (pol X), a member of the X family of DNA polymerases, is thought to be involved in base excision repair. Kinetics data indicate that pol X catalyzes DNA polymerization with low fidelity, suggesting a role in viral mutagenesis. Though pol X lacks the fingers domain that binds the DNA in other members of the X family, it binds DNA tightly. To help interpret details of this interaction, molecular dynamics simulations of free pol X at different salt concentrations and of pol X bound to gapped DNA, in the presence and in the absence of the incoming nucleotide, are performed. Anchors for the simulations are two NMR structures of pol X without DNA and a model of one NMR structure plus DNA and incoming nucleotide. Our results show that, in its free form, pol X can exist in two stable conformations that interconvert to one another depending on the salt concentration. When gapped double stranded DNA is introduced near the active site, pol X prefers an open conformation, regardless of the salt concentration. Finally, under physiological conditions, in the presence of both gapped DNA and correct incoming nucleotide, and two divalent ions, the thumb subdomain of pol X undergoes a large conformational change, closing upon the DNA. These results predict for pol X a substrate-induced conformational change triggered by the presence of DNA and the correct incoming nucleotide in the active site, as in DNA polymerase beta. The simulations also suggest specific experiments (e.g., for mutants Phe-102Ala, Val-120Gly, and Lys-85Val that may reveal crucial DNA binding and active-site organization roles) to further elucidate the fidelity mechanism of pol X.

Fidelity discrimination in DNA polymerase beta: differing closing profiles for a mismatched (G:A) versus matched (G:C) base pair

Understanding fidelity-the faithful replication or repair of DNA by polymerases-requires tracking of the structural and energetic changes involved, including the elusive transient intermediates, for nucleotide incorporation at the template/primer DNA junction. We report, using path sampling simulations and a reaction network model, strikingly different transition states in DNA polymerase beta's conformational closing for correct dCTP versus incorrect dATP incoming nucleotide opposite a template G. The cascade of transition states leads to differing active-site assembly processes toward the "two-metal-ion catalysis" geometry. We demonstrate that these context-specific pathways imply different selection processes: while active-site assembly occurs more rapidly with the correct nucleotide and leads to primer extension, the enzyme remains open longer, has a more transient closed state, and forms product more slowly when an incorrect nucleotide is present. Our results also suggest that the rate-limiting step in pol beta's conformational closing is not identical to that for overall nucleotide insertion and that the rate-limiting step in the overall nucleotide incorporation process for matched as well as mismatched systems occurs after the closing conformational change.

Frequent mutation related with overexpression of DNA polymerase beta in primary tumors and precancerous lesions of human stomach

To explore whether DNA polymerase beta (pol beta) contributes to the malignant transformation of gastric mucosa, we examined pol beta in gastric tumor cell lines, primary tumors and precancerous lesions. Point mutations of pol beta were detected in 6 of 13 cell lines and 23 of 104 tissues including 35.0% (14/40) of gastric cancer (GC), 30.0% (3/10) of dysplasia (Dys), 28.6% (4/14) of intestinal metaplasia (IM) and 10.5% (2/19) of chronic atrophic gastritis (CAG), respectively. A frequent mutation was a T to C transition at nucleotide 889, which was observed in 4 GC cell lines, 7 GC, 2 Dys, and 2 IM. The level of pol beta expression in tumors was higher than that of their matched normal tissues and gradual changes from GC, Dys, CAG to IM. These results indicate that the mutation and overexpression of pol beta may influence the progression during gastric carcinogenesis.