Effect of 8-oxoguanine on DNA structure and deformability

DNA damage and repair in the nucleosome: insights from computational methods

Cellular DNA is constantly exposed to endogenous or exogenous factors that can induce lesions. Several types of lesions have been described that can result from UV/ionizing irradiations, oxidative stress, or free radicals, among others. In order to overcome the deleterious effects of such damages, i.e., mutagenicity or cytotoxicity, cells possess a highly complex DNA repair machinery, involving repair enzymes targeting specific types of lesions through dedicated cellular pathways. In addition, DNA is highly compacted in the nucleus, the first level of compaction consisting of ~ 147 DNA base pairs wrapped around a core of histones, the so-called nucleosome core particle. In this complex environment, the DNA structure is highly constrained, and fine-tuned mechanisms involving remodeling processes are required to expose the DNA to repair enzymes and to facilitate the damage removal. However, these nucleosome-specific mechanisms remain poorly understood, and computational methods emerged only recently as powerful tools to investigate DNA damages in such complex systems as the nucleosome. In this mini-review, we summarize the latest advances brought out by computational approaches in the field, opening new exciting perspectives for the study of DNA damage and repair in the nucleosome context.

The potential for OGG1 inhibition to be a therapeutic strategy for pulmonary diseases

INTRODUCTION

Pulmonary diseases impose a daunting burden on healthcare systems and societies. Current treatment approaches primarily address symptoms, underscoring the urgency for the development of innovative pharmaceutical solutions. A noteworthy focus lies in targeting enzymes recognizing oxidatively modified DNA bases within gene regulatory elements, given their pivotal role in governing gene expression.

AREAS COVERED

This review delves into the intricate interplay between the substrate-specific binding of 8-oxoguanine DNA glycosylase 1 (OGG1) and epigenetic regulation, with a focal point on elucidating the molecular underpinnings and their biological implications. The absence of OGG1 distinctly attenuates the binding of transcription factors to cis elements, thereby modulating pro-inflammatory or pro-fibrotic transcriptional activity. Through a synergy of experimental insights gained from cell culture studies and murine models, utilizing prototype OGG1 inhibitors (O8, TH5487, and SU0268), a promising panorama emerges. These investigations underscore the absence of cytotoxicity and the establishment of a favorable tolerance profile for these OGG1 inhibitors.

EXPERT OPINION

Thus, the strategic targeting of the active site pocket of OGG1 through the application of small molecules introduces an innovative trajectory for advancing redox medicine. This approach holds particular significance in the context of pulmonary diseases, offering a refined avenue for their management.

How Clustered DNA Damage Can Change the Electronic Properties of ds-DNA—Differences between GAG, GAOXOG, and OXOGAOXOG

Every 24 h, roughly 3 × 1017 incidences of DNA damage are generated in the human body as a result of intra- or extra-cellular factors. The structure of the formed lesions is identical to that formed during radio- or chemotherapy. Increases in the clustered DNA damage (CDL) level during anticancer treatment have been observed compared to those found in untreated normal tissues. 7,8-dihydro-8-oxo-2′-deoxyguanosine (OXOG) has been recognized as the most common lesion. In these studies, the influence of OXOG, as an isolated (oligo-OG) or clustered DNA lesion (oligo-OGOG), on charge transfer has been analyzed in comparison to native oligo-G. DNA lesion repair depends on the damage recognition step, probably via charge transfer. Here the electronic properties of short ds-oligonucleotides were calculated and analyzed at the M062x/6-31++G** level of theory in a non-equilibrated and equilibrated solvent state. The rate constant of hole and electron transfer according to Marcus’ theory was also discussed. These studies elucidated that OXOG constitutes the sink for migrated radical cations. However, in the case of oligo-OGOG containing a 5′-OXOGAXOXG-3′ sequence, the 3′-End OXOG becomes predisposed to electron-hole accumulation contrary to the undamaged GAG fragment. Moreover, it was found that the 5′-End OXOG present in an OXOGAOXOG fragment adopts a higher adiabatic ionization potential than the 2′-deoxyguanosine of an undamaged analog if both ds-oligos are present in a cationic form. Because increases in CDL formation have been observed during radio- or chemotherapy, understanding their role in the above processes can be crucial for the efficiency and safety of medical cancer treatment.

Dynamics and energetics of PCBP1 binding to severely oxidized RNA

