Ionization Patterns and Chemical Reactivity of Cytosine-Guanine Watson-Crick Pairs.

The front cover artwork is provided by Prof. Papadantonakis' group. The image shows a Watson-Crick Guanine-Cytosine pair, and the difference between vertical and adiabatic ionization potentials. Read the full text of the Research Article at 10.1002/cphc.202300946.

A DFTB study on the electronic response of encapsulated DNA nucleobases onto chiral CNTs as a sequencer.

Sequencing the DNA nucleobases is essential in the diagnosis and treatment of many diseases related to human genes. In this article, the encapsulation of DNA nucleobases with some of the important synthesized chiral (7, 6), (8, 6), and (10, 8) carbon nanotubes were investigated. The structures were modeled by applying density functional theory based on tight binding method (DFTB) by considering semi-empirical basis sets. Encapsulating DNA nucleobases on the inside of CNTs caused changes in the electronic properties of the selected chiral CNTs. The results confirmed that van der Waals (vdW) interactions, π-orbitals interactions, non-bonded electron pairs, and the presence of high electronegative atoms are the key factors for these changes. The result of electronic parameters showed that among the CNTs, CNT (8, 6) is a suitable choice in sequencing guanine (G) and cytosine (C) DNA nucleobases. However, they are not able to sequence adenine (A) and thymine (T). According to the band gap energy engineering approach and absorption energy, the presence of G and C DNA nucleobases decreased the band gap energy of CNTs. Hence selected CNTs suggested as biosensor substrates for sequencing G and C DNA nucleobases.

Direct Measurement of 8OG Syn-Anti Flips in Mutagenic 8OG·A and Long-Range Damage-Dependent Hoogsteen Breathing Dynamics Using 1H CEST NMR.

Elucidating how damage impacts DNA dynamics is essential for understanding the mechanisms of damage recognition and repair. Many DNA lesions alter their propensities to form low-populated and short-lived conformational states. However, NMR methods to measure these dynamics require isotopic enrichment, which is difficult for damaged nucleotides. Here, we demonstrate the utility of the 1H chemical exchange saturation transfer (CEST) NMR experiment in measuring the dynamics of oxidatively damaged 8-oxoguanine (8OG) in the mutagenic 8OGsyn·Aanti mismatch. Using 8OG-H7 as an NMR probe of the damaged base, we directly measured 8OG syn-anti flips to form a lowly populated (pop. ∼ 5%) and short-lived (lifetime ∼50 ms) nonmutagenic 8OGanti·Aanti. These exchange parameters were in quantitative agreement with values from 13C off-resonance R1ρ and CEST on the labeled partner adenine. The Watson-Crick-like 8OGsyn·Aanti mismatch also rescued the kinetics of Hoogsteen motions at distant A-T base pairs, which the G·A mismatch had slowed down. The results lend further support for 8OGanti·Aanti as a minor conformational state of 8OG·A, reveal that 8OG damage can impact Hoogsteen dynamics at a distance, and demonstrate the utility of 1H CEST for measuring damage-dependent dynamics in unlabeled DNA.

Identification of bona fide RNA G-quadruplex binding proteins.

RNAs often accomplish their diverse functions through direct interactions with RNA-binding proteins (RBPs) in a sequence- and/or structure-dependent manner. RNA G-quadruplexes (rG4s) are unique secondary structures formed by guanine-rich RNA sequences which impact RNA function independently and in combination with RBPs. Efforts from several labs have identified dozens of rG4 specific RBPs (rG4BPs), although the research is still in the growing phase. Here we present methods for the systematic identification of rG4BPs using a pull-down approach that takes advantage of the chemical modification of guanine bases. This allows abolishing the rG4 structures while still maintaining the base composition intact, and hence helps in recognizing true rG4BPS (in contrast to G-rich motif binders). In combination with other biochemical assays, such an approach can be efficiently used for the identification and characterization of bona fide rG4BPs.

Comprehensive analysis of intramolecular G-quadruplex structures: furthering the understanding of their formalism.

