Local deformations revealed by dynamics simulations of DNA polymerase Beta with DNA mismatches at the primer terminus

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.

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.

Interactions of the DNA Repair Enzyme Human Thymine DNA Glycosylase with Cognate and Noncognate DNA

Glycosylases specifically recognize and flip their target base out of the DNA helix into the enzyme's active site. Our simulations show that a partially flipped state, already present in free DNA carrying a T:G mispair, becomes the more probable state compared to the closed state after binding of thymine DNA glycosylase (TDG). Paired thymine (T:A) or methyl-cytosine (mC:G) does not exhibit a partially flipped state in free or complexed DNA. Important enzyme-DNA interactions exhibit significant strength in the intrahelical and extrahelical TDG-DNA complexes. The computed binding free energy differences suggest these interactions account for the stabilization of the partially flipped state, thereby driving the T:G mispair toward base flip. In the fully flipped state, the cognate base thymine is significantly better accommodated in the enzyme's active site than noncognate bases are, suggesting the hydrolysis step as the last of several stages at which base recognition can be achieved.

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.

How does the long G·G* Watson-Crick DNA base mispair comprising keto and enol tautomers of the guanine tautomerise? The results of a QM/QTAIM investigation

The double proton transfer (DPT) in the long G·G* Watson-Crick base mispair (|C6N1(G*)N1C6(G)| = 36.4°; C1 symmetry), involving keto and enol tautomers of the guanine (G) nucleobase, along two intermolecular neighboring O6H···O6 (8.39) and N1···HN1 (6.14 kcal mol(-1)) H-bonds that were established to be slightly anti-cooperative, leads to its transformation into the G*·G base mispair through a single transition state (|C6N1N1C6| = 37.1°; C1), namely to the interconversion into itself. It was shown that the G·G* ↔ G*·G tautomerisation via the DPT is assisted by the third specific contact, that sequentially switches along the intrinsic reaction coordinate (IRC) in an original way: (G)N2H···N2(G*) H-bond (-25.13 to -10.37) → N2···N2 van der Waals contact (-10.37 to -9.23) → (G)N2···HN2(G*) H-bond (-9.23 to 0.79) → (G*)N2···HN2(G) H-bond (0.79 to 7.35 Bohr). The DPT tautomerisation was found to proceed through the asynchronous concerted mechanism by employing the QM/QTAIM approach and the methodology of the scans of the geometric, electron-topological, energetic, polar and NBO properties along the IRC. Nine key points, that can be considered as part of the tautomerisation repertoire, have been established and analyzed in detail. Furthermore, it was shown that the G·G* or G*·G base mispair is a thermodynamically and dynamically stable structure with a lifetime of 8.22 × 10(-10) s and all 6 low-frequency intermolecular vibrations are able to develop during this time span. Lastly, our results highlight the importance of the G·G* ↔ G*·G DPT tautomerisation, which can have implications for biological and chemical sensing applications.

Effects of stability of base pairs containing an oxazolone on DNA elongation

The nucleoside 2,2,4-triamino-5(2H)-oxazolone (Oz) can result from oxidative damage to guanine residues in DNA. Despite differences among the three polymerases (Pol β, KF exo⁻, and Pol η) regarding nucleotide incorporation patterns opposite Oz, all three polymerases can incorporate guanine opposite Oz. Based on ab initio calculations, we proposed a structure for a stable Oz:G base pair. Here, to assess the stability of each Oz-containing base pair (Oz:G, Oz:A, Oz:C, and Oz:T) upon DNA replication, we determined the efficiency of Pol β-, KF exo⁻-, or Pol η-catalyzed primer extension beyond each base pair. With each polymerase, extension beyond Oz:G was more efficient than that beyond Oz:A, Oz:C, or Oz:T. Moreover, thermal denaturation studies revealed that the T_m value for the duplex containing Oz:G was significantly higher than those obtained for duplexes containing Oz:A, Oz:C, or Oz:T. Therefore, the results from ab initio calculations along with those from DNA replication assays and thermal denaturation experiments supported the conclusion that Oz:G is the most stable of the Oz-containing base pairs.

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.

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.

