Looking for Waldo: a potential thermodynamic signature to DNA damage

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.

Impact of 3-deazapurine nucleobases on RNA properties

Deazapurine nucleosides such as 3-deazaadenosine (c3A) are crucial for atomic mutagenesis studies of functional RNAs. They were the key for our current mechanistic understanding of ribosomal peptide bond formation and of phosphodiester cleavage in recently discovered small ribozymes, such as twister and pistol RNAs. Here, we present a comprehensive study on the impact of c3A and the thus far underinvestigated 3-deazaguanosine (c3G) on RNA properties. We found that these nucleosides can decrease thermodynamic stability of base pairing to a significant extent. The effects are much more pronounced for 3-deazapurine nucleosides compared to their constitutional isomers of 7-deazapurine nucleosides (c7G, c7A). We furthermore investigated base pair opening dynamics by solution NMR spectroscopy and revealed significantly enhanced imino proton exchange rates. Additionally, we solved the X-ray structure of a c3A-modified RNA and visualized the hydration pattern of the minor groove. Importantly, the characteristic water molecule that is hydrogen-bonded to the purine N3 atom and always observed in a natural double helix is lacking in the 3-deazapurine-modified counterpart. Both, the findings by NMR and X-ray crystallographic methods hence provide a rationale for the reduced pairing strength. Taken together, our comparative study is a first major step towards a comprehensive understanding of this important class of nucleoside modifications.

Toxicity and repair of DNA adducts produced by the natural product yatakemycin

Yatakemycin (YTM) is an extraordinarily toxic DNA alkylating agent with potent antimicrobial and antitumor properties and the most recent addition to the CC-1065 and duocarmycin family of natural products. While bulky DNA lesions the size of those produced by YTM are normally removed from the genome by the nucleotide excision repair (NER) pathway, YTM adducts are also a substrate for the bacterial DNA glycosylases AlkD and YtkR2, unexpectedly implicating base excision repair (BER) in their elimination. The reason for the extreme toxicity of these lesions and the molecular basis for how they are eliminated by BER have been unclear. Here, we describe the structural and biochemical properties of YTM adducts responsible for their toxicity, and define the mechanism by which they are excised by AlkD. These findings delineate an alternative strategy for repair of bulky DNA damage and establish the cellular utility of this pathway relative to that of NER.

Probing the Dynamic Interaction between Damaged DNA and a Cellular Responsive Protein Using a Piezoelectric Mass Biosensor

The binding events between damaged DNA and recognition biomolecules are of great interest for understanding the activity of DNA-damaging drugs and the related DNA repair networks. Herein, a simple and sensitive sensor system was tailored for real-time probing of the dynamic molecular recognition between cisplatin-damaged-DNA (cisPt-DNA) and a cellular responsive protein, high-mobility-group box 1 (HMGB1). By integration of flow injection analysis (FIA) with quartz crystal microbalance (QCM), the interaction time-course of cisPt-DNA and HMGB1 domain A (HMGB1a) was investigated. The highly specific sensing interface was carefully designed and fabricated using cisPt-DNA as recognition element. A hybrid self-assembled monolayer consisting of cysteamine and mercaptohexanol was introduced to resist nonspecific adsorption. The calculated kinetic parameters (k_ass and k_diss) and the dissociation constant (K_D) demonstrated the rapid recognition and tight binding of HMGB1a toward cisPt-DNA. Molecular docking was employed to simulate the complex formed by cisPt-DNA and HMGB1a. The tight binding of such a DNA-damage responsive complex is appealing for the downstream molecular recognition event related to the resistance to DNA repair. This continuous-flow QCM biosensor is an ideal tool for studying specific interactions between drug-damaged-DNAs and their recognition proteins in a physiological-relevant environment, and will provide a potential sensor platform for rapid screening and evaluating metal anticancer drugs.

Using Wavelet Analysis To Assist in Identification of Significant Events in Molecular Dynamics Simulations

Long time scale molecular dynamics (MD) simulations of biological systems are becoming increasingly commonplace due to the availability of both large-scale computational resources and significant advances in the underlying simulation methodologies. Therefore, it is useful to investigate and develop data mining and analysis techniques to quickly and efficiently extract the biologically relevant information from the incredible amount of generated data. Wavelet analysis (WA) is a technique that can quickly reveal significant motions during an MD simulation. Here, the application of WA on well-converged long time scale (tens of μs) simulations of a DNA helix is described. We show how WA combined with a simple clustering method can be used to identify both the physical and temporal locations of events with significant motion in MD trajectories. We also show that WA can not only distinguish and quantify the locations and time scales of significant motions, but by changing the maximum time scale of WA a more complete characterization of these motions can be obtained. This allows motions of different time scales to be identified or ignored as desired.