G-quadruplexes (G4) are helical structures found in guanine-rich DNA or RNA sequences. Generally, their formalism is based on a few dozen structures, which can produce some inconsistencies or incompleteness. Using the website ASC-G4, we analyzed the structures of 333 intramolecular G4s, of all types, which allowed us to clarify some key concepts and present new information. To each of the eight distinguishable topologies corresponds a groove-width signature and a predominant glycosidic configuration (gc) pattern governed by the directions of the strands. The relative orientations of the stacking guanines within the strands, which we quantified and related to their vertical gc successions, determine the twist and tilt of the helices. The latter impact the minimum groove widths, which represent the space available for lateral ligand binding. The G4 four helices have similar twists, even when these twists are irregular, meaning that they have various angles along the strands. Despite its importance, the vertical gc succession has no strict one-to-one relationship with the topology, which explains the discrepancy between some topologies and their corresponding circular dichroism spectra. This study allowed us to introduce the new concept of platypus G4s, which are structures with properties corresponding to several topologies.

Structure of native four-repeat satellite III sequence with non-canonical base interactions.

Tandem-repetitive DNA (where two or more DNA bases are repeated numerous times) can adopt non-canonical secondary structures. Many of these structures are implicated in important biological processes. Human Satellite III (HSat3) is enriched for tandem repeats of the sequence ATGGA and is located in pericentromeric heterochromatin in many human chromosomes. Here, we investigate the secondary structure of the four-repeat HSat3 sequence 5'-ATGGA ATGGA ATGGA ATGGA-3' using X-ray crystallography, NMR, and biophysical methods. Circular dichroism spectroscopy, thermal stability, native PAGE, and analytical ultracentrifugation indicate that this sequence folds into a monomolecular hairpin with non-canonical base pairing and B-DNA characteristics at concentrations below 0.9 mM. NMR studies at 0.05-0.5 mM indicate that the hairpin is likely folded-over into a compact structure with high dynamics. Crystallographic studies at 2.5 mM reveal an antiparallel self-complementary duplex with the same base pairing as in the hairpin, extended into an infinite polymer. The non-canonical base pairing includes a G-G intercalation sandwiched by sheared A-G base pairs, leading to a cross-strand four guanine stack, so called guanine zipper. The guanine zippers are spaced throughout the structure by A-T/T-A base pairs. Our findings lend further insight into recurring structural motifs associated with the HSat3 and their potential biological functions.

Ionization Patterns and Chemical Reactivity of Cytosine-Guanine Watson-Crick Pairs.

Gas-phase and aqueous vertical ionization potentials, vIPgas and vIPaq respectively and measurements of the molecular electrostatic and local ionization maps calculated at the DFT/B3LYP-D3/ 6-311+G** level of theory and the C-PCM reaction field model for single- and double-stranded CpG and 5MeCpG pairs show that the vIPaq for single- and double-stranded pairs of C-G and 5MeC-G are practically the same, in the range of 5.79 to 5.81 eV. The aqueous adiabatic ionization potentials for single-stranded CpG and 5MeCpG are 5.52 eV and 5.51 eV respectively and they reflect the nuclear reorganization that takes place after the abstraction of the electron. The aqueous adiabatic ionization energy values that correspond to the CpG+. radical cation and the hydrated electron, e-,, being at infinite distance, adIPaq+Vo, are 3.92 eV and 3.91 eV respectively with (Vo=-1.6 eV) Analysis of data suggest that the HOMO-LUMO energy gap in the hard/soft-acid/base (HSAB) concept cannot be used a priori to determine the effect of cytosine methylation on the guanine enhanced oxidative damage in DNA.

Formation and Nanomechanical Properties of Silver-Mediated Guanine DNA Duplexes in Aqueous Solution.

