Crystal structures of human DNA polymerase beta complexed with gapped and nicked DNA: evidence for an induced fit mechanism

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

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

Characterization of the RNA binding energetics of the Candida albicans poly(A) polymerase

The 3' ends of eukaryotic mRNAs are characterized by the presence of a poly(A) tail, which plays a critical role in stability, transport, and translation of the mRNAs. In the present study, we report the expression, purification and enzymatic characterization of the poly(A) polymerase of Candida albicans, an important human pathogen. As a first step toward elucidating the nature of the interaction between RNA and the enzyme, fluorescence spectroscopy assays were also performed to monitor the binding of RNA to the protein. Our assays revealed that the initial interaction between RNA and the enzyme is characterized by a high enthalpy of association and that the minimal RNA binding site of the enzyme is eight nucleotides. Moreover, both the kinetics of real-time RNA binding and the contribution of electrostatic interactions to the overall binding energy were investigated. Finally, we also correlated the effect of RNA binding on protein structure, using both circular dichroism and guanidium hydrochloride-induced denaturation studies as structural indicators. Our data indicate that the protein undergoes structural modifications upon RNA binding, although the interaction does not significantly modify the stability of the protein. In addition to the determination of the energetics of RNA binding, our study provides a better understanding of the molecular basis of RNA binding by poly(A) polymerases.

Molecular dynamic simulations of cisplatin- and oxaliplatin-d(GG) intrastand cross-links reveal differences in their conformational dynamics

Mismatch repair proteins, DNA damage-recognition proteins and translesion DNA polymerases discriminate between Pt-GG adducts containing cis-diammine ligands (formed by cisplatin (CP) and carboplatin) and trans-RR-diaminocyclohexane ligands (formed by oxaliplatin (OX)) and this discrimination is thought to be important in determining differences in the efficacy, toxicity and mutagenicity of these platinum anticancer agents. We have postulated that these proteins recognize differences in conformation and/or conformational dynamics of the DNA containing the adducts. We have previously determined the NMR solution structure of OX-DNA, CP-DNA and undamaged duplex DNA in the 5'-d(CCTCAGGCCTCC)-3' sequence context and have shown the existence of several conformational differences in the vicinity of the Pt-GG adduct. Here we have used molecular dynamics simulations to explore differences in the conformational dynamics between OX-DNA, CP-DNA and undamaged DNA in the same sequence context. Twenty-five 10 ns unrestrained fully solvated molecular dynamics simulations were performed starting from two different DNA conformations using AMBER v8.0. All 25 simulations reached equilibrium within 4 ns, were independent of the starting structure and were in close agreement with previous crystal and NMR structures. Our data show that the cis-diammine (CP) ligand preferentially forms hydrogen bonds on the 5' side of the Pt-GG adduct, while the trans-RR-diaminocyclohexane (OX) ligand preferentially forms hydrogen bonds on the 3' side of the adduct. In addition, our data show that these differences in hydrogen bond formation are strongly correlated with differences in conformational dynamics, specifically the fraction of time spent in different DNA conformations in the vicinity of the adduct, for CP- and OX-DNA adducts. We postulate that differential recognition of CP- and OX-GG adducts by mismatch repair proteins, DNA damage-recognition proteins and DNA polymerases may be due, in part, to differences in the fraction of time that the adducts spend in a conformation favorable for protein binding.

DNA polymerase beta catalysis: are different mechanisms possible?

DNA polymerases are crucial constituents of the complex cellular machinery for replicating and repairing DNA. Discerning mechanistic pathways of DNA polymerase on the atomic level is important for revealing the origin of fidelity discrimination. Mammalian DNA polymerase beta (pol beta), a small (39 kDa) member of the X-family, represents an excellent model system to investigate polymerase mechanisms. Here, we explore several feasible low-energy pathways of the nucleotide transfer reaction of pol beta for correct (according to Watson-Crick hydrogen bonding) G:C basepairing versus the incorrect G:G case within a consistent theoretical framework. We use mixed quantum mechanics/molecular mechanics (QM/MM) techniques in a constrained energy minimization protocol to effectively model not only the reactive core but also the influence of the rest of the enzymatic environment and explicit solvent on the reaction. The postulated pathways involve initial proton abstraction from the terminal DNA primer O3'H group, nucleophilic attack that extends the DNA primer chain, and elimination of pyrophosphate. In particular, we analyze several possible routes for the initial deprotonation step: (i) direct transfer to a phosphate oxygen O(Palpha) of the incoming nucleotide, (ii) direct transfer to an active site Asp group, and (iii) transfer to explicit water molecules. We find that the most probable initial step corresponds to step (iii), involving initial deprotonation to water, which is followed by proton migration to active site Asp residues, and finally to the leaving pyrophosphate group, with an activation energy of about 15 kcal/mol. We argue that initial deprotonation steps (i) and (ii) are less likely as they are at least 7 and 11 kcal/mol, respectively, higher in energy. Overall, the rate-determining step for both the correct and the incorrect nucleotide cases is the initial deprotonation in concert with nucleophilic attack at the phosphate center; however, the activation energy we obtain for the mismatched G:G case is 5 kcal/mol higher than that of the matched G:C complex, due to active site structural distortions. Taken together, our results support other reported mechanisms and help define a framework for interpreting nucleotide specificity differences across polymerase families, in terms of the concept of active site preorganization or the so-called "pre-chemistry avenue".

Characterization of terminal deoxynucleotidyl transferase and polymerase mu in zebrafish

Terminal deoxynucleotidyl transferase (TdT) contributes to the junctional diversity of immunoglobulin and T-cell receptors by incorporating nucleotides in a template-independent manner. A closely related enzyme, polymerase mu (polmu), a template-directed polymerase, plays a role in general end-joining double-strand break repair. We cloned zebrafish TdT and polmu and found them to be 43% identical in amino acid sequence. Comparisons with sequences of other species revealed conserved residues typical for TdT in the zebrafish sequence that support the template independence of this enzyme. Some but not all of these features were identified in zebrafish polmu. In adult fish, TdT expression was most prominent in thymus, pro- and mesonephros, the primary lymphoid organs in teleost fish and in spleen, intestine, and the tissue around the intestine. Polmu expression was detected not only in pro- and mesonephros, the major sites for B-lymphocyte development, but also in ovary and testis and in all tissue preparations to a low extent. TdT expression starts at 4 dpf and increases thereafter. Polmu is expressed at all times to a similar extent. In situ studies showed a strong expression of TdT and polmicro in the thymic cortex of 8-week-old fish. The characterization of zebrafish TdT and polmu provide new insights in fish lymphopoiesis and addresses the importance and evolution of TdT and polmu themselves.

