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.     

Validating the vitality strategy for fighting drug resistance

The current challenge in designing effective drugs against HIV-1 is to find novel candidates with high potency, but with a lower susceptibility to mutations associated with drug resistance. Trying to address this challenge, we developed in our previous study (Ishikita and Warshel, Angew Chem Int Ed Engl 2008; 47:697-700) a novel computational strategy for fighting drug resistance by predicting the likely moves of the virus through constraints on binding and catalysis. This has been based on calculating the ratio between the vitality values ((K(i) k(cat)/K(M))(mutant)/(K(i) k(cat)/K(M))(wild-type)) and using it as a guide for predicting the moves of the virus. The corresponding calculations of the binding affinity, K(i), were carried out using the semi-macroscopic version of the protein dipole Langevin dipole (PDLD/S) in its linear response approximation (LRA) in its β version (PDLD/S-LRA/β). We also calculate the proteolytic efficiency, k(cat)/K(M), by evaluating the transition state (TS) binding free energies using the PDLD/S-LRA/β method. Here we provide an extensive validation of our strategy by calculating the vitality of six existing clinical and experimental drug candidates. It is found that the computationally determined vitalities correlate reasonably well with those derived from the corresponding experimental data. This indicates that the calculated vitality may be used to identify mutations that would be most effective for the survival of the virus. Thus, it should be possible to use our approach in screening for mutations that would provide the most effective resistance to any proposed antiviral drug. This ability should be very useful in guiding the design of drug molecules that will lead to the slowest resistance.     

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

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

phi 29 DNA polymerase residue Leu384, highly conserved in motif B of eukaryotic type DNA replicases,is involved in nucleotide insertion fidelity

Replicative DNA polymerases achieve insertion fidelity by geometric selection of a complementary nucleotide followed by induced fit: movement of the fingers subdomain toward the active site to enclose the incoming and templating nucleotides generating a binding pocket for the nascent base pair. Several residues of motif B of DNA polymerases from families A and B, localized in the fingers subdomain, have been described to be involved in template/primer binding and dNTP selection. Here we complete the analysis of this motif, which has the consensus "KLX2NSXYG" in DNA polymerases from family B, characterized by mutational analysis of conserved leucine, Leu384 of phi 29 DNA polymerase. Mutation of Leu384 into Arg resulted in a phi 29 DNA polymerase with reduced nucleotide insertion fidelity during DNA-primed polymerization and protein-primed initiation reactions. However, the mutation did not alter the intrinsic affinity for the different dNTPs, as shown in the template-independent terminal protein-deoxynucleotidylation reaction. We conclude that Leu384 of phi 29 DNA polymerase plays an important role in positioning the templating nucleotide at the polymerization active site and in controlling nucleotide insertion fidelity. This agrees with the localization of the corresponding residue in the closed ternary complexes of family A and family B DNA polymerases, contributing to form the binding pocket for the nascent base pair. As an additional effect, mutant polymerase L384R was strongly reduced in DNA binding, resulting in reduced processivity during polymerization.     

Identification of the magnesium ion binding site in the catalytic center of Escherichia coli primase by iron cleavage

Magnesium is essential for the catalysis reaction of Escherichia coli primase, the enzyme synthesizing primer RNA chains for initiation of DNA replication. To map the Mg(2+) binding site in the catalytic center of primase, we have employed the iron cleavage method in which the native bound Mg(2+) ions were replaced with Fe(2+) ions and the protein was then cleaved in the vicinity of the metal binding site by adding DTT which generated free hydroxyl radicals from the bound iron. Three Fe(2+) cleavages were generated at sites designated I, II, and III. Adding Mg(2+) or Mn(2+) ions to the reaction strongly inhibited Fe(2+) cleavage; however, adding Ca(2+) or Ba(2+) ions had much less effect. Mapping by chemical cleavage and subsequent site-directed mutagensis demonstrated that three acidic residues, Asp345 and Asp347 of a conserved DPD sequence and Asp269 of a conserved EGYMD sequence, were the amino acid residues that chelated Mg(2+) ions in the catalytic center of primase. Cleavage data suggested that binding to D345 is significantly stronger than to D347 and somewhat stronger than to D269.     

