Monte Carlo Simulation of Hydration of the Guanine-Uracil Pairs with Guanine in Two Tautomeric Forms: Contribution of Water Bridging to Relative Stability of Mispairs

Abstract

Results are presented from Monte Carlo simulation of hydration of guanine-uracil mispairs by 25 and 50 water molecules. The hydration shells of three mispairs formed between “normal” dioxo form of uracil (U) and three forms of guanine (“normal” amino-oxo tautomer G and two rotamers of the “rare” amino-hydroxy tautomer G*) depend on the tautomeric forms of the guanine molecule. The simulation shows the important role of hydration effects on the relative stability of the mispairs.     

A molecular dynamics study of the effect of G.T mispairs on the conformation of DNA in solution

The effect of G.T mispair incorporation into a double-helical environment was examined by molecular dynamics simulation. The 60-ps simulations performed on the two hexanucleotide duplexes d (G3C3)2 and d(G3TC2)2 included 10 Na+ counterions and first hydration shell waters. The resulting backbone torsional angle trajectories were analyzed to select time spans representative of conformational domains. The average backbone angles and helical parameters of the last time span for both duplexes are reported. During the simulation the hexamers retained B-type DNA structures that differed from typical A- or B-DNA forms. The overall helical structures for the two duplexes are vary similar. The presence of G.T mispairs did not alter the overall helical structure of the oligonucleotide duplex. Large propeller twist and buckle angles were obtained for both duplexes. The purine/pyrimidine crossover step showed a large decrease in propeller twist in the normal duplex but not in the mismatch duplex. Upon the formation of wobble mispairs in the mismatched duplex, the guanines moved into the minor groove and the thymines moved into the major groove. This helped prevent purine/purine clash and created a deformation in the relative orientation of the glycosidic bonds. It also exposed the free O4 of the thymines in the major groove and N2 of the guanines in the minor groove to interactions with solvent and counterions. These factors seemed to contribute to the apparently higher rigidity of the mismatched duplex during the simulation.     

The structure of guanosine-thymidine mismatches in B-DNA at 2.5-A resolution

The structure of the deoxyoligomer d(C-G-C-G-A-A-T-T-T-G-C-G) was determined at 2.5-A resolution by single crystal x-ray diffraction techniques. The final R factor is 18% with the location of 71 water molecules. The oligomer crystallizes in a B-DNA-type conformation, with two strands interacting to form a dodecamer duplex. The double helix consists of four A X T and six G X C Watson-Crick base pairs and two G X T mismatches. The G X T pairs adopt a "wobble" structure with the thymine projecting into the major groove and the guanine into the minor groove. The mispairs are accommodated in the normal double helix by small adjustments in the conformation of the sugar phosphate backbone. A comparison with the isomorphous parent compound containing only Watson-Crick base pairs shows that any changes in the structure induced by the presence of G X T mispairs are highly localized. The global conformation of the duplex is conserved. The G X T mismatch has already been studied by x-ray techniques in A and Z helices where similar results were found. The geometry of the mispair is essentially identical in all structures so far examined, irrespective of the DNA conformation. The hydration is also similar with solvent molecules bridging the functional groups of the bases via hydrogen bonds. Hydration may be an important factor in stabilizing G X T mismatches. A characteristic of Watson-Crick paired A X T and G X C bases is the pseudo 2-fold symmetry axis in the plane of the base pairs. The G X T wobble base pair is pronouncedly asymmetric. This asymmetry, coupled with the disposition of functional groups in the major and minor grooves, provides a number of features which may contribute to the recognition of the mismatch by repair enzymes.     

Crystal structure of human methyl-binding domain IV glycosylase bound to abasic DNA