Silver cations can mediate base pairing of guanine (G) DNA oligomers, yielding linear parallel G-Ag+-G duplexes with enhanced stabilities compared to those of canonical DNA duplexes. To enable their use in programmable DNA nanotechnologies, it is critical to understand solution-state formation and the nanomechanical stiffness of G-Ag+-G duplexes. Using temperature-controlled circular dichroism (CD) spectroscopy, we find that heating mixtures of G oligomers and silver salt above 50 °C fully destabilizes G-quadruplex structures and converts oligomers to G-Ag+-G duplexes. Electrospray ionization mass spectrometry supports that G-Ag+-G duplexes form at stoichiometries of 1 Ag+ per base pair, and CD spectroscopy suggests that as the Ag+/base stoichiometry increases further, G-Ag+-G duplexes undergo additional morphological changes. Using liquid-phase atomic force microscopy, we find that this excess Ag+ enables assembly of long fiberlike structures with ∼2.5 nm heights equivalent to a single DNA duplex but with lengths that far exceed a single duplex. Finally, using the conditions established to form single G-Ag+-G duplexes, we use a surface forces apparatus (SFA) to compare the solution-phase stiffness of single G-Ag+-G duplexes with dG-dC Watson-Crick-Franklin duplexes. SFA shows that G-Ag+-G duplexes are 1.3 times stiffer than dG-dC duplexes, confirming gas-phase ion mobility spectrometry measurements and computational predictions. These findings may guide the development of structural DNA nanotechnologies that rely on silver-mediated base pairing.

Structure of a DNA G-quadruplex that Modulates SP1 Binding Sites Architecture in HIV-1 Promoter.

Nucleic acid sequences containing guanine tracts are able to form non-canonical DNA or RNA structures known as G-quadruplexes (or G4s). These structures, based on the stacking of G-tetrads, are involved in various biological processes such as gene expression regulation. Here, we investigated a G4 forming sequence, HIVpro2, derived from the HIV-1 promoter. This motif is located 60 nucleotides upstream of the proviral Transcription Starting Site (TSS) and overlaps with two SP1 transcription factor binding sites. Using NMR spectroscopy, we determined that HIVpro2 forms a hybrid type G4 structure with a core that is interrupted by a single nucleotide bulge. An additional reverse-Hoogsteen AT base pair is stacked on top of the tetrad. SP1 transcription factor is known to regulate transcription activity of many genes through the recognition of Guanine-rich duplex motifs. Here, the formation of HIVpro2 G4 may modulate SP1 binding sites architecture by competing with the formation of the canonical duplex structure. Such DNA structural switch potentially participates to the regulation of viral transcription and may also interfere with HIV-1 reactivation or viral latency.

Best method to determine DNA G-quadruplex folding: The 1H-13C HSQC NMR experiment.

NMR spectroscopy is the major method for G-quadruplex structure determination under physiologically relevant solution conditions. Unlike duplex B-DNA, in which all nucleotides adopt an anti glycosidic conformation, the core tetrad-guanines in a G-quadruplex can adopt anti or syn glycosidic conformation depending on the folding structure. An experimental method that can clearly and unambiguously determine syn and anti tetrad-Gs in a G-quadruplex is highly desirable and necessary. In the present study, we exploit the advantages of the 1H-13C HSQC experiment to determine tetrad-G's glycosidic conformation and thus folding topology of G-quadruplexes. We use several examples to demonstrate the clear and straightforward determination of the guanine glycosidic conformations and G-quadruplex folding structures. Moreover, 1H-13C HSQC data can readily identify adenine H2 resonances as well as determine unusual syn conformation in loop and flanking sequences, a challenging task by standard 2D NOESY.

Sensing of DNA modifications by pAgo proteins in vitro.

Many prokaryotic Argonaute (pAgo) proteins act as programmable nucleases that use small guide DNAs for recognition and cleavage of complementary target DNA. Recent studies suggested that pAgos participate in cell defense against invader DNA and may also be involved in other genetic processes, including DNA replication and repair. The ability of pAgos to recognize specific targets potentially make them an invaluable tool for DNA manipulations. Here, we demonstrate that DNA-guided DNA-targeting pAgo nucleases from three bacterial species, DloAgo from Dorea longicatena, CbAgo from Clostridium butyricum and KmAgo from Kurthia massiliensis, can sense site-specific modifications in the target DNA, including 8-oxoguanine, thymine glycol, ethenoadenine and pyrimidine dimers. The effects of DNA modifications on the activity of pAgos strongly depend on their positions relative to the site of cleavage and are comparable to or exceed the effects of guide-target mismatches at corresponding positions. For all tested pAgos, the strongest effects are observed when DNA lesions are located at the cleavage position. The results demonstrate that DNA cleavage by pAgos is strongly affected by DNA modifications, thus making possible their use as sensors of DNA damage.

