Sequence-dependent conformational perturbation in DNA duplexes containing an epsilonA.T mismatch using molecular dynamics simulation

In What Ways Do Synthetic Nucleotides and Natural Base Lesions Alter the Structural Stability of G-Quadruplex Nucleic Acids?

Synthetic analogs of natural nucleotides have long been utilized for structural studies of canonical and noncanonical nucleic acids, including the extensively investigated polymorphic G-quadruplexes (GQs). Dependence on the sequence and nucleotide modifications of the folding landscape of GQs has been reviewed by several recent studies. Here, an overview is compiled on the thermodynamic stability of the modified GQ folds and on how the stereochemical preferences of more than 70 synthetic and natural derivatives of nucleotides substituting for natural ones determine the stability as well as the conformation. Groups of nucleotide analogs only stabilize or only destabilize the GQ, while the majority of analogs alter the GQ stability in both ways. This depends on the preferred syn or anti N-glycosidic linkage of the modified building blocks, the position of substitution, and the folding architecture of the native GQ. Natural base lesions and epigenetic modifications of GQs explored so far also stabilize or destabilize the GQ assemblies. Learning the effect of synthetic nucleotide analogs on the stability of GQs can assist in engineering a required stable GQ topology, and exploring the in vitro action of the single and clustered natural base damage on GQ architectures may provide indications for the cellular events.

Structural studies of malonaldehyde–glyoxal and malonaldehyde–methylglyoxal etheno adducts of adenine nucleosides based on spectroscopic methods and DFT-GIAO calculations

The substitution position in the etheno rings of M₁Gx-A and M₁MGx-dA was determined based on the comparison of data derived from NMR spectra with results obtained from computational calculations.

Helical structure of xylose-DNA

Synthetic biology and systems chemistry demonstrate a growing interest in modified nucleotides to achieve an enzymatically stable artificial nucleic acid. A potential candidate system is xylose-DNA, in which the 2'-deoxy-beta-D-ribo-furanose is substituted by 2'-deoxy-beta-D-xylo-furanose. We present here the helical structure and conformational analysis of xylose-DNA on the basis of 35 ns MD simulations of a 29-base-pair DNA duplex. Starting from a right-handed xylose-DNA helix, we observe a remarkable conformational transition from right- to left-handed helix. The left-handed xylose-DNA is highly dynamic, involving screwing and unscrewing motion of the helix. The sugar pucker induced helical changes influence the backbone to adopt the backbone angles for xylose-DNA while retaining the Watson-Crick base pairing and stacking interactions. The results demonstrate the chiral orthogonality of the ribose and xylose based episomes. As far as stability and compactness of information storage is concerned, the ribose based natural DNA is unsurpassed.

Substrate specificity and sequence-dependent activity of the Saccharomyces cerevisiae 3-methyladenine DNA glycosylase (Mag)

DNA glycosylases initiate base excision repair by first binding, then excising aberrant DNA bases. Saccharomyces cerevisiae encodes a 3-methyladenine (3MeA) DNA glycosylase, Mag, that recognizes 3MeA and various other DNA lesions including 1,N6-ethenoadenine (epsilon A), hypoxanthine (Hx) and abasic (AP) sites. In the present study, we explore the relative substrate specificity of Mag for these lesions and in addition, show that Mag also recognizes cisplatin cross-linked adducts, but does not catalyze their excision. Through competition binding and activity studies, we show that in the context of a random DNA sequence Mag binds epsilon A and AP-sites the most tightly, followed by the cross-linked 1,2-d(ApG) cisplatin adduct. While epsilon A binding and excision by Mag was robust in this sequence context, binding and excision of Hx was extremely poor. We further studied the recognition of epsilon A and Hx by Mag, when these lesions are present at different positions within A:T and G:C tracts. Overall, epsilon A was slightly less well excised from each position within the A:T and G:C tracts compared to excision from the random sequence, whereas Hx excision was greatly increased in these sequence contexts (by up to 7-fold) compared to the random sequence. However, given most sequence contexts, Mag had a clear preference for epsilon A relative to Hx, except in the TTXTT (X=epsilon A or Hx) sequence context from which Mag removed both lesions with almost equal efficiency. We discuss how DNA sequence context affects base excision by various 3MeA DNA glycosylases.

Alkylpurine-DNA-N-glycosylase excision of 7-(hydroxymethyl)-1,N6-ethenoadenine, a glycidaldehyde-derived DNA adduct