Free energy simulations of uncatalyzed DNA replication fidelity: structure and stability of T.G and dTTP.G terminal DNA mismatches flanked by a single dangling nucleotide

A reference system for DNA replication fidelity was studied by free energy perturbation (FEP) and linear interaction energy (LIE) methods. The studied system included a hydrated duplex DNA with the 5'-CG dangling end of the templating strand, and dCTP4-.Mg2+ or dTTP4-.Mg2+ inserted opposite the dangling G to form a correct (i.e., Watson-Crick) or incorrect (i.e., wobble) base pair, respectively. The average distance between the 3'-terminal oxygen of the primer strand and the alpha-phosphorus of dNTP was found to be 0.2 A shorter for the correct base pair than for the incorrect base pair. Binding of the incorrect dNTP was found to be disfavored by 0.4 kcal/mol relative to the correct dNTP. We estimated that improved binding and more near-attack configurations sampled by the correct base pair should translate in aqueous solution and in the absence of DNA polymerase into a six times faster rate for the incorporation of the correct dNTP into DNA. The accuracy of the calculated binding free energy difference was verified by examining the relative free energy for melting duplex DNA containing GC and GT terminal base pairs flanked by a 5' dangling C. The calculated LIE and FEP free energies of 1.7 and 1.1 kcal/mol, respectively, compared favorably with the experimental estimate of 1.4 kcal/mol obtained using the nearest neighbor parameters. To decompose the calculated free energies into additive electrostatic and van der Waals contributions and to provide a set of rigorous theoretical data for the parametrization of the LIE method, we suggested a variant of the FEP approach, for which we coined a binding-relevant free energy (BRFE) acronym. BRFE approach is characterized by its unique perturbation pathway and by its exclusion of the intramolecular energy of a rigid part of the ligand from the total potential energy.

Decomposition of the solvation free energies of deoxyribonucleoside triphosphates using the free energy perturbation method

Free energy perturbation (FEP) calculations using the Amber 95 force field and the TIP3P water model were carried out to evaluate the solvation free energy of deoxyribonucleoside triphosphates in aqueous solution. Solvation free energies of -307.5, -311.5, -314.1, and -317.0 kcal/mol were calculated for the (Mg x dTTP)2-, (Mg x dATP)2-, (Mg x dCTP)2-, and (Mg x dGTP)2- complexes, respectively. Structural origins of the relative solvation free energies of deoxyribonucleoside phosphates were examined by calculating the contribution of the interaction of the base moiety with its surroundings. We showed that for each nucleobase the magnitude of this contribution is unaffected by substituting the 5'-OH group of the corresponding nucleoside with the charged mono- or triphosphate groups. This free energy contribution was further decomposed into the sum of free energies originating from the interactions of the base with itself, its substituent, water, and Na+ ions. Although the sum of these components was nearly constant over a wide range of solutes the individual free energy constituents varied significantly. Furthermore, this decomposition showed a high degree of additivity. Computational conditions necessary for obtaining additive free energy decomposition for the systems studied here within the framework of the FEP method included the use of a single mutation pathway and a subdivision of the FEP protocol into 51 or more windows.

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

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

Following an environmental carcinogen N2-dG adduct through replication: elucidating blockage and bypass in a high-fidelity DNA polymerase

We have investigated how a benzo[a]pyrene-derived N²-dG adduct, 10S(+)-trans-anti-[BP]-N²-dG ([BP]G*), is processed in a well-characterized Pol I family model replicative DNA polymerase, Bacillus fragment (BF). Experimental results are presented that reveal relatively facile nucleotide incorporation opposite the lesion, but very inefficient further extension. Computational studies follow the possible bypass of [BP]G* through the pre-insertion, insertion and post-insertion sites as BF alternates between open and closed conformations. With dG* in the normal B-DNA anti conformation, BP seriously disturbs the polymerase structure, positioning itself either deeply in the pre-insertion site or on the crowded evolving minor groove side of the modified template, consistent with a polymerase-blocking conformation. With dG* in the less prevalent syn conformation, BP causes less distortion: it is either out of the pre-insertion site or in the major groove open pocket of the polymerase. Thus, the syn conformation can account for the observed relatively easy incorporation of nucleotides, with mutagenic purines favored, opposite the [BP]G* adduct. However, with the lesion in the BF post-insertion site, more serious distortions caused by the adduct even in the syn conformation explain the very inefficient extension observed experimentally. In vivo, a switch to a potentially error-prone bypass polymerase likely dominates translesion bypass.

Different structural states in oligonucleosomes are required for early versus late steps of base excision repair

Chromatin in eukaryotic cells is folded into higher order structures of folded nucleosome filaments, and DNA damage occurs at all levels of this structural hierarchy. However, little is known about the impact of higher order folding on DNA repair enzymes. We examined the catalytic activities of purified human base excision repair (BER) enzymes on uracil-containing oligonucleosome arrays, which are folded primarily into 30 nm structures when incubated in repair reaction buffers. The catalytic activities of uracil DNA glycosylase (UDG) and apyrimidinic/apurinic endonuclease (APE) digest G:U mismatches to completion in the folded oligonucleosomes without requiring significant disruption. In contrast, DNA polymerase β (Pol β) synthesis is inhibited in a major fraction (∼80%) of the oligonucleosome array, suggesting that single strand nicks in linker DNA are far more accessible to Pol β in highly folded oligonucleosomes. Importantly, this barrier in folded oligonucleosomes is removed by purified chromatin remodeling complexes. Both ISW1 and ISW2 from yeast significantly enhance Pol β accessibility to the refractory nicked sites in oligonucleosomes. These results indicate that the initial steps of BER (UDG and APE) act efficiently on highly folded oligonucleosome arrays, and chromatin remodeling may be required for the latter steps of BER in intact chromatin.

Conformational proofreading: the impact of conformational changes on the specificity of molecular recognition

To perform recognition, molecules must locate and specifically bind their targets within a noisy biochemical environment with many look-alikes. Molecular recognition processes, especially the induced-fit mechanism, are known to involve conformational changes. This raises a basic question: Does molecular recognition gain any advantage by such conformational changes? By introducing a simple statistical-mechanics approach, we study the effect of conformation and flexibility on the quality of recognition processes. Our model relates specificity to the conformation of the participant molecules and thus suggests a possible answer: Optimal specificity is achieved when the ligand is slightly off target; that is, a conformational mismatch between the ligand and its main target improves the selectivity of the process. This indicates that deformations upon binding serve as a conformational proofreading mechanism, which may be selected for via evolution.