The mammalian repair protein MBD4 (methyl-CpG-binding domain IV) excises thymine from mutagenic G·T mispairs generated by deamination of 5-methylcytosine (mC), and downstream base excision repair proteins restore a G·C pair. MBD4 is also implicated in active DNA demethylation by initiating base excision repair of G·T mispairs generated by a deaminase enzyme. The question of how mismatch glycosylases attain specificity for excising thymine from G·T, but not A·T, pairs remains largely unresolved. Here, we report a crystal structure of the glycosylase domain of human MBD4 (residues 427-580) bound to DNA containing an abasic nucleotide paired with guanine, providing a glimpse of the enzyme-product complex. The mismatched guanine remains intrahelical, nestled into a recognition pocket. MBD4 provides selective interactions with the mismatched guanine (N1H, N2H(2)) that are not compatible with adenine, which likely confer mismatch specificity. The structure reveals no interactions that would be expected to provide the MBD4 glycosylase domain with specificity for acting at CpG sites. Accordingly, we find modest 1.5- to 2.7-fold reductions in G·T activity upon altering the CpG context. In contrast, 37- to 580-fold effects were observed previously for thymine DNA glycosylase. These findings suggest that specificity of MBD4 for acting at CpG sites depends largely on its methyl-CpG-binding domain, which binds preferably to G·T mispairs in a methylated CpG site. MBD4 glycosylase cannot excise 5-formylcytosine (fC) or 5-carboxylcytosine (caC), intermediates in a Tet (ten eleven translocation)-initiated DNA demethylation pathway. Our structure suggests that MBD4 does not provide the electrostatic interactions needed to excise these oxidized forms of mC.     

Fidelity of mispair formation and mispair extension is dependent on the interaction between the minor groove of the primer terminus and Arg668 of DNA polymerase I of Escherichia coli

The hydrogen bonding interactions between the Klenow fragment of Escherichia coli DNA polymerase I with the proofreading exonuclease inactivated (KF(-)) and the minor groove of DNA were examined with modified oligodeoxynucleotides in which 3-deazaguanine (3DG) replaced guanine. This substitution would prevent a hydrogen bond from forming between the polymerase and that one site on the DNA. If the hydrogen bonding interaction were important, then we should observe a decrease in the rate of reaction. The steady-state and pre-steady-state kinetics of DNA replication were measured with 10 different oligodeoxynucleotide duplexes in which 3DG was placed at different positions. The largest decrease in the rate of replication was observed when 3DG replaced guanine at the 3'-terminus of the primer. The effect of this substitution on mispair extension and formation was then probed. The G to 3DG substitution at the primer terminus decreased the k(pol) for the extension past G/C, G/A, and G/G base pairs but not the G/T base pair. The G to 3DG substitution at the primer terminus also decreased the formation of correct base pairs as well as incorrect base pairs. However, in all but two mispairs, the effect on correct base pairs was much greater than that of mispairs. These results indicate that the hydrogen bond between Arg668 and the minor groove of the primer terminus is important in the fidelity of both formation and extension of mispairs. These experiments support a mechanism in which Arg668 forms a hydrogen bonding fork between the minor groove of the primer terminus and the ring oxygen of the deoxyribose moiety of the incoming dNTP to align the 3'-hydroxyl group with the alpha-phosphate of the dNTP. This is one mechanism by which the polymerase can use the geometry of the base pairs to modulate the rate of formation and extension of mispairs.     

Refined crystal structure of an octanucleotide duplex with G . T mismatched base-pairs

Single crystal X-ray diffraction techniques have been used to determine the structure of the DNA octamer d(G-G-G-G-C-T-C-C) at a resolution of 2.25 A. The asymmetric unit consists of two strands coiled about each other to produce an A-type DNA helix. The double helix contains six G . C Watson-Crick base-pairs and two G . T mismatched base-pairs. The mismatches adopt a "wobble" type structure in which both bases retain their major tautomer forms. The double helix is able to accommodate this G . T pairing with little distortion of the overall helical conformation. Crystals of this octamer melt at a substantially lower temperature than do those of a related octamer also containing two G . T base-pairs. We attribute this destabilization to disruption of the hydration network around the mismatch site combined with changes in intermolecular packing. Full details are given of conformational parameters, base stacking, intermolecular contacts and hydration involving 52 solvent molecules.     

Modified SV40 for analysis of mismatch repair in simian and human cells

