Molecular conformations of DNA and RNA subunits

Issued as Report of expenditures, Project no. G-41-64

A new anti conformation for N-(deoxyguanosin-8-yl)-2-acetylaminofluorene (AAF-dG) allows Watson–Crick pairing in the Sulfolobus solfataricus P2 DNA polymerase IV (Dpo4)

Primer extension studies have shown that the Y-family DNA polymerase IV (Dpo4) from Sulfolobus solfataricus P2 can preferentially insert C opposite N-(deoxyguanosin-8-yl)-2-acetylaminofluorene (AAF-dG) [F. Boudsocq, S. Iwai, F. Hanaoka and R. Woodgate (2001) Nucleic Acids Res., 29, 4607–4616]. Our goal is to elucidate on a structural level how AAF-dG can be harbored in the Dpo4 active site opposite an incoming dCTP, using molecular modeling and molecular dynamics simulations, since AAF-dG prefers the syn glycosidic torsion. Both anti and syn conformations of the templating AAF-dG in a Dpo4 ternary complex were investigated. All four dNTPs were studied. We found that an anti glycosidic torsion with C1′-exo deoxyribose conformation allows AAF-dG to be Watson–Crick hydrogen-bonded with dCTP with modest polymerase perturbation, but other nucleotides are more distorting. The AAF is situated in the Dpo4 major groove open pocket with fluorenyl rings 3′- and acetyl 5′-directed along the modified strand, irrespective of dNTP. With AAF-dG syn, the fluorenyl rings are in the small minor groove pocket and the active site region is highly distorted. The anti-AAF-dG conformation with C1′-exo sugar pucker can explain the preferential incorporation of dC by Dpo4. Possible relevance of our new major groove structure for AAF-dG to other polymerases, lesion repair and solution conformations are discussed

Influence of ribose 2′-O-methylation on GpC conformation by classical potential energy calculations

Potential energy calculations were employed to examine the effect of ribose 2′-O-methylation on the conformation of GpC. Minimum energy conformations and allowed conformational regions were calculated for 2′MeGpC and Gp2′MeC. The two lowest energy conformations of 2′MeGpC and Gp2′MeC are similar to those of GpC itself. The helical RNA conformation (sugar pucker-C(3′)-endo, ω′ and ω,g−g−, bases-anti) is the global minimum, and a helix-reversing conformation with ω′, ω in the vicinity of 20°, 80° is next in energy. However, subtle differences between the three molecules are noted. When the substitution is on the 5′ ribose (Gp2′MeC), the energy of the helical conformation is less than that of GpC, due to favorable interactions of the added methyl group. When the substitution is at the 3′ ribose (2′MeGpC) these stabilizing interactions are outweighed by steric restrictions, and the helical conformation is of higher energy than for GpC. Furthermore, the statistical weight of the 2′MeGpC g− g− helical region is substantially less than the corresponding weight for Gp2′MeC. In addition, 2′MeGpC′s methoxy group is conformationally restricted to a narrow range centered at 76°. This group has a broadly allowed region between 50 and 175° in Gp2′MeC. These differences occur because the appended methyl group in 2′MeGpC is located in the interior of the helix cylinder, as it would be in polynucleotide, while it hangs unimpeded in Gp2′MeC. These findings suggest that 2′-O-methylation has both stabilizing and destabilizing influences on the helical conformation of RNA. For 2′MeGpC the destabilizing steric hindrance imposed by the nature of the guanine base dominates

Conformational stability in dinucleoside phosphate crystals. Semiempirical potential energy calculations for uridylyl-3'-5'-adenosine monophosphate (UpA) and guanylyl-3',5'-cytidine monophosphate (GpC)

Classical potential energy calculations were performed for the dinucleoside phosphates UpA and GpC. Two widely accessible low-energy regions of conformation space were found for the w', w pair. That of lowest energy contains conformations similar to helical RNA, with w' and w in the vicinity of 300° and 280°, respectively. All five experimental observations of crystalline GpC, two of ApU, and the helical fragment of ApApA fall in this range. The second lowest region has w' and w at about 20° and 80°, respectively, which is in the general region of one experimentally observed crystalline conformer of UpA, and the nonhelical region of ApApA. It is concluded that GpC and ApU, which were crystallized as either sodium or calcium salts, are shielded from each other in the crystal by the water of hydration and are therefore free to adopt their predicted in vacuo minimum energy helical conformations. By contrast, crystalline UpA had only 1/2 water per molecule, and was forced into higher energy conformations in order to maximize intermolecular hydrogen bonding

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

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

Mutagenic nucleotide incorporation and hindered translocation by a food carcinogen C8-dG adduct in Sulfolobus solfataricus P2 DNA polymerase IV (Dpo4): modeling and dynamics studies

Bulky carcinogen-DNA adducts commonly cause replicative polymerases to stall, leading to a switch to bypass polymerases. We have investigated nucleotide incorporation opposite the major adduct of 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP) in the DinB family polymerase, Dpo4, using molecular modeling and molecular dynamics (MD) simulations. PhIP, the most prevalent heterocyclic aromatic amine formed by cooking of proteinaceous food, is mutagenic in mammalian cells and is implicated in mammary and colon tumors. Our results show that the dG-C8-PhIP adduct can be accommodated in the spacious major groove Dpo4 open pocket, with Dpo4 capable of incorporating dCTP, dTTP or dATP opposite the adduct reasonably well. However, the PhIP ring system on the minor groove side would seriously disturb the active site, regardless of the presence and identity of dNTP. Furthermore, the simulations indicate that dATP and dTTP are better incorporated in the damaged system than in their respective mismatched but unmodified controls, suggesting that the PhIP adduct enhances incorporation of these mismatches. Finally, bulky C8-dG adducts, situated in the major groove, are likely to impede translocation in this polymerase (Rechkoblit et al. (2006), PLoS Biol., 4, e11). However, N(2)-dG adducts, which can reside on the minor groove side, appear to cause less hindrance when in this position

Classical potential energy calculations for ApA, CpC, GpG, and UpU. The influence of the bases on RNA subunit conformations

Classical potential energy calculations have been made for the ribodinucleoside monophosphates ApA, CpC, GpG, and UpU. Van der Waal's, electrostatic, and torsional contributions to the energy were calculated, and the energy was minimized with the seven backbone conformational angles as simultaneously variable parameters. At the global minimum, ApA and CpC have conformations like double helical RNA: the angles ω′ and ω are g−g−, the sugar pucker is C3′-endo, and the bases are anti. GpG and UpU, on the other hand, have the ω′,ω angle pair g−t at the global minimum, and for GpG the bases are syn. Energy contour maps for ω′ and ω show two broad, low energy regions for ApA, CpC, and UpU: one is g−g−, and the second encompasses g−t and g+g+ within a single low energy contour. The two regions are connected by a path at 10–13 kcal./mole. For GpG, with bases syn, however, only a small low-energy region at g−t is found. The helical ‘A’ RNA conformation is 8.5 kcal/mole higher for this molecule. Thus, the base composition is shown to influence the conformations adopted by dinucleoside phosphates. Comparison of calculations with experimental data, where available, show good agreement