Solution structures of a DNA dodecamer duplex with and without a cisplatin 1,2-d(GG) intrastrand cross-link: comparison with the same DNA duplex containing an oxaliplatin 1,2-d(GG) intrastrand cross-link

Proteins that discriminate between cisplatin-DNA adducts and oxaliplatin-DNA adducts are thought to be responsible for the differences in tumor range, toxicity, and mutagenicity of these two important chemotherapeutic agents. However, the structural basis for differential protein recognition of these adducts has not been determined and could be important for the design of more effective platinum anticancer agents. We have determined high-resolution NMR structures for cisplatin-GG and undamaged DNA dodecamers in the AGGC sequence context and have compared these structures with the oxaliplatin-GG structure in the same sequence context determined previously in our laboratory. This structural study allows the first direct comparison of cisplatin-GG DNA and oxaliplatin-GG DNA solution structures referenced to undamaged DNA in the same sequence context. Non-hydrogen atom rmsds of 0.81 and 1.21 were determined for the 15 lowest-energy structures for cisplatin-GG DNA and undamaged DNA, respectively, indicating good structural convergence. The theoretical NOESY spectra obtained by back-calculation from the final average structures showed excellent agreement with the experimental data, indicating that the final structures are consistent with the NMR data. Several significant conformational differences were observed between the cisplatin-GG adduct and the oxaliplatin-GG adduct, including buckle at the 5' G6.C19 base pair, opening at the 3' G7.C18 base pair, twist at the A5G6.T20C19 base pair step, slide, twist, and roll at the G6G7.C19C18 base pair step, slide at the G7C8.C18G17 base pair step, G6G7 dihedral angle, and overall bend angle. We hypothesize that these conformational differences may be related to the ability of various DNA repair proteins, DNA binding proteins, and DNA polymerases to discriminate between cisplatin-GG and oxaliplatin-GG adducts.

Loop II of DNA polymerase beta is important for polymerization activity and fidelity

The accurate replication and transmission of genetic information is critical in the life of an organism. During its entire lifespan, the genetic information is constantly under attack from endogenous and exogenous sources of damage. To ensure that the content of its genetic information is faithfully preserved for synthesis and transmission, eukaryotic cells have developed a complex system of genomic quality control. Key players in this process are DNA polymerases, the enzymes responsible for synthesizing the DNA, because errors introduced into the genome by polymerase can result in mutations. We use DNA polymerase beta (pol β) as a model system to investigate mechanisms of preserving fidelity during nucleotide incorporation. In the study described here, we characterized the role that loop II of pol β plays in maintaining the activity and fidelity of pol β. We report here that the absence or shortening of loop II compromises the catalytic activity of pol β. Our data also show that loop variants of a specific length have a lower fidelity when compared to the wild-type polymerase. Taken together, our results indicate that loop II is important for the catalytic activity and fidelity of pol β.

Participation of the Fingers Subdomain of Escherichia coli DNA Polymerase I in the Strand Displacement Synthesis of DNA

The replication of the genome requires the removal of RNA primers from the Okazaki fragments and their replacement by DNA. In prokaryotes, this process is completed by DNA polymerase I by means of strand displacement DNA synthesis and 5 '-nuclease activity. Here, we demonstrate that the strand displacement DNA synthesis is facilitated by the collective participation of Ser(769), Phe(771), and Arg(841) present in the fingers subdomain of DNA polymerase I. The steady and presteady state kinetic analysis of the properties of appropriate mutant enzymes suggest that: (a) Ser(769) and Phe(771) together are involved in the strand separation via the formation of a flap structure, and (b) Arg(841) interacts with the template strand to achieve the optimal strand separation and DNA synthesis. The amino acid residues Ser(769) and Phe(771) are constituents of the O1-helix, which together with O and O2 helices form a 3-helix bundle structure. We note that this 3-helix bundle motif also exists in prokaryotic RNA polymerase. Thus in both DNA and RNA polymerases, this motif may have been adopted to achieve the strand separation function.

A water-mediated and substrate-assisted catalytic mechanism for Sulfolobus solfataricus DNA polymerase IV

DNA polymerases are enzymes responsible for the synthesis of DNA from nucleotides. Understanding their molecular fundamentals is a prerequisite for elucidating their aberrant activities in diseases such as cancer. Here we have carried out ab initio quantum mechanical/molecular mechanical (QM/MM) studies on the nucleotidyl-transfer reaction catalyzed by the lesion-bypass DNA polymerase IV (Dpo4) from Sulfolobus solfataricus, with template guanine and Watson-Crick paired dCTP as the nascent base pair. The results suggested a novel water-mediated and substrate-assisted (WMSA) mechanism: the initial proton transfer to the alpha-phosphate of the substrate via a bridging crystal water molecule is the rate-limiting step, the nucleotidyl-transfer step is associative with a metastable pentacovalent phosphorane intermediate, and the pyrophosphate leaving is facilitated by a highly coordinated proton relay mechanism through mediation of water which neutralizes the evolving negative charge. The conserved carboxylates, which retain their liganding to the two Mg2+ ions during the reaction process, are found to be essential in stabilizing transition states. This WMSA mechanism takes specific advantage of the unique structural features of this low-fidelity lesion-bypass Y-family polymerase, which has a more spacious and solvent-exposed active site than replicative and repair polymerases.

Stable complexes formed by HIV-1 reverse transcriptase at distinct positions on the primer-template controlled by binding deoxynucleoside triphosphates or foscarnet

Binding of the next complementary dNTP by the binary complex containing HIV-1 reverse transcriptase (RT) and primer-template induces conformational changes that have been implicated in catalytic function of RT. We have used DNase I footprinting, gel electrophoretic mobility shift, and exonuclease protection assays to characterize the interactions between HIV-1 RT and chain-terminated primer-template in the absence and presence of various ligands. Distinguishable stable complexes were formed in the presence of foscarnet (an analog of pyrophosphate), the dNTP complementary to the first (+1) templating nucleotide or the dNTP complementary to the second (+2) templating nucleotide. The position of HIV-1 RT on the primer-template in each of these complexes is different. RT is located upstream in the foscarnet complex, relative to the +1 complex, and downstream in the +2 complex. These results suggest that HIV-1 RT can translocate along the primer-template in the absence of phosphodiester bond formation. The ability to form a specific foscarnet complex might explain the inhibitory properties of this compound. The ability to recognize the second templating nucleotide has implications for nucleotide misincorporation.