We have developed a way of introducing specific mispairs into the genome of simian virus 40 and of determining the fate of the mispaired bases in simian and human cells. Mispairs are introduced into viral DNA within the intron of the gene coding for the large T antigen. Each DNA molecule harbors a single mispair in a defined orientation. Transfection of mismatch-containing SV40 DNA into host cells yields plaques, each corresponding to a productive infection initiated by a single viral DNA molecule. Isolation of DNA derived from individual plaques and determination of the DNA sequence at the site of the mispair reveals whether correction occurred and what the repair products are. Here we describe repair patterns for G/T and A/C mispairs in CV-1 African green monkey kidney cells, and for G/T mispairs in human fibroblasts derived from 3 normal individuals, 1 patient with xeroderma pigmentosum (complementation group A), and 3 patients with Bloom's syndrome. G/T mispairs, which arise in resting DNA through the deamination of 5-methylcytosine (mC) to form thymine, are corrected in all cases with extremely high efficiency and nearly always in favor of guanine. In contrast, A/C mispairs are corrected randomly and relatively inefficiently in simian cells.     

A specific mismatch repair event protects mammalian cells from loss of 5-methylcytosine

5-Methylcytosine spontaneously deaminates to form thymine, thus generating G/T mispairs in DNA. We investigated the way in which these lesions are addressed in mammalian cells by introducing specific G/T mispairs into the genome of SV40 and determining the fate of the mismatched bases in simian cells. Mispairs were incorporated in 12 bp synthetic duplexes ligated into SV40 DNA between the BstXI and TaqI restriction sites. Analysis of 347 plaques obtained after transfection of this modified DNA indicated that mispairs were corrected in 343 cases (99%), revealing 314 repair events in favor of guanine (90%) and 29 in favor of thymine (8%). Correction in favor of guanine occurred regardless of the orientation of the mispair in DNA and regardless of whether the mispair was in the commonly methylated CpG dinucleotide. These results attest to a specific mismatch repair pathway that restores G/C pairs lost through deamination of 5-methylcytosine residues.     

Oxidatively generated damage to DNA at 5-methylcytosine mispairs

Oxidatively generated damage to DNA has been implicated as causing mutations that lead to aging and disease. The one-electron oxidation of normal DNA leads to formation of a nucleobase radical cation that hops through the DNA until it is trapped irreversibly, primarily by reaction at guanine. It has been observed that 5-methylcytosine (C(m)) is a mutational "hot-spot". However, C(m) in a Watson-Crick base pair with G is not especially susceptible to oxidatively induced damage. Radical cation hopping is inhibited in duplexes that contain C-A or C-T mispairs, but no reaction is detected at cytosine. In contrast, we find that the one-electron oxidation of DNA that contains C(m)-A or C(m)-T mispairs results primarily in reaction at C(m) even in the presence of GG steps. The reaction at C(m) is attributed to proton coupled electron transfer, which provides a relatively low activation barrier path for reaction at 5-methylcytosine. This enhanced reactivity of C(m) in mispairs may contribute to the formation of mutational hot spots at C(m).     

Structural features and hydration of a dodecamer duplex containing two C.A mispairs.

X-ray diffraction techniques have been used to characterise the crystal and molecular structure of the deoxyoligomer d(C-G-C-A-A-A-T-T-C-G-C-G) at 2.5A resolution. The final R factor is 0.19 with the location of 78 solvent molecules. The oligomer crystallises in a B-DNA type conformation with two strands coiled about each other to produce a duplex. This double helix consists of four A.T and six G.C Watson-Crick base pairs and two C.A mispairs. The mismatched base pairs adopt a "wobble" type structure with the cytosine displaced laterally into the major groove, the adenine into the minor groove. We have proposed that the two close contacts observed in the C.A pairing represent two hydrogen bonds one of which results from protonation of adenine. The mispairs are accommodated in the double helix with small adjustments in the conformation of the sugar-phosphate backbone. Details of the backbone conformation, base stacking interactions, thermal parameters and the hydration are now presented and compared with those of the native oligomer d(C-G-C-G-A-A-T-T-C-G-C-G) and with variations of this sequence containing G.T and G.A mispairs.     