Partial base flipping is sufficient for strand slippage near DNA duplex termini

Strand slippage is a structural mechanism by which insertion-deletion (indel) mutations are introduced during replication by polymerases. Three-dimensional atomic-resolution structural pathways are still not known for the decades-old template slippage description. The dynamic nature of the process and the higher energy intermediates involved increase the difficulty of studying these processes experimentally. In the present study, restrained and unrestrained molecular dynamics simulations, carried out using multiple nucleic acid force fields, are used to demonstrate that partial base-flipping can be sufficient for strand slippage at DNA duplex termini. Such strand slippage can occur in either strand, i.e. near either the 3' or the 5' terminus of a DNA strand, which suggests that similar structural flipping mechanisms can cause both primer and template slippage. In the repetitive mutation hot-spot sequence studied, non-canonical base-pairing with exposed DNA groove atoms of a neighboring G:C base-pair stabilizes a partially flipped state of the cytosine. For its base-pair partner guanine, a similar partially flipped metastable intermediate was not detected, and the propensity for sustained slippage was also found to be lower. This illustrates that a relatively small metastable DNA structural distortion in polymerase active sites could allow single base insertion or deletion mutations to occur, and stringent DNA groove molecular recognition may be required to maintain intrinsic DNA polymerase fidelity. The implications of a close relationship between base-pair dissociation, base unstacking, and strand slippage are discussed in the context of sequence dependence of indel mutations.

The effect of a G:T mispair on the dynamics of DNA

Distortions in the DNA sequence such as damages or mispairs are specifically recognized and processed by DNA repair enzymes. A particular challenge for the enzymatic specificity is the recognition of a wrongly-placed native nucleotide such as thymine in T:G mispairs. An important step of substrate binding which is observed in many repair proteins is the flipping of the target base out of the DNA helix into the enzyme's active site. In this work we investigate how much the intrinsic dynamics of mispaired DNA is changed compared to canonical DNA. Our molecular dynamics simulations of DNA with and without T:G mispairs show significant differences in the conformation of paired and mispaired DNA. The wobble pair T:G shows local distortions such as twist, shear and stretch which deviate from canonical B form values. Moreover, the T:G mispair is found to be kinetically less stable, exhibiting two states with respect to base opening: a closed state comparable to the canonical base pairs, and a more open state, indicating a proneness for base flip. In addition, we observe that the thymine base in a T:G mispair is significantly more probable to be flipped than thymine in a T:A pair or cytosine in a C:G pair. Such local deformations and in particular the existence of a second, more-open state can be speculated to help the target-site recognition by repair enzymes.

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.

Base-pairing preferences, physicochemical properties and mutational behaviour of the DNA lesion 8-nitroguanine

8-Nitro-2'-deoxyguanosine (8-nitrodG) is a relatively unstable, mutagenic lesion of DNA that is increasingly believed to be associated with tissue inflammation. Due to the lability of the glycosidic bond, 8-nitrodG cannot be incorporated into oligodeoxynucleotides (ODNs) by chemical DNA synthesis and thus very little is known about its physicochemical properties and base-pairing preferences. Here we describe the synthesis of 8-nitro-2'-O-methylguanosine, a ribonucleoside analogue of this lesion, which is sufficiently stable to be incorporated into ODNs. Physicochemical studies demonstrated that 8-nitro-2'-O-methylguanosine adopts a syn conformation about the glycosidic bond; thermal melting studies and molecular modelling suggest a relatively stable syn-8-nitroG·anti-G base pair. Interestingly, when this lesion analogue was placed in a primer-template system, extension of the primer by either avian myeloblastosis virus reverse transcriptase (AMV-RT) or human DNA polymerase β (pol β), was significantly impaired, but where incorporation opposite 8-nitroguanine did occur, pol β showed a 2:1 preference to insert dA over dC, while AMV-RT incorporated predominantly dC. The fact that no 8-nitroG·G base pairing is seen in the primer extension products suggests that the polymerases may discriminate against this pairing system on the basis of its poor geometric match to a Watson-Crick pair.

