Structure of a Stable Interstrand DNA Cross-Link Involving a β-N-Glycosyl Linkage Between an N6-dA Amino Group and an Abasic Site

NEIL3: A unique DNA glycosylase involved in interstrand DNA crosslink repair

Endonuclease VIII-like 3 (NEIL3) is a versatile DNA glycosylase that repairs a diverse array of chemical modifications to DNA. Unlike other glycosylases, NEIL3 has a preference for lesions within single-strand DNA and at single/double-strand DNA junctions. Beyond its canonical role in base excision repair of oxidized DNA, NEIL3 initiates replication-dependent interstrand DNA crosslink repair as an alternative to the Fanconi Anemia pathway. This review outlines our current understanding of NEIL3's biological functions, role in disease, and three-dimensional structure as it pertains to substrate specificity and catalytic mechanism.

Stable Interstrand Cross-Links Generated from the Repair of 1,N6-Ethenoadenine in DNA by α-Ketoglutarate/Fe(II)-Dependent Dioxygenase ALKBH2

DNA cross-links severely challenge replication and transcription in cells, promoting senescence and cell death. In this paper, we report a novel type of DNA interstrand cross-link (ICL) produced as a side product during the attempted repair of 1,N6-ethenoadenine (εA) by human α-ketoglutarate/Fe(II)-dependent enzyme ALKBH2. This stable/nonreversible ICL was characterized by denaturing polyacrylamide gel electrophoresis analysis and quantified by high-resolution LC-MS in well-matched and mismatched DNA duplexes, yielding 5.7% as the highest level for cross-link formation. The binary lesion is proposed to be generated through covalent bond formation between the epoxide intermediate of εA repair and the exocyclic N6-amino group of adenine or the N4-amino group of cytosine residues in the complementary strand under physiological conditions. The cross-links occur in diverse sequence contexts, and molecular dynamics simulations rationalize the context specificity of cross-link formation. In addition, the cross-link generated from attempted εA repair was detected in cells by highly sensitive LC-MS techniques, giving biological relevance to the cross-link adducts. Overall, a combination of biochemical, computational, and mass spectrometric methods was used to discover and characterize this new type of stable cross-link both in vitro and in human cells, thereby uniquely demonstrating the existence of a potentially harmful ICL during DNA repair by human ALKBH2.

Interaction of mitoxantrone with abasic sites - DNA strand cleavage and inhibition of apurinic/apyrimidinic endonuclease 1, APE1

Mitoxantrone (1,4-dihydroxy-5,8-bis[2-(2-hydroxyethylamino)ethylamino]-anthracene-9,10-dione) is a clinically-relevant synthetic anthracenedione that functions as a topoisomerase II poison by trapping DNA double-strand break intermediates. Mitoxantrone binds to DNA via both stacking interactions with DNA bases and hydrogen bonding with the sugar-phosphate backbone. It has been shown that mitoxantrone inhibits apurinic/apyrimidinic (AP) endonuclease 1 (APE1)-catalyzed incision of DNA containing a tetrahydrofuran (THF) moiety and more recently, that mitoxantrone forms Schiff base conjugates at AP sites in DNA. In this study, mitoxantrone-mediated inhibition of APE1 at THF sites was shown to be consistent with preferential binding to, and thermal stabilization of DNA containing a THF site as compared to non-damaged DNA. Investigations into the properties of mitoxantrone at AP and 3' α,β-unsaturated aldehyde sites demonstrated that in addition to being a potent inhibitor of APE1 at these biologically-relevant substrates (∼ 0.5 μM IC50 on AP site-containing DNA), mitoxantrone also incised AP site-containing DNA by catalyzing β- and β/δ-elimination reactions. The efficiency of these reactions to generate the 3' α,β-unsaturated aldehyde and 3' phosphate products was modulated by DNA structure. Although these cell-free reactions revealed that mitoxantrone can generate 3' phosphates, cells lacking polynucleotide kinase phosphatase did not show increased sensitivity to mitoxantrone treatment. Consistent with its ability to inhibit APE1 activity on DNAs containing either an AP site or a 3' α,β-unsaturated aldehyde, combined exposures to clinically-relevant concentrations of mitoxantrone and a small molecule APE1 inhibitor revealed additive cytotoxicity. These data suggest that in a cellular context, mitoxantrone may interfere with APE1 DNA repair functions.