Base flipping mechanism and binding strength of methyl-damaged DNA during the interaction with AGT.

Methyl damage to DNA bases is common in the cell nucleus. O6-alkylguanine-DNA alkyl transferase (AGT) may be a promising candidate for direct damage reversal in methylated DNA (mDNA) at the O6 point of the guanine. Indeed, atomic-level investigations in the contact region of AGT-DNA complex can provide an in-depth understanding of their binding mechanism, allowing to evaluate the silico-drug nature of AGT and its utility in removing methyl damage in DNA. In this study, molecular dynamics (MD) simulation was utilized to examine the flipping of methylated nucleotide, the binding mechanism between mDNA and AGT, and the comparison of binding strength prior and post methyl transfer to AGT. The study reveals that methylation at the O6 atom of guanine weakens the hydrogen bond (H-bond) between guanine and cytosine, permitting for the flipping of such nucleotide. The formation of a H-bond between the base pair of methylated nucleotide (i.e., cytosine) and the intercalated arginine of AGT also forces the nucleotide to rotate. Following that, electrostatics and van der Waals contacts as well as hydrogen bonding contribute to form the complex of DNA and protein. The stronger binding of AGT with DNA before methyl transfer creates the suitable condition to transfer methyl adduct from DNA to AGT.

Multiscale simulations reveal the role of PcrA helicase in protecting against spontaneous point mutations in DNA.

Proton transfer across hydrogen bonds in DNA can produce non-canonical nucleobase dimers and is a possible source of single-point mutations when these forms mismatch under replication. Previous computational studies have revealed this process to be energetically feasible for the guanine-cytosine (GC) base pair, but the tautomeric product (G[Formula: see text]C[Formula: see text]) is short-lived. In this work we reveal, for the first time, the direct effect of the replisome enzymes on proton transfer, rectifying the shortcomings of existing models. Multi-scale quantum mechanical/molecular dynamics (QM/MM) simulations reveal the effect of the bacterial PcrA Helicase on the double proton transfer in the GC base pair. It is shown that the local protein environment drastically increases the activation and reaction energies for the double proton transfer, modifying the tautomeric equilibrium. We propose a regime in which the proton transfer is dominated by tunnelling, taking place instantaneously and without atomic rearrangement of the local environment. In this paradigm, we can reconcile the metastable nature of the tautomer and show that ensemble averaging methods obscure detail in the reaction profile. Our results highlight the importance of explicit environmental models and suggest that asparagine N624 serves a secondary function of reducing spontaneous mutations in PcrA Helicase.

Quantitative Structure-Activity Relationship in the Series of 5-Ethyluridine, N2-Guanine, and 6-Oxopurine Derivatives with Pronounced Anti-Herpetic Activity.

A quantitative analysis of the relationship between the structure and inhibitory activity against the herpes simplex virus thymidine kinase (HSV-TK) was performed for the series of 5-ethyluridine, N2-guanine, and 6-oxopurines derivatives with pronounced anti-herpetic activity (IC50 = 0.09 ÷ 160,000 µmol/L) using the GUSAR 2019 software. On the basis of the MNA and QNA descriptors and whole-molecule descriptors using the self-consistent regression, 12 statistically significant consensus models for predicting numerical pIC50 values were constructed. These models demonstrated high predictive accuracy for the training and test sets. Molecular fragments of HSV-1 and HSV-2 TK inhibitors that enhance or diminish the anti-herpetic activity are considered. Virtual screening of the ChEMBL database using the developed QSAR models revealed 42 new effective HSV-1 and HSV-2 TK inhibitors. These compounds are promising for further research. The obtained data open up new opportunities for developing novel effective inhibitors of TK.

Strengthened cooperativity of DNA-based cyclic hydrogen-bonded rosettes by subtle functionalization.

Cooperative effects cause extra stabilization of hydrogen-bonded supramolecular systems. In this work we have designed hydrogen-bonded rosettes derived from a guanine-cytosine Janus-type motif with the aim of finding a monomer that enhances the synergy of supramolecular systems. For this, relativistic dispersion-corrected density functional theory computations have been performed. Our proposal involves a monomer with three hydrogen-bonds pointing in the same direction, which translates into shorter bonds, stronger donor-acceptor interactions, and more attractive electrostatic interactions, thus giving rise to rosettes with strengthened cooperativity. This newly designed rosette has triple the cooperativity found for the naturally occurring guanine quadruplex.