DNA polymerase family X: function, structure, and cellular roles

The X family of DNA polymerases in eukaryotic cells consists of terminal transferase and DNA polymerases beta, lambda, and mu. These enzymes have similar structural portraits, yet different biochemical properties, especially in their interactions with DNA. None of these enzymes possesses a proofreading subdomain, and their intrinsic fidelity of DNA synthesis is much lower than that of a polymerase that functions in cellular DNA replication. In this review, we discuss the similarities and differences of three members of Family X: polymerases beta, lambda, and mu. We focus on biochemical mechanisms, structural variation, fidelity and lesion bypass mechanisms, and cellular roles. Remarkably, although these enzymes have similar three-dimensional structures, their biochemical properties and cellular functions differ in important ways that impact cellular function.

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.

Mismatched base-pair simulations for ASFV Pol X/DNA complexes help interpret frequent G*G misincorporation

DNA polymerase X (pol X) from the African swine fever virus is a 174-amino-acid repair polymerase that likely participates in a viral base excision repair mechanism, characterized by low fidelity. Surprisingly, pol X's insertion rate of the G*G mispair is comparable to that of the four Watson-Crick base pairs. This behavior is in contrast with another X-family polymerase, DNA polymerase beta (pol beta), which inserts G*G mismatches poorly, and has higher DNA repair fidelity. Using molecular dynamics simulations, we previously provided support for an induced-fit mechanism for pol X in the presence of the correct incoming nucleotide. Here, we perform molecular dynamics simulations of pol X/DNA complexes with different incoming incorrect nucleotides in various orientations [C*C, A*G, and G*G (anti) and A*G and G*G (syn)] and compare the results to available kinetic data and prior modeling. Intriguingly, the simulations reveal that the G*G mispair with the incoming nucleotide in the syn configuration undergoes large-scale conformational changes similar to that observed in the presence of correct base pair (G*C). The base pairing in the G*G mispair is achieved via Hoogsteen hydrogen bonding with an overall geometry that is well poised for catalysis. Simulations for other mismatched base pairs show that an intermediate closed state is achieved for the A*G and G*G mispair with the incoming dGTP in anti conformation, while the protein remains near the open conformation for the C*C and the A*G syn mismatches. In addition, catalytic site geometry and base pairing at the nascent template-incoming nucleotide interaction reveal distortions and misalignments that range from moderate for A*G anti to worst for the C*C complex. These results agree well with kinetic data for pol X and provide a structural/dynamic basis to explain, at atomic level, the fidelity of this polymerase compared with other members of the X family. In particular, the more open and pliant active site of pol X, compared to pol beta, allows pol X to accommodate bulkier mismatches such as guanine opposite guanine, while the more structured and organized pol beta active site imposes higher discrimination, which results in higher fidelity. The possibility of syn conformers resonates with other low-fidelity enzymes such as Dpo4 (from the Y family), which readily accommodate oxidative lesions.

The Asp285 variant of DNA polymerase beta extends mispaired primer termini via increased nucleotide binding

Endogenous DNA damage occurs at a rate of at least 20,000 lesions per cell per day. Base excision repair (BER) is a key pathway for maintaining genome stability. Several pol beta variants were identified as conferring resistance to 3'-azido-3'-deoxythymidine (AZT) in Escherichia coli (Kosa et al. (1999) J. Biol. Chem. 274, 3851-3858). Detailed biochemical studies on one of these AZT-resistant variants, His285 to Asp, have shown that the H285D variant of pol beta possesses pre-steady-state kinetics that are similar to the wild-type polymerase. In gap filling assays with 5-bp gapped DNA, H285D showed a slight mutator phenotype. In depth single turnover kinetic analysis revealed that H285D is much more efficient than wild-type pol beta at extending mispaired primer termini. This mispair extension property of H285D is attributed to a greatly increased binding to the next correct nucleotide in the presence of a mispair. This change in K d(dNTP),app is not accompanied by a change in k pol; values for k pol are the same for both H285D and wild-type. Close examination of available structural data, as well as molecular modeling, has shown that residue 285 is able to make several stabilizing contacts in the fingers domain of the polymerase, and the introduction of a negatively charged side chain could have important effects on the enzyme. It is postulated that the loss of the contact between His285, Lys289, and Ile323 is responsible for the ability of H285D to extend mispairs through disruption of contacts near the C-terminal end of pol beta and propagation into the nucleotide binding pocket.