Multiple functions of DNA polymerases

The primary role of DNA polymerases is to accurately and efficiently replicate the genome in order to ensure the maintenance of the genetic information and its faithful transmission through generations. This is not a simple task considering the size of the genome and its constant exposure to endogenous and environmental DNA damaging agents. Thus, a number of DNA repair pathways operate in cells to protect the integrity of the genome. In addition to their role in replication, DNA polymerases play a central role in most of these pathways. Given the multitude and the complexity of DNA transactions that depend on DNA polymerase activity, it is not surprising that cells in all organisms contain multiple highly specialized DNA polymerases, the majority of which have only recently been discovered. Five DNA polymerases are now recognized in Escherichia coli, 8 in Saccharomyces cerevisiae, and at least 15 in humans. While polymerases in bacteria, yeast and mammalian cells have been extensively studied much less is known about their counterparts in plants. For example, the plant model organism Arabidopsis thaliana is thought to contain 12 DNA polymerases, whose functions are mostly unknown. Here we review the properties and functions of DNA polymerases focusing on yeast and mammalian cells but paying special attention to the plant enzymes and the special circumstances of replication and repair in plant cells.

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

Background

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

Results

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

Conclusion

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

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

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

Single-turnover kinetic analysis of the mutagenic potential of 8-oxo-7,8-dihydro-2'-deoxyguanosine during gap-filling synthesis catalyzed by human DNA polymerases lambda and beta

In the presence of 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxodG) damage, many DNA polymerases exhibit a dual coding potential which facilitates efficient incorporation of matched dCTP or mismatched dATP. This also holds true for the insertion of 8-oxodGTP opposite template bases dC and dA. Employing single-turnover kinetic methods, we examined human DNA polymerase beta and its novel X-family homolog, human DNA polymerase lambda, to determine which nucleotide and template base was preferred when encountering 8-oxodG and 8-oxodGTP, respectively. While DNA polymerase beta preferentially incorporated dCTP over dATP, DNA polymerase lambda did not modulate a preference for either dCTP or dATP when opposite 8-oxodG in single-nucleotide gapped DNA, as incorporation proceeded with essentially equal efficiency and probability. Moreover, DNA polymerase lambda is more efficient than DNA polymerase beta to fill this oxidized single-nucleotide gap. Insertion of 8-oxodGTP by both DNA polymerases lambda and beta occurred predominantly against template dA, thereby reiterating how the asymmetrical design of the polymerase active site differentially accommodated the anti and syn conformations of 8-oxodG and 8-oxodGTP. Although the electronegative oxygen at the C8 position of 8-oxodG may induce DNA structural perturbations, human DNA ligase I was found to effectively ligate the incorporated 8-oxodGMP to a downstream strand, which sealed the nicked DNA. Consequently, the erroneous nucleotide incorporations catalyzed by DNA polymerases lambda and beta as well as the subsequent ligation catalyzed by a DNA ligase during base excision repair are a threat to genomic integrity.

Inhibition of hepatitis B virus polymerase by entecavir

Entecavir (ETV; Baraclude) is a novel deoxyguanosine analog with activity against hepatitis B virus (HBV). ETV differs from the other nucleoside/tide reverse transcriptase inhibitors approved for HBV therapy, lamivudine (LVD) and adefovir (ADV), in several ways: ETV is >100-fold more potent against HBV in culture and, at concentrations below 1 microM, displays no significant activity against human immunodeficiency virus (HIV). Additionally, while LVD and ADV are obligate DNA chain terminators, ETV halts HBV DNA elongation after incorporating a few additional bases. Three-dimensional homology models of the catalytic center of the HBV reverse transcriptase (RT)-DNA-deoxynucleoside triphosphate (dNTP) complex, based on the HIV RT-DNA structure, were used with in vitro enzyme kinetic studies to examine the mechanism of action of ETV against HBV RT. A novel hydrophobic pocket in the rear of the RT dNTP binding site that accommodates the exocyclic alkene moiety of ETV was predicted, establishing a basis for the superior potency observed experimentally. HBV DNA chain termination by ETV was accomplished through disfavored energy requirements as well as steric constraints during subsequent nucleotide addition. Validation of the model was accomplished through modeling of LVD resistance substitutions, which caused an eightfold decrease in ETV susceptibility and were predicted to reduce, but not eliminate, the ETV-binding pocket, in agreement with experimental observations. ADV resistance changes did not affect the ETV docking model, also agreeing with experimental results. Overall, these studies explain the potency, mechanism, and cross-resistance profile of ETV against HBV and account for the successful treatment of naive and LVD- or ADV-experienced chronic HBV patients.

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.

X-ray crystallographic and steady state fluorescence characterization of the protein dynamics of yeast polyadenylate polymerase

Polyadenylate polymerase (PAP) catalyzes the synthesis of poly(A) tails on the 3'-end of pre-mRNA. PAP is composed of three domains: an N-terminal nucleotide-binding domain (homologous to the palm domain of DNA and RNA polymerases), a middle domain (containing other conserved, catalytically important residues), and a unique C-terminal domain (involved in protein-protein interactions required for 3'-end formation). Previous X-ray crystallographic studies have shown that the domains are arranged in a V-shape such that they form a central cleft with the active site located at the base of the cleft at the interface between the N-terminal and middle domains. In the previous studies, the nucleotides were bound directly to the N-terminal domain and exhibited a conspicuous lack of adenine-specific interactions that would constitute nucleotide recognition. Furthermore, it was postulated that base-specific contacts with residues in the middle domain could occur either as a result of a change in the conformation of the nucleotide or domain movement. To address these issues and to better characterize the structural basis of substrate recognition and catalysis, we report two new crystal structures of yeast PAP. A comparison of these structures reveals that the N-terminal and C-terminal domains of PAP move independently as rigid bodies along two well defined axes of rotation. Modeling of the nucleotide into the most closed state allows us to deduce specific nucleotide interactions involving residues in the middle domain (K215, Y224 and N226) that are proposed to be involved in substrate binding and specificity. To further investigate the nature of PAP domain flexibility, 2-aminopurine labeled molecular probes were employed in steady state fluorescence and acrylamide quenching experiments. The results suggest that the closed domain conformation is stabilized upon recognition of the correct subtrate, MgATP, in an enzyme-substrate ternary complex. The implications of these results on the enzyme mechanism of PAP and the possible role for domain motion in an induced fit mechanism are discussed.

