On molecular mechanisms of nucleic acid synthesis fidelity aspects. 1. Contribution of base interactions

This paper is concerned with the molecular mechanisms of spontaneous replacements of base pairs in the processes of template synthesis of nucleic acids. The method of atom-atom potential functions was used to calculate the energies of interaction in non-complementary base pairs arranged in a common plane so that the mutual position of the glycosidic bonds does not differ appreciably from their position for Watson-Crick pairs in the DNA double helix. A number of local minima of this energy have been found, which could occur in template synthesis and result in insertion of incorrect bases into the double helix. The calculation results are indicative of formation of purine-purine pairs with one of the purine nucleotides in syn-conformation, which can be regarded as a typical pathway of transversion, and that of wobble-pairs TG and AC, which can be regarded as a typical pathway of transition. The contribution of intramolecular interactions of nucleic acids as well as interactions of polynucleotide chains with an enzyme to the fidelity of template synthesis of nucleic acids is discussed. The calculation results are compared with the experimental data on the frequency of spontaneous mutations and the frequency of errors in template synthesis of nucleic acids in vitro.     

Isostericity and tautomerism of base pairs in nucleic acids

The natural bases of nucleic acids form a great variety of base pairs with at least two hydrogen bonds between them. They are classified in twelve main families, with the Watson–Crick family being one of them. In a given family, some of the base pairs are isosteric between them, meaning that the positions and the distances between the C1′ carbon atoms are very similar. The isostericity of Watson–Crick pairs between the complementary bases forms the basis of RNA helices and of the resulting RNA secondary structure. Several defined suites of non-Watson–Crick base pairs assemble into RNA modules that form recurrent, rather regular, building blocks of the tertiary architecture of folded RNAs. RNA modules are intrinsic to RNA architecture are therefore disconnected from a biological function specifically attached to a RNA sequence. RNA modules occur in all kingdoms of life and in structured RNAs with diverse functions. Because of chemical and geometrical constraints, isostericity between non-Watson–Crick pairs is restricted and this leads to higher sequence conservation in RNA modules with, consequently, greater difficulties in extracting 3D information from sequence analysis. Nucleic acid helices have to be recognised in several biological processes like replication or translational decoding. In polymerases and the ribosomal decoding site, the recognition occurs on the minor groove sides of the helical fragments. With the use of alternative conformations, protonated or tautomeric forms of the bases, some base pairs with Watson–Crick-like geometries can form and be stabilized. Several of these pairs with Watson–Crick-like geometries extend the concept of isostericity beyond the number of isosteric pairs formed between complementary bases. These observations set therefore limits and constraints to geometric selection in molecular recognition of complementary Watson–Crick pairs for fidelity in replication and translation processes.     

Interaction of Drosophila DNA polymerase alpha holoenzyme with synthetic template-primers containing mismatched primer bases or propanodeoxyguanosine adducts at various positions in template and primer regions.

