Energy barriers and rates of tautomeric transitions in DNA bases: ab initio quantum chemical study

Tautomeric transitions of DNA bases are proton transfer reactions, which are important in biology. These reactions are involved in spontaneous point mutations of the genetic material. In the present study, intrinsic reaction coordinates (IRC) analyses through ab initio quantum chemical calculations have been carried out for the individual DNA bases A, T, G, C and also A:T and G:C base pairs to estimate the kinetic and thermodynamic barriers using MP2/6-31G** method for tautomeric transitions. Relatively higher values of kinetic barriers (about 50-60 kcal/mol) have been observed for the single bases, indicating that tautomeric alterations of isolated single bases are quite unlikely. On the other hand, relatively lower values of the kinetic barriers (about 20-25 kcal/mol) for the DNA base pairs A:T and G:C clearly suggest that the tautomeric shifts are much more favorable in DNA base pairs than in isolated single bases. The unusual base pairing A':C, T':G, C':A or G':T in the daughter DNA molecule, resulting from a parent DNA molecule with tautomeric shifts, is found to be stable enough to result in a mutation. The transition rate constants for the single DNA bases in addition to the base pairs are also calculated by computing the free energy differences between the transition states and the reactants.     

X-ray analysis of a DNA dodecamer containing 2'-deoxy-N4-methoxycytidine.

The crystal structure of DNA dodecamer with the sequence of d(CGCAAATTXGCG), where X is 2'-deoxy-N4-methoxycytidine, has been determined by X-ray analysis. The dodecamers form a double helix with B-form conformation. The electron density indicates that the two modified cytosine bases respectively make a pair with the adenine bases on the opposite strand in a manner of Watson-Crick geometry and that the methoxy groups are in anti conformation to the N3 atom.     

Modeling of Escherichia coli Endonuclease V structure in complex with DNA

Endonuclease V (EndoV) is a metal-dependent DNA repair enzyme involved in removal of deaminated bases (e.g., deoxyuridine, deoxyinosine, and deoxyxanthosine), with pairing specificities different from the original bases. Homologs of EndoV are present in all major phyla from bacteria to humans and their function is quite well analyzed. EndoV has been combined with DNA ligase to develop an enzymatic method for mutation scanning and has been engineered to obtain variants with different substrate specificities that serve as improved tools in mutation recognition and cancer mutation scanning. However, little is known about the structure and mechanism of substrate DNA binding by EndoV. Here, we present the results of a bioinformatic analysis and a structural model of EndoV from Escherichia coli in complex with DNA. The structure was obtained by a combination of fold-recognition, comparative modeling, de novo modeling and docking methods. The modeled structure provides a convenient tool to study protein sequence-structure-function relationships in EndoV and to engineer its further variants.     

X-ray analyses of two DNA dodecamers containing N4-methoxy-cytosine paired with adenine or guanine.

To investigate mismatch of base-pairings in relation to mutagenesis by oxyamines, crystal structures of two DNA dodecamers with the sequence d(CGCZAA TTmo4CGCG) (Z = A or G), containing N4-methoxy-cytosine (mo4C), have been determined by X-ray analysis. These dodecamers essentially form right-handed B-form duplexes, respectively. In the dodecamer with Z = A, the two mo4C residues are adapted in imino form with the anti methoxyl group to form pairs with A on the opposite strand in a manner of Watson-Crick fashion. While in the dodecamer with Z = G, one mo4C in amino form with the anti methoxyl group forms a normal Watson-Crick pair with G, but the other one in imino form with syn methoxyl conformation wobbles with G. Based on these results, possible mutation mechanism has been proposed.     

Spectrum of spontaneous frameshift mutations. Sequences of bacteriophage T4 rII gene frameshifts