Structural insight into the substrate specificity of DNA Polymerase mu

DNA polymerase mu (Pol mu) is a family X enzyme with unique substrate specificity that contributes to its specialized role in nonhomologous DNA end joining (NHEJ). To investigate Pol mu's unusual substrate specificity, we describe the 2.4 A crystal structure of the polymerase domain of murine Pol mu bound to gapped DNA with a correct dNTP at the active site. This structure reveals substrate interactions with side chains in Pol mu that differ from other family X members. For example, a single amino acid substitution, H329A, has little effect on template-dependent synthesis by Pol mu from a paired primer terminus, but it reduces both template-independent and template-dependent synthesis during NHEJ of intermediates whose 3' ends lack complementary template strand nucleotides. These results provide insight into the substrate specificity and differing functions of four closely related mammalian family X DNA polymerases.

Methylation of DNA polymerase beta by protein arginine methyltransferase 1 regulates its binding to proliferating cell nuclear antigen

DNA polymerase beta (pol beta) is a key player in DNA base excision repair (BER). Here, we describe the complex formation of pol beta with the protein arginine methyltransferase 1 (PRMT1). PRMT1 specifically methylated pol beta in vitro and in vivo. Arginine 137 was identified in pol beta as an important target for methylation by PRMT1. Neither the polymerase nor the dRP-lyase activities of pol beta were affected by PRMT1 methylation. However, this modification abolished the interaction of pol beta with proliferating cell nuclear antigen (PCNA). Together, our results provide evidence that PRMT1 methylation of pol beta might play a regulatory role in BER by preventing the involvement of pol beta in PCNA-dependent DNA metabolic events.

Increased flexibility enhances misincorporation: temperature effects on nucleotide incorporation opposite a bulky carcinogen-DNA adduct by a Y-family DNA polymerase

The Y-family DNA polymerase Dpo4, from the thermophilic crenarchaeon Sulfolobus solfataricus P2, offers a valuable opportunity to investigate the effect of conformational flexibility on the bypass of bulky lesions because of its ability to function efficiently at a wide range of temperatures. Combined molecular modeling and experimental kinetic studies have been carried out for 10S-(+)-trans-anti-[BP]-N2-dG ((+)-ta-[BP]G), a lesion derived from the covalent reaction of a benzo[a]pyrene metabolite with guanine in DNA, at 55 degrees C and results compared with an earlier study at 37 degrees C (Perlow-Poehnelt, R. A., Likhterov, I., Scicchitano, D. A., Geacintov, N. E., and Broyde, S. (2004) J. Biol. Chem. 279, 36951-36961). The experimental results show that there is more overall nucleotide insertion opposite (+)-ta-[BP]G due to particularly enhanced mismatch incorporation at 55 degrees C compared with 37 degrees C. The molecular dynamics simulations suggest that mismatched nucleotide insertion opposite (+)-ta-[BP]G is increased at 55 degrees C compared with 37 degrees C because the higher temperature shifts the preference of the damaged base from the anti to the syn conformation, with the carcinogen on the more open major groove side. The mismatched dNTP structures are less distorted when the damaged base is syn than when it is anti, at the higher temperature. However, with the normal partner dCTP, the anti conformation with close to Watson-Crick alignment remains more favorable. The molecular dynamics simulations are consistent with the kcat values for nucleotide incorporation opposite the lesion studied, providing structural interpretation of the experimental observations. The observed temperature effect suggests that conformational flexibility plays a role in nucleotide incorporation and bypass fidelity opposite (+)-ta-[BP]G by Dpo4.

Magnesium-cationic dummy atom molecules enhance representation of DNA polymerase beta in molecular dynamics simulations: improved accuracy in studies of structural features and mutational effects

Human DNA polymerase beta (pol beta) fills gaps in DNA as part of base excision DNA repair. Due to its small size it is a convenient model enzyme for other DNA polymerases. Its active site contains two Mg(2+) ions, of which one binds an incoming dNTP and one catalyzes its condensation with the DNA primer strand. Simulating such binuclear metalloenzymes accurately but computationally efficiently is a challenging task. Here, we present a magnesium-cationic dummy atom approach that can easily be implemented in molecular mechanical force fields such as the ENZYMIX or the AMBER force fields. All properties investigated here, namely, structure and energetics of both Michaelis complexes and transition state (TS) complexes were represented more accurately using the magnesium-cationic dummy atom model than using the traditional one-atom representation for Mg(2+) ions. The improved agreement between calculated free energies of binding of TS models to different pol beta variants and the experimentally determined activation free energies indicates that this model will be useful in studying mutational effects on catalytic efficiency and fidelity of DNA polymerases. The model should also have broad applicability to the modeling of other magnesium-containing proteins.

Structure of DNA polymerase beta with a benzo[c]phenanthrene diol epoxide-adducted template exhibits mutagenic features

We have determined the crystal structure of the human base excision repair enzyme DNA polymerase beta (Pol beta) in complex with a 1-nt gapped DNA substrate containing a template N2-guanine adduct of the tumorigenic (-)-benzo[c]phenanthrene 4R,3S-diol 2S,1R-epoxide in the gap. Nucleotide insertion opposite this adduct favors incorrect purine nucleotides over the correct dCMP and hence can be mutagenic. The structure reveals that the phenanthrene ring system is stacked with the base pair immediately 3' to the modified guanine, thereby occluding the normal binding site for the correct incoming nucleoside triphosphate. The modified guanine base is displaced downstream and prevents the polymerase from achieving the catalytically competent closed conformation. The incoming nucleotide binding pocket is distorted, and the adducted deoxyguanosine is in a syn conformation, exposing its Hoogsteen edge, which can hydrogen-bond with dATP or dGTP. In a reconstituted base excision repair system, repair of a deaminated cytosine (i.e., uracil) opposite the adducted guanine was dramatically decreased at the Pol beta insertion step, but not blocked. The efficiency of gap-filling dCMP insertion opposite the adduct was diminished by >6 orders of magnitude compared with an unadducted templating guanine. In contrast, significant misinsertion of purine nucleotides (but not dTMP) opposite the adducted guanine was observed. Pol beta also misinserts a purine nucleotide opposite the adduct with ungapped DNA and exhibits limited bypass DNA synthesis. These results indicate that Pol beta-dependent base excision repair of uracil opposite, or replication through, this bulky DNA adduct can be mutagenic.

