Incorrect nucleotide insertion at the active site of a G:A mismatch catalyzed by DNA polymerase beta

Polγ coordinates DNA synthesis and proofreading to ensure mitochondrial genome integrity

Accurate replication of mitochondrial DNA (mtDNA) by DNA polymerase γ (Polγ) is essential for maintaining cellular energy supplies, metabolism, and cell cycle control. To illustrate the structural mechanism for Polγ coordinating polymerase (pol) and exonuclease (exo) activities to ensure rapid and accurate DNA synthesis, we determined four cryo-EM structures of Polγ captured after accurate or erroneous incorporation to a resolution of 2.4-3.0 Å. The structures show that Polγ employs a dual-checkpoint mechanism to sense nucleotide misincorporation and initiate proofreading. The transition from replication to error editing is accompanied by increased dynamics in both DNA and enzyme, in which the polymerase relaxes its processivity and the primer-template DNA unwinds, rotates, and backtracks to shuttle the mismatch-containing primer terminus 32 Å to the exo site for editing. Our structural and functional studies also provide a foundation for analyses of Polγ mutation-induced human diseases and aging.

Probing the general base for DNA polymerization in telomerase: a molecular dynamics investigation

Telomerase is an RNA-dependent DNA polymerase that plays a role in the maintenance of the 3' end of the eukaryotic chromosome, known as a telomere, by catalyzing the DNA polymerization reaction in cancer and embryonic stem cells. The detailed molecular details of the DNA polymerization by telomerase, especially the general base for deprotonating the terminal 3'-hydroxyl, which triggers the chemical reaction, remain elusive. We conducted a computational investigation using hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations to probe the detailed mechanism of the reaction. Our simulations started with the telomerase:RNA:DNA:dNTP ternary complex, and by using enhanced sampling QM/MM MD simulations, we probed the general base involved directly in the polymerization. We report the participation of an aspartate (Asp344) coordinated to Mg and an active site water molecule, jointly acting as a base during nucleic acid addition. The Asp344 residue remains transiently protonated during the course of the reaction, and later it deprotonates by transferring its proton to the water at the end of the reaction.

Structural insights into the promutagenic bypass of the major cisplatin-induced DNA lesion

The cisplatin-1,2-d(GpG) (Pt-GG) intrastrand cross-link is the predominant DNA lesion generated by cisplatin. Cisplatin has been shown to predominantly induce G to T mutations and Pt-GG permits significant misincorporation of dATP by human DNA polymerase β (polβ). In agreement, polβ overexpression, which is frequently observed in cancer cells, is linked to cisplatin resistance and a mutator phenotype. However, the structural basis for the misincorporation of dATP opposite Pt-GG is unknown. Here, we report the first structures of a DNA polymerase inaccurately bypassing Pt-GG. We solved two structures of polβ misincorporating dATP opposite the 5'-dG of Pt-GG in the presence of Mg2+ or Mn2+. The Mg2+-bound structure exhibits a sub-optimal conformation for catalysis, while the Mn2+-bound structure is in a catalytically more favorable semi-closed conformation. In both structures, dATP does not form a coplanar base pairing with Pt-GG. In the polβ active site, the syn-dATP opposite Pt-GG appears to be stabilized by protein templating and pi stacking interactions, which resembles the polβ-mediated dATP incorporation opposite an abasic site. Overall, our results suggest that the templating Pt-GG in the polβ active site behaves like an abasic site, promoting the insertion of dATP in a non-instructional manner.

Bypass of the Major Alkylative DNA Lesion by Human DNA Polymerase η

A wide range of endogenous and exogenous alkylating agents attack DNA to generate various alkylation adducts. N7-methyl-2-deoxyguanosine (Fm7dG) is the most abundant alkylative DNA lesion. If not repaired, Fm7dG can undergo spontaneous depurination, imidazole ring-opening, or bypass by translesion synthesis DNA polymerases. Human DNA polymerase η (polη) efficiently catalyzes across Fm7dG in vitro, but its structural basis is unknown. Herein, we report a crystal structure of polη in complex with templating Fm7dG and an incoming nonhydrolyzable dCTP analog, where a 2'-fluorine-mediated transition destabilization approach was used to prevent the spontaneous depurination of Fm7dG. The structure showed that polη readily accommodated the Fm7dG:dCTP base pair with little conformational change of protein and DNA. In the catalytic site, Fm7dG and dCTP formed three hydrogen bonds with a Watson-Crick geometry, indicating that the major keto tautomer of Fm7dG is involved in base pairing. The polη-Fm7dG:dCTP structure was essentially identical to the corresponding undamaged structure, which explained the efficient bypass of the major methylated lesion. Overall, the first structure of translesion synthesis DNA polymerase bypassing Fm7dG suggests that in the catalytic site of Y-family DNA polymerases, small N7-alkylguanine adducts may be well tolerated and form the canonical Watson-Crick base pair with dCTP through their keto tautomers.

Uncovering a unique approach for damaged DNA replication: A computational investigation of a mutagenic tobacco-derived thymine lesion