A lysine residue in the fingers subdomain of T7 DNA polymerase modulates the miscoding potential of 8-oxo-7,8-dihydroguanosine

8-oxo-7,8-dihydroguanosine (8oG) is a highly mutagenic DNA lesion that stably pairs with adenosine, forming 8oG(syn).dA(anti) Hoogsteen base pairs. DNA polymerases show different propensities to insert dCMP or dAMP opposite 8oG, but the molecular mechanisms that determine faithful or mutagenic bypass are poorly understood. Here, we report kinetic and structural data providing evidence that, in T7 DNA polymerase, residue Lys536 is responsible for attenuating the miscoding potential of 8oG. The Lys536Ala polymerase shows a significant increase in mutagenic 8oG bypass versus wild-type polymerase, and a crystal structure of the Lys536Ala mutant reveals a closed complex with an 8oG(syn).dATP mismatch in the polymerase active site, in contrast to the unproductive, open complex previously obtained by using wild-type polymerase. We propose that Lys536 acts as a steric and/or electrostatic filter that attenuates the miscoding potential of 8oG by normally interfering with the binding of 8oG in a syn conformation that pairs with dATP.

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.

Exploration of the conformational space of myosin recovery stroke via molecular dynamics

Muscle contractions are driven by cyclic conformational changes of myosin, whose molecular mechanisms of operation are being elucidated by recent advances in crystallographic studies and single molecule experiments. To complement such structural studies and consider the energetics of the conformational changes of myosin head, umbrella sampling molecular dynamics (MD) simulations were performed with the all-atom model of the scallop myosin sub-fragment 1 (S1) with a bound ATP in solution in explicit water using the crystallographic near-rigor and transition state conformations as two references. The constraints on RMSD reaction coordinates used for the umbrella sampling were found to steer the conformational changes efficiently, and relatively close correlations have been observed between the set of characteristic structural changes including the lever arm rotation and the closing of the nucleotide binding pocket. The lever arm angle and key residue interaction distances in the nucleotide binding pocket and the relay helix show gradual changes along the recovery stroke reaction coordinate, consistent with previous crystallographic and computational minimum energy studies. Thermal fluctuations, however, appear to make the switch-2 coordination of ATP more flexible than suggested by crystal structures. The local solvation environment of the fluorescence probe, Trp 507 (scallop numbering), also appears highly mobile in the presence of thermal fluctuations.

The Werner syndrome protein operates in base excision repair and cooperates with DNA polymerase beta

Genome instability is a characteristic of cancer and aging, and is a hallmark of the premature aging disorder Werner syndrome (WS). Evidence suggests that the Werner syndrome protein (WRN) contributes to the maintenance of genome integrity through its involvement in DNA repair. In particular, biochemical evidence indicates a role for WRN in base excision repair (BER). We have previously reported that WRN helicase activity stimulates DNA polymerase beta (pol β) strand displacement synthesis in vitro. In this report we demonstrate that WRN exonuclease activity can act cooperatively with pol β, a polymerase lacking 3′–5′ proofreading activity. Furthermore, using small interference RNA technology, we demonstrate that WRN knockdown cells are hypersensitive to the alkylating agent methyl methanesulfonate, which creates DNA damage that is primarily repaired by the BER pathway. In addition, repair assays using whole cell extracts from WRN knockdown cells indicate a defect in long patch (LP) BER. These findings demonstrate that WRN plays a direct role in the repair of methylation-induced DNA damage, and suggest a role for both WRN helicase and exonuclease activities together with pol β during LP BER.

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.

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.