The DNA sequences of 185 independent spontaneous frameshift mutations in the rIIB gene of bacteriophage T4 are described. Approximately half of the frameshifts, including those at hot spot sites, are fully consistent with classical proposals that frameshift mutations are produced by a mechanism involving the misaligned pairing of repeated DNA sequences. However, the remaining frameshifts are inconsistent with this model. Correlations between the positions of two base-pair frameshifts and the bases of DNA hairpins suggest that local DNA topology might influence frameshift mutation. Warm spots for larger deletions share the property of having endpoints adjacent to DNA sequences whose complementarity to sequences a few base-pairs away suggest that non-classical DNA misalignments may participate in deletion mutation. A model for duplication mutation as a consequence of strand displacement synthesis is discussed. In all, 15 frameshifts were complex combinations of frameshifts and base substitutions. Three of these were identical, and have extended homology to a sequence 256 base-pairs away that is likely to participate in the mutational event; the remainder are unique combinations of frameshifts and transversions. The frequency and diversity of complex mutants suggest a challenge to the assumption that the molecular evolution of DNA must depend primarily upon the accumulation of single nucleotide changes.     

Theoretical study on the binding mechanism between N6-methyladenine and natural DNA bases

N6-methyladenine (m⁶A) is a rare base naturally occurring in DNA. It is different from the base adenine due to its N-CH₃. Therefore, the base not only pairs with thymine, but also with other DNA bases (cytosine, adenine and guanine). In this work, Møller-Plesset second-order (MP2) method has been used to investigate the binding mechanism between m⁶A and natural DNA bases in gas phase and in aqueous solution. The results show that N-CH₃ changed the way of N6-methyladenine binding to natural DNA bases. The binding style significantly influences the stability of base pairs. The trans-m⁶A:G and trans-m⁶A:C conformers are the most stable among all the base pairs. The existence of solvent can remarkably reduce the stability of the base pairs, and the DNA bases prefer pairing with trans-m⁶A to cis-m⁶A. Besides, the properties of these hydrogen bonds have been analyzed by atom in molecules (AIM) theory, natural bond orbital (NBO) analysis and Wiberg bond indexes (WBI). In addition, pairing with m⁶A decreases the binding energies compared to the normal Watson-Crick base pairs, it may explain the instability of the N6 site methylated DNA in theory.

Proofreading activity of DNA polymerase Pol2 mediates 3'-end processing during nonhomologous end joining in yeast

Genotoxic agents that cause double-strand breaks (DSBs) often generate damage at the break termini. Processing enzymes, including nucleases and polymerases, must remove damaged bases and/or add new bases before completion of repair. Artemis is a nuclease involved in mammalian nonhomologous end joining (NHEJ), but in Saccharomyces cerevisiae the nucleases and polymerases involved in NHEJ pathways are poorly understood. Only Pol4 has been shown to fill the gap that may form by imprecise pairing of overhanging 3' DNA ends. We previously developed a chromosomal DSB assay in yeast to study factors involved in NHEJ. Here, we use this system to examine DNA polymerases required for NHEJ in yeast. We demonstrate that Pol2 is another major DNA polymerase involved in imprecise end joining. Pol1 modulates both imprecise end joining and more complex chromosomal rearrangements, and Pol3 is primarily involved in NHEJ-mediated chromosomal rearrangements. While Pol4 is the major polymerase to fill the gap that may form by imprecise pairing of overhanging 3' DNA ends, Pol2 is important for the recession of 3' flaps that can form during imprecise pairing. Indeed, a mutation in the 3'-5' exonuclease domain of Pol2 dramatically reduces the frequency of end joins formed with initial 3' flaps. Thus, Pol2 performs a key 3' end-processing step in NHEJ.     

The intrinsic mechanics of B-DNA in solution characterized by NMR

Experimental characterization of the structural couplings in free B-DNA in solution has been elusive, because of subtle effects that are challenging to tackle. Here, the exploitation of the NMR measurements collected on four dodecamers containing a substantial set of dinucleotide sequences provides new, consistent correlations revealing the DNA intrinsic mechanics. The difference between two successive residual dipolar couplings (ΔRDCs) involving C6/8-H6/8, C3′-H3′ and C4′-H4′ vectors are correlated to the ³¹P chemical shifts (δP), which reflect the populations of the BI and BII backbone states. The δPs are also correlated to the internucleotide distances (D_inter) involving H6/8, H2′ and H2″ protons. Calculations of NMR quantities on high resolution X-ray structures and controlled models of DNA enable to interpret these couplings: the studied ΔRDCs depend mostly on roll, while D_inter are mainly sensitive to twist or slide. Overall, these relations demonstrate how δP measurements inform on key inter base parameters, in addition to probe the BI↔BII backbone equilibrium, and shed new light into coordinated motions of phosphate groups and bases in free B-DNA in solution. Inspection of the 5′ and 3′ ends of the dodecamers also supplies new information on the fraying events, otherwise neglected.     