Crystal structure of an RNA 16-mer duplex R(GCAGAGUUAAAUCUGC)2 with nonadjacent G(syn).A+(anti) mispairs

G.A mispairs are one of the most common noncanonical structural motifs of RNA. The 1.9 A resolution crystal structure of the RNA 16-mer r(GCAGAGUUAAAUCUGC)2 has been determined with two isolated or nonadjacent G.A mispairs. The molecule crystallizes with one duplex in the asymmetric unit in space group R3 and unit cell dimensions a = b = c = 49.24 A and alpha = beta = gamma = 51.2 degrees. It is the longest known oligonucleotide duplex at this resolution and isomorphous to the 16-mer duplex with the C.A+ mispairs [Pan, et al., (1998) J. Mol. Biol. 283, 977-984]. The C.A+ mispair behaves like a wobble pair while the G.A+ does not. The G.A mispairs are protonated at N1 of the adenines as in the C.A+ mispairs, and two hydrogen bonds in the G(syn).A+(anti) conformation are formed. The syn guanine is stabilized by an intranucleotide hydrogen bond between the 2-amino and the 5'-phosphate groups. The G(syn).A+(anti) conformation can provide a different surface for recognition in the grooves compared to other G.A hydrogen bonding schemes. The major groove is widened between the two mispairs allowing access to ligands. One of the 3-fold axes is occupied by a sodium ion and a water molecule, while a second is occupied by another water molecule.     

Why do O-alkylguanine and O-alkylthymine miscode? The relationship between the structure of DNA containing O-alkylguanine and O-alkylthymine and the mutagenic properties of these bases

The carcinogenic and mutagenic N-nitroso compounds produce GC to AT and TA to GC transition mutations because they alkylate O⁶ of guanine and O⁴ of thymine. It has been generally assumed that these mutations occur because O⁶-alkylguanine forms a stable mispair with thymine and O⁴-alkylthymine forms a mispair with guanine. Recent studies have shown that this view is mistaken and that the alkylG·T and alkylT·G mispairs are not more stable than their alkylG·C or alkylT·A counterparts. Two possible explanations based on recent structural studies are put forward to account for the miscoding. The first possibility is that the DNA polymerase might mistake O⁶-alkylguanine for adenine, and O⁴-alkylthymine for cytosine, because of the physical similarity of these bases. O⁶-Methylguanine and adenine are similarly lipophilic and X-ray crystallography of the nucleosides has shown a close similarity in bond angles and lengths between O⁶-methylguanine and adenine, and between O⁴-methylthymine and cytosine. The second possible explanation is that the important factor in the miscoding is that the alkylG·T and alkylT·G mispairs retain the Watson-Crick alignment with N1 of the purine juxtaposed to N3 of the pyrimidine while the alkylG·C and alkylT·A pairs adopt a wobble conformation. ³¹P NMR of DNA duplexes show that the phosphodiester links both 3′ and 5′ to the C have to be distorted to accomodate the O⁶-ethylguanine:C pair, whereas there is less distortion of the phosphodiesters 3′ and 5′ to the T in an ethylG·T pair. Recent kinetic measurements show that the essential aspect of base selection in DNA synthesis is the ease of formation of the phosphodiester links on both the 3′ and 5′ side of the incoming base. The Watson-Crick alignment of the alkylG·T and alkylT·G mispairs may facilitate formation of these phosphodiester links, and this alignment rather than the strength of the base pairs and the extent of hydrogen bonding between them may be the crucial factor in the miscoding. If either hypothesis is correct it suggests that previously too much emphasis has been placed on the stability of the normal pairs in the replication of DNA.     

Possible conformations of double-helical polynucleotides containing incorrect base pairs