Glycidaldehyde (GDA) is a bifunctional alkylating agent that has been shown to be mutagenic in vitro and carcinogenic in rodents. However, the molecular mechanism by which it exerts these effects is not established. GDA is capable of forming exocyclic hydroxymethyl-substituted etheno adducts on base residues in vitro. One of them, 7-(hydroxymethyl)-1,N6-ethenoadenine (7-hm-epsilonA), was identified as the principal adduct in mouse skin treated with GDA or a glycidyl ether. In this work, using defined oligonucleotides containing a site-specific 7-hm-epsilonA, the human and mouse alkylpurine-DNA-N-glycosylases (APNGs), responsible for the removal of the analogous 1,N6-ethenoadenine (epsilonA) adduct, are shown to recognize and excise 7-hm-epsilonA. Such an activity can be significantly modulated by both 5' neighboring and opposite sequence contexts. The efficiency of human or mouse APNG for excision of 7-hm-epsilonA is about half that, or similar to the excision of epsilonA, respectively. When human or mouse cell-free extracts were tested, however, the extent of 7-hm-epsilonA excision is dramatically lower than that for epsilonA, suggesting that, in the crude extracts, the APNG activities toward these two adducts are differentially affected. Using cell-free extracts from APNG deficient mice, this enzyme is shown to be the primary glycosylase excising 7-hm-epsilonA. A structural approach, using molecular modeling, was employed to examine how the structure of the 7-hm-epsilonA adduct affects DNA conformation, as compared to the epsilonA adduct. These novel substrate specificities could have both biological and structural implications.

The p-benzoquinone DNA adducts derived from benzene are highly mutagenic

Benzene is a human leukemia carcinogen, resulting from its cellular metabolism. A major benzene metabolite is p-benzoquinone (pBQ), which can damage DNA by forming the exocyclic base adducts pBQ-dC, pBQ-dA, and pBQ-dG in vitro. To gain insights into the role of pBQ in benzene genotoxicity, we examined in vitro translesion synthesis and in vivo mutagenesis of these pBQ adducts. Purified REV1 and Polkappa were essentially incapable of translesion synthesis in response to the pBQ adducts. Opposite pBQ-dA and pBQ-dC, purified human Poliota was capable of error-prone nucleotide insertion, but was unable to perform extension synthesis. Error-prone translesion synthesis was observed with Poleta. However, DNA synthesis largely stopped opposite the lesion. Consistent with in vitro results, replication of site-specifically damaged plasmids was strongly inhibited by pBQ adducts in yeast cells, which depended on both Polzeta and Poleta. In wild-type cells, the majority of translesion products were deletions at the site of damage, accounting for 91%, 90%, and 76% for pBQ-dA, pBQ-dG, and pBQ-dC, respectively. These results show that the pBQ-dC, pBQ-dA, and pBQ-dG adducts are strong blocking lesions, and are highly mutagenic by predominantly inducing deletion mutations. These results are consistent with the lesion structures predicted by molecular dynamics simulation. Our results led to the following model. Translesion synthesis normally occurs by directly copying the lesion site through base insertion and extension synthesis. When the lesion becomes incompatible in accommodating a base opposite the lesion in DNA, translesion synthesis occurs by a less efficient lesion loop-out mechanism, resulting in avoiding copying the damaged base and leading to deletion.

Chloroethylnitrosourea-derived ethano cytosine and adenine adducts are substrates for Escherichia coli glycosylases excising analogous etheno adducts

Exocyclic ethano DNA adducts are saturated etheno ring derivatives formed mainly by therapeutic chloroethylnitrosoureas (CNUs), which are also mutagenic and carcinogenic. In this work, we report that two of the ethano adducts, 3,N4-ethanocytosine (EC) and 1,N6-ethanoadenine (EA), are novel substrates for the Escherichia coli mismatch-specific uracil-DNA glycosylase (Mug) and 3-methyladenine DNA glycosylase II (AlkA), respectively. It has been shown previously that Mug excises 3,N4-ethenocytosine (epsilonC) and AlkA releases 1,N6-ethenoadenine (epsilonA). Using synthetic oligonucleotides containing a single ethano or etheno adduct, we found that both glycosylases had a approximately 20-fold lower excision activity toward EC or EA than that toward their structurally analogous epsilonC or epsilonA adduct. Both enzymes were capable of excising the ethano base paired with any of the four natural bases, but with varying efficiencies. The Mug activity toward EC could be stimulated by E. coli endonuclease IV and, more efficiently, by exonuclease III. Molecular dynamics (MD) simulations showed similar structural features of the etheno and ethano derivatives when present in DNA duplexes. However, also as shown by MD, the stacking interaction between the EC base and Phe 30 in the Mug active site is reduced as compared to the epsilonC base, which could account for the lower EC activity observed in this study.