Molecular dynamics studies on free and bound targets of the bovine papillomavirus type I e2 protein: the protein binding effect on DNA and the recognition mechanism

Molecular dynamics simulations of a total duration of 30 ns in explicit solvent were carried out on the BPV-1-E2 protein complexed to a high-affinity DNA target containing the two hydrogen-bonded ACCG.CGGT half-sites separated by the noncontacted ACGT sequence. The analysis of the trajectories focuses on the DNA structure and on the dynamics. The data are compared to those issued from recent simulations made on three free targets that recognize E2 with different affinities. E2 does not drastically perturb the mechanic properties of the free DNA: the structural relationships between the BI/BII backbone substates and some helical parameters are preserved in the complex despite a severe slowing down of the phosphate group motions. The structures of both free and bound half-sites are very close to each other although the conformational space explored by these regions is narrowed when they are contacted by the protein. The enhanced plasticity found in the best free target spacers, mainly manifested through the backbone motions, allows a clear overlap between several free and bound global DNA features such as the base displacement. Furthermore, this flexibility is preserved in the complex. Our results support the hypothesis that E2 takes advantage of free predistorted structures that may minimize the DNA deformation cost. In addition, we observe that E2 is far from totally stiffening the DNA, suggesting that the entropic penalty inherent in the complex formation could be limited.

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.

Highly organized but pliant active site of DNA polymerase beta: compensatory mechanisms in mutant enzymes revealed by dynamics simulations and energy analyses

To link conformational transitions noted for DNA polymerases with kinetic results describing catalytic efficiency and fidelity, we investigate the role of key DNA polymerase beta residues on subdomain motion through simulations of five single-residue mutants: Arg-283-Ala, Tyr-271-Ala, Asp-276-Val, Arg-258-Lys, and Arg-258-Ala. Since a movement toward a closed state was only observed for R258A, we suggest that Arg(258) is crucial in modulating motion preceding chemistry. Analyses of protein/DNA interactions in the mutant active site indicate distinctive hydrogen bonding and van der Waals patterns arising from compensatory structural adjustments. By comparing closed mutant complexes with the wild-type enzyme, we interpret experimentally derived nucleotide binding affinities in molecular terms: R283A (decreased), Y271A (increased), D276V (increased), and R258A (decreased). Thus, compensatory interactions (e.g., in Y271A with adjacent residues Phe(272), Asn(279), and Arg(283)) increase the overall binding affinity for the incoming nucleotide although direct interactions may decrease. Together with energetic analyses, we predict that R258G might increase the rate of nucleotide insertion and maintain enzyme fidelity as R258A; D276L might increase the nucleotide binding affinity more than D276V; and R283A/K280A might decrease the nucleotide binding affinity and increase misinsertion more than R283A. The combined observations regarding key roles of specific residues (e.g., Arg(258)) and compensatory interactions echo the dual nature of polymerase active site, namely versatility (to accommodate various basepairs) and specificity (for preserving fidelity) and underscore an organized but pliant active site essential to enzyme function.

In silico evidence for DNA polymerase-beta's substrate-induced conformational change

Structural information for mammalian DNA pol-beta combined with molecular and essential dynamics studies have provided atomistically detailed views of functionally important conformational rearrangements that occur during DNA repair and replication. This conformational closing before the chemical reaction is explored in this work as a function of the bound substrate. Anchors for our study are available in crystallographic structures of the DNA pol-beta in "open" (polymerase bound to gapped DNA) and "closed" (polymerase bound to gapped DNA and substrate, dCTP) forms; these different states have long been used to deduce that a large-scale conformational change may help the polymerase choose the correct nucleotide, and hence monitor DNA synthesis fidelity, through an "induced-fit" mechanism. However, the existence of open states with bound substrate and closed states without substrates suggest that substrate-induced conformational closing may be more subtle. Our dynamics simulations of two pol-beta/DNA systems (with/without substrates at the active site) reveal the large-scale closing motions of the thumb and 8-kDa subdomains in the presence of the correct substrate--leading to nearly perfect rearrangement of residues in the active site for the subsequent chemical step of nucleotidyl transfer--in contrast to an opening trend when the substrate is absent, leading to complete disassembly of the active site residues. These studies thus provide in silico evidence for the substrate-induced conformational rearrangements, as widely assumed based on a variety of crystallographic open and closed complexes. Further details gleaned from essential dynamics analyses clarify functionally relevant global motions of the polymerase-beta/DNA complex as required to prepare the system for the chemical reaction of nucleotide extension.