Kinetic effect of a downstream strand and its 5'-terminal moieties on single nucleotide gap-filling synthesis catalyzed by human DNA polymerase lambda

During short-patch base excision repair, the excision of a 5'-terminal 2-deoxyribose-5-phosphate moiety of the downstream strand by the 5'-2-deoxyribose-5-phosphate lyase activity of either DNA polymerase beta or lambda is believed to occur after each respective enzyme catalyzes gap-filling DNA synthesis. Yet the effects of this 5'-terminal 2-deoxyribose-5-phosphate moiety on the polymerase activities of these two enzymes have never been quantitatively determined. Moreover, x-ray crystal structures of truncated polymerase lambda have revealed that the downstream strand and its 5'-phosphate group of gapped DNA interact intensely with the dRPase domain, but the kinetic effect of these interactions is unclear. Here, we utilized pre-steady state kinetic methods to systematically investigate the effect of a downstream strand and its 5'-moieties on the polymerase activity of the full-length human polymerase lambda. The downstream strand and its 5'-phosphate were both found to increase nucleotide incorporation efficiency (kp/Kd) by 15 and 11-fold, respectively, with the increase procured by the effect on the nucleotide incorporation rate constant kp rather than the ground state nucleotide binding affinity Kd. With 4 single nucleotide-gapped DNA substrates containing a 1,2-dideoxyribose-5-phosphate moiety, a 2-deoxyribose-5-phosphate mimic, we measured the incorporation efficiencies of 16 possible nucleotides. Our results demonstrate that although this 5'-terminal 2-deoxyribose-5-phosphate mimic does not affect the fidelity of polymerase lambda, it moderately decreased the polymerase efficiency by 3.4-fold. Moreover, this decrease in polymerase efficiency is due to a drop of similar magnitude in kp rather than Kd. The implication of the downstream strand and its 5'-moieties on the kinetics of gap-filling synthesis is discussed.

Site-directed mutagenesis in the fingers subdomain of HIV-1 reverse transcriptase reveals a specific role for the beta3-beta4 hairpin loop in dNTP selection

HIV-1 reverse transcriptase shares the key features of high fidelity polymerases, such as a closed architecture of the active site, but displays a level of fidelity that is intermediate to that of high fidelity, replicative polymerases and low fidelity translesion synthesis (TLS) polymerases. The beta3-beta4 loop of the HIV-1 RT fingers subdomain makes transient contacts with the dNTP and template base. To investigate the role of active site architecture in HIV-1 RT fidelity, we truncated the beta3-beta4 loop, eliminating contact between Lys65 and the gamma-phosphate of dNTP. The mutant, in a manner reminiscent of TLS polymerases, was only able to incorporate a nucleotide that was capable of base-pairing with the template nucleotide, but not a nucleotide shape-analog incapable of Watson-Crick hydrogen bonding. Unexpectedly, however, the deletion mutant differed from the TLS polymerases in that it displayed an increased fidelity. The increased fidelity was associated with reduced dNTP binding affinity as measured using the dead end complex formation. In an effort to delineate the specific amino acid residue in the deleted segment responsible for this phenotype, we examined the K65 residue. Two substitution mutants, K65R and K65A were studied. The K65A mutant behaved similarly to the deletion mutant displaying dependence on Watson-Crick hydrogen bonding, increased fidelity and reduced dNTP-binding, while the K65R was more akin to wild-type enzyme. These results underscore the key role of the K65 residue in the phenotype observed in the deletion mutant. Based on the well-known electrostatic interaction between K65 and the gamma-phosphate moiety of incoming dNTP substrate in the ternary complex structure of HIV-1 RT, we conclude that non-discriminatory interactions between beta3-beta4 loop and the dNTP in wild-type HIV-1 RT help lower dNTP selectivity. Our results show that the fidelity of dNTP insertion is influenced by protein interactions with the triphosphate moiety.

Correct and incorrect nucleotide incorporation pathways in DNA polymerase beta

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

Site-directed mutational analysis of structural interactions of low molecule compounds binding to the N-terminal 8 kDa domain of DNA polymerase beta

We previously reported the mode of inhibition of DNA polymerase beta (pol. beta) by long chain fatty acids and a bile acid, involving binding analyses to the N-terminal 8-kDa DNA binding domain. Here we describe a site-directed mutational analysis in which the key amino acids (L11, K35, H51, K60, L77, and T79), which are direct interaction sites in the domain, were substituted with K, A, A, A, K, and A, respectively. And their pol. beta interactions with a C24-long chain fatty acid, nervonic acid (NA), and a bile acid, lithocholic acid (LCA), were investigated by gel mobility shift assay and NMR spectroscopy. In the case of K35A, there was complete loss of DNA binding activity while K60A hardly has any activity. In contrast the other mutations had no appreciable effects. Thus, K35 and K60 are key amino acid sites for binding to template DNA. The DNA binding activities of L11K, H51A, and T79A as well as the wild type were inhibited by NA to the same extent. T79A demonstrated a disturbed interaction with LCA. 1H-15N HSQC NMR analysis indicated that despite their many similarities, the wild-type and the mutant proteins displayed some significant chemical shift differences. Not only were the substituted amino acid residues three-dimensionally shifted, but some amino acids which are positioned far distant from the key amino acids showed a shift. These results suggest that the interaction surface was significantly distorted with the result that LCA could not bind to the domain. These findings confirm our previous biochemical and 3D structural proposals concerning inhibition by NA and LCA.

Crystal structures of clinically relevant Lys103Asn/Tyr181Cys double mutant HIV-1 reverse transcriptase in complexes with ATP and non-nucleoside inhibitor HBY 097