4-(Methylnitrosamino)-1-(3-pyridyl)-1-butanone is a potent nicotine carcinogen that leads to many DNA lesions, the most persistent being the O2-[4-oxo-4-(3-pyridyl)butyl]thymine adduct (POB-T). Although the experimental mutagenic profile for the minor groove POB-T lesion has been previously reported, the findings are puzzling in terms of the human polymerases involved. Specifically, while pol κ typically replicates minor groove adducts, in vivo studies indicate pol η replicates POB-T despite being known for processing major groove adducts. Our multiscale modeling approach reveals that the canonical (anti) glycosidic orientation of POB-T can fit in the pol κ active site, but only a unique (syn) POB-T conformation is accommodated by pol η. These distinct binding orientations rationalize the differential in vitro mutagenic spectra based on the preferential stabilization of dGTP and dTTP opposite the lesion for pol κ and η, respectively. Overall, by uncovering the first evidence for the replication of a damaged pyrimidine in the syn glycosidic orientation, the current work provides the insight necessary to clarify a discrepancy in the DNA replication literature, expand the biological role of the critical human pol η, and understand the mutational signature in human cancers associated with tobacco exposure.

Insights into the effect of minor groove interactions and metal cofactors on mutagenic replication by human DNA polymerase β

DNA polymerases accommodate various base-pair conformations in the event of incorrect insertions. In particular, Watson-Crick-like dG:dTTP base pair has been observed at the insertion site of human DNA polymerase β (pol β). A potential factor contributing to the diverse conformations of base-pair mismatches is minor groove interactions. To gain insights into the effect of minor groove interactions on base-pair conformations, we generated an Asn279Ala polβ mutant that cannot make minor groove contacts with an incoming nucleotide. We conducted structural and kinetic studies of Asn279Ala polβ in complex with incoming dTTP and templating dG or O6-methyl-dG. The crystal structure of the Asn279Ala polβ-G:T complex showed a wobble dG:dTTP base pair, indicating that the previously observed Watson-Crick-like dG:dTTP conformation was induced by the minor groove contact. In contrast, O6-methyl-dG, an analog of the enol tautomer of guanine, formed a Watson-Crick-like base pair with dTTP in the absence of the minor groove contact. These results suggest that the Watson-Crick-like G:T base pair at the insertion site is formed by the rare enol tautomers of G or T, whose population is increased by the minor groove hydrogen bond with Asn279. Kinetic studies showed that Asn279Ala mutation decreased dG:dTTP misincorporation rate six-fold in the presence of Mg²⁺ but increased the rate three-fold in the presence of Mn²⁺, highlighting the effect of minor groove interactions and metal ions on promutagenic replication by polβ.

Revealing the role of the product metal in DNA polymerase β catalysis

DNA polymerases catalyze a metal-dependent nucleotidyl transferase reaction during extension of a DNA strand using the complementary strand as a template. The reaction has long been considered to require two magnesium ions. Recently, a third active site magnesium ion was identified in some DNA polymerase product crystallographic structures, but its role is not known. Using quantum mechanical/ molecular mechanical calculations of polymerase β, we find that a third magnesium ion positioned near the newly identified product metal site does not alter the activation barrier for the chemical reaction indicating that it does not have a role in the forward reaction. This is consistent with time-lapse crystallographic structures following insertion of Sp-dCTPαS. Although sulfur substitution deters product metal binding, this has only a minimal effect on the rate of the forward reaction. Surprisingly, monovalent sodium or ammonium ions, positioned in the product metal site, lowered the activation barrier. These calculations highlight the impact that an active site water network can have on the energetics of the forward reaction and how metals or enzyme side chains may interact with the network to modulate the reaction barrier. These results also are discussed in the context of earlier findings indicating that magnesium at the product metal position blocks the reverse pyrophosphorolysis reaction.

Metal Fluorides: Tools for Structural and Computational Analysis of Phosphoryl Transfer Enzymes

The phosphoryl group, PO₃^-, is the dynamic structural unit in the biological chemistry of phosphorus. Its transfer from a donor to an acceptor atom, with oxygen much more prevalent than nitrogen, carbon, or sulfur, is at the core of a great majority of enzyme-catalyzed reactions involving phosphate esters, anhydrides, amidates, and phosphorothioates. The serendipitous discovery that the phosphoryl group could be labeled by "nuclear mutation," by substitution of PO₃^- by MgF₃^- or AlF₄^-, has underpinned the application of metal fluoride (MF _x ) complexes to mimic transition states for enzymatic phosphoryl transfer reactions, with sufficient stability for experimental analysis. Protein crystallography in the solid state and ¹⁹F NMR in solution have enabled direct observation of ternary and quaternary protein complexes embracing MF _x transition state models with precision. These studies have underpinned a radically new mechanistic approach to enzyme catalysis for a huge range of phosphoryl transfer processes, as varied as kinases, phosphatases, phosphomutases, and phosphohydrolases. The results, without exception, have endorsed trigonal bipyramidal geometry (tbp) for concerted, "in-line" stereochemistry of phosphoryl transfer. QM computations have established the validity of tbp MF _x complexes as reliable models for true transition states, delivering similar bond lengths, coordination to essential metal ions, and virtually identical hydrogen bond networks. The emergence of protein control of reactant orbital overlap between bond-forming species within enzyme transition states is a new challenging theme for wider exploration.

Uniform Free-Energy Profiles of the P-O Bond Formation and Cleavage Reactions Catalyzed by DNA Polymerases β and λ