Oxidatively generated lesions such as 8-oxo-7, 8-dihydroguanine (8-oxoG) on RNA strands constitute a hallmark marker of the oxidative stress in the cell. Poly-C binding protein 1 (PCBP1) is able to specifically recognize severely damaged RNA strands containing two 8-oxoG lesions separated by five nucleobases, which trigger a signaling pathway leading to cell apoptosis. We apply an in silico protocol based on microsecond timescale all-atom classical molecular dynamics simulations associated with conformational and energy analyses to unveil the specific recognition mechanism at a molecular level. By comparing the RNA and protein behavior for sequences with six different damage profiles, our results highlight an allosteric mechanism, allowing a stronger binding of the oxidized guanine at position 9 only if another 8-oxoG lesion is present at position 15, in full agreement with experiments. We assess the role of lysine K23 and the additional ketone group of the oxidized guanine, thanks to computational site-directed mutagenesis.

7,8-Dihydro-8-oxoguanosine Lesions Inhibit the Theophylline Aptamer or Change Its Selectivity

Aptamers are attractive constructs due to their high affinity/selectivity towards a target. Here 7,8-dihydro-8-oxoguanosine (8-oxoG) has been used, due in part to its unique H-bonding capabilities (Watson-Crick or Hoogsteen), to expand the "RNA alphabet". Its impact on the theophylline RNA aptamer was explored by modifying its binding pocket at positions G11, G25, or G26. Structural probing, with RNases A and T1 , showed that modification at G11 leads to a drastic structural change, whereas the G25-/G26-modified analogues exhibited cleavage patterns similar to that of the canonical construct. The recognition properties towards three xanthine derivatives were then explored through thermophoresis. Modifying the aptamer at position G11 led to binding inhibition. Modification at G25, however, changed the selectivity towards theobromine (Kd ≈160 μm), with a poor affinity for theophylline (Kd >1.5 mm) being observed. Overall, 8-oxoG can have an impact on the structures of aptamers in a position-dependent manner, leading to altered target selectivity.

Telomere uncapping by common oxidative guanine lesions: Insights from atomistic models

Oxidative damage to DNA is widely known to contribute to aging and disease. This relationship has been extensively studied for telomeres - structures that cap chromosome ends - due to their role in cell proliferation and senescence, and exceptional susceptibility to oxidation. Indeed, the repetitive telomeric DNA sequence contains the 5'-GGG-3' motif that has the lowest ionization potential of all trinucleotides. Accordingly, experiments consistently show that telomeric oxidative lesions are more abundant and persistent than elsewhere in the genome. This led to a hypothesis that telomeres act as sensors of prolonged oxidative stress and prevent carcinogenesis, as disruption of telomeric integrity triggers senescence or apoptosis. Here, we use atomistic alchemical Molecular Dynamics simulations to perform a combinatorial assessment of changes in DNA binding affinity of telomeric proteins induced by oxidative guanine lesions. We rank lesions by their effect on telomere integrity, as well as telomeric proteins by their sensitivity to DNA oxidation. While the binding of most proteins is abolished by DNA oxidation, HOT1 emerges as a notable exception, suggesting its potential role in sensing of oxidative damage. Through statistical analysis and free energy decomposition, we also identify common trends in structural responses of protein-DNA complexes that contribute to decreased binding affinity.

How methyl-sugar interactions determine DNA structure and flexibility

The sequence dependent structure and flexibility of the DNA double helix is of key importance for gene expression and DNA packing and it can be modulated by DNA modifications. The presence of a C5′-methyl group in thymine or the frequent C5′-methylated-cytosine affects the DNA fine structure, however, the underlying mechanism and steric origins have remained largely unexplained. Employing Molecular Dynamics free energy simulations that allow switching on or off interactions with the methyl groups in several DNA sequences, we systematically identified the physical origin of the coupling between methyl groups and DNA backbone fine structure. Whereas methyl-solvent and methyl–nucleobase interactions were found to be of minor importance, the methyl group interaction with the 5′ neighboring sugar was identified as main cause for influencing the population of backbone substates. The sterical methyl sugar clash prevents the formation of unconventional stabilizing hydrogen bonds between nucleobase and backbone. The technique was also used to study the contribution of methyl groups to DNA flexibility and served to explain why the presence of methyl sugar clashes in thymine and methyl-cytosine can result in an overall local increase of DNA flexibility.

Oxidative damage to epigenetically methylated sites affects DNA stability, dynamics and enzymatic demethylation

