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

All-or-None Selectivity in Probing Polarity-Determined Trinucleotide Repeat Foldings with a Parity Resolution by a Beyond-Size-Matching Ligand

Abnormal amplification of trinucleotide repeats (TNRs) is associated with neurodegenerative diseases by forming a particular hairpin bulge. It is well known that the polarity and parity of TNRs can regulate the formed hairpin structures. Therefore, there is a great challenge to efficiently discriminate the hairpin structures of TNRs with substantial selectivity. Herein, we developed a fluorescent ligand of pseudohypericin (Pse) with a beyond-size-matching (BSM) geometry to selectively sense hairpin structures of GTC and CTG TNRs. The GTC hairpin structures can bind with Pse dominantly at extreme T-T mismatches by the virtue of their most extrahelical conformations, while there is no binding event to occur with the polarity-inverted counterpart CTG hairpin structures because of the limited space provided by their intrahelical T-T mismatches. In addition, this all-or-none response with the polarity-dependent folding (PoDF) is independent of the length of these TNRs. Interestingly, the parity-dependent folding (PaDF) of GTC hairpin structures can also be resolved. Besides pure TNRs, the competency of this BSM ligand to sense the PoDF and PaDF effects was also generalized to DNAs with TNRs occurring at loop and stem end regions. To our knowledge, this is the first experimental observation with the state-of-the-art performance over the fluorescence measurement of PoDF and PaDF in TNRs. Our work provides an expedient way to elucidate the TNR folding by designing ligands having BSM features.

Transient Hoogsteen Base Pairs Observed in Unbiased Molecular Dynamics Simulations of DNA

Duplex DNA is modeled as canonical B-DNA displaying the characteristic Watson-Crick base pairs. A less common and short-lived pairing of the nucleobases is the Hoogsteen (HG) conformation. The low population of the HG base pairs (<1%) necessitates extended sampling times in order to analyze through unbiased molecular dynamics (MD) simulations. We have discovered that with extended sampling times using multiple independent copies of an 18-mer sequence, the MD trajectories reproduce the expected and transient HG base pairing. Consistent with experiment, the percentage of the HG events are within the range of ∼0.1-1.0% over the combined aggregate sampling time of more than 3.6 ms. We present the reliability of the current AMBER set of nucleic acid force fields and tools to accurately simulate naturally occurring base-pairing and opening events without any bias or restraints. The mechanism consists of base pair fraying, flipping of the purine, and reformation with HG base pairs.

Atomic-Scale Tailoring and Molecular-Level Tracking of Oxygen-Containing Tungsten Single-Atom Catalysts with Enhanced Singlet Oxygen Generation

The local coordination structure of metal atoms in single-atom catalysts (SACs) greatly influences their catalytic performance. And for most SACs, single metal atoms were anchored on carbon materials with N or C coordination. However, the rational design of oxygen-containing SACs and analyzing its structure-performance relationship remain challenging. Herein, we used amino-rich compounds to tailor the metatungstate and fix the W atoms and finally obtained the oxygen-containing W-SACs. The structural evolution of tungsten and its coordination atoms were tracked by electrospray ionization high-definition mass spectrometry. Furthermore, aberration-corrected transmission electron microscopy, X-ray absorption fine-structure spectroscopy, and first-principles calculation results revealed that different from the traditional SACs, the WO₂N₂ moiety (W coordinated with two O atoms and two N atoms) may be the favored structure for W species. This special structure promoted the energy transfer for enhancing singlet oxygen generation. This work presents an efficient way to prepare more high-efficiency SACs by atomic-scale tailoring and structural evolution tracking at the molecular level.

Recognition of a tandem lesion by DNA bacterial formamidopyrimidine glycosylases explored combining molecular dynamics and machine learning