Structural insights by molecular dynamics simulations into specificity of the major human AP endonuclease toward the benzene-derived DNA adduct, pBQ-C

The benzetheno exocyclic adduct of the cytosine (C) base (pBQ-C) is a product of reaction between DNA and a stable metabolite of the human carcinogen benzene, p-benzoquinone (pBQ). We reported previously that the pBQ-C-containing duplex is a substrate for the human AP endonuclease (APE1), an enzyme that cleaves an apurinic/apyrimidinic (AP) site from double stranded DNA. In this work, using molecular dynamics simulation (MD), we provided a structural explanation for the recognition of the pBQ-C adduct by APE1. Molecular modeling of the DNA duplex containing pBQ-C revealed significant displacement of this adduct toward the major groove with pronounced kinking of the DNA at the lesion site, which could serve as a structural element recognized by the APE1 enzyme. Using 3 ns MD it was shown that the position of the pBQ-C adduct is stabilized by two hydrogen bonds formed between the adduct and the active site amino acids Asp 189 and Ala 175. The pBQ-C/APE1 complex, generated by MD, has a similar hydrogen bond network between target phosphodiester bond at the pBQ-C site and key amino acids at the active site, as in the crystallographically determined APE1 complexed with an AP site-containing DNA duplex. The position of the adduct at the enzyme active site, together with the hydrogen bond network, suggests a similar reaction mechanism for phosphodiester bond cleavage of oligonucleotide containing pBQ-C as reported for the AP site.

Structural insights by molecular dynamics simulations into differential repair efficiency for ethano-A versus etheno-A adducts by the human alkylpurine-DNA N-glycosylase

1,N6-ethenoadenine adducts (epsilonA) are formed by known environmental carcinogens and found to be removed by human alkylpurine-DNA N-glycosylase (APNG). 1,N6-ethanoadenine (EA) adducts differ from epsilonA by change of a double bond to a single bond in the 5-member exocyclic ring and are formed by chloroethyl nitrosoureas, which are used in cancer therapy. In this work, using purified recombinant human APNG, we show that EA is a substrate for the enzyme. However, the excision efficiency of EA was 65-fold lower than that of epsilonA. Molecular dynamics simulation produced similar structural motifs for epsilonA and EA when incorporated into a DNA duplex, suggesting that there are no specific conformational features in the DNA duplex which can account for the differences in repair efficiency. However, when EA was modeled into the APNG active site, based on the APNG/epsilonA-DNA crystallographic coordinates, in structures produced by 2 ns molecular dynamics simulation, we observed weakening in the stacking interaction between EA and aromatic side chains of the key amino acids in the active site. In contrast, the planar epsilonA is better stacked at the enzyme active site. We propose that the observed destabilization of the EA adduct at the active site, such as reduced stacking interactions, could account for the biochemically observed weaker recognition of EA by APNG as compared to epsilonA.

Molecular dynamics simulation of nucleic acids: successes, limitations, and promise

In the last five years we have witnessed a significant increase in the number publications describing accurate and reliable all-atom molecular dynamics simulations of nucleic acids. This increase has been facilitated by the development of fast and efficient methods for treating the long-range electrostatic interactions, the availability of faster parallel computers, and the development of well-validated empirical molecular mechanical force fields. With these technologies, it has been demonstrated that simulation is not only capable of consistently reproducing experimental observations of sequence specific fine structure of DNA, but also can give detailed insight into prevalent problems in nucleic acid structure, ion association and specific hydration of nucleic acids, polyadenine tract bending, and the subtle environmental dependence of the A-DNA-B-DNA duplex equilibrium. Despite the advances, there are still issues with the methods that need to be resolved through rigorous controlled testing. In general, these relate to deficiencies of the underlying molecular mechanical potentials or applied methods (such as the imposition of true periodicity in Ewald simulations and the need for energy conservation), and significant limits in effective conformational sampling. In this perspective, we provide an overview of our experiences, provide some cautionary notes, and provide recommendations for further study in molecular dynamics simulation of nucleic acids.