Human X-family DNA polymerases β (Polβ) and λ (Polλ) catalyze the nucleotidyl-transfer reaction in the base excision repair pathway of the cellular DNA damage response. Using empirical valence bond and free-energy perturbation simulations, we explore the feasibility of various mechanisms for the deprotonation of the 3′-OH group of the primer DNA strand, and the subsequent formation and cleavage of P–O bonds in four Polβ, two truncated Polλ (tPolλ), and two tPolλ Loop1 mutant (tPolλΔL1) systems differing in the initial X-ray crystal structure and nascent base pair. The average calculated activation free energies of 14, 18, and 22 kcal mol^–1 for Polβ, tPolλ, and tPolλΔL1, respectively, reproduce the trend in the observed catalytic rate constants. The most feasible reaction pathway consists of two successive steps: specific base (SB) proton transfer followed by rate-limiting concerted formation and cleavage of the P–O bonds. We identify linear free-energy relationships (LFERs) which show that the differences in the overall activation and reaction free energies among the eight studied systems are determined by the reaction free energy of the SB proton transfer. We discuss the implications of the LFERs and suggest pK_a of the 3′-OH group as a predictor of the catalytic rate of X-family DNA polymerases.

Structures of DNA Polymerase Mispaired DNA Termini Transitioning to Pre-catalytic Complexes Support an Induced-Fit Fidelity Mechanism

High-fidelity DNA synthesis requires that polymerases display a strong preference for right nucleotide insertion. When the wrong nucleotide is inserted, the polymerase deters extension from the mismatched DNA terminus. Twenty-three crystallographic structures of DNA polymerase β with terminal template-primer mismatches were determined as binary DNA and ternary pre-catalytic substrate complexes. These structures indicate that the mismatched termini adopt various distorted conformations that attempt to satisfy stacking and hydrogen-bonding interactions. The binary complex structures indicate an induced strain in the mismatched template nucleotide. Addition of a non-hydrolyzable incoming nucleotide stabilizes the templating nucleotide with concomitant strain in the primer terminus. Several dead-end ternary complex structures suggest that DNA synthesis might occur as the enzyme transitions from an open to a closed complex. The structures are consistent with an induced-fit mechanism where a mismatched terminus is misaligned relative to the correct incoming nucleotide to deter or delay further DNA synthesis.

Exploring the mechanism of DNA polymerases by analyzing the effect of mutations of active site acidic groups in Polymerase β

Elucidating the catalytic mechanism of DNA polymerase is crucial for a progress in the understanding of the control of replication fidelity. This work tries to advance the mechanistic understanding by analyzing the observed effect of mutations of the acidic groups in the active site of Polymerase β as well as the pH effect on the rate constant. The analysis involves both empirical valence bond (EVB) free energy calculations and considerations of the observed pH dependence of the reaction. The combined analysis indicates that the proton transfer (PT) from the nucleophilic O3' has two possible pathways, one to D256 and the second to the bulk. We concluded based on calculations and the experimental pH profile that the most likely path for the wild-type (WT) and the D256E and D256A mutants is a PT to the bulk, although the WT may also use a PT to Asp 256. Our analysis highlights the need for very extensive sampling in the calculations of the activation barrier and also clearly shows that ab initio QM/MM calculations that do not involve extensive sampling are unlikely to give a clear quantitative picture of the reaction mechanism. Proteins 2016; 84:1644-1657. © 2016 Wiley Periodicals, Inc.

Requirement for transient metal ions revealed through computational analysis for DNA polymerase going in reverse

DNA polymerases facilitate faithful insertion of nucleotides, a central reaction occurring during DNA replication and repair. DNA synthesis (forward reaction) is "balanced," as dictated by the chemical equilibrium by the reverse reaction of pyrophosphorolysis. Two closely spaced divalent metal ions (catalytic and nucleotide-binding metals) provide the scaffold for these reactions. The catalytic metal lowers the pKa of O3' of the growing primer terminus, and the nucleotide-binding metal facilitates substrate binding. Recent time-lapse crystallographic studies of DNA polymerases have identified an additional metal ion (product metal) associated with pyrophosphate formation, leading to the suggestion of its possible involvement in the reverse reaction. Here, we establish a rationale for a role of the product metal using quantum mechanical/molecular mechanical calculations of the reverse reaction in the confines of the DNA polymerase β active site. Additionally, site-directed mutagenesis identifies essential residues and metal-binding sites necessary for pyrophosphorolysis. The results indicate that the catalytic metal site must be occupied by a magnesium ion for pyrophosphorolysis to occur. Critically, the product metal site is occupied by a magnesium ion early in the pyrophosphorolysis reaction path but must be removed later. The proposed dynamic nature of the active site metal ions is consistent with crystallographic structures. The transition barrier for pyrophosphorolysis was estimated to be significantly higher than that for the forward reaction, consistent with kinetic activity measurements of the respective reactions. These observations provide a framework to understand how ions and active site changes could modulate the internal chemical equilibrium of a reaction that is central to genome stability.

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

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

The spontaneous replication error and the mismatch discrimination mechanisms of human DNA polymerase β