Lys103Asn and Tyr181Cys are the two mutations frequently observed in patients exposed to various non-nucleoside reverse transcriptase inhibitor drugs (NNRTIs). Human immunodeficiency virus (HIV) strains containing both reverse transcriptase (RT) mutations are resistant to all of the approved NNRTI drugs. We have determined crystal structures of Lys103Asn/Tyr181Cys mutant HIV-1 RT with and without a bound non-nucleoside inhibitor (HBY 097, (S)-4-isopropoxycarbonyl-6-methoxy-3-(methylthio-methyl)-3,4-dihydroquinoxalin-2(1H)-thione) at 3.0 A and 2.5 A resolution, respectively. The structure of the double mutant RT/HBY 097 complex shows a rearrangement of the isopropoxycarbonyl group of HBY 097 compared to its binding with wild-type RT. HBY 097 makes a hydrogen bond with the thiol group of Cys181 that helps the drug retain potency against the Tyr181Cys mutation. The structure of the unliganded double mutant HIV-1 RT showed that Lys103Asn mutation facilitates coordination of a sodium ion with Lys101 O, Asn103 N and O(delta1), Tyr188 O(eta), and two water molecules. The formation of the binding pocket requires the removal of the sodium ion. Although the RT alone and the RT/HBY 097 complex were crystallized in the presence of ATP, only the RT has an ATP coordinated with two Mn(2+) at the polymerase active site. The metal coordination mimics a reaction intermediate state in which complete octahedral coordination was observed for both metal ions. Asp186 coordinates at an axial position whereas the carboxylates of Asp110 and Asp185 are in the planes of coordination of both metal ions. The structures provide evidence that NNRTIs restrict the flexibility of the YMDD loop and prevent the catalytic aspartate residues from adopting their metal-binding conformations.

Active site dynamics and combined quantum mechanics/molecular mechanics (QM/MM) modelling of a HIV-1 reverse transcriptase/DNA/dTTP complex

We have investigated the structure and dynamics of the HIV-1 reverse transcriptase (HIV-RT) active site, by modelling the active conformation of the HIV-1 RT/DNA/deoxythymidine triphosphate (dTTP) ternary complex. This has included molecular dynamics simulations with the CHARMM27 force field, and combined quantum mechanics/molecular mechanics (QM/MM) calculations. Three different ternary systems were studied to investigate the effects of different protonation states of the dTTP substrate (a deprotonated and two different mono-protonated triphosphate forms of dTTP at the active site), and the effects of different possible protonation state of potentially catalytic aspartate residues (Asp185 and Asp186) were tested. Several potentially important hydrogen-bonding interactions with amino acids and bound water molecules in the deoxyribonucleoside triphosphate (dNTP) binding pocket were examined. The model of the deprotonated form of the dTTP substrate with the two aspartates in their charged (basic) form seemed to be the most stable and its orientation was in good agreement with X-ray crystallographic structure. In addition, two different semiempirical (AM1/CHARMM and PM3/CHARMM) QM/MM methods were tested for the HIV-RT system, in structural optimizations. Both methods provided conformations of the triphosphate moiety in either fully deprotonated or mono-protonated forms, which agreed well with the experimental structure of dTTP. The only significant difference between the AM1/CHARMM and PM3/CHARMM minimized structures is that the PM3/CHARMM Palpha-O3' optimized distance (important for nucleotide addition) is longer by 0.66 A in the deprotonated system but shorter by 0.37 A in the mono-protonated triphosphate system as compared with those obtained from AM1/CHARMM minimized structure. The obtained results suggest that both of these QM/MM methods, and the stochastic boundary molecular dynamics approach applied in this work, can give reasonable results for modelling the catalytically active complex of this important enzyme. The results provide insight into the structure and interactions of the active site of this important enzyme, with implications for its mechanism, which may be useful in inhibitor design.

A specific loop in human DNA polymerase mu allows switching between creative and DNA-instructed synthesis

Human DNA polymerase mu (Polμ) is a family X member that has terminal transferase activity but, in spite of a non-orthodox selection of the template information, displays its maximal catalytic efficiency in DNA-templated reactions. As terminal deoxynucleotidyl transferase (TdT), Polμ has a specific loop (loop1) that could provide this enzyme with its terminal transferase activity. When loop1 was deleted, human Polμ lacked TdT activity but improved DNA-binding and DNA template-dependent polymerization. Interestingly, when loop1 from TdT was inserted in Polμ (substituting its cognate loop1), the resulting chimaera displayed TdT activity, preferentially inserting dGTP residues, but had a strongly reduced template-dependent polymerization activity. Therefore, a specialized loop in Polμ, that could adopt alternative conformations, appears to provide this enzyme with a dual capacity: (i) template independency to create new DNA information, in which loop1 would have an active role by acting as a ‘pseudotemplate’; (ii) template-dependent polymerization, in which loop1 must allow binding of the template strand. Recent in vivo and in vitro data suggest that such a dual capacity could be advantageous to resolve microhomology-mediated end-joining reactions.

The Structure of T. aquaticus DNA Polymerase III Is Distinct from Eukaryotic Replicative DNA Polymerases

The crystal structure of Thermus aquaticus DNA polymerase III alpha subunit reveals that the structure of the catalytic domain of the eubacterial replicative polymerase is unrelated to that of the eukaryotic replicative polymerase but rather belongs to the Polbeta-like nucleotidyltransferase superfamily. A model of the polymerase complexed with both DNA and beta-sliding clamp interacting with a reoriented binding domain and internal beta binding site was constructed that is consistent with existing biochemical data. Within the crystal, two C-terminal domains are interacting through a surface that is larger than many dimer interfaces. Since replicative polymerases of eubacteria and eukaryotes/archaea are not homologous, the nature of the replicative polymerase in the last common ancestor is unknown. Although other possibilities have been proposed, the plausibility of a ribozyme DNA polymerase should be considered.

Crystal Structure of the Catalytic α Subunit of E. coli Replicative DNA Polymerase III

Bacterial replicative DNA polymerases such as Polymerase III (Pol III) share no sequence similarity with other polymerases. The crystal structure, determined at 2.3 A resolution, of a large fragment of Pol III (residues 1-917), reveals a unique chain fold with localized similarity in the catalytic domain to DNA polymerase beta and related nucleotidyltransferases. The structure of Pol III is strikingly different from those of members of the canonical DNA polymerase families, which include eukaryotic replicative polymerases, suggesting that the DNA replication machinery in bacteria arose independently. A structural element near the active site in Pol III that is not present in nucleotidyltransferases but which resembles an element at the active sites of some canonical DNA polymerases suggests that, at a more distant level, all DNA polymerases may share a common ancestor. The structure also suggests a model for interaction of Pol III with the sliding clamp and DNA.

Energy analysis of chemistry for correct insertion by DNA polymerase beta

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

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

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

A new paradigm for DNA polymerase specificity