Effect of metal ion binding on the structural stability of the hepatitis C virus RNA polymerase

The RNA polymerase activity of the hepatitis C virus, a major human pathogen, has previously been shown to be supported by metal ions. In the present study, we report a systematic analysis of the effect of metal ion binding on the structural stability of the hepatitis C virus RNA polymerase. Chemical and thermal denaturation assays revealed that the stability of the protein is increased significantly in the presence of metal ions. Structural analyses clearly established that metal ion binding increases hydrophobic exposure on the RNA polymerase surface. Furthermore, our denaturation studies, coupled with polymerization assays, demonstrate that the active site region of the polymerase is more sensitive to chemical denaturant than other structural scaffolds. We also report the first detailed study of the thermodynamic parameters involved in the interaction between the hepatitis C virus RNA polymerase and metal ions. Finally, a mutational analysis was also performed to investigate the importance of Asp(220), Asp(318), and Asp(319) for metal ion binding. This mutational study underscores a strict requirement for each of the residues for metal binding, indicating that the active center of the HCV RNA polymerase is intolerant to virtually any perturbations of the metal coordination sphere, thereby highlighting the critical role of the enzyme-bound metal ions. Overall, our results indicate that metal ions play a dual modulatory role in the RNA polymerase reaction by promoting both a favorable geometry of the active site for catalysis and by increasing the structural stability of the enzyme.     

The roles of Tyr391 and Tyr619 in RB69 DNA polymerase replication fidelity

In the family-B DNA polymerase of bacteriophage RB69, the conserved aromatic palm-subdomain residues Tyr391 and Tyr619 interact with the last primer-template base-pair. Tyr619 interacts via a water-mediated hydrogen bond with the phosphate of the terminal primer nucleotide. The main-chain amide of Tyr391 interacts with the corresponding template nucleotide. A hydrogen bond has been postulated between Tyr391 and the hydroxyl group of Tyr567, a residue that plays a key role in base discrimination. This hydrogen bond may be crucial for forcing an infrequent Tyr567 rotamer conformation and, when the bond is removed, may influence fidelity. We investigated the roles of these residues in replication fidelity in vivo employing phage T4 rII reversion assays and an rI forward assay. Tyr391 was replaced by Phe, Met and Ala, and Tyr619 by Phe. The Y391A mutant, reported previously to decrease polymerase affinity for incoming nucleotides, was unable to support DNA replication in vivo, so we used an in vitro fidelity assay. Tyr391F/M replacements affect fidelity only slightly, implying that the bond with Tyr567 is not essential for fidelity. The Y391A enzyme has no mutator phenotype in vitro. The Y619F mutant displays a complex profile of impacts on fidelity but has almost the same mutational spectrum as the parental enzyme. The Y619F mutant displays reduced DNA binding, processivity, and exonuclease activity on single-stranded DNA and double-stranded DNA substrates. The Y619F substitution would disrupt the hydrogen bond network at the primer terminus and may affect the alignment of the 3' primer terminus at the polymerase active site, slowing chemistry and overall DNA synthesis.     

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.     

Identification of residues critical for the polymerase activity of the Klenow fragment of DNA polymerase I from Escherichia coli

The Klenow fragment structure, together with many biochemical experiments, has suggested a region of the protein that may contain the polymerase active site. We have changed 7 amino acid residues within this region by site-directed mutagenesis, yielding 12 mutant proteins which have been purified and analyzed in vitro. The results of steady-state kinetic determinations of Km(dNTP) and kcat for the polymerase reaction, together with measurements of DNA binding affinity, suggest strongly that this study has succeeded in targeting important active site residues. Moreover, the in vitro data allow dissection of the proposed active site region into two clusters of residues that are spatially, as well as functionally, fairly distinct. Mutations in Tyr766, Arg841, and Asn845 cause an increase in Km(dNTP), suggesting that contacts with the incoming dNTP are made in this region. Mutations in the second cluster of residues, Gln849, Arg668, and Asp882, cause a large decrease in kcat, suggesting a role for these residues in catalysis of the polymerase reaction. The DNA-binding properties of mutations at positions 849 and 668 may indicate that the catalytic role of these side chains is associated with their interaction with the DNA substrate. Screening of the mutations in vivo for the classical polA-defective phenotype (sensitivity to DNA damage) demonstrated that a genetic screen of this type may be a reasonable predictor or kcat or of DNA binding affinity in future mutational studies.     