To provide molecular-level insights into the spontaneous replication error and the mismatch discrimination mechanisms of human DNA polymerase β (polβ), we report four crystal structures of polβ complexed with dG•dTTP and dA•dCTP mismatches in the presence of Mg2+ or Mn2+. The Mg(2+)-bound ground-state structures show that the dA•dCTP-Mg2+ complex adopts an 'intermediate' protein conformation while the dG•dTTP-Mg2+ complex adopts an open protein conformation. The Mn(2+)-bound 'pre-chemistry-state' structures show that the dA•dCTP-Mn2+ complex is structurally very similar to the dA•dCTP-Mg2+ complex, whereas the dG•dTTP-Mn2+ complex undergoes a large-scale conformational change to adopt a Watson-Crick-like dG•dTTP base pair and a closed protein conformation. These structural differences, together with our molecular dynamics simulation studies, suggest that polβ increases replication fidelity via a two-stage mismatch discrimination mechanism, where one is in the ground state and the other in the closed conformation state. In the closed conformation state, polβ appears to allow only a Watson-Crick-like conformation for purine•pyrimidine base pairs, thereby discriminating the mismatched base pairs based on their ability to form the Watson-Crick-like conformation. Overall, the present studies provide new insights into the spontaneous replication error and the replication fidelity mechanisms of polβ.

Transition-state destabilization reveals how human DNA polymerase β proceeds across the chemically unstable lesion N7-methylguanine

N7-Methyl-2′-deoxyguanosine (m7dG) is the predominant lesion formed by methylating agents. A systematic investigation on the effect of m7dG on DNA replication has been difficult due to the chemical instability of m7dG. To gain insights into the m7dG effect, we employed a 2′-fluorine-mediated transition-state destabilzation strategy. Specifically, we determined kinetic parameters for dCTP insertion opposite a chemically stable m7dG analogue, 2′-fluoro-m7dG (Fm7dG), by human DNA polymerase β (polβ) and solved three X-ray structures of polβ in complex with the templating Fm7dG paired with incoming dCTP or dTTP analogues. The kinetic studies reveal that the templating Fm7dG slows polβ catalysis ∼300-fold, suggesting that m7dG in genomic DNA may impede replication by some DNA polymerases. The structural analysis reveals that Fm7dG forms a canonical Watson–Crick base pair with dCTP, but metal ion coordination is suboptimal for catalysis in the polβ-Fm7dG:dCTP complex, which partially explains the slow insertion of dCTP opposite Fm7dG by polβ. In addition, the polβ-Fm7dG:dTTP structure shows open protein conformations and staggered base pair conformations, indicating that N7-methylation of dG does not promote a promutagenic replication. Overall, the first systematic studies on the effect of m7dG on DNA replication reveal that polβ catalysis across m7dG is slow, yet highly accurate.

Metal-dependent conformational activation explains highly promutagenic replication across O6-methylguanine by human DNA polymerase β

Human DNA polymerase β (polβ) inserts, albeit slowly, T opposite the carcinogenic lesion O6-methylguanine (O6MeG) ∼30-fold more frequently than C. To gain insight into this promutagenic process, we solved four ternary structures of polβ with an incoming dCTP or dTTP analogue base-paired with O6MeG in the presence of active-site Mg²⁺ or Mn²⁺. The Mg²⁺-bound structures show that both the O6MeG·dCTP/dTTP–Mg²⁺ complexes adopt an open protein conformation, staggered base pair, and one active-site metal ion. The Mn²⁺-bound structures reveal that, whereas the O6Me·dCTP–Mn²⁺ complex assumes the similar altered conformation, the O6MeG·dTTP–Mn²⁺ complex adopts a catalytically competent state with a closed protein conformation and pseudo-Watson–Crick base pair. On the basis of these observations, we conclude that polβ slows nucleotide incorporation opposite O6MeG by inducing an altered conformation suboptimal for catalysis and promotes mutagenic replication by allowing Watson–Crick-mode for O6MeG·T but not for O6MeG·C in the enzyme active site. The O6MeG·dTTP–Mn²⁺ ternary structure, which represents the first structure of mismatched polβ ternary complex with a closed protein conformation and coplanar base pair, the first structure of pseudo-Watson–Crick O6MeG·T formed in the active site of a DNA polymerase, and a rare, if not the first, example of metal-dependent conformational activation of a DNA polymerase, indicate that catalytic metal-ion coordination is utilized as a kinetic checkpoint by polβ and is crucial for the conformational activation of polβ. Overall, our structural studies not only explain the promutagenic polβ catalysis across O6MeG but also provide new insights into the replication fidelity of polβ.

Structural basis for promutagenicity of 8-halogenated guanine

8-Halogenated guanine (haloG), a major DNA adduct formed by reactive halogen species during inflammation, is a promutagenic lesion that promotes misincorporation of G opposite the lesion by various DNA polymerases. Currently, the structural basis for such misincorporation is unknown. To gain insights into the mechanism of misincorporation across haloG by polymerase, we determined seven x-ray structures of human DNA polymerase β (polβ) bound to DNA bearing 8-bromoguanine (BrG). We determined two pre-catalytic ternary complex structures of polβ with an incoming nonhydrolyzable dGTP or dCTP analog paired with templating BrG. We also determined five binary complex structures of polβ in complex with DNA containing BrG·C/T at post-insertion and post-extension sites. In the BrG·dGTP ternary structure, BrG adopts syn conformation and forms Hoogsteen base pairing with the incoming dGTP analog. In the BrG·dCTP ternary structure, BrG adopts anti conformation and forms Watson-Crick base pairing with the incoming dCTP analog. In addition, our polβ binary post-extension structures show Hoogsteen BrG·G base pair and Watson-Crick BrG·C base pair. Taken together, the first structures of haloG-containing DNA bound to a protein indicate that both BrG·G and BrG·C base pairs are accommodated in the active site of polβ. Our structures suggest that Hoogsteen-type base pairing between G and C8-modified G could be accommodated in the active site of a DNA polymerase, promoting G to C mutation.