The role of water in protein-DNA recognition

Is it by design or by default that water molecules are observed at the interfaces of some protein-DNA complexes? Both experimental and theoretical studies on the thermodynamics of protein-DNA binding overwhelmingly support the extended hydrophobic view that water release from interfaces favors binding. Structural and energy analyses indicate that the waters that remain at the interfaces of protein-DNA complexes ensure liquid-state packing densities, screen the electrostatic repulsions between like charges (which seems to be by design), and in a few cases act as linkers between complementary charges on the biomolecules (which may well be by default). This review presents a survey of the current literature on water in protein-DNA complexes and a critique of various interpretations of the data in the context of the role of water in protein-DNA binding and principles of protein-DNA recognition in general.

Gene structure, purification and characterization of DNA polymerase beta from Xiphophorus maculatus

Cloning of the Xiphophorus maculatus Polbeta gene and overexpression of the recombinant Polbeta protein has been performed. The organization of the XiphPolbeta introns and exons, including intron-exon boundaries, have been assigned and were found to be similar to that for human Polbeta with identical exon sizes except for exon XII coding for an additional two amino acid residues in Xiphophorus. The cDNA sequence encoding the 337-amino acid X. maculatus DNA polymerase beta (Polbeta) protein was subcloned into the Escherichia coli expression plasmid pET. Induction of transformed E. coli cells resulted in the high-level expression of soluble recombinant Polbeta, which catalyzed DNA synthesis on template-primer substrates. The steady-state Michaelis constants (Km) and catalytic efficiencies (kcat/Km) of the recombinant XiphPolbeta for nucleotide insertion opposite single-nucleotide gap DNA substrates were measured and compared with previously published values for recombinant human Polbeta. Steady-state in vitro Km and kcat/Km values for correct nucleotide insertion by XiphPolbeta and human Polbeta were similar, although the recombinant Xiphophorus protein exhibited 2.5-7-fold higher catalytic efficiencies for dGTP and dCTP insertion versus human Polbeta. In contrast, the recombinant XiphPolbeta displayed significantly lower fidelities than human Polbeta for dNTP insertion opposite a single-nucleotide gap at 37 degrees C.

Influence of DNA structure on DNA polymerase beta active site function: extension of mutagenic DNA intermediates

In the ternary substrate complex of DNA polymerase (pol) beta, the nascent base pair (templating and incoming nucleotides) is sandwiched between the duplex DNA terminus and polymerase. To probe molecular interactions in the dNTP-binding pocket, we analyzed the kinetic behavior of wild-type pol beta on modified DNA substrates that alter the structure of the DNA terminus and represent mutagenic intermediates. The DNA substrates were modified to 1) alter the sequence of the duplex terminus (matched and mismatched), 2) introduce abasic sites near the nascent base pair, and 3) insert extra bases in the primer or template strands to mimic frameshift intermediates. The results indicate that the nucleotide insertion efficiency (k(cat)/K(m), dGTP-dC) is highly dependent on the sequence identity of the matched (i.e. Watson-Crick base pair) DNA terminus (template/primer, G/C approximately A/T > T/A approximately C/G). Mismatches at the primer terminus strongly diminish correct nucleotide insertion efficiency but do not affect DNA binding affinity. Transition intermediates are generally extended more easily than transversions. Most mismatched primer termini decrease the rate of insertion and binding affinity of the incoming nucleotide. In contrast, the loss of catalytic efficiency with homopurine mismatches at the duplex DNA terminus is entirely due to the inability to insert the incoming nucleotide, since K(d)((dGTP)) is not affected. Abasic sites and extra nucleotides in and around the duplex terminus decrease catalytic efficiency and are more detrimental to the nascent base pair binding pocket when situated in the primer strand than the equivalent position in the template strand.