Unprecedented reactivity of polyamines with aldehydic DNA modifications: structural determinants of reactivity, characterization and enzymatic stability of adducts

Apurinic/apyrimidinic (AP) sites, 5-formyluracil (fU) and 5-formylcytosine (fC) are abundant DNA modifications that share aldehyde-type reactivity. Here, we demonstrate that polyamines featuring at least one secondary 1,2-diamine fragment in combination with aromatic units form covalent DNA adducts upon reaction with AP sites (with concomitant cleavage of the AP strand), fU and, to a lesser extent, fC residues. Using small-molecule mimics of AP site and fU, we show that reaction of secondary 1,2-diamines with AP sites leads to the formation of unprecedented 3'-tetrahydrofuro[2,3,4-ef]-1,4-diazepane ('ribodiazepane') scaffold, whereas the reaction with fU produces cationic 2,3-dihydro-1,4-diazepinium adducts via uracil ring opening. The reactivity of polyamines towards AP sites versus fU and fC can be tuned by modulating their chemical structure and pH of the reaction medium, enabling up to 20-fold chemoselectivity for AP sites with respect to fU and fC. This reaction is efficient in near-physiological conditions at low-micromolar concentration of polyamines and tolerant to the presence of a large excess of unmodified DNA. Remarkably, 3'-ribodiazepane adducts are chemically stable and resistant to the action of apurinic/apyrimidinic endonuclease 1 (APE1) and tyrosyl-DNA phosphoesterase 1 (TDP1), two DNA repair enzymes known to cleanse a variety of 3' end-blocking DNA lesions.

MIM-ML: A Novel Quantum Chemical Fragment-Based Random Forest Model for Accurate Prediction of NMR Chemical Shifts of Nucleic Acids

We developed a random forest machine learning (ML) model for the prediction of 1H and 13C NMR chemical shifts of nucleic acids. Our ML model is trained entirely on reproducing computed chemical shifts obtained previously on 10 nucleic acids using a Molecules-in-Molecules (MIM) fragment-based density functional theory (DFT) protocol including microsolvation effects. Our ML model includes structural descriptors as well as electronic descriptors from an inexpensive low-level semiempirical calculation (GFN2-xTB) and trained on a relatively small number of DFT chemical shifts (2080 1H chemical shifts and 1780 13C chemical shifts on the 10 nucleic acids). The ML model is then used to make chemical shift predictions on 8 new nucleic acids ranging in size from 600 to 900 atoms and compared directly to experimental data. Though no experimental data was used in the training, the performance of our model is excellent (mean absolute deviation of 0.34 ppm for 1H chemical shifts and 2.52 ppm for 13C chemical shifts for the test set), despite having some nonstandard structures. A simple analysis suggests that both structural and electronic descriptors are critical for achieving reliable predictions. This is the first attempt to combine ML from fragment-based DFT calculations to predict experimental chemical shifts accurately, making the MIM-ML model a valuable tool for NMR predictions of nucleic acids.

Accurate and Cost-Effective NMR Chemical Shift Predictions for Nucleic Acids Using a Molecules-in-Molecules Fragmentation-Based Method