Clustered DNA lesions containing 5-formyluracil and AP site: repair via the BER system

Lesions in the DNA arise under ionizing irradiation conditions or various chemical oxidants as a single damage or as part of a multiply damaged site within 1–2 helical turns (clustered lesion). Here, we explored the repair opportunity of the apurinic/apyrimidinic site (AP site) composed of the clustered lesion with 5-formyluracil (5-foU) by the base excision repair (BER) proteins. We found, that if the AP site is shifted relative to the 5-foU of the opposite strand, it could be repaired primarily via the short-patch BER pathway. In this case, the cleavage efficiency of the AP site-containing DNA strand catalyzed by human apurinic/apyrimidinic endonuclease 1 (hAPE1) decreased under AP site excursion to the 3'-side relative to the lesion in the other DNA strand. DNA synthesis catalyzed by DNA polymerase lambda was more accurate in comparison to the one catalyzed by DNA polymerase beta. If the AP site was located exactly opposite 5-foU it was expected to switch the repair to the long-patch BER pathway. In this situation, human processivity factor hPCNA stimulates the process.

Amino acid substitution in the active site of DNA polymerase β explains the energy barrier of the nucleotidyl transfer reaction

DNA polymerase β (pol β) is a bifunctional enzyme widely studied for its roles in base excision DNA repair, where one key function is gap-filling DNA synthesis. In spite of significant progress in recent years, the atomic level mechanism of the DNA synthesis reaction has remained poorly understood. Based on crystal structures of pol β in complex with its substrates and theoretical considerations of amino acids and metals in the active site, we have proposed that a nearby carboxylate group of Asp256 enables the reaction by accepting a proton from the primer O3'group, thus activating O3'as the nucleophile in the reaction path. Here, we tested this proposal by altering the side chain of Asp256 to Glu and then exploring the impact of this conservative change on the reaction. The D256E enzyme is more than 1000-fold less active than the wild-type enzyme, and the crystal structures are subtly different in the active sites of the D256E and wild-type enzymes. Theoretical analysis of DNA synthesis by the D256E enzyme shows that the O3'proton still transfers to the nearby carboxylate of residue 256. However, the electrostatic stabilization and location of the O3' proton transfer during the reaction path are dramatically altered compared with wild-type. Surprisingly, this is due to repositioning of the Arg254 side chain in the Glu256 enzyme active site, such that Arg254 is not in position to stabilize the proton transfer from O3'. The theoretical results with the wild-type enzyme indicate an early charge reorganization associated with the O3' proton transfer, and this does not occur in the D256E enzyme. The charge reorganization is mediated by the catalytic magnesium ion in the active site.

Substrate recognition by norovirus polymerase: microsecond molecular dynamics study