We studied recognition and binding of synthetic template-primers by Drosophila DNA polymerase alpha (pol alpha) holoenzyme. The template-primers used contained either mismatched base pairs at various positions in the primer region or exocyclic propanodeoxyguanosine (PdG) adducts at various positions in both template and primer.pol alpha requires primer-terminal complementarity of greater than or equal to 4 base pairs for efficient binding and incorporation. When a mismatched base pair is at the -4 position relative to the 3'-primer terminus, minimal but detectable binding occurs. This is consistent with the ability of pol alpha to incorporate a single nucleotide on a template-primer containing a mismatch at this position, but at a rate of only 7% relative to incorporation on a perfectly matched template-primer. No binding or incorporation (less than 1% of incorporation on a perfectly matched template-primer) was evident when a mismatched base pair was at the -3 position or closer, relative to the 3'-primer terminus. Similar results were obtained when PdG was placed at various positions in the primer region. When a PdG residue was located in the template region (+ 3 position relative to the 3'-primer terminus), single-nucleotide incorporation was stimulated 3-4-fold. These observations suggest that there are intrinsic aspects to the mechanism of nucleotide incorporation by pol alpha which ensure the fidelity of DNA synthesis by this enzyme and may provide novel insights into the fundamental mechanism of polymerase translocation along templates.     

The fidelity of base selection by the polymerase subunit of DNA polymerase III holoenzyme.

In common with other DNA polymerases, DNA polymerase III holoenzyme of E. coli selects the biologically correct base pair with remarkable accuracy. DNA polymerase III is particularly useful for mechanistic studies because the polymerase and editing activities reside on separate subunits. To investigate the biochemical mechanism for base insertion fidelity, we have used a gel electrophoresis assay to measure kinetic parameters for the incorporation of correct and incorrect nucleotides by the polymerase (alpha) subunit of DNA polymerase III. As judged by this assay, base selection contributes a factor of roughly 10(4)-10(5) to the overall fidelity of genome duplication. The accuracy of base selection is determined mainly by the differential KM of the enzyme for correct vs. incorrect deoxynucleoside triphosphate. The misinsertion of G opposite template A is relatively efficient, comparable to that found for G opposite T. Based on a variety of other work, the G:A pair may require a special correction mechanism, possibly because of a syn-anti pairing approximating Watson-Crick geometry. We suggest that precise recognition of the equivalent geometry of the Watson-Crick base pairs may be the most critical feature for base selection.     

Experimental and theoretical study of energetics of complex formation between nucleic acid bases and the bases with amide group.

A number of nucleic acid base pairs and complexes between the bases and the amide group of acrylamide have been studied experimentally by using mass spectrometry and theoretically by the method of atom-atom potential function calculations. It has been found from temperature dependencies of peak intensities in mass spectra of m2.2.9(3) Gua.m1Ura, m9 Ade.m1Cyt, m2.2.9(3) Gua.m1Gua.m1Cyt pairs that enthalpy values, delta H, of the complex formation are equal to 14.2 +/- 1.1, 13.5 +/- 1.3 and 16.4 +/- 1.4 kcal/M, respectively, and those of acrylamide with m1.3(2) Ura and m1Thy corresponds to 9.7 +/- 1.0 and 6.8 +/- 0.6 kcal/M. There is a good agreement of the experimental data with calculations when taking into account both the amino-oxo and the amino-hydroxy tautomeric forms of guanine. A combined use of the data allows us to determine the energy, the modes of interaction and the structure of the complexes. The results are discussed in connection with the modelling of molecular structure of biopolymers by the method of classical potential functions, protein-nucleic acids recognition and fidelity of nucleic acids biosynthesis.     

Inefficient proofreading and biased error rates during inaccurate DNA synthesis by a mutant derivative of Saccharomyces cerevisiae DNA polymerase delta

DNA polymerase delta (pol delta) is a high fidelity eukaryotic enzyme that participates in DNA repair and is essential for DNA replication. Toward the goal of dissecting its multiple biological functions, here we describe the biochemical properties of Saccharomyces cerevisiae pol delta with a methionine replacing conserved leucine 612 at the polymerase active site. Compared with wild type pol delta, L612M pol delta has normal processivity and slightly higher polymerase specific activity. L612M pol delta also has normal 3' exonuclease activity, yet it is impaired in partitioning mismatches to the exonuclease active site, thereby reducing DNA synthesis fidelity. Error rates in vitro for L612M pol delta are elevated for both base substitutions and single base deletions but in a highly biased manner. For each of the six possible pairs of reciprocal mismatches that could arise during replication of complementary DNA strands to account for any particular base substitution in vivo (e.g. T-dGMP or A-dCMP for T to C transitions), L612M pol delta error rates are substantially higher for one mismatch than the other. These results provide a biochemical explanation for our observation, which confirms earlier genetic studies, that a haploid pol3-L612M S. cerevisiae strain has an elevated spontaneous mutation rate that is likely due to reduced replication fidelity in vivo.     

DNA polymerase insertion fidelity. Gel assay for site-specific kinetics

A quantitative assay based on gel electrophoresis is described to measure nucleotide insertion kinetics at an arbitrary DNA template site. The assay is used to investigate kinetic mechanisms governing the fidelity of DNA synthesis using highly purified Drosophila DNA polymerase alpha holoenzyme complex and M13 primer-template DNA. Km and Vmax values are reported for correct insertion of A and misinsertion of G, C, and T opposite a single template T site. The misinsertion frequencies are 2 X 10(-4) for G-T and 5 X 10(-5) for both C-T and T-T relative to normal A-T base pairs. The dissociation constant of the polymerase-DNA-dNTP complex, as measured by Km, plays a dominant role in determining the rates of forming right and wrong base pairs. Compared with Km for insertion of A opposite T (3.7 +/- 0.7 microM), the Km value is 1100-fold greater for misinsertion of G opposite T (4.2 +/- 0.4 mM), and 2600-fold greater for misinsertion of either C or T opposite T (9.8 +/- 4.2 mM). These Km differences indicate that in the enzyme binding site the stability of A-T base pairs is 4.3 kcal/mol greater than G-T pairs and 4.9 kcal/mol greater than C-T or T-T pairs. In contrast to the large differences in Km, differences in Vmax are relatively small. There is only a 4-fold reduction in Vmax for insertion of G opposite T and an 8-fold reduction for C or T opposite T, compared with the correct insertion of A. For the specific template T site investigated, the nucleotide insertion fidelity for Drosophila polymerase alpha seems to be governed primarily by a Km discrimination mechanism.     

Structure of human tryptophanyl-tRNA synthetase in complex with tRNATrp reveals the molecular basis of tRNA recognition and specificity

Aminoacyl-tRNA synthetases (aaRSs) are a family of enzymes responsible for the covalent link of amino acids to their cognate tRNAs. The selectivity and species-specificity in the recognitions of both amino acid and tRNA by aaRSs play a vital role in maintaining the fidelity of protein synthesis. We report here the first crystal structure of human tryptophanyl-tRNA synthetase (hTrpRS) in complex with tRNA^Trp and Trp which, together with biochemical data, reveals the molecular basis of a novel tRNA binding and recognition mechanism. hTrpRS recognizes the tRNA acceptor arm from the major groove; however, the 3′ end CCA of the tRNA makes a sharp turn to bind at the active site with a deformed conformation. The discriminator base A73 is specifically recognized by an α-helix of the unique N-terminal domain and the anticodon loop by an α-helix insertion of the C-terminal domain. The N-terminal domain appears to be involved in Trp activation, but not essential for tRNA binding and acylation. Structural and sequence comparisons suggest that this novel tRNA binding and recognition mechanism is very likely shared by other archaeal and eukaryotic TrpRSs, but not by bacterial TrpRSs. Our findings provide insights into the molecular basis of tRNA specificity and species-specificity.     

Hinge residue Ile260 of DNA polymerase beta is important for enzyme activity and fidelity

DNA polymerases ensure efficient insertion of the correct dNTP into the DNA substrate. They have evolved mechanisms for discriminating among very similar dNTP substrates. DNA polymerase beta is a repair polymerase that provides a model system for a direct study of insertion fidelity. In this study, we examined the role of hinge residue Ile260 of the rat Polbeta on enzyme activity and accuracy. We changed residue I260 to every other amino acid residue and used genetic screens to assess the activity and fidelity of the resulting mutants. The I260D, -E, -K, -N, and -R mutants are significantly less active than wild-type Polbeta. Interestingly, I260H and I260Q are active but exhibit mutator activity. This suggests that the nonpolar nature of residue 260 is important for maintaining the activity and fidelity of Polbeta. We employ molecular modeling as an aid in explaining the observed phenotypes and propose a mechanism whereby the positioning of the DNA substrate in the enzyme and within the surface of the hinge may be a key player in forming an optimal active site for phosphodiester bond formation between Watson-Crick base pairs.     

Influence of enzyme-substrate contacts located outside the EcoRI recognition site on cleavage of duplex oligodeoxyribonucleotide substrates by EcoRI endonuclease.

A complete understanding of the sequence-specific interaction between the EcoRI restriction endonuclease and its DNA substrate requires identification of all contacts between the enzyme and substrate, and evaluation of their significance. We have searched for possible contacts adjacent to the recognition site, GAATTC, by using a series of substrates with differing lengths of flanking sequence. Each substrate is a duplex of non-self-complementary oligodeoxyribonucleotides in which the recognition site is flanked by six base pairs on one side and from zero to three base pairs on the other. Steady-state kinetic values were determined for the cleavage of each strand of these duplexes. A series of substrates in which the length of flanking sequence was varied on both sides of the hexamer was also examined. The enzyme cleaved both strands of each of the substrates. Decreasing the flanking sequence to fewer than three base pairs on one side of the recognition site induced an asymmetry in the rates of cleavage of the two strands. The scissile bond nearest the shortening sequence was hydrolyzed with increasing rapidity as base pairs were successively removed. Taken together, the KM and kcat values obtained may be interpreted to indicate the relative importance of several likely enzyme-substrate contacts located outside the canonical hexameric recognition site.     

The fidelity of DNA synthesis catalyzed by derivatives of Escherichia coli DNA polymerase I

The fidelity of DNA synthesis by an exonuclease-proficient DNA polymerase results from the selectivity of the polymerization reaction and from exonucleolytic proofreading. We have examined the contribution of these two steps to the fidelity of DNA synthesis catalyzed by the large Klenow fragment of Escherichia coli DNA polymerase I, using enzymes engineered by site-directed mutagenesis to inactivate the proofreading exonuclease. Measurements with two mutant Klenow polymerases lacking exonuclease activity but retaining normal polymerase activity and protein structure demonstrate that the base substitution fidelity of polymerization averages one error for each 10,000 to 40,000 bases polymerized, and can vary more than 30-fold depending on the mispair and its position. Steady-state enzyme kinetic measurements of selectivity at the initial insertion step by the exonuclease-deficient polymerase demonstrate differences in both the Km and the Vmax for incorrect versus correct nucleotides. Exonucleolytic proofreading by the wild-type enzyme improves the average base substitution fidelity by 4- to 7-fold, reflecting efficient proofreading of some mispairs and less efficient proofreading of others. The wild-type polymerase is highly accurate for -1 base frameshift errors, with an error rate of less than or equal to 10(-6). The exonuclease-deficient polymerase is less accurate, suggesting that proofreading also enhances frameshift fidelity. Even without a proofreading exonuclease, Klenow polymerase has high frameshift fidelity relative to several other DNA polymerases, including eucaryotic DNA polymerase-alpha, an exonuclease-deficient, 4-subunit complex whose catalytic subunit is almost three times larger. The Klenow polymerase has a large (46 kDa) domain containing the polymerase active site and a smaller (22 kDa) domain containing the active site for the 3'----5' exonuclease. Upon removal of the small domain, the large polymerase domain has altered base substitution error specificity when compared to the two-domain but exonuclease-deficient enzyme. It is also less accurate for -1 base errors at reiterated template nucleotides and for a 276-nucleotide deletion error. Thus, removal of a protein domain of a DNA polymerase can affect its fidelity.     

Enzymatic and nucleotide sequence studies of a kanamycin-inactivating enzyme encoded by a plasmid from thermophilic bacilli in comparison with that encoded by plasmid pUB110

The product of a kanamycin resistance gene encoded by plasmid pTB913 isolated from a thermophilic bacillus was identified as a kanamycin nucleotidyltransferase which is similar to that encoded by plasmid pUB110 from a mesophile, Staphylococcus aureus. The enzyme encoded by pTB913 was more thermostable than that encoded by pUB110. In view of a close resemblance of restriction endonuclease cleavage maps around the BglII site in the structural genes of both enzymes, ca. 1,200 base pairs were sequenced, followed by amino-terminal amino acid sequencing of the enzyme. The two nucleotide sequences were found to be identical to each other except for only one base in the midst of the structural gene. Each structural gene, initiating from a GUG codon as methionine, was composed of 759 base pairs and 253 amino acid residues (molecular weight, ca. 29,000). The sole difference was transversion from a cytosine (pUB110) to an adenine (pTB913) at a position + 389, counting the first base of the initiation codon as + 1. That is, a threonine at position 130 for the pUB110-coded kanamycin nucleotidyltransferase was replaced by a lysine for the pTB913-coded enzyme. The difference in thermostability between the two enzymes caused by a single amino acid replacement is discussed in light of electrostatic effects.     

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

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

Multiple cleavage activities of endonuclease V from Thermotoga maritima: recognition and strand nicking mechanism

Endonuclease V is a deoxyinosine 3'-endonuclease which initiates removal of inosine from damaged DNA. A thermostable endonuclease V from the hyperthermophilic bacterium Thermotoga maritima has been cloned and expressed in Escherichia coli. The DNA recognition and reaction mechanisms were probed with both double-stranded and single-stranded oligonucleotide substrates which contained inosine, abasic site (AP site), uracil, or mismatches. Gel mobility shift and kinetic analyses indicate that the enzyme remains bound to the cleaved inosine product. This slow product release may be required in vivo to ensure an orderly process of repairing deaminated DNA. When the enzyme is in excess, the primary nicked products experience a second nicking event on the complementary strand, leading to a double-stranded break. Cleavage at AP sites suggests that the enzyme may use a combination of base contacts and local distortion for recognition. The weak binding to uracil sites may preclude the enzyme from playing a significant role in repair of such sites, which may be occupied by uracil-specific DNA glycosylases. Analysis of cleavage patterns of all 12 natural mismatched base pairs suggests that purine bases are preferrentially cleaved, showing a general hierarchy of A = G > T > C. A model accounting for the recognition and strand nicking mechanism of endonuclease V is presented.     

Clamping down on weak terminal base pairs: oligonucleotides with molecular caps as fidelity-enhancing elements at the 5'- and 3'-terminal residues

The base-pairing fidelity of oligonucleotides depends on the identity of the nucleobases involved and the position of matched or mismatched base pairs in the duplex. Nucleobases forming weak base pairs, as well as a terminal position favor mispairing. We have searched for 5′-appended acylamido caps that enhance the stability and base-pairing fidelity of oligonucleotides with a 5′-terminal 2′-deoxyadenosine residue using combinatorial synthesis and MALDI-monitored nuclease selections. This provided the residue of 4-(pyren-1-yl)butyric acid as a lead. Lead optimization gave (S)-N-(pyren-1-ylmethyl)pyrrolidine-3-phosphate as a cap that increases duplex stability and base-pairing fidelity. For the duplex of 5′-AGGTTGAC-3′ with its fully complementary target, this cap gives an increase in the UV melting point T_m of +10.9°C. The T_m is 6.3–8.3°C lower when a mismatched nucleobase faces the 5′-terminal dA residue. The optimized cap can be introduced via automated DNA synthesis. It was combined with an anthraquinone carboxylic acid residue as a cap for the 3′-terminal residue. A doubly capped dodecamer thus prepared gives a melting point decrease for double-terminal mismatches that is 5.7–5.9°C greater than that for the unmodified control duplex.     

Initiation signals for complementary strand DNA synthesis in the region of the replication origin of the Escherichia coli chromosome.

We have used an in vivo plasmid-phi X174 packaging system to detect replication initiation signals in the region of the replication origin (oriC) of the Escherichia coli chromosome. The results obtained are summarized as follows: (i) Neither within nor close to oriC effective signals for initiating complementary strand synthesis could be detected. We conclude that initiation mechanisms for leading and lagging strand synthesis at oriC are not identical to any known priming mechanism of DNA synthesis. (ii) At least five signals that can function as complementary strand origins on ss-plasmid DNA are located in a region about 2000-3300 base pairs away from oriC in the clockwise direction on the chromosome. We suggest that these signals are protein n' like recognition sequences since they are dependent for their activity on dnaB protein and show sequence similarities to other putative n' recognition sequences. Surprisingly, some of the signals are located on the template for leading strand synthesis.