We have developed, implemented, and assessed an efficient protocol for the prediction of NMR chemical shifts of large nucleic acids using our molecules-in-molecules (MIM) fragment-based quantum chemical approach. To assess the performance of our approach, MIM-NMR calculations are calibrated on a test set of three nucleic acids, where the structure is derived from solution-phase NMR studies. For DNA systems with multiple conformers, the one-layer MIM method with trimer fragments (MIM1trimer) is benchmarked to get the lowest energy structure, with an average error of only 0.80 kcal/mol with respect to unfragmented full molecule calculations. The MIMI-NMRdimer calibration with respect to unfragmented full molecule calculations shows a mean absolute deviation (MAD) of 0.06 and 0.11 ppm, respectively, for 1H and 13C nuclei, but the performance with respect to experimental NMR chemical shifts is comparable to the more expensive MIM1-NMR and MIM2-NMR methods with trimer subsystems. To compare with the experimental chemical shifts, a standard protocol is derived using DNA systems with Protein Data Bank (PDB) IDs 1SY8, 1K2K, and 1KR8. The effect of structural minimizations is employed using a hybrid mechanics/semiempirical approach and used for computations in solution with implicit and explicit-implicit solvation models in our MIM1-NMRdimer methodology. To demonstrate the applicability of our protocol, we tested it on seven nucleic acids, including structures with nonstandard residues, heteroatom substitutions (F and B atoms), and side chain mutations with a size ranging from ∼300 to 1100 atoms. The major improvement for predicted MIM1-NMRdimer calculations is obtained from structural minimizations and implicit solvation effects. A significant improvement with the explicit-implicit solvation model is observed only for two smaller nucleic acid systems (1KR8 and 7NBK), where the expensive first solvation shell is replaced by the microsolvation model, in which a single water molecule is added for each solvent-exposed amino and imino protons, along with the implicit solvation. Overall, our target accuracy of ∼0.2-0.3 ppm for 1H and ∼2-3 ppm for 13C has been achieved for large nucleic acids. The proposed MIM-NMR approach is accurate and cost-effective (linear scaling with system size), and it can aid in the structural assignments of a wide range of complex biomolecules.

Effects of Local Sequence, Reaction Conditions, and Various Additives on the Formation and Stability of Interstrand Cross-Links Derived from the Reaction of an Abasic Site with an Adenine Residue in Duplex DNA

The experiments described here examined the effects of reaction conditions, various additives, and local sequence on the formation and stability interstrand cross-links (ICLs) derived from the reaction of an apurinic/apyrimidinic (AP) site with the exocyclic amino group of an adenine residue on the opposing strand in duplex DNA. Cross-link formation was observed in a range of different buffers, with faster formation rates observed at pH 5. Inclusion of the base excision repair enzyme alkyladenine DNA glycosylase (hAAG) which binds tightly to AP-containing duplexes decreased, but did not completely prevent, formation of the dA-AP ICL. Formation of the dA-AP ICL was not altered by the presence of the biological metal ion Mg²⁺ or the biological thiol, glutathione. Several organocatalysts of imine formation did not enhance the rate of dA-AP ICL formation. Duplex length did not have a large effect on dA-AP yield, so long as the melting temperature of the duplex was not significantly below the reaction temperature (the duplex must remain hybridized for efficient ICL formation). Formation of the dA-AP ICL was examined in over 40 different sequences that varied the neighboring and opposing bases at the cross-linking site. The results indicate that ICL formation can occur in a wide variety of sequence contexts under physiological conditions. Formation of the dA-AP ICL was strongly inhibited by the aldehyde-trapping agents methoxyamine and hydralazine, by NaBH₃CN, by the intercalator ethidium bromide, and by the minor groove-binding agent netropsin. ICL formation was inhibited to some extent in bicarbonate and Tris buffers. The dA-AP ICL showed substantial inherent stability under a variety of conditions and was not a substrate for AP-processing enzymes APE1 or Endo IV. Finally, we characterized cross-link formation in a small (11 bp) stem-loop (hairpin) structure and in DNA-RNA hybrid duplexes.

Establishing Linkages Among DNA Damage, Mutagenesis, and Genetic Diseases