Theoretical conformational analysis using classical potential functions has shown the possibility of incorporation of nucleotide mispairs with the bases in normal tautomeric forms into the DNA double helix. Incorrect purine-pyrimidine, purine-purine and pyrimidine-pyrimidine pairs can be incorporated into the double helix existing both in A- and B-conformations. The most energy favourable conformations of fragments containing a mispair have all the dihedral angles of the sugar-phosphate backbone within the limits characteristic of double helices consisting of Watson-Crick nucleotide pairs. Incorporation of mispairs is possible practically without the appearance of reduced interatomic contacts. Mutual position of bases in the incorporated mispair does not differ much from their position at the energy minimum of the corresponding isolated base pairs. Conformational parameters of irregular regions of double-stranded polynucleotides containing G:U, I:A, I:A* (syn) and U:C pairs are presented. Distortion of the sugar-phosphate backbone is the least upon incorporation of the G:U pair. Formation of mispairs in the processes of nucleic acid biosynthesis and spontaneous mutagenesis is discussed.     

Long-range oxidative damage in duplex DNA: the effect of bulged G in a G–C tract and tandem G/A mispairs

Two series of duplex DNA oligomers were prepared having an anthraquinone derivative (AQ) covalently linked at a 5′-terminus. Irradiation of the AQ at 350 nm leads to injection of an electron hole (radical cation) into the DNA. The radical cation migrates through the DNA causing reaction primarily at G_n sequences. In one series, GA tandem mispairs are inserted between GG steps to assess the effect of the mispair on the transport of the radical cation, reaction (damage) caused by the radical cation at the mispair, and repair of the resulting damage by formamidopyrimidine DNA glycosylase (Fpg). In the second series, a bulged guanine in a G₃C₂ sequence is interposed between the GG steps. These experiments reveal that neither G/A tandem mispairs nor bulged guanines are significant barriers to long-range charge migration in DNA. The radical cation does not cause reaction at guanines in the G/A tandem mispair. Reaction does occur at the bulged guanine, but it is repaired by Fpg.     

Structure of purine-purine mispairs during misincorporation and extension by Escherichia coli DNA polymerase I

The mechanism by which purine-purine mispairs are formed and extended was examined with the high-fidelity Klenow fragment of Escherichia coli DNA polymerase I with the proofreading exonuclease activity inactivated. The structures of the purine-purine mispairs were examined by comparing the kinetics of mispair formation with adenine versus 7-deazaadenine and guanine versus 7-deazaguanine at four positions in the DNA, the incoming dNTP, the template base, and both positions of the terminal base pair. A decrease in rate associated with a 7-deazapurine substitution would suggest that the nucleotide is in a syn conformation in a Hoogsteen base pair with the opposite base. During mispair formation, the k(pol)/K(d) values for the insertion of dATP opposite A (dATP/A) as well as dATP/G and dGTP/G were decreased greater than 10-fold with the deazapurine in the dNTP. These results suggest that during mispair formation the newly forming base pair is in a Hoogsteen geometry with the incoming dNTP in the syn conformation and the template base in the anti conformation. During mispair extension, the only decrease in k(pol)/K(d) was associated with the G/G base pair in which 7-deazaguanine was in the template strand. These results as well as previous results [McCain et al. (2005) Biochemistry 44, 5647-5659] in which a hydrogen bond was found between the 3-position of guanine at the primer terminus and Arg668 during G/A and G/G mispair extension indicate that the conformation of the purine at the primer terminus is in the anti conformation during mispair extension. These results suggest that purine-purine mispairs are formed via a Hoogsteen geometry in which the dNTP is in the syn conformation and the template is in the anti conformation. During extension, however, the conformation of the primer terminus changes to an anti configuration while the template base may be in either the syn or anti conformations.     

Unique misinsertion specificity of poliota may decrease the mutagenic potential of deaminated cytosines