Orchestration of cooperative events in DNA synthesis and repair mechanism unraveled by transition path sampling of DNA polymerase beta's closing

Our application of transition path sampling to a complex biomolecular system in explicit solvent, the closing transition of DNA polymerase beta, unravels atomic and energetic details of the conformational change that precedes the chemical reaction of nucleotide incorporation. The computed reaction profile offers detailed mechanistic insights into, as well as kinetic information on, the complex process essential for DNA synthesis and repair. The five identified transition states extend available experimental and modeling data by revealing highly cooperative dynamics and critical roles of key residues (Arg-258, Phe-272, Asp-192, and Tyr-271) in the enzyme's function. The collective cascade of these sequential conformational changes brings the DNA/DNA polymerase beta system to a state nearly competent for the chemical reaction and suggests how subtle residue motions and conformational rate-limiting steps affect reaction efficiency and fidelity; this complex system of checks and balances directs the system to the chemical reaction and likely helps the enzyme discriminate the correct from the incorrect incoming nucleotide. Together with the chemical reaction, these conformational features may be central to the dual nature of polymerases, requiring specificity (for correct nucleotide selection) as well as versatility (to accommodate different templates at every step) to maintain overall fidelity. Besides leading to these biological findings, our developed protocols open the door to other applications of transition path sampling to long-time, large-scale biomolecular reactions.

Molecular modeling of the major benzo[a]pyrene N2-dG adduct in cases where mutagenesis results are known in double stranded DNA

The potent mutagen/carcinogen benzo[a]pyrene (B[a]P) is metabolically activated to (+)-anti-B[a]PDE, which induces a full spectrum of mutations (e.g. GC-->TA, GC-->AT, etc.). One hypothesis for this complexity is that different mutations are induced by different conformations of its major adduct [+ta]-B[a]P-N2-dG when bypassed during DNA replication (probably by different DNA polymerases). Previous molecular modeling studies suggested that B[a]P-N2-dG adducts can in principle adopt at least 16 potential conformational classes in ds-DNA. Herein we report on molecular modeling studies with the eight conformations most likely to be relevant to base substitution mutagenesis in 10 cases where mutagenesis has been studied in ds-DNA plasmids in E. coli with B[a]P-N2-dG adducts of differing stereoisomers and DNA sequence contexts, as well as in five cases where the conformation is known by NMR. Of the approximately 11,000 structures generated in this study, the computed lowest energy structures are reported for 120 cases (i.e. eight conformations and 15 examples), and their conformations compared. Of the eight conformations, four are virtually always computed to be high in energy. The remaining four lower energy conformations include two with the BP moiety in the minor groove (designated: BPmi5 and BPmi3), and two base-displaced conformations, one with the dG moiety in the major groove (designated: Gma5) and one with the dG in the minor groove (designated: Gmi3). Interestingly, these four are the only conformations that have been observed for B[a]P-N2-dG adducts in NMR studies. Independent of sequence contexts and adduct stereochemistry, BPmi5 structures tend to look reasonably similar, as do BPmi3 structures, while the base-displaced structures Gma5 and BPmi3 tend to show greater variability in structure. A correlation was sought between modeling and mutagenesis results in the case of the low energy conformations BPmi5, BPmi3, Gma5 and Gma3. Plots of log[(G-->T)/(G-->A)] versus energy[(conformation X)-(conformation Y)] were constructed for all six pairwise combinations of these four conformations, and the only plot giving a straight line involved Gma5 and Gmi3. While this finding is striking, its significance is unclear (as discussed).