Crystallographic studies on damaged DNAs: III. N(4)-methoxycytosine can form both Watson-Crick type and wobbled base pairs in a B-form duplex.

To investigate the mutation mechanism of purine transition in DNA damaged with methoxyamine, a DNA dodecamer with the sequence d(CGCGAATTmo(4)CGCG), where mo(4)C is 2'-deoxy-N(4)-methoxycytidine, has been synthesized and its crystal structure determined. Two dodecamers form a B-form duplex. Electron density maps clearly show that one of the two mo(4)C residues forms a pair with a guanine residue of the opposite strand, the geometry being the canonical Watson-Crick type, and that the other mo(4)C residue forms a wobble pair with the opposite guanine residue. These two pairings are ascribed to the tautomerization of the methoxylated cytosine moieties between the amino and imino forms.     

Probing the conformations of eight cloned DNA dodecamers; CGCGAATTCGCG, CGCGTTAACGCG, CGCGTATACGCG, CGCGATATCGCG, CGCAAATTTGCG, CGCTTTAAAGCG, CGCGGATCCGCG and CGCGGTACCGCG.

The self complementary DNA dodecamers d(CGCGAATTCGCG), d(CGCGTTAACGCG), d(CGCGTATACGCG), d(CGCGATATCGCG), d(CGCAAATTTGCG), d(CGCTTTAAAGCG), d(CGCGGATCCGCG) and d(CGCGGTACCGCG) have been cloned into the Smal site of plasmid pUC19. Radiolabelled polylinker fragments containing these inserts have been digested with nucleases and chemical agents, probing the structure of the central AT base pairs. The sequences AATT and AAATTT are relatively resistant to digestion by DNase I, micrococcal nuclease and hydroxyl radicals, consistent with the suggestion that they possess a narrow minor groove. Nuclease digestion of TTAA is much more even, and comparable to that at mixed sequence DNA. TpA steps in ATAT, TATA and GTAC are cut less well by DNAse I than in TTAA. DNasel cleavage of surrounding bases, especially CpG is strongly influenced by the nature of the central sequence.     

Proofreading Activity of DNA Polymerase Pol2 Mediates 3′-End Processing during Nonhomologous End Joining in Yeast

Genotoxic agents that cause double-strand breaks (DSBs) often generate damage at the break termini. Processing enzymes, including nucleases and polymerases, must remove damaged bases and/or add new bases before completion of repair. Artemis is a nuclease involved in mammalian nonhomologous end joining (NHEJ), but in Saccharomyces cerevisiae the nucleases and polymerases involved in NHEJ pathways are poorly understood. Only Pol4 has been shown to fill the gap that may form by imprecise pairing of overhanging 3′ DNA ends. We previously developed a chromosomal DSB assay in yeast to study factors involved in NHEJ. Here, we use this system to examine DNA polymerases required for NHEJ in yeast. We demonstrate that Pol2 is another major DNA polymerase involved in imprecise end joining. Pol1 modulates both imprecise end joining and more complex chromosomal rearrangements, and Pol3 is primarily involved in NHEJ-mediated chromosomal rearrangements. While Pol4 is the major polymerase to fill the gap that may form by imprecise pairing of overhanging 3′ DNA ends, Pol2 is important for the recession of 3′ flaps that can form during imprecise pairing. Indeed, a mutation in the 3′-5′ exonuclease domain of Pol2 dramatically reduces the frequency of end joins formed with initial 3′ flaps. Thus, Pol2 performs a key 3′ end-processing step in NHEJ.     

A generalized conformational energy function of DNA derived from molecular dynamics simulations