The combination of several closely spaced DNA lesions, which can be induced by a single radical hit, constitutes a hallmark in the DNA damage landscape and radiation chemistry. The occurrence of such a tandem base lesion gives rise to a strong coupling with the double helix degrees of freedom and induces important structural deformations, in contrast to DNA strands containing a single oxidized nucleobase. Although such complex lesions are known to be refractory to repair by DNA glycosylases, there is still a lack of structural evidence to rationalize these phenomena. In this contribution, we explore, by numerical modeling and molecular simulations, the behavior of the bacterial glycosylase responsible for base excision repair (MutM), specialized in excising oxidatively-damaged defects such as 7,8-dihydro-8-oxoguanine (8-oxoG). The difference in lesion recognition between a simple damage and a tandem lesion featuring an additional abasic site is assessed at atomistic resolution owing to microsecond molecular dynamics simulations and machine learning postprocessing, allowing to extensively pinpoint crucial differences in the interaction patterns of the damaged bases. Our results reveal substantial changes in the interaction network surrounding the 8-oxoG upon addition of an adjacent abasic site, leading to the perturbation of the intercalation triad which is crucial for lesion recognition and processing. The recognition process might also be impacted by a more constrained MutM-DNA binding upon tandem damage, as shown by the machine learning post-processing. This work advocates for the use of such high throughput numerical simulations for exploring the complex combinatorial chemistry of tandem DNA lesions repair and more generally local multiple damaged sites of the utmost significance in radiation chemistry.

Automated AFM analysis of DNA bending reveals initial lesion sensing strategies of DNA glycosylases

Base excision repair is the dominant DNA repair pathway of chemical modifications such as deamination, oxidation, or alkylation of DNA bases, which endanger genome integrity due to their high mutagenic potential. Detection and excision of these base lesions is achieved by DNA glycosylases. To investigate the remarkably high efficiency in target site search and recognition by these enzymes, we applied single molecule atomic force microscopy (AFM) imaging to a range of glycosylases with structurally different target lesions. Using a novel, automated, unbiased, high-throughput analysis approach, we were able to resolve subtly different conformational states of these glycosylases during DNA lesion search. Our results lend support to a model of enhanced lesion search efficiency through initial lesion detection based on altered mechanical properties at lesions. Furthermore, its enhanced sensitivity and easy applicability also to other systems recommend our novel analysis tool for investigations of diverse, fundamental biological interactions.

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.

Bending of DNA duplexes with mutation motifs

Mutations can be induced by environmental factors but also arise spontaneously during DNA replication or due to deamination of methylated cytosines at CpG dinucleotides. Sites where mutations occur with higher frequency than would be expected by chance are termed hotspots while sites that contain mutations rarely are termed coldspots. Mutations are permanently scanned and repaired by repair systems. Among them, the mismatch repair targets base pair mismatches, which are discriminated from canonical base pairs by probing altered elasticity of DNA. Using biased molecular dynamics simulations, we investigated the elasticity of coldspots and hotspots motifs detected in human genes associated with inherited disorders, and also of motifs with Czech population hotspots and de novo mutations. Main attention was paid to mutations leading to G/T and A+/C pairs. We observed that hotspots without CpG/CpHpG sequences are less flexible than coldspots, which indicates that flexible sequences are more effectively repaired. In contrary, hotspots with CpG/CpHpG sequences exhibited increased flexibility as coldspots. Their mutability is more likely related to spontaneous deamination of methylated cytosines leading to C > T mutations, which are primarily targeted by base excision repair. We corroborated conclusions based on computer simulations by measuring melting curves of hotspots and coldspots containing G/T mismatch.

Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.

Unwinding Induced Melting of Double-Stranded DNA Studied by Free Energy Simulations

DNA unwinding plays a major role in many biological processes, such as replication, transcription, and repair. It can lead to local melting and strand separation and can serve as a key mechanism to promote access to the separate strands of a double-stranded DNA. While DNA unwinding has been investigated extensively by DNA cyclization and single-molecule studies on a length-scale of kilo base pairs, it is neither fully understood at the base pair level nor at the level of molecular interactions. By employing a torque acting on the termini of DNA oligonucleotides during molecular dynamics free energy simulations, we locally unwind the central part of a DNA beyond an elastic (harmonic) regime. The simulations reproduce experimental results on the twist elasticity in the harmonic regime (characterized by a mostly quadratic free energy change with respect to changes in twist) and a deformation up to 7° was found as a limit of the harmonic response. Beyond this limit the free energy increase per twist change dropped dramatically coupled to local base pair disruptions and significant deformation of the nucleic acid backbone structure. Restriction of the DNA bending flexibility resulted in a stiffer harmonic response and an earlier onset of the anharmonic response. Whereas local melting with a complete disruption of base pairing and flipping of nucleotides was observed in case of an AT rich central segment strong backbone changes and changes in the stacking arrangements were observed in case of a GC rich segment. Unrestrained MD simulations starting from locally melted DNA reformed regular B-DNA after 50-300 ns simulation time. The simulations may have important implications for understanding DNA recognition processes coupled with significant structural alterations.