DNA damage by chemicals, radiation, or oxidative stress leads to a mutational spectrum, which is complex because it is determined in part by lesion structure, the DNA sequence context of the lesion, lesion repair kinetics, and the type of cells in which the lesion is replicated. Accumulation of mutations may give rise to genetic diseases such as cancer and therefore understanding the process underlying mutagenesis is of immense importance to preserve human health. Chemical or physical agents that cause cancer often leave their mutational fingerprints, which can be used to back-calculate the molecular events that led to disease. To make a clear link between DNA lesion structure and the mutations a given lesion induces, the field of single-lesion mutagenesis was developed. In the last three decades this area of research has seen much growth in several directions, which we attempt to describe in this Perspective.

Model of abasic site DNA cross-link repair; from the architecture of NEIL3 DNA binding domains to the X-structure model

Covalent DNA interstrand crosslinks are toxic DNA damage lesions that block the replication machinery that can cause a genomic instability. Ubiquitous abasic DNA sites are particularly susceptible to spontaneous cross-linking with a base from the opposite DNA strand. Detection of a crosslink induces the DNA helicase ubiquitination that recruits NEIL3, a DNA glycosylase responsible for the lesion removal. NEIL3 utilizes several zinc finger domains indispensable for its catalytic NEI domain repairing activity. They recruit NEIL3 to the repair site and bind the single-stranded DNA. However, the molecular mechanism underlying their roles in the repair process is unknown. Here, we report the structure of the tandem zinc-finger GRF domain of NEIL3 and reveal the molecular details of its interaction with DNA. Our biochemical data indicate the preferential binding of the GRF domain to the replication fork. In addition, we obtained a structure for the catalytic NEI domain in complex with the DNA reaction intermediate that allowed us to construct and validate a model for the interplay between the NEI and GRF domains in the recognition of an interstrand cross-link. Our results suggest a mechanism for recognition of the DNA replication X-structure by NEIL3, a key step in the interstrand cross-link repair.

Hybridization-specific chemical reactions to create interstrand crosslinking and threaded structures of nucleic acids

The interstrand crosslinking and threaded structures of nucleic acids have high potential in oligonucleotide therapeutics, chemical biology, and nanotechnology. For example, properly designed crosslinking structures provide high activity and nuclease resistance for anti-miRNAs. The noncovalent labeling and modification by the threaded structures are useful as new chemical biology tools. Photoreversible crosslinking creates smart materials, such as reversible photoresponsive gels and DNA origami objects. This review introduces the creation of interstrand crosslinking and threaded structures, such as catenanes and rotaxanes, based on hybridization-specific chemical reactions and their functions and perspectives.

The rate of formation and stability of abasic site interstrand crosslinks in the DNA duplex

DNA interstrand crosslinks (ICLs) strands pose an impenetrable barrier for DNA replication. Different ICLs are known to recruit distinct DNA repair pathways. NEIL3 glycosylase has been known to remove an abasic (Ap) site derived DNA crosslink (Ap-ICL). An Ap-ICL forms spontaneously from the Ap site with an adjacent adenine in the opposite strand. Lack of genetic models and a poor understanding of the fate of these lesions leads to many questions about the occurrence and the toxicity of Ap-ICL in cells. Here, we investigate the circumstances of Ap-ICL formation. With an array of different oligos, we have investigated the rates of formation, the yields, and the stability of Ap-ICL. Our findings point out how different bases in the vicinity of the Ap site change crosslink formation in vitro. We reveal that AT-rich rather than GC-rich regions in the surrounding Ap site lead to higher rates of Ap-ICL formation. Overall, our data reveal that Ap-ICL can be formed in virtually any DNA sequence context surrounding a hot spot of a 5'-Ap-dT pair, albeit with significantly different rates and yields. Based on Ap-ICL formation in vitro, we attempt to predict the number of Ap-ICLs in the cell.