DNA damage can affect various regulatory elements of the genome, with the consequences for DNA structure, dynamics, and interaction with proteins remaining largely unexplored. We used solution NMR spectroscopy, restrained and free molecular dynamics to obtain the structures and investigate dominant motions for a set of DNA duplexes containing CpG sites permuted with combinations of 5-methylcytosine (mC), the primary epigenetic base, and 8-oxoguanine (oxoG), an abundant DNA lesion. Guanine oxidation significantly changed the motion in both hemimethylated and fully methylated DNA, increased base pair breathing, induced BI→BII transition in the backbone 3′ to the oxoG and reduced the variability of shift and tilt helical parameters. UV melting experiments corroborated the NMR and molecular dynamics results, showing significant destabilization of all methylated contexts by oxoG. Notably, some dynamic and thermodynamic effects were not additive in the fully methylated oxidized CpG, indicating that the introduced modifications interact with each other. Finally, we show that the presence of oxoG biases the recognition of methylated CpG dinucleotides by ROS1, a plant enzyme involved in epigenetic DNA demethylation, in favor of the oxidized DNA strand. Thus, the conformational and dynamic effects of spurious DNA oxidation in the regulatory CpG dinucleotide can have far-reaching biological consequences.

8-Oxoguanine Affects DNA Backbone Conformation in the EcoRI Recognition Site and Inhibits Its Cleavage by the Enzyme

8-oxoguanine is one of the most abundant and impactful oxidative DNA lesions. However, the reasons underlying its effects, especially those not directly explained by the altered base pairing ability, are poorly understood. We report the effect of the lesion on the action of EcoRI, a widely used restriction endonuclease. Introduction of 8-oxoguanine inside, or adjacent to, the GAATTC recognition site embedded within the Drew—Dickerson dodecamer sequence notably reduced the EcoRI activity. Solution NMR revealed that 8-oxoguanine in the DNA duplex causes substantial alterations in the sugar—phosphate backbone conformation, inducing a BI→BII transition. Moreover, molecular dynamics of the complex suggested that 8-oxoguanine, although does not disrupt the sequence-specific contacts formed by the enzyme with DNA, shifts the distribution of BI/BII backbone conformers. Based on our data, we propose that the disruption of enzymatic cleavage can be linked with the altered backbone conformation and dynamics in the free oxidized DNA substrate and, possibly, at the protein—DNA interface.

Global deformation facilitates flipping of damaged 8-oxo-guanine and guanine in DNA

Oxidation of guanine (Gua) to form 7,8-dihydro-8-oxoguanine (8oxoG) is a frequent mutagenic DNA lesion. DNA repair glycosylases such as the bacterial MutM can effciently recognize and eliminate the 8oxoG damage by base excision. The base excision requires a 8oxoG looping out (flipping) from an intrahelical base paired to an extrahelical state where the damaged base is in the enzyme active site. It is still unclear how the damage is identified and flipped from an energetically stable stacked and paired state without any external energy source. Free energy simulations have been employed to study the flipping process for globally deformed DNA conformational states. DNA deformations were generated by systematically untwisting the DNA to mimic its conformation in repair enzyme encounter complex. The simulations indicate that global DNA untwisting deformation toward the enzyme bound form alone (without protein) significantly reduces the penalty for damage flipping to about half of the penalty observed in regular DNA. The finding offers a mechanistic explanation how binding free energy that is transformed to binding induced DNA deformation facilitates flipping and helps to rapidly detect a damaged base. It is likely of general relevance since repair enzyme binding frequently results in significant deformation of the target DNA.

Adenine Glycosylase MutY of Corynebacterium pseudotuberculosis presents the antimutator phenotype and evidences of glycosylase/AP lyase activity in vitro