Proteins recognize DNA sequences by two different mechanisms. The first is direct readout, in which recognition is mediated by direct interactions between the protein and the DNA bases. The second is indirect readout, which is caused by the dependence of conformation and the deformability of the DNA structure on the sequence. Various energy functions have been proposed to evaluate the contribution of indirect readout to the free-energy changes in complex formations. We developed a new generalized energy function to estimate the dependence of the deformability of DNA on the sequence. This function was derived from molecular dynamics simulations previously conducted on B-DNA dodecamers, each of which had one possible tetramer sequence embedded at its center. By taking the logarithm of the probability distribution function (PDF) for the base-step parameters of the central base-pair step of the tetramer, its ability to distinguish the native sequence from random ones was superior to that with the previous method that approximated the energy function in harmonic form. From a comparison of the energy profiles calculated with these two methods, we found that the harmonic approximation caused significant errors in the conformational energies of the tetramers that adopted multiple stable conformations.     

Cognate base‐pair selectivity of hydrophobic unnatural bases in DNA ligation by T4 DNA ligase

We present cognate base pair selectivity in template‐dependent ligation by T4 DNA ligase using a hydrophobic unnatural base pair (UBP), Ds‐Pa. T4 DNA ligase efficiently recognizes the Ds‐Pa pairing at the conjugation position, and Ds excludes the noncognate pairings with the natural bases. Our results indicate that the hydrophobic base pairing is allowed in enzymatic ligation with higher cognate base‐pair selectivity, relative to the hydrogen‐bond interactions between pairing bases. The efficient ligation using Ds‐Pa can be employed in recombinant DNA technology using genetic alphabet expansion, toward the creation of semi‐synthetic organisms containing UBPs.     

Structure and axis curvature in two dA6 x dT6 DNA oligonucleotides: comparison of molecular dynamics simulations with results from crystallography and NMR spectroscopy

Molecular dynamics (MD) simulations have been performed on the A6 containing DNA dodecamers d(GGCAAAAAACGG) solved by NMR and d(CGCAAAAAAGCG) solved by crystallography. The experimental structures differ in the direction of axis bending and in other small but important aspects relevant to the DNA curvature problem. Five nanosecond MD simulations of each sequence have been performed, beginning with both the NMR and crystal forms as well as canonical B-form DNA. The results show that all simulations converge to a common form in close proximity to the observed NMR structure, indicating that the structure obtained in the crystal is likely a strained form due to packing effects. A-tracts in the MD model are essentially straight. The origin of axis curvature is found at pyrimidine-purine steps in the flanking sequences.     

The mutational specificity of DNA polymerase-beta during in vitro DNA synthesis. Production of frameshift, base substitution, and deletion mutations

The frequency and specificity of mutations produced in vitro by eucaryotic DNA polymerase-beta have been determined in a forward mutation assay using a 250-base target sequence in M13mp2 DNA. Homogeneous DNA polymerase-beta, isolated from four different sources, produces mutations at a frequency of 4-6%/single round of gap-filling DNA synthesis. DNA sequence analyses of 460 independent mutants resulting from this error-prone DNA synthesis demonstrate a wide variety of mutational events. Frameshift and base substitutions are made at approximately equal frequency and together comprise about 90% of all mutations. Two mutational "hot spots" for frameshift and base substitution mutations were observed. The characteristics of the mutations at these sites suggest that certain base substitution errors result from dislocation of template bases rather than from direct mispair formation by DNA polymerase-beta. When considering the entire target sequence, single-base frameshift mutations occur primarily in runs of identical bases, usually pyrimidines. The loss of a single base occurs 20-80 times more frequently than single-base additions and much more frequently than the loss of two or more bases. Base substitutions occur at many sites throughout the target, representing a wide spectrum of mispair formations. Averaged over a large number of phenotypically detectable sites, the base substitution error frequency is greater than one mistake for every 5000 bases polymerized. Large deletion mutations are also observed, at a frequency more than 10-fold over background, indicating that purified DNA polymerases alone are capable of producing such deletions. These data are discussed in relation to the physical and kinetic properties of the purified enzymes and with respect to the proposed role for this DNA polymerase in vivo.     

A conformational rationale for the replication of nucleic acids under prebiotic conditions