We show that T7 DNA polymerase exists in three distinct structural states, as reported by a conformationally sensitive fluorophore attached to the recognition (fingers) domain. The conformational change induced by a correct nucleotide commits the substrate to the forward reaction, and the slow reversal of the conformational change eliminates the rate of the chemistry step from any contribution toward enzyme specificity. Discrimination against mismatches is enhanced by the rapid release of mismatched nucleotides from the ternary E.DNA.deoxynucleoside triphosphate complex and by the use of substrate-binding energy to actively misalign catalytic residues to reduce the rate of misincorporation. Our refined model for enzyme selectivity extends traditional thermodynamic formalism by including substrate-induced structural alignment or misalignment of catalytic residues as a third dimension on the free-energy profile and by including the rate of substrate dissociation as a key kinetic parameter.

Molecular flexibility in protein–DNA interactions

In living cells protein-DNA interactions are fundamental processes. Here, we compare the 3D structures of several DNA-binding proteins frequently determined with and without attached DNA. We studied the global structure (backbone-traces) as well as the local structure (binding sites) by comparing pair-wise the related atoms. The DNA-interaction sites of uncomplexed proteins show conspicuously high local structural flexibility. Binding to DNA results in specific local conformations, which are clearly distinct from the unbound states. The adaptation of the protein's binding site to DNA can never be described by the lock and key model but in all cases by the induced fit model. Conformational changes in the seven protein backbone traces take place in different ways. Two of them dock onto DNA without a significant change, while the other five proteins are characterized by a backbone conformation change caused by DNA docking. In the case of three proteins of the latter group the DNA-complexed conformation also occurs in a few uncomplexed structures. This behavior can be described by a conformational ensemble, which is narrowed down by DNA docking until only one single DNA-complexed conformation occurs. Different docking models are discussed and each of the seven proteins is assigned to one of them.

The properties of steric gate mutants reveal different constraints within the active sites of Y-family and A-family DNA polymerases

Y-family (lesion-bypass) DNA polymerases show the same overall structural features seen in other members of the polymerase superfamily, yet their active sites are more open, with fewer contacts to the DNA and nucleotide substrates. This raises the question of whether analogous active-site side chains play equivalent roles in the bypass polymerases and their classical DNA polymerase counterparts. In Klenow fragment, an A-family DNA polymerase, the steric gate side chain (Glu710) not only prevents ribonucleotide incorporation but also plays an important role in discrimination against purine-pyrimidine mispairs. In this work we show that the steric gate (Phe12) of the Y-family polymerase Dbh plays a very minor role in fidelity, despite its analogous role in sugar selection. Using ribonucleotide discrimination to report on the positioning of a mispaired dNTP, we found that the pyrimidine of a Pu-dPyTP nascent mispair occupies a similar position to that of a correctly paired dNTP in the Dbh active site, whereas in Klenow fragment the mispaired dNTP sits higher in the active site pocket. If purine-pyrimidine mispairs adopt the expected wobble geometry, the difference between the two polymerases can be attributed to the binding of the templating base, with the looser binding site of Dbh permitting a variety of template conformations with only minimal adjustment at the incoming dNTP. In Klenow fragment the templating base is more rigidly held, so that changes in base pair geometry would affect the dNTP position, allowing the Glu710 side chain to serve as a sensor of nascent mispairs.

Magnesium-induced assembly of a complete DNA polymerase catalytic complex

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

Crystal structures of the vaccinia virus polyadenylate polymerase heterodimer: insights into ATP selectivity and processivity

Polyadenylation of mRNAs in poxviruses, crucial for virion maturation, is carried out by a poly(A) polymerase heterodimer composed of a catalytic component, VP55, and a processivity factor, VP39. The ATP-gamma-S bound and unbound crystal structures of the vaccinia polymerase reveal an unusual architecture for VP55 that comprises of N-terminal, central or catalytic, and C-terminal domains with different topologies and that differs from many polymerases, including the eukaryotic poly(A) polymerases. Residues in the active site of VP55, located between the catalytic and C-terminal domains, make specific interactions with the adenine of the ATP analog, establishing the molecular basis of ATP recognition. VP55's concave surface docks the globular VP39. A model for RNA primer binding that involves all three VP55 domains and VP39 is proposed. The model supports biochemical evidence that VP39 functions as a processivity factor by partially enclosing the RNA primer at the heterodimer interface.

RNA aptamers selected against DNA polymerase beta inhibit the polymerase activities of DNA polymerases beta and kappa

DNA polymerase β (polβ), a member of the X family of DNA polymerases, is the major polymerase in the base excision repair pathway. Using in vitro selection, we obtained RNA aptamers for polβ from a variable pool of 8 × 10¹² individual RNA sequences containing 30 random nucleotides. A total of 60 individual clones selected after seven rounds were screened for the ability to inhibit polβ activity. All of the inhibitory aptamers analyzed have a predicted tri-lobed structure. Gel mobility shift assays demonstrate that the aptamers can displace the DNA substrate from the polβ active site. Inhibition by the aptamers is not polymerase specific; inhibitors of polβ also inhibited DNA polymerase κ, a Y-family DNA polymerase. However, the RNA aptamers did not inhibit the Klenow fragment of DNA polymerase I and only had a minor effect on RB69 DNA polymerase activity. Polβ and κ, despite sharing little sequence similarity and belonging to different DNA polymerase families, have similarly open active sites and relatively few interactions with their DNA substrates. This may allow the aptamers to bind and inhibit polymerase activity. RNA aptamers with inhibitory properties may be useful in modulating DNA polymerase actvity in cells.

Up-regulation of the fidelity of human DNA polymerase lambda by its non-enzymatic proline-rich domain

DNA repair pathways are essential for maintaining genome stability. DNA polymerase beta plays a critical role in base-excision repair in vivo. DNA polymerase lambda, a recently identified X-family homolog of DNA polymerase beta, is hypothesized to be a second polymerase involved in base-excision repair. The full-length DNA polymerase lambda is comprised of three domains: a C-terminal DNA polymerase beta-like domain, an N-terminal BRCA1 C-terminal domain, and a previously uncharacterized proline-rich domain. Strikingly, pre-steady-state kinetic analyses reveal that, although human DNA polymerase lambda has almost identical fidelity to human DNA polymerase beta, the C-terminal DNA polymerase beta-like domain alone displays a dramatic, up to 100-fold loss in fidelity. We further demonstrate that the non-enzymatic proline-rich domain confers the increase in fidelity of DNA polymerase lambda by significantly lowering incorporation rate constants of incorrect nucleotides. Our studies illustrate a novel mechanism, in which the DNA polymerase fidelity is controlled not by an accessory protein or a proofreading exonuclease domain but by an internal regulatory domain.