Corynebacterium pseudotuberculosis is the etiological agent of caseous lymphadenitis, a disease that predominantly affects small ruminants, causing significant economic losses worldwide. As a facultative intracellular pathogen, this bacterium is exposed to an environment rich in reactive oxygen species (ROS) within macrophages. To ensure its genetic stability, C. pseudotuberculosis relies on efficient DNA repair pathways for excision of oxidative damage such as 8-oxoguanine, a highly mutagenic lesion. MutY is an adenine glycosylase involved in adenine excision from 8-oxoG:A mismatches avoiding genome mutation incorporation. The purpose of this study was to characterize MutY protein from C. pseudotuberculosis and determine its involvement with DNA repair. In vivo functional complementation assay employing mutY gene deficient Escherichia coli transformed with CpmutY showed a 13.5-fold reduction in the rate of spontaneous mutation, compared to cells transformed with empty vector. Also, under oxidative stress conditions, CpMutY protein favored the growth of mutY deficient E. coli, relative to the same strain in the absence of CpMutY. To demonstrate the involvement of this enzyme in recognition and excision of 8-oxoguanine lesion, an in vitro assay was performed. CpMutY protein was capable of recognizing and excising 8-oxoG:A but not 8-oxoG:C presenting evidences of glycosylase/AP lyase activity in vitro. In silico structural characterization revealed the presence of preserved motifs related to the MutY activity on DNA repair, such as catalytic residues involved in glycosylase/AP lyase activity and structural DNA-binding elements, such as the HhH motif and the [4Fe-4S] cluster. The three-dimensional structure of CpMutY, generated by comparative modeling, exhibits a catalytic domain very similar to that of E. coli MutY. Taken together, these results indicate that the CpmutY encodes a functional protein homologous to MutY from E. coli and is involved in the prevention of mutations and the repair of oxidative DNA lesions.

A dynamic checkpoint in oxidative lesion discrimination by formamidopyrimidine-DNA glycosylase

In contrast to proteins recognizing small-molecule ligands, DNA-dependent enzymes cannot rely solely on interactions in the substrate-binding centre to achieve their exquisite specificity. It is widely believed that substrate recognition by such enzymes involves a series of conformational changes in the enzyme-DNA complex with sequential gates favoring cognate DNA and rejecting nonsubstrates. However, direct evidence for such mechanism is limited to a few systems. We report that discrimination between the oxidative DNA lesion, 8-oxoguanine (oxoG) and its normal counterpart, guanine, by the repair enzyme, formamidopyrimidine-DNA glycosylase (Fpg), likely involves multiple gates. Fpg uses an aromatic wedge to open the Watson-Crick base pair and everts the lesion into its active site. We used molecular dynamics simulations to explore the eversion free energy landscapes of oxoG and G by Fpg, focusing on structural and energetic details of oxoG recognition. The resulting energy profiles, supported by biochemical analysis of site-directed mutants disturbing the interactions along the proposed path, show that Fpg selectively facilitates eversion of oxoG by stabilizing several intermediate states, helping the rapidly sliding enzyme avoid full extrusion of every encountered base for interrogation. Lesion recognition through multiple gating intermediates may be a common theme in DNA repair enzymes.

Enhanced sampling simulations of DNA step parameters

A novel approach for the selection of step parameters as reaction coordinates in enhanced sampling simulations of DNA is presented. The method uses three atoms per base and does not require coordinate overlays or idealized base pairs. This allowed for a highly efficient implementation of the calculation of all step parameters and their Cartesian derivatives in molecular dynamics simulations. Good correlation between the calculated and actual twist, roll, tilt, shift, and slide parameters is obtained, while the correlation with rise is modest. The method is illustrated by its application to the methylated and unmethylated 5'-CATGTGACGTCACATG-3' double stranded DNA sequence. One-dimensional umbrella simulations indicate that the flexibility of the central CG step is only marginally affected by methylation.

Coupled binding-bending-folding: The complex conformational dynamics of protein-DNA binding studied by atomistic molecular dynamics simulations

BACKGROUND

Protein-DNA binding often involves dramatic conformational changes such as protein folding and DNA bending. While thermodynamic aspects of this behavior are understood, and its biological function is often known, the mechanism by which the conformational changes occur is generally unclear. By providing detailed structural and energetic data, molecular dynamics simulations have been helpful in elucidating and rationalizing protein-DNA binding.

SCOPE OF REVIEW

This review will summarize recent atomistic molecular dynamics simulations of the conformational dynamics of DNA and protein-DNA binding. A brief overview of recent developments in DNA force fields is given as well.

MAJOR CONCLUSIONS

Simulations have been crucial in rationalizing the intrinsic flexibility of DNA, and have been instrumental in identifying the sequence of binding events, the triggers for the conformational motion, and the mechanism of binding for a number of important DNA-binding proteins.

GENERAL SIGNIFICANCE

Molecular dynamics simulations are an important tool for understanding the complex binding behavior of DNA-binding proteins. With recent advances in force fields and rapid increases in simulation time scales, simulations will become even more important for future studies. This article is part of a Special Issue entitled Recent developments of molecular dynamics.