Effect of Guanine Adduct Size, Shape, and Linker Type on the Conformation of Adducted DNA: A DFT and Molecular Dynamics Study.

DNA is damaged through various exogenous sources (e.g., automobile exhaust, tobacco smoke, and processed foods), which can yield diverse C8-dG bulky aryl adducts. Adducts are known to induce structural changes to DNA that can lead to various biological outcomes, ranging from cell death to diseases such as cancer. Unfortunately, the relationship between the chemical composition of the damaged product, the adducted DNA structure, and the biological consequences is not well understood, which limits the development of disease detection and prevention strategies. The present study uses density functional theory (DFT) calculations and quintuplicate 1 µs molecular dynamics (MD) simulations to characterize the structure of DNA containing 21 model C8-dG adducts that systematically differ in size (phenyl to pyrenyl), shape (α (2,3), ß (3,4) fusion, or ring substitution), and nucleobase-aryl group linkage (N, O, and C-linked). DFT calculations reveal that the inherent structural features of the G nucleobase adducts are impacted by linker type and bulky moiety shape, but not size, with the conformational flexibility reducing with α-ring fusion and linker composition as N > O > C. These structural properties are maintained in nucleoside models, which also reveal an increased propensity for anti-to-syn rotation about the glycosidic bond with N < O < C linker type. Although these diverse chemical features do not influence the global structure of adducted DNA, the adducts differentially impact the conformation local to the adducted site, including the relative populations of structures with the bulky moiety in the major groove (B conformer) and intercalated (stacked) into the helix (S conformer). Specifically, while the smallest phenyl adducts favor the B conformation and the largest pyrenyl-derived adducts stabilize the S conformation, the B/S ratio decreases with an increase in ring size and N > O > C linker composition. The shape and size (length) of the adduct can further finetune the B/S ratio, with ß-fused naphthyl or α-fused phenanthryl N-linked adducts and O or C-linked adducts containing ring substitution increasing the prevalence of the S adducted DNA conformation. Overall, this work uncovers the significant effect of bulky moiety size and linker type, as well as the lesser impact of aryl group shape, on adducted DNA structure, which suggests differential replication and repair outcomes, and thereby represents an important step toward rationalizing connections between the structure and biological consequences of diverse DNA adducts.

8-oxoguanine riboswitches in bacteria detect and respond to oxidative DNA damage.

Riboswitches rely on structured aptamer domains to selectively sense their target ligands and regulate gene expression. However, some riboswitch aptamers in bacteria carry mutations in their otherwise strictly conserved binding pockets that change ligand specificities. The aptamer domain of a riboswitch class originally found to selectively sense guanine forms a three-stem junction that has since been observed to exploit numerous alterations in its ligand-binding pocket. These rare variants have modified their ligand specificities to sense other purines or purine derivatives, including adenine, 2'-deoxyguanosine (three classes), and xanthine. Herein, we report the characteristics of a rare variant that is narrowly distributed in the Paenibacillaceae family of bacteria. Known representatives are always associated with genes encoding 8-oxoguanine deaminase. As predicted from this gene association, these variant riboswitches tightly bind 8-oxoguanine (8-oxoG), strongly discriminate against other purine derivatives, and function as genetic "ON" switches. Following exposure of cells to certain oxidative stresses, a representative 8-oxoG riboswitch activates gene expression, likely caused by the accumulation of 8-oxoG due to oxidative damage to G nucleobases in DNA, RNA, and the nucleotide pool. Furthermore, an engineered version of the variant aptamer was prepared that exhibits specificity for 8-oxoadenine, further demonstrating that RNA aptamers can acquire mutations that expand their ability to detect and respond to oxidative damage.

Investigating the Spectroscopy of the Gas Phase Guanine-Cytosine Pair: Keto versus Enol Configurations.