Metal-induced DNA translocation leads to DNA polymerase conformational activation

Binding of the catalytic divalent ion to the ternary DNA polymerase β/gapped DNA/dNTP complex is thought to represent the final step in the assembly of the catalytic complex and is consequently a critical determinant of replicative fidelity. We have analyzed the effects of Mg(2+) and Zn(2+) on the conformational activation process based on NMR measurements of [methyl-(13)C]methionine DNA polymerase β. Unexpectedly, both divalent metals were able to produce a template base-dependent conformational activation of the polymerase/1-nt gapped DNA complex in the absence of a complementary incoming nucleotide, albeit with different temperature thresholds. This conformational activation is abolished by substituting Glu295 with lysine, thereby interrupting key hydrogen bonds necessary to stabilize the closed conformation. These and other results indicate that metal-binding can promote: translocation of the primer terminus base pair into the active site; expulsion of an unpaired pyrimidine, but not purine, base from the template-binding pocket; and motions of polymerase subdomains that close the active site. We also have performed pyrophosphorolysis studies that are consistent with predictions based on these results. These findings provide new insight into the relationships between conformational activation, enzyme activity and polymerase fidelity.     

Insight into the catalytic mechanism of DNA polymerase beta: structures of intermediate complexes

The catalytic reaction mediated by DNA polymerases is known to require two Mg(II) ions, one associated with dNTP binding and the other involved in metal ion catalysis of the chemical step. Here we report a functional intermediate structure of a DNA polymerase with only one metal ion bound, the DNA polymerase beta-DNA template-primer-chromium(III).2'-deoxythymidine 5'-beta,gamma-methylenetriphosphate [Cr(III).dTMPPCP] complex, at 2.6 A resolution. The complex is distinct from the structures of other polymerase-DNA-ddNTP complexes in that the 3'-terminus of the primer has a free hydroxyl group. Hence, this structure represents a fully functional intermediate state. Support for this contention is provided by the observation of turnover in biochemical assays of crystallized protein as well as from the determination that soaking Pol beta crystals with Mn(II) ions leads to formation of the product complex, Pol beta-DNA-Cr(III).PCP, whose structure is also reported. An important feature of both structures is that the fingers subdomain is closed, similar to structures of other ternary complexes in which both metal ion sites are occupied. These results suggest that closing of the fingers subdomain is induced specifically by binding of the metal-dNTP complex prior to binding of the catalytic Mg(2+) ion. This has led us to reevaluate our previous evidence regarding the existence of a rate-limiting conformational change in Pol beta's reaction pathway. The results of stopped-flow studies suggest that there is no detectable rate-limiting conformational change step.     

Y265H mutator mutant of DNA polymerase beta. Proper teometric alignment is critical for fidelity

DNA polymerases have the unique ability to select a specific deoxynucleoside triphosphate from a pool of similarly structured substrates. One of these enzymes, DNA polymerase beta, offers a simple system to relate polymerase structure to the fidelity of DNA synthesis. In this study, a mutator DNA polymerase beta, Y265H, was identified using an in vivo genetic screen. Purified Y265H produced errors at a 40-fold higher frequency than the wild-type protein in a forward mutation assay. At 37 degrees C, transient kinetic analysis demonstrated that the alteration caused a 111-fold decrease in the maximum rate of polymerization and a 117-fold loss in fidelity for G misincorporation opposite template A. Our data suggest that the maximum rate of polymerization was reduced, because Y265H was dramatically impaired in its ability to perform nucleotidyl transfer in the presence of the correct nucleotide substrate. In contrast, at 20 degrees C, the mutant protein had a fidelity similar to wild-type enzyme. Both proteins at 20 degrees C demonstrate a rapid change in protein conformation, followed by a slow chemical step. These data suggest that proper geometric alignment of template, 3'-OH of the primer, magnesium ions, dNTP substrates, and the active site residues of DNA polymerase beta are important factors in polymerase fidelity and provide the first evidence that Tyr-265 is important for this alignment to occur properly in DNA polymerase beta.     

