Experimental and theoretical studies of the conformational perturbations induced by an abasic site

Examining Ty3 polypurine tract structure and function by nucleoside analog interference

We have combined nucleoside analog interference with chemical footprinting, thermal denaturation, NMR spectroscopy, and biochemical studies to understand recognition of the polypurine tract (PPT) primer of the Saccharomyces cerevisiae long terminal repeat-containing retrotransposon Ty3 by its cognate reverse transcriptase. Locked nucleic acid analogs, which constrain sugar ring geometry, were introduced pairwise throughout the PPT (-)-DNA template, whereas abasic tetrahydrofuran linkages, which lack the nucleobase but preserve the sugar phosphate backbone, were introduced throughout the (-)-strand DNA template and (+)-strand RNA primer. Collectively, our data suggest that both the 5'- and 3'-portions of the PPT-containing RNA/DNA hybrid are sensitive to nucleoside analog substitution, whereas the intervening region can be modified without altering cleavage specificity. These two regions most likely correspond to portions of the PPT that make close contact with the Ty3 reverse transcriptase thumb subdomain and RNase H catalytic center, respectively. Achieving a similar phenotype with nucleoside analogs that have different effects on duplex geometry reveals structural features that are important mediators of Ty3 PPT recognition. Finally, the results from introducing tetrahydrofuran lesions around the scissile PPT/unique 3'-sequence junction indicate that template nucleobase -1 is dispensable for catalysis, whereas a primer nucleobase on either side of the junction is necessary.

Conformational heterogeneity observed in simulations of a pyrene-substituted DNA

NMR studies previously carried out for a DNA system with a pyrene-substituted base did not observe NOEs involving the adenine located 5' to the pyrene, and thus the conformation of the adenine was poorly defined in the resulting family of refined structures. However, chemical shift data suggested that an AT base pair may be present. We have carried out fully unrestrained molecular dynamics simulations starting from several members of the family of structures, and these simulations support the existence of an AT base pair for this region. Simulations in both explicit and implicit solvent were carried out, with each converging to either anti or syn conformation for adenine and base pairing in all cases. During these simulations, large and dramatic conformational changes are observed that suggest pathways for complex conformational changes in the highly packed DNA interior. Our analysis reveals little difference in the energies of these syn and anti conformations, in contrast to control calculations carried out for standard DNA (in the absence of a neighboring pyrene). While no interconversion between the conformations was observed in standard simulations, reversible anti/syn exchange was directly simulated using the locally enhanced sampling approach. No exchange was seen in the non-pyrene control sequence. Together, these results suggest that an increased flexibility is introduced as a consequence of the pyrene substitution, offering an explanation that is consistent with the available experimental data. These results increase our optimism that simulations in atomic detail may provide accurate models for experimental observations in complex systems.

The major human abasic endonuclease: formation, consequences and repair of abasic lesions in DNA

DNA continuously suffers the loss of its constituent bases, and thereby, a loss of potentially vital genetic information. Sites of missing bases--termed abasic or apurinic/apyrimidinic (AP) sites--form spontaneously, through damage-induced hydrolytic base release, or by enzyme-catalyzed removal of modified or mismatched bases during base excision repair (BER). In this review, we discuss the structural and biological consequences of abasic lesions in DNA, as well as the multiple repair pathways for such damage, while emphasizing the mechanistic operation of the multi-functional human abasic endonuclease APE1 (or REF-1) and its potential relationship to disease.

Abasic DNA structure, reactivity, and recognition

Loss of a base in DNA, i.e., creation of an abasic site leaving a deoxyribose residue in the strand, is a frequent lesion that may occur spontaneously, or under the action of radiations and alkylating agents, or enzymatically as an intermediate in the repair of modified or abnormal bases. The abasic site lesion is mutagenic or lethal if not repaired. From a chemical point of view,the abasic site is an alkali-labile residue that leads to strand breakage through beta- and delta- elimination. Progress in the understanding of the chemistry and enzymology of abasic DNA largely relies upon the study of synthetic abasic duplexes. Several efficient synthetic methods have thus been developed to introduce the lesion (or a stable analogue) at defined position in the sequence. Physicochemical and spectroscopic examination of such duplexes, including calorimetry, melting temperature, high-field nmr and molecular modeling indicate that the lesion strongly destabilizes the duplex, although remaining in the canonical B-form with structural modifications strictly located at the site of the lesion. Probes have been developed to titrate the damage in DNA in vitro. Series of molecules have been devised to recognize specifically the abasic site, exhibiting a cleavage activity and mimicking the AP nucleases. Others have been prepared that bind strongly to the abasic site and show promise in potentiating the cytotoxic and antitumor activity of the clinically used nitrosourea (bis-chloroethylnitrosurea).