We report on a vibrational study of the guanine-cytosine dimer tautomers using state-of-the-art quasiclassical trajectory and semiclassical vibrational spectroscopy. The latter includes possible quantum mechanical effects. Through an accurate comparison to the experimental spectra, we are able to shine a light on the hydrogen bond network of one of the main subunits of DNA and put the experimental assignment on a solid footing. Our calculations corroborate the experimental conclusion that the global minimum Watson-and-Crick structure is not detected in the spectra, and there is no evidence of tunnel-effect-based double proton hopping. Our accurate assignment of the spectral features may also serve as a basis for the development of precise force fields to study the guanine-cytosine dimer.

Reaction of human telomeric unit TTAGGG and a photoactivatable Pt(IV) anticancer prodrug.

The interaction of a photoactivatable diazidodihydroxido Pt(IV) prodrug, trans,trans,trans-[Pt(N3)2(OH)2(py)2] (py = pyridine; 1), with a hexamer straight human telomeric DNA unit sequence (5'-T1T2A3G4G5G6-3', I) upon light irradiation was investigated by electrospray ionization mass spectroscopy (ESI-MS). In the primary mass spectrum, two major mono-platinated I adducts with the bound Pt moieties, trans-[PtII(N3)(py)2]+ (1') and trans-[PtII(py)2]2+ (1''), respectively, were detected. It is rare to observe such high abundance and nearly equal intensity platinated DNA adducts formed by these two PtII species because 1' is usually the only major reduced Pt(II) species produced by the photodecomposition of complex 1 in the presence of DNA while 1'' was rarely detected as the major reduced PtII species reported previously. Subsequent tandem mass spectrometric analysis by collision-induced dissociation (CID) showed that in the former adduct {I + 1'}2+, G6 and A3 were the platination sites. While in the latter adduct {I + 1''}2+, a potential intrastrand crosslink was speculated after G4 and G6 sites were identified. Additionally, other minor platinated adducts like di-platinated I adduct by 1' with platination sites at G4 and G6 and mono-platinated I adducts containing base oxidation were also detected by mass spectrometry. Due to the rich guanines and their sensitivity to oxidation, the oxidation induced by 1 most probably occurred at guanine. The oxidation adducts were proposed as 8-hydroxyl guanine, spiroiminodihydantoin (Sp), 2,6-diamino-4-hydroxy-5-formamidopyrimidine (FapyG), 5-guanidinohydantoin (Gh), and/or dehydroguanidinohydantoin (DGh) referring to previous reports. The obtained results provide useful chemical information about the photoreaction between photoactivatable Pt(IV) anticancer prodrugs and human telomeric DNA. Such special damages of Pt(IV) prodrugs on human telomeric DNA implicate its active role in the mechanism of Pt(IV) prodrugs and further support the unique sequence-dependent photointeraction profile of complex 1 reacting with DNA.

RNA structural probing of guanine and uracil nucleotides in yeast.

RNA structure can be essential for its cellular function. Therefore, methods to investigate the structure of RNA in vivo are of great importance for understanding the role of cellular RNAs. RNA structure probing is an indirect method to asess the three-dimensional structure of RNA by analyzing the reactivity of different nucleotides to chemical modifications. Dimethyl sulfate (DMS) is a well-established compound that reports on base pairing context of adenine (A) and cytidine (C) in-vitro and in-vivo, but is not reactive to guanine (G) or uracil (U). Recently, new compounds were used to modify Gs and Us in plant, bacteria, and human cells. To complement the scope of RNA structural probing by chemical modifications in the model organism yeast, we analyze the effectiveness of guanine modification by the glyoxal family in Saccharomyces cerevisiae and Candida albicans. We show that within glyoxal family of compounds, phenylglyoxal (PGO) is the best guanine probe for structural probing in S. cerevisiae and C. albicans. Further, we show that PGO treatment does not affect the processing of different RNA species in the cell and is not toxic for the cells under the conditions we have established for RNA structural probing. We also explore the effectiveness of uracil modification by Cyclohexyl-3-(2-Morpholinoethyl) Carbodiimide metho-p-Toluenesulfonate (CMCT) in vivo and demonstrate that uracils can be modified by CMCT in S. cerevisiae in vivo. Our results provide the conditions for in vivo probing the reactivity of guanine and uracil nucleotides in RNA structures in yeast and offer a valuable tool for studying RNA structure and function in two widely used yeast model systems.