Interactions of the 8-kDa domain of rat DNA polymerase beta with DNA

Interactions between the isolated 8-kDa domain of the rat DNA polymerase beta and DNA have been studied, using the quantitative fluorescence titration technique. The obtained results show that the number of nucleotide residues occluded in the native 8-kDa domain complex with the ssDNA (the site size) is strongly affected by Mg2+ cations. In the absence of Mg2+, the domain occludes 13 +/- 0.7 nucleotide residues, while in the presence of Mg2+ the site size decreases to 9 +/- 0.6 nucleotides. The high affinity of the magnesium cation binding, as well as the dramatic changes in the monovalent salt effect on the protein-ssDNA interactions in the presence of Mg2+, indicates that the site size decrease results from the Mg2+ binding to the domain. The site size of the isolated domain-ssDNA complex is significantly larger than the 5 +/- 2 site size determined for the (pol beta)5 binding mode formed by an intact polymerase, indicating that the intact enzyme, but not the isolated domain, has the ability to use only part of the domain DNA-binding site in its interactions with the nucleic acid. Salt effect on the intrinsic interactions of the domain with the ssDNA indicates that a net release of m approximately 5 ions accompanies the complex formation. Independence of the number of ions released upon the type of anion in solution strongly suggests that the domain forms as many as seven ionic contacts with the ssDNA. Experiments with different ssDNA oligomers show that the affinity decreases gradually with the decreasing number of nucleotide residues in the oligomer. The data indicate a continuous, energetically homogeneous structure of the DNA-binding site of the domain, with crucial, nonspecific contacts between the protein and the DNA evenly distributed over the entire binding site. The DNA-binding site shows little base specificity. Moreover, the domain has an intrinsic affinity and site size of its complex with the dsDNA conformation, similar to the affinity and site size with the ssDNA. The significance of these results for the mechanistic role of the 8-kDa domain in the functioning of rat pol beta is discussed.     

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

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

Subunit protein-affinity isolation of Drosophila DNA polymerase catalytic subunit

gfLittle is known at present about the biochemical properties of very large-sized Drosophila DNA polymerases. In a previous study, we tried to purify Drosophila pol. catalytic subunit from embryos through seven column chromatographies and study its biochemical properties. However, we failed to characterize it precisely because an insufficient amount of the enzyme was generated. In this report, we describe direct purification from Drosophila embryos to near homogeneity using Drosophila DNA polymerase second subunit (Drosophila pol. 2) protein-conjugated affinity column chromatography and characterization of the enzyme in detail. To our knowledge this is the first demonstration of native DNA polymerase purification with activity using a subunit protein-affinity column. We observed new characteristics of Drosophila pol. catalytic subunit as follows: Drosophila pol. catalytic subunit synthesized DNA processively in the presence of both Mn(2+) and Mg(2+) ions, but Mn(2+) inhibited the 3'-5' proofreading activity, thereby decreasing the fidelity of DNA replication by 50%.     

Structural insights into single-stranded DNA binding and cleavage by F factor TraI

Conjugative plasmid transfer between bacteria disseminates antibiotic resistance and diversifies prokaryotic genomes. Relaxases, proteins essential for conjugation, cleave one plasmid strand sequence specifically prior to transfer. Cleavage occurs through a Mg(2+)-dependent transesterification involving a tyrosyl hydroxyl and a DNA phosphate. The structure of the F plasmid TraI relaxase domain, described here, is a five-strand beta sheet flanked by alpha helices. The protein resembles replication initiator protein AAV-5 Rep but is circularly permuted, yielding a different topology. The beta sheet forms a binding cleft lined with neutral, nonaromatic residues, unlike most single-stranded DNA binding proteins which use aromatic and charged residues. The cleft contains depressions, suggesting base recognition occurs in a knob-into-hole fashion. Unlike most nucleases, three histidines but no acidic residues coordinate a Mg(2+) located near the catalytic tyrosine. The full positive charge on the Mg(2+) and the architecture of the active site suggest multiple roles for Mg(2+) in DNA cleavage.