DNA polymerase iota (poliota) is a distributive error-prone enzyme that can incorporate nucleotides opposite a variety of DNA lesions. Further elongation is, however, either substantially inhibited or completely abolished. Here, we provide evidence that poliota can facilitate the efficient bypass of uracil and its derivatives as well as oxidized cytosine and guanine residues. The fidelity of translesion replication depends upon the lesion encountered. Correct nucleotides were inserted preferentially opposite 7,8-dihydro-8-oxoguanine (8-oxoG) and 5-hydroxycytosine (5-OHC). However, when bypassing uracil, 5-hydroxyuracil (5-OHU) or 5,6-dihydrouracil (5,6-DHU), poliota inserted T and G with a 4- to 26-fold preference over the Watson-Crick base, A. While the T:U, T:5-OHU and T:5,6-DHU mispairs were extended poorly, the G:U, G:5-OHU and G:5,6-DHU mispairs were extended with equal or greater efficiency than the correctly paired primer termini. Thus, poliota-dependent misinsertion of G opposite uracil and its derivatives may actually provide a mechanism whereby mammalian cells can decrease the mutagenic potential of lesions formed via the deamination of cytosine.     

Opposite base-dependent reactions of a human base excision repair enzyme on DNA containing 7,8-dihydro-8-oxoguanine and abasic sites.

The guanine modification 7,8-dihydro-8-oxoguanine (8-oxoG) is a potent premutagenic lesion formed spontaneously at high frequencies in the genomes of aerobic organisms. We have characterized a human DNA repair glycosylase for 8-oxoG removal, hOGH1 (human yeast OGG1 homologue), by molecular cloning and functional analysis. Expression of the human cDNA in a repair deficient mutator strain of Escherichia coli (fpg mutY) suppressed the spontaneous mutation frequency to almost normal levels. The hOGH1 enzyme was localized to the nucleus in cells transfected by constructs of hOGH1 fused to green fluorescent protein. Enzyme purification yielded a protein of 38 kDa removing both formamidopyrimidines and 8-oxoG from DNA. The enzymatic activities of hOGH1 was analysed on DNA containing single residues of 8-oxoG or abasic sites opposite each of the four normal bases in DNA. Excision of 8-oxoG opposite C was the most efficient and was followed by strand cleavage via beta-elimination. However, significant removal of 8-oxoG from mispairs (8-oxoG: T >G >A) was also demonstrated, but essentially without an associated strand cleavage reaction. Assays with abasic site DNA showed that strand cleavage was indeed dependent on the presence of C in the opposite strand, irrespective of the prior removal of an 8-oxoG residue. It thus appears that strand incisions are made only if repair completion results in correct base insertion, whereas excision from mispairs preserves strand continuity and hence allows for error-free correction by a postreplicational repair mechanism.     

A repair system for 8-oxo-7,8-dihydrodeoxyguanine.

Active oxygen species can damage DNA and may play a role in aging and carcinogenesis. We have tested MutY glycosylase for activity on undamaged mispairs as well as mispairs formed with the oxidatively damaged substrates, 8-oxo-7,8-dihydrodeoxyguanine (GO) or 8-oxo-7,8-dihydrodeoxyadenine (AO). MutY acts as a glycosylase on four of the heteroduplexes tested, A/G, A/GO, A/C, and A/AO, removing the undamaged adenine from each substrate. Genetic data suggest that the primary substrate for MutY glycosylase in vivo is the A/GO mispair. We present biochemical evidence demonstrating that MutY glycosylase is an important part of a repair system that includes the MutM and MutT proteins. The GO repair system is dedicated to the repair of the oxidatively damaged guanine and the mutations it can induce.     

Aqueous Hydration of Nucleic Acid Constituents: Monte Carlo Computer Simulation Studies

Abstract

Monte Carlo computer simulations were performed on dilute aqueous solutions of thymine, cytosine, uracil, adenine, guanine, the dimethyl phosphate anion in the gauche-gauche conformation and a ribose and deoxyribose derivative. The aqueous hydration of each molecule was analysed in terms of quasi-component distribution functions based on the Proximity Criterion, and partitioned into hydrophobic, hydrophilic and ionic contributions. Color stereo views of selected hydration complexes are also presented. A preliminary discussion of the transferability of functional group coordination numbers is given. The results enable to comment on two current problems related to the hydration of nucleic acids: a) the theory of Dickerson and coworkers on the role of water in the relative stability of the A and B forms of DNA and b) the idea of water bridges and filaments emerging from the computer simulation results on the hydration of DNA fragments by Clementi.