Molecular Dynamics Simulations of 5-Hydroxycytosine Damaged DNA

Oxidation of cytosine is a leading cause of mutations and can lead to cancer. Here we report molecular dynamics simulations that characterized the structure and flexibility of 5-hydroxycytosine damaged DNA. A total of four systems were studied: undamaged DNA, damaged DNA base paired to a matching guanine, damaged DNA base paired to a mismatching adenine, and the corresponding undamaged mismatched strand. The simulations showed high spatial similarity between undamaged and damaged DNA; however, the matched damaged strand had greater overtwisting flexibility, and for both the matched and unmatched strands sugar puckering was much more flexible at the damaged site. The mismatch introduced larger changes, notably a loss in hydrogen bonding and a gain in stacking interactions, as well as effects on base pair and step geometry and solvation. Implications for damage recognition are discussed.

Hydrogen abstraction by photoexcited benzophenone: consequences for DNA photosensitization

We report a computational investigation of the hydrogen abstraction (H-abstraction) induced by triplet benzophenone (³BP) on thymine nucleobase and backbone sugar.

Recognition of Damaged DNA for Nucleotide Excision Repair: A Correlated Motion Mechanism with a Mismatched cis-syn Thymine Dimer Lesion

Mammalian global genomic nucleotide excision repair requires lesion recognition by XPC, whose detailed binding mechanism remains to be elucidated. Here we have delineated the dynamic molecular pathway and energetics of lesion-specific and productive binding by the Rad4/yeast XPC lesion recognition factor, as it forms the open complex [Min, J. H., and Pavletich, N. P. (2007) Nature 449, 570–575; Chen, X., et al. (2015) Nat. Commun. 6, 5849] that is required for excision. We investigated extensively a cis-syn cyclobutane pyrimidine dimer in mismatched duplex DNA, using high-level computational approaches. Our results delineate a preferred correlated motion mechanism, which provides for the first time an atomistic description of the sequence of events as Rad4 productively binds to the damaged DNA.

Understanding DNA under oxidative stress and sensitization: the role of molecular modeling

DNA is constantly exposed to damaging threats coming from oxidative stress, i.e., from the presence of free radicals and reactive oxygen species. Sensitization from exogenous and endogenous compounds that strongly enhance the frequency of light-induced lesions also plays an important role. The experimental determination of DNA lesions, though a difficult subject, is somehow well established and allows to elucidate even extremely rare DNA lesions. In parallel, molecular modeling has become fundamental to clearly understand the fine mechanisms related to DNA defects induction. Indeed, it offers an unprecedented possibility to get access to an atomistic or even electronic resolution. Ab initio molecular dynamics may also describe the time-evolution of the molecular system and its reactivity. Yet the modeling of DNA (photo-)reactions does necessitate elaborate multi-scale methodologies to tackle a damage induction reactivity that takes place in a complex environment. The double-stranded DNA environment is first characterized by a very high flexibility, but also a strongly inhomogeneous electrostatic embedding. Additionally, one aims at capturing more subtle effects, such as the sequence selectivity which is of critical important for DNA damage. The structure and dynamics of the DNA/sensitizers complexes, as well as the photo-induced electron- and energy-transfer phenomena taking place upon sensitization, should be carefully modeled. Finally the factors inducing different repair ratios for different lesions should also be rationalized. In this review we will critically analyze the different computational strategies used to model DNA lesions. A clear picture of the complex interplay between reactivity and structural factors will be sketched. The use of proper multi-scale modeling leads to the in-depth comprehension of DNA lesions mechanisms and also to the rational design of new chemo-therapeutic agents.

On the absence of intrahelical DNA dynamics on the μs to ms timescale

DNA helices display a rich tapestry of motion on both short (<100 ns) and long (>1 ms) timescales. However, with the exception of mismatched or damaged DNA, experimental measures indicate that motions in the 1 μs to 1 ms range are effectively absent, which is often attributed to difficulties in measuring motions in this time range. We hypothesized that these motions have not been measured because there is effectively no motion on this timescale, as this provides a means to distinguish faithful Watson-Crick base-paired DNA from damaged DNA. The absence of motion on this timescale would present a 'static' DNA sequence-specific structure that matches the encounter timescales of proteins, thereby facilitating recognition. Here we report long-timescale (~10-44 μs) molecular dynamics simulations of a B-DNA duplex structure that addresses this hypothesis using both an 'Anton' machine and large ensembles of AMBER GPU simulations.