Molecular dynamics simulations of complexes between Norwalk virus RNA dependent RNA polymerase and its natural CTP and 2dCTP (both containing the O5'-C5'-C4'-O4' sequence of atoms bridging the triphosphate and sugar moiety) or modified coCTP (C5'-O5'-C4'-O4'), cocCTP (C5'-O5'-C4'-C4'') substrates were produced by means of CUDA programmable graphical processing units and the ACEMD software package. It enabled us to gain microsecond MD trajectories clearly showing that similar nucleoside triphosphates can bind surprisingly differently into the active site of the Norwalk virus RNA dependent RNA polymerase. It corresponds to their different modes of action (CTP-substrate, 2dCTP-poor substrate, coCTP-chain terminator, cocCTP-inhibitor). Moreover, extremely rare events-as repetitive pervasion of Arg182 into a potentially reaction promoting arrangement-were captured.

Why nature really chose phosphate

Phosphoryl transfer plays key roles in signaling, energy transduction, protein synthesis, and maintaining the integrity of the genetic material. On the surface, it would appear to be a simple nucleophile displacement reaction. However, this simplicity is deceptive, as, even in aqueous solution, the low-lying d-orbitals on the phosphorus atom allow for eight distinct mechanistic possibilities, before even introducing the complexities of the enzyme catalyzed reactions. To further complicate matters, while powerful, traditional experimental techniques such as the use of linear free-energy relationships (LFER) or measuring isotope effects cannot make unique distinctions between different potential mechanisms. A quarter of a century has passed since Westheimer wrote his seminal review, 'Why Nature Chose Phosphate' (Science 235 (1987), 1173), and a lot has changed in the field since then. The present review revisits this biologically crucial issue, exploring both relevant enzymatic systems as well as the corresponding chemistry in aqueous solution, and demonstrating that the only way key questions in this field are likely to be resolved is through careful theoretical studies (which of course should be able to reproduce all relevant experimental data). Finally, we demonstrate that the reason that nature really chose phosphate is due to interplay between two counteracting effects: on the one hand, phosphates are negatively charged and the resulting charge-charge repulsion with the attacking nucleophile contributes to the very high barrier for hydrolysis, making phosphate esters among the most inert compounds known. However, biology is not only about reducing the barrier to unfavorable chemical reactions. That is, the same charge-charge repulsion that makes phosphate ester hydrolysis so unfavorable also makes it possible to regulate, by exploiting the electrostatics. This means that phosphate ester hydrolysis can not only be turned on, but also be turned off, by fine tuning the electrostatic environment and the present review demonstrates numerous examples where this is the case. Without this capacity for regulation, it would be impossible to have for instance a signaling or metabolic cascade, where the action of each participant is determined by the fine-tuned activity of the previous piece in the production line. This makes phosphate esters the ideal compounds to facilitate life as we know it.

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.

Effect of β,γ-CHF- and β,γ-CHCl-dGTP halogen atom stereochemistry on the transition state of DNA polymerase β

Recently, we synthesized the first individual β,γ-CHX-dGTP diastereomers [(R)- or (S)-CHX, where X is F or Cl] and determined their structures in ternary complexes with DNA polymerase β (pol β). We now report stereospecificity by pol β on the mixed β,γ-CHX diastereomer pairs using nuclear magnetic resonance and on the separate diastereomers using transient kinetics. For both the F and Cl diastereomers, the R isomer is favored over the S isomer for G·C correct incorporation, with stereospecificities [(k(pol)/K(d))(R)/(k(pol)/K(d))(S)] of 3.8 and 6.3, respectively, and also for G·T misincorporation, with stereospecificities of 11 and 7.8, respectively. Stereopreference for the (R)-CHF-dGTP diastereomer was abolished for k(pol) but not K(d) with mutant pol β (R183A). These compounds constitute a new class of stereochemical probes for active site interactions involving halogen atoms. As Arg183 is unique in family X pols, the design of CXY deoxyribonucleotide analogues to enhance interaction is a possible strategy for inhibiting BER selectively in cancer cells.

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.

β,γ-CHF- and β,γ-CHCl-dGTP diastereomers: synthesis, discrete 31P NMR signatures, and absolute configurations of new stereochemical probes for DNA polymerases

Deoxynucleoside 5'-triphosphate analogues in which the β,γ-bridging oxygen has been replaced with a CXY group are useful chemical probes to investigate DNA polymerase catalytic and base-selection mechanisms. A limitation of such probes has been that conventional synthetic methods generate a mixture of diastereomers when the bridging carbon substitution is nonequivalent (X ≠ Y). We report here a general solution to this long-standing problem with four examples of β,γ-CXY dNTP diastereomers: (S)- and (R)-β,γ-CHCl-dGTP (12a-1/12a-2) and (S)- and (R)-β,γ-CHF-dGTP (12b-1/12b-2). Central to their preparation was conversion of the prochiral parent bisphosphonic acids to the P,C-dimorpholinamide derivatives 7 of their (R)-mandelic acid monoesters, which provided access to the individual diastereomers 7a-1, 7a-2, 7b-1, and 7b-2 by preparative HPLC. Selective acidic hydrolysis of the P-N bond then afforded "portal" diastereomers, which were readily coupled to morpholine-activated dGMP. Removal of the chiral auxiliary by H(2) (Pd/C) gave the four individual diastereomeric nucleotides 12, which were characterized by (31)P, (1)H, and (19)F NMR spectroscopy and by mass spectrometry. After treatment with Chelex-100 to remove traces of paramagnetic ions, at pH ~10 the diastereomer pairs 12a,b exhibit discrete P(α) and P(β)(31)P resonances. The more upfield P(α) and more downfield P(β) resonances (and also the more upfield (19)F NMR resonance in 12b) are assigned to the R configuration at the P(β)-CHX-P(γ) carbons on the basis of the absolute configurations of the individual diastereomers as determined from the X-ray crystallographic structures of their ternary complexes with DNA and polymerase β.

Bidentate and tridentate metal-ion coordination states within ternary complexes of RB69 DNA polymerase

Two divalent metal ions are required for primer-extension catalyzed by DNA polymerases. One metal ion brings the 3'-hydroxyl of the primer terminus and the α-phosphorus atom of incoming dNTP together for bond formation so that the catalytically relevant conformation of the triphosphate tail of the dNTP is in an α,β,γ-tridentate coordination complex with the second metal ion required for proper substrate alignment. A probable base selectivity mechanism derived from structural studies on Dpo4 suggests that the inability of mispaired dNTPs to form a substrate-aligned, tridentate coordination complex could effectively cause the mispaired dNTPs to be rejected before catalysis. Nevertheless, we found that mispaired dNTPs can actually form a properly aligned tridentate coordination complex. However, complementary dNTPs occasionally form misaligned complexes with mutant RB69 DNA polymerases (RB69pols) that are not in a tridentate coordination state. Here, we report finding a β,γ-bidentate coordination complex that contained the complementary dUpNpp opposite dA in the structure of a ternary complex formed by the wild type RB69pol at 1.88 Å resolution. Our observations suggest that several distinct metal-ion coordination states can exist at the ground state in the polymerase active site and that base selectivity is unlikely to be based on metal-ion coordination alone.

Mechanism for the autophosphorylation of CheA histidine kinase: QM/MM calculations

The CheA histidine kinase, a model of TCS (the two-component system), mediates the signal transduction pathway of bacterial chemotaxis via autophosphorylation. Since the TCSs are rarely found in mammalians, they have become attractive targets for the development of new antibiotics. To characterize the autophosphoryl-transfer mechanism of CheA histidine kinase, molecular dynamics simulations combined with quantum mechanics/molecular mechanics calculations were employed on the constructed 3D model of P1-P4-ATP complex. A two-step reaction mechanism was proposed and confirmed by our computations: the autophosphoryl-transfer reaction takes place followed by a rapid and reversible conformational change from ground state to prechemistry state. In addition, a two-dimensional potential energy surface was calculated for autophosphorylation, and the transition state displays an associative character. Moreover, we found Lys48 serves as the catalytic acid to stabilize transition state through a water-mediated proton-transfer pathway, and Glu67 acts as not only a hydrogen bond acceptor but also a structure anchor to modulate the imidazole ring of His45 in the active site. Our findings clearly provide a detailed autophosphoryl-transfer mechanism of CheA histidine kinase and thus are important for discovering new antibiotics.

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.

An abridged transition state model to derive structure, dynamics, and energy components of DNA polymerase β fidelity

We show how a restricted reaction surface can be used to facilitate the calculation of biologically important contributions of active site geometries and dynamics to DNA polymerase fidelity. Our analysis, using human DNA polymerase beta (pol β), is performed within the framework of an electrostatic linear free energy response (EFER) model. The structure, dynamics, and energetics of pol β-DNA-dNTP interactions are computed between two points on the multidimensional reaction free energy surface. "Point 1" represents a ground state activation intermediate (GSA), which is obtained by deprotonating the terminal 3'OH group of the primer DNA strand. "Point 2" is the transition state (PTS) for the attack of the 3'O(-) (O(nuc)) on the P(α) atom of dNTP substrate, having the electron density of a dianionic phosphorane intermediate. Classical molecular dynamics simulations are used to compute the geometric and dynamic contributions to the formation of right and wrong O(nuc)-P chemical bonds. Matched dCTP·G and mismatched dATP·G base pairs are used to illustrate the analysis. Compared to the dCTP·G base pair, the dATP·G mismatch has fewer GSA configurations with short distances between O(nuc) and P(α) atoms and between the oxygen in the scissile P-O bond (O(lg)) and the nearest structural water. The thumb subdomain conformation of the GSA complex is more open for the mismatch, and the H-bonds in the mispair become more extended during the nucleophilic attack than in the correct pair. The electrostatic contributions of pol β and DNA residues to catalysis of the right and wrong P-O(nuc) bond formation are 5.3 and 3.1 kcal/mol, respectively, resulting in an 80-fold contribution to fidelity. The EFER calculations illustrate the considerable importance of Arg183 and an O(lg)-proximal water molecule to pol β fidelity.

Lee Pedersen’s work in theoretical and computational chemistry and biochemistry

Nature at the lab level in biology and chemistry can be described by the application of quantum mechanics. In many cases, a reasonable approximation to quantum mechanics is classical mechanics realized through Newton’s equations of motion. Dr. Pedersen began his career using quantum mechanics to describe the properties of small molecular complexes that could serve as models for biochemical systems. To describe large molecular systems required a drop-back to classical means and this led surprisingly to a major improvement in the classical treatment of electrostatics for all molecules, not just biological molecules. Recent work has involved the application of quantum mechanics for the putative active sites of enzymes to gain greater insight into the key steps in enzyme catalysis.

Replication infidelity via a mismatch with Watson-Crick geometry

In describing the DNA double helix, Watson and Crick suggested that "spontaneous mutation may be due to a base occasionally occurring in one of its less likely tautomeric forms." Indeed, among many mispairing possibilities, either tautomerization or ionization of bases might allow a DNA polymerase to insert a mismatch with correct Watson-Crick geometry. However, despite substantial progress in understanding the structural basis of error prevention during polymerization, no DNA polymerase has yet been shown to form a natural base-base mismatch with Watson-Crick-like geometry. Here we provide such evidence, in the form of a crystal structure of a human DNA polymerase λ variant poised to misinsert dGTP opposite a template T. All atoms needed for catalysis are present at the active site and in positions that overlay with those for a correct base pair. The mismatch has Watson-Crick geometry consistent with a tautomeric or ionized base pair, with the pH dependence of misinsertion consistent with the latter. The results support the original idea that a base substitution can originate from a mismatch having Watson-Crick geometry, and they suggest a common catalytic mechanism for inserting a correct and an incorrect nucleotide. A second structure indicates that after misinsertion, the now primer-terminal G • T mismatch is also poised for catalysis but in the wobble conformation seen in other studies, indicating the dynamic nature of the pathway required to create a mismatch in fully duplex DNA.

Computational delineation of the catalytic step of a high-fidelity DNA polymerase

The Bacillus fragment, belonging to a class of high-fidelity polymerases, demonstrates high processivity (adding approximately 115 bases per DNA binding event) and exceptional accuracy (1 error in 10(6) nucleotide incorporations) during DNA replication. We present analysis of structural rearrangements and energetics just before and during the chemical step (phosphodiester bond formation) using a combination of classical molecular dynamics, mixed quantum mechanics molecular mechanics simulations, and free energy computations. We find that the reaction is associative, proceeding via the two-metal-ion mechanism, and requiring the proton on the terminal primer O3' to transfer to the pyrophosphate tail of the incoming nucleotide before the formation of the pentacovalent transition state. Different protonation states for key active site residues direct the system to alternative pathways of catalysis and we estimate a free energy barrier of approximately 12 kcal/mol for the chemical step. We propose that the protonation of a highly conserved catalytic aspartic acid residue is essential for the high processivity demonstrated by the enzyme and suggest that global motions could be part of the reaction free energy landscape.

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.

A computational study of the hydrolysis of dGTP analogues with halomethylene-modified leaving groups in solution: implications for the mechanism of DNA polymerases

DNA polymerases make up a family of enzymes responsible for regulating DNA replication and repair, which in turn maintains the integrity of the genome. However, despite intensive kinetic, crystallographic, and computational studies, elucidation of the detailed enzymatic mechanism still presents a significant challenge. We recently developed an alternative strategy for exploring the fidelity and mechanism of DNA polymerases, by probing leaving group effects on nucleotidyl transfer using a series of dGTP bisphosphonate analogues in which the beta,gamma-bridging oxygen was replaced by a series of substituted methylene groups (X = CYZ, where Y and Z = H, halogen, or another substituent). Pre-steady state kinetic measurements of DNA polymerase-catalyzed incorporation of correctly base paired (R) and mispaired (W) analogues demonstrated a strong linear free energy relationship (LFER) between the polymerase rate constant (k(pol)) and the highest pK(a) of the free bisphosphonic acid corresponding to the leaving group. However, unexpectedly, the data segregated into two distinctly different linear correlations depending on the nature of the substituent. The discrepancy between the two lines was considerably greater when the dGTP analogue formed an incorrect (G.T) rather than a correct (G.C) base pair, although the reason for this phenomenon remains unexplained. Here, we have evaluated the complete free energy surfaces for bisphosphonate hydrolysis in aqueous solution and evaluated the corresponding LFER. Our study, which employs several alternative solvation models, finds a split of the calculated LFER for the mono- and dihalogen compounds into two parallel lines, reflecting their behavior in the polymerase-catalyzed condensation reaction. We suggest that the division into two linear subsets may be a generalized solvation phenomenon involving the overall electrostatic interaction between the substrates and their surroundings and would also be observed in polar solvents in the absence of the enzyme, if the reaction in solvent is in fact identical to that of the enzyme. However, the amplified differences between the LFER lines for the incorporation of matched and mismatched deoxynucleotides probably reflects the differences in the electrostatic interaction between the TS charges in the polymerase active site. An understanding of the mechanism of this reaction in solution could thereby provide a steppingstone for understanding the factors governing the fidelity of DNA 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.

Polymerase-tailored variations in the water-mediated and substrate-assisted mechanism for nucleotidyl transfer: insights from a study of T7 DNA polymerase

The nucleotidyl transfer reaction catalyzed by DNA polymerases is the critical step governing the accurate transfer of genetic information during DNA replication, and its malfunctioning can cause mutations leading to human diseases, including cancer. Here, utilizing ab initio quantum mechanical/molecular mechanical calculations with free-energy perturbation, we carried out an extensive investigation of the nucleotidyl transfer reaction mechanism in the well-characterized high-fidelity replicative DNA polymerase from phage T7. Our defined mechanism entails an initial concerted deprotonation of a conserved crystal water molecule with protonation of the gamma-phosphate of the deoxynucleotide triphosphate(dNTP) via a solvent water molecule, and then the proton on the primer 3'-terminus is transferred to the resulting hydroxide ion. Subsequently, the nucleophilic attack takes place, with the formation of a metastable pentacovalent phosphorane intermediate. Finally, the pyrophosphate leaves, facilitated by the relay of the proton on the gamma-phosphate to the alpha-beta bridging oxygen via solvent water. The computed activation free-energy barrier is consistent with kinetic data for the chemistry step with correct nucleotide incorporation in T7 DNA polymerase. This variant of the water-mediated and substrate-assisted mechanism has features tailored to the structure of the T7 DNA polymerase. However, a unifying theme in the water-mediated and substrate-assisted mechanism is the cycling through crystal and solvent water molecules of the proton originating from the primer 3'-terminus to the alpha-beta bridging oxygen of the deoxynucleotide triphosphate; this neutralizes the evolving negative charge as pyrophosphate leaves and restores the polymerase to its pre-chemistry state. These unifying features are likely requisite elements for nucleotidyl transfer reactions.

DNA nucleobase properties and photoreactivity: Modeling environmental effects

The accurate ab initio quantum chemical (QM) method multiconfigurational second-order perturbation (CASSPT2)/complete active space self-consistent field (CASSCF) has been used in conjunction with molecular mechanics (MM) procedures to compute molecular properties and photoinduced reactivity of DNA/RNA nucleobases (NABs) in isolation and within a realistic environment, in which the double helix strand, the aqueous media, and the external counterions are included. It is illustrated that the use of an MM model is helpful both to account for short- and long-range effects of the system surrounding the QM molecular core and to provide the proper structural constraints that allow more accurate QM geometry determinations.

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.