The Watson-Crick model at once gave an explanation for the mechanism of replication of DNA. But the hydrogen-bonding forces between the bases alone are not enough for the specificity of base-pairing mechanisms, since any pair of bases can be positioned to have at least two hydrogen bonds. In the present-day biological organisms, sophisticated enzymatic machinery is supposed to constrain the ribose-phosphate backbone to have regular structure, aiding the self-templating duplication. For the prebiotic stage, whence sophisticated enzymes would not have been evolved, we propose a novel double helical conformation of DNA wherein the two sugar-phosphate backbones are pulled towards each other by (CH  O) hydrogen bonds conferring stereospecificity for the formation of (A : T)- and (G : C)-pairs, in the self-templating chains of DNA. Our model-building efforts and computer calculations endorse the stereochemical feasibility of the conformation. The pairing of homologous sequences of two double helices of DNA is explained by direct hydrogen-bonding interactions in our model and it is thus relevant to the present-day biological functions also, at least in some stages of the cell-cycle.     

Effect of base-pair sequence on the conformations and thermally induced transitions in oligodeoxyribonucleotides containing only AT base pairs

Tm curves, CD spectra, and kinetics results of the self-complementary DNA dodecamers d(A6T6), d(A3T3A3T3), d(A2T2A2T2A2T2), d(ATATATATATAT), and d(T6A6) demonstrate that the thermal transitions of these oligomers at low salt concentration involve a hairpin intermediate. At high salt concentrations (greater than 0.1 M Na+) only a duplex to denatured-strand transition appears to occur. The temperature and salt-concentration regions of the transitions are very sequence dependent. Alternating-type AT sequences have a lower duplex stability and a greater tendency to form hairpins than sequences containing more nonalternating AT base pairs. Of the two nonalternating sequences, d(T6A6) is significantly less stable than d(A6T6). Both oligomers have CD curves that are very similar to the unusual CD spectrum of poly(dA).poly(dT). The Raman spectra of these two oligomers are also quite similar, but at low temperature, small intensity differences in two backbone modes and three nucleoside vibrations are obtained. The hairpin to duplex transition for the AT dodecamers was examined by salt-jump kinetics measurements. The transition is faster than transitions for palindromic-sequence oligomers containing terminal GC base pairs. Stopped-flow kinetics studies indicate that the transition is second order and has a relatively low activation energy. The reaction rate increases with increasing ionic strength. These results are consistent with a three-step mechanism for the hairpin to duplex reaction: (i) fraying of the hairpin oligomers' terminal base pairs, (ii) a rate-determining bimolecular step involving formation of a cruciform-type intermediate from two hairpin oligomers with open terminal base pairs, and (iii) base-pair migration and formation in the intermediate to give the duplex.     

Alternative pathways for the in vivo repair of O6-alkylguanine and O4-alkylthymine in Escherichia coli: the adaptive response and nucleotide excision repair.

The in vivo removal of three different O-alkylated bases from DNA was measured in Escherichia coli. Using monoclonal antibodies specific for O6-methylguanine, O6-ethylguanine and O4-ethylthymine we have monitored the removal of these lesions from six different strains to assess the relative contributions of the adaptive response and of nucleotide excision repair. During the first hour after DNA alkylation, O6-methylguanine, O6-ethylguanine and O4-ethylthymine lesions were repaired almost exclusively by nucleotide excision, except when the adaptive response was being constitutively expressed. In wild-type E. coli the adaptive response began to contribute to O6-methylguanine repair about one hour after alkylation, the time required for the full induction of the ada DNA methyltransferase. In contrast, the adaptive response did not play such a large role in the repair of O6-ethylguanine and O4-ethylthymine in wild-type E. coli, presumably because DNA ethylation damage is a poor inducer of the adaptive response; possible reasons for this poor induction are discussed. The repair of all three O-alkylated lesions was virtually absent in ada- uvr- bacteria suggesting that no alternative pathway is available for their repair, at least during the first two hours after alkylation. When the repair of O-alkylated bases was compromised by an ada- or by a uvr- mutation, the bacteria became more sensitive to alkylation induced killing and mutation.     

DNA models for A, B, C and D conformations related to fiber X-ray, infrared and NMR measurements

A conformational analysis of the A, B, C and D DNA forms was made in order to establish molecular models presenting a good agreement with experimental data obtained from fiber X-ray, infrared linear dichroism and 31P NMR. The proposed models have been refined and do present good stereochemistry and optimized H-bond distances between bases associated with the Watson-Crick pairing. The DNA conformations proposed are a left handed double helix for the C form and right handed helices for A, B and D. Relations to conformational transitions between these forms are discussed.