New insights into the structure of abasic DNA from molecular dynamics simulations

Abasic (AP) sites constitute a common form of DNA damage, arising from the spontaneous or enzymatic breakage of the N-glycosyl bond and the loss of a nucleotide base. To examine the effects of such damage on DNA structure, especially in the vicinity of the abasic sugar, four 1.5 ns molecular dynamics simulations of double-helical DNA dodecamers with and without a single abasic (tetrahydrofuran, X) lesion in a 5'-d(CXT) context have been performed and analyzed. The results indicate that the abasic site does not maintain a hole or gap in the DNA, but instead perturbs the canonical structure and induces additional flexibility close to the abasic site. In the apurinic simulations (i.e., when a pyrimidine is opposite the AP site), the abasic sugar flipped in and out of the minor groove, and the gap was water filled, except during the occurrence of a novel non-Watson-Crick C-T base pair across the abasic site. The apyrimidinic gap was not penetrated by water until the abasic sugar flipped out and remained extrahelical. Both AP helices showed kinks of 20-30 degrees at the abasic site. The Watson-Crick hydrogen bonds are more transient throughout the DNA double helices containing an abasic site. The abasic sugar displayed an unusually broad range of sugar puckers centered around the northern pucker. The increased motion of the bases and backbone near the abasic site appear to correlate with sequence-dependent helical stability. The data indicate that abasic DNA contorts more easily and in specific ways relative to unmodified DNA, an aspect likely to be important in abasic site recognition and hydrolysis.

Mapping the protein-DNA interface and the metal-binding site of the major human apurinic/apyrimidinic endonuclease

Apurinic/apyrimidinic (AP) endonuclease Ape1 is a key enzyme in the mammalian base excision repair pathway that corrects AP sites in the genome. Ape1 cleaves the phosphodiester bond immediately 5' to AP sites through a hydrolytic reaction involving a divalent metal co-factor. Here, site-directed mutagenesis, chemical footprinting techniques, and molecular dynamics simulations were employed to gain insights into how Ape1 interacts with its metal cation and AP DNA. It was found that Ape1 binds predominantly to the minor groove of AP DNA, and that residues R156 and Y128 contribute to protein-DNA complex stability. Furthermore, the Ape1-AP DNA footprint does not change along its reaction pathway upon active-site coordination of Mg(2+) or in the presence of DNA polymerase beta (polbeta), an interactive protein partner in AP site repair. The DNA region immediately 5' to the abasic residue was determined to be in close proximity to the Ape1 metal-binding site. Experimental evidence is provided that amino acid residues E96, D70, and D308 of Ape1 are involved in metal coordination. Molecular dynamics simulations, starting from the active site of the Ape1 crystal structure, suggest that D70 and E96 bind directly to the metal, while D308 coordinates the cation through the first hydration shell. These studies define the Ape1-AP DNA interface, determine the effect of polbeta on the Ape1-DNA interaction, and reveal new insights into the Ape1 active site and overall protein dynamics.

The impact of abasic sites on DNA flexibility

We use internal coordinate molecular mechanics calculations to study the impact of abasic sites on the conformation and the mechanics of the DNA double helix. Abasic sites, which are common mutagenic lesions, are shown to locally modify both the groove geometry and the curvature of DNA in a sequence dependent manner. By controlled twisting and bending, it is also shown that these lesions modify the deformability of the duplex, generally increasing its flexibility, but again to an extent which depends on the nature of the abasic site and on the surrounding base sequence. Both the conformational and mechanical influence of this type of DNA damage may be significant for recognition and repair mechanisms.

Abasic sites in duplex DNA: molecular modeling of sequence-dependent effects on conformation

Molecular modeling calculations using JUnction Minimization of Nucleic Acids (JUMNA) have been used to study sequence effects on the conformation of abasic sites within duplex DNA. We have considered lesions leading to all possible unpaired bases (X), adenine, guanine, cytosine, or thymine contained within two distinct sequence contexts, CXC and GXG. Calculations were carried out on DNA 11-mers using extensive conformational search techniques to locate the most stable abasic conformations and using Poisson-Boltzmann corrected electrostatics to account for solvation effects. The results, which are in very good agreement with available experimental data, point to strong sequence effects on both the position of the unpaired base (intra or extrahelical) and on the overall curvature induced by the abasic lesion. For CXC, unpaired purines are found to lie within the helix, while unpaired pyrimidines are either extrahelical or in equilibrium between the intra and extrahelical forms. For GXG, all unpaired bases lead to intrahelical forms, but with marked, sequence-dependent differences in induced curvature.