Triple helical DNA in a duplex context and base pair opening

Toxicity of rare earth elements europium and samarium on zebrafish development and locomotor performance

Rare earth elements are integral to modern technology, but their increasing environmental distribution due to anthropogenic activities poses potential health risks to humans. This study utilized zebrafish as a model to assess developmental and locomotor performance effects of europium and samarium. Exposure to Eu or Sm induced a reduction in heart rate, growth inhibition, and morphological deformities. RNA-Seq revealed gene expression alterations linked to critical biological processes and functions following Eu or Sm exposure. Impaired organogenesis in liver and exocrine pancreas, evident through fluorescence imaging, was confirmed transcriptionally. Exposure to Eu or Sm significantly impaired the burst and spontaneous swimming behaviors of zebrafish larvae, characterized by pronounced reductions in movement distance, frequency, and velocity. These observations indicate severe locomotor dysfunction in zebrafish exposed to Eu and Sm. The comprehensive downregulation of the oxidative phosphorylation pathway is likely a primary factor contributing to these motor impairments. Apoptosis induced by Eu and Sm, confirmed through acridine orange staining, was accompanied by the upregulation of the intrinsic apoptosis pathway. Our findings contribute critical insights into the health risks of rare earth elements, informing risk assessment and management strategies.

A comprehensive study of Z-DNA density and its evolutionary implications in birds

BACKGROUND

Z-DNA, a left-handed helical form of DNA, plays a significant role in genomic stability and gene regulation. Its formation, associated with high GC content and repetitive sequences, is linked to genomic instability, potentially leading to large-scale deletions and contributing to phenotypic diversity and evolutionary adaptation.

RESULTS

In this study, we analyzed the density of Z-DNA-prone motifs of 154 avian genomes using the non-B DNA Motif Search Tool (nBMST). Our findings indicate a higher prevalence of Z-DNA motifs in promoter regions across all avian species compared to other genomic regions. A negative correlation was observed between Z-DNA density and developmental time in birds, suggesting that species with shorter developmental periods tend to have higher Z-DNA densities. This relationship implies that Z-DNA may influence the timing and regulation of development in avian species. Furthermore, Z-DNA density showed associations with traits such as body mass, egg mass, and genome size, highlighting the complex interactions between genome architecture and phenotypic characteristics. Gene Ontology (GO) analysis revealed that Z-DNA motifs are enriched in genes involved in nucleic acid binding, kinase activity, and translation regulation, suggesting a role in fine-tuning gene expression essential for cellular functions and responses to environmental changes. Additionally, the potential of Z-DNA to drive genomic instability and facilitate adaptive evolution underscores its importance in shaping phenotypic diversity.

CONCLUSIONS

This study emphasizes the role of Z-DNA as a dynamic genomic element contributing to gene regulation, genomic stability, and phenotypic diversity in avian species. Future research should experimentally validate these associations and explore the molecular mechanisms by which Z-DNA influences avian biology.

RNA:DNA triplexes: a mechanism for epigenetic communication between hosts and microbes?

Molecular communication between host and microbe is mediated by the transfer of many different classes of macromolecules. Recently, the trafficking of RNA molecules between organisms has gained prominence as an efficient way to manipulate gene expression via RNA interference (RNAi). Here, we posit a new epigenetic control mechanism based on triple helix (triplex) structures comprising nucleic acids from both host and microbe. Indeed, RNA:DNA triplexes are known to regulate gene expression in humans, but it is unknown whether interkingdom triplexes are formed either to manipulate host processes during pathogenesis or as a host defense response. We hypothesize that a fraction of the extracellular RNAs commonly released by microbes (e.g., bacteria, fungi, and protists) and their hosts form triplexes with the genome of the other species, thereby impacting chromatin conformation and gene expression. We invite the field to consider interkingdom triplexes as unexplored weaponry in the arms race between host and microbe.

Structure and stability of different triplets involving artificial nucleobases: clues for the formation of semisynthetic triple helical DNA

A triple helical DNA can control gene expression, help in homologous recombination, induce mutations to facilitate DNA repair mechanisms, suppress oncogene formations, etc. However, the structure and function of semisynthetic triple helical DNA are not known. To understand this, various triplets formed between eight artificial nucleobases (P, Z, J, V, B, S, X, and K) and four natural DNA bases (G, C, A, and T) are studied herein by employing a reliable density functional theoretic (DFT) method. Initially, the triple helix-forming artificial nucleobases interacted with the duplex DNA containing GC and AT base pairs, and subsequently, triple helix-forming natural bases (G and C) interacted with artificial duplex DNA containing PZ, JV, BS, and XK base pairs. Among the different triplets formed in the first category, the C-JV triplet is found to be the most stable with a binding energy of about - 31 kcal/mol. Similarly, among the second category of triplets, the Z-GC and V-GC triplets are the most stable. Interestingly, Z-GC and V-GC are found to be isoenergetic with a binding energy of about - 30 kcal/mol. The C-JV, and Z-GC or V-GC triplets are about 12-14 kcal/mol more stable than the JV and GC base pairs respectively. Microsolvation of these triplets in 5 explicit water molecules further enhanced their stability by 16-21 kcal/mol. These results along with the consecutive stacking of the C-JV triplet (C-JV/C-JV) data indicate that the synthetic nucleobases can form stable semisynthetic triple helical DNA. However, consideration of a full-length DNA containing one or more semisynthetic bases or base pairs is necessary to understand the formation of semisynthetic DNA in living cells.

The chromatin - triple helix connection

Mammalian genomes are extensively transcribed, producing a large number of coding and non-coding transcripts. A large fraction of the nuclear RNAs is physically associated with chromatin, functioning in gene activation and silencing, shaping higher-order genome organisation, such as involvement in long-range enhancer-promoter interactions, transcription hubs, heterochromatin, nuclear bodies and phase transitions. Different mechanisms allow the tethering of these chromatin-associated RNAs (caRNA) to chromosomes, including RNA binding proteins, the RNA polymerases and R-loops. In this review, we focus on the sequence-specific targeting of RNA to DNA by forming triple helical structures and describe its interplay with chromatin. It turns out that nucleosome positioning at triple helix target sites and the nucleosome itself are essential factors in determining the formation and stability of triple helices. The histone H3-tail plays a critical role in triple helix stabilisation, and the role of its epigenetic modifications in this process is discussed.

DNA hydrogels and nanogels for diagnostics, therapeutics, and theragnostics of various cancers

As an efficient class of hydrogel-based therapeutic drug delivery systems, deoxyribonucleic acid (DNA) hydrogels (particularly DNA nanogels) have attracted massive attention in the last five years. The main contributor to this is the programmability of these 3-dimensional (3D) scaffolds that creates fundamental effects, especially in treating cancer diseases. Like other active biological ingredients (ABIs), DNA hydrogels can be functionalized with other active agents that play a role in targeting drug delivery and modifying the half-life of the therapeutic cargoes in the body's internal environment. Considering the brilliant advantages of DNA hydrogels, in this survey, we intend to submit an informative collection of feasible methods for the design and preparation of DNA hydrogels and nanogels, and the responsivity of the immune system to these therapeutic cargoes. Moreover, the interactions of DNA hydrogels with cancer biomarkers are discussed in this account. Theragnostic DNA nanogels as an advanced species for both detection and therapeutic purposes are also briefly reviewed.

The effects of RNA.DNA-DNA triple helices on nucleosome structures and dynamics

Noncoding RNAs (ncRNAs) are an emerging epigenetic factor and have been recognized as playing a key role in many gene expression pathways. Structurally, binding of ncRNAs to isolated DNA is strongly dependent on sequence complementary and results in the formation of an RNA.DNA-DNA (RDD) triple helix. However, in vivo DNA is not isolated but is rather packed in chromatin fibers, the fundamental unit of which is the nucleosome. Biochemical experiments have shown that ncRNA binding to nucleosomal DNA is elevated at DNA entry and exit sites and is dependent on the presence of the H3 N-terminal tails. However, the structural and dynamical bases for these mechanisms remain unknown. Here, we have examined the mechanisms and effects of RDD formation in the context of the nucleosome using a series of all-atom molecular dynamics simulations. Results highlight the importance of DNA sequence on complex stability, elucidate the effects of the H3 tails on RDD structures, show how RDD formation impacts the structure and dynamics of the H3 tails, and show how RNA alters the local and global DNA double-helical structure. Together, our results suggest ncRNAs can modify nucleosome, and potentially higher-order chromatin, structures and dynamics as a means of exerting epigenetic control.

Three's a crowd - stabilisation, structure, and applications of DNA triplexes

DNA is a strikingly flexible molecule and can form a variety of secondary structures, including the triple helix, which is the subject of this review. The DNA triplex may be formed naturally, during homologous recombination, or can be formed by the introduction of a synthetic triplex forming oligonucleotide (TFO) to a DNA duplex. As the TFO will bind to the duplex with sequence specificity, there is significant interest in developing TFOs with potential therapeutic applications, including using TFOs as a delivery mechanism for compounds able to modify or damage DNA. However, to combine triplexes with functionalised compounds, a full understanding of triplex structure and chemical modification strategies, which may increase triplex stability or in vivo degradation, is essential - these areas will be discussed in this review. Ruthenium polypyridyl complexes, which are able to photooxidise DNA and act as luminescent DNA probes, may serve as a suitable photophysical payload for a TFO system and the developments in this area in the context of DNA triplexes will also be reviewed.

Recognition of ATT Triplex and DNA:RNA Hybrid Structures by Benzothiazole Ligands

Interactions of an array of nucleic acid structures with a small series of benzothiazole ligands (bis-benzothiazolyl-pyridines—group 1, 2-thienyl/2-benzothienyl-substituted 6-(2-imidazolinyl)benzothiazoles—group 2, and three 2-aryl/heteroaryl-substituted 6-(2-imidazolinyl)benzothiazoles—group 3) were screened by competition dialysis. Due to the involvement of DNA:RNA hybrids and triplex helices in many essential functions in cells, this study’s main aim is to detect benzothiazole-based moieties with selective binding or spectroscopic response to these nucleic structures compared to regular (non-hybrid) DNA and RNA duplexes and single-stranded forms. Complexes of nucleic acids and benzothiazoles, selected by this method, were characterized by UV/Vis, fluorescence and circular dichroism (CD) spectroscopy, isothermal titration calorimetry, and molecular modeling. Two compounds (1 and 6) from groups 1 and 2 demonstrated the highest affinities against 13 nucleic acid structures, while another compound (5) from group 2, despite lower affinities, yielded higher selectivity among studied compounds. Compound 1 significantly inhibited RNase H. Compound 6 could differentiate between B- (binding of 6 dimers inside minor groove) and A-type (intercalation) helices by an induced CD signal, while both 5 and 6 selectively stabilized ATT triplex in regard to AT duplex. Compound 3 induced strong condensation-like changes in CD spectra of AT-rich DNA sequences.

Dynamics Studies of DNA with Non-canonical Structure Using NMR Spectroscopy

The non-canonical structures of nucleic acids are essential for their diverse functions during various biological processes. These non-canonical structures can undergo conformational exchange among multiple structural states. Data on their dynamics can illustrate conformational transitions that play important roles in folding, stability, and biological function. Here, we discuss several examples of the non-canonical structures of DNA focusing on their dynamic characterization by NMR spectroscopy: (1) G-quadruplex structures and their complexes with target proteins; (2) i-motif structures and their complexes with proteins; (3) triplex structures; (4) left-handed Z-DNAs and their complexes with various Z-DNA binding proteins. This review provides insight into how the dynamic features of non-canonical DNA structures contribute to essential biological processes.

Nucleosomes Stabilize ssRNA-dsDNA Triple Helices in Human Cells

Chromatin-associated non-coding RNAs modulate the epigenetic landscape and its associated gene expression program. The formation of triple helices is one mechanism of sequence-specific targeting of RNA to chromatin. With this study, we show an important role of the nucleosome and its relative positioning to the triplex targeting site (TTS) in stabilizing RNA-DNA triplexes in vitro and in vivo. Triplex stabilization depends on the histone H3 tail and the location of the TTS close to the nucleosomal DNA entry-exit site. Genome-wide analysis of TTS-nucleosome arrangements revealed a defined chromatin organization with an enrichment of arrangements that allow triplex formation at active regulatory sites and accessible chromatin. We further developed a method to monitor nucleosome-RNA triplexes in vivo (TRIP-seq), revealing RNA binding to TTS sites adjacent to nucleosomes. Our data strongly support an activating role for RNA triplex-nucleosome complexes, pinpointing triplex-mediated epigenetic regulation in vivo.

The ability of locked nucleic acid oligonucleotides to pre-structure the double helix: A molecular simulation and binding study

Locked nucleic acid (LNA) oligonucleotides bind DNA target sequences forming Watson-Crick and Hoogsteen base pairs, and are therefore of interest for medical applications. To be biologically active, such an oligonucleotide has to efficiently bind the target sequence. Here we used molecular dynamics simulations and electrophoresis mobility shift assays to elucidate the relation between helical structure and affinity for LNA-containing oligonucleotides. In particular, we have studied how LNA substitutions in the polypyrimidine strand of a duplex (thus forming a hetero duplex, i.e. a duplex with a DNA polypurine strand and an LNA/DNA polypyrimidine strand) enhance triplex formation. Based on seven polypyrimidine single strand oligonucleotides, having LNAs in different positions and quantities, we show that alternating LNA with one or more non-modified DNA nucleotides pre-organizes the hetero duplex toward a triple-helical-like conformation. This in turn promotes triplex formation, while consecutive LNAs distort the duplex structure disfavoring triplex formation. The results support the hypothesis that a pre-organization in the hetero duplex structure enhances the binding of triplex forming oligonucleotides. Our findings may serve as a criterion in the design of new tools for efficient oligonucleotide hybridization.

Modeling p K Shift in DNA Triplexes Containing Locked Nucleic Acids

The protonation states for nucleic acid bases are difficult to assess experimentally. In the context of DNA triplex, the protonation state of cytidine in the third strand is particularly important, because it needs to be protonated in order to form Hoogsteen hydrogen bonds. A sugar modification, locked nucleic acid (LNA), is widely used in triplex forming oligonucleotides to target sites in the human genome. In this study, the parameters for LNA are developed in line with the CHARMM nucleic acid force field and validated toward the available structural experimental data. In conjunction, two computational methods were used to calculate the protonation state of the third strand cytidine in various DNA triplex environments: λ-dynamics and multiple pH regime. Both approaches predict p K of this cytidine shifted above physiological pH when cytidine is in the third strand in a triplex environment. Both methods show an upshift due to cytidine methylation, and a small downshift when the sugar configuration is locked. The predicted p K values for cytidine in DNA triplex environment can inform the design of better-binding oligonucleotides.

Purine- and pyrimidine-triple-helix-forming oligonucleotides recognize qualitatively different target sites at the ribosomal DNA locus

Triplexes are noncanonical DNA structures, which are functionally associated with regulation of gene expression through ncRNA targeting to chromatin. Based on the rules of Hoogsteen base-pairing, polypurine sequences of a duplex can potentially form triplex structures with single-stranded oligonucleotides. Prediction of triplex-forming sequences by bioinformatics analyses have revealed enrichment of potential triplex targeting sites (TTS) at regulatory elements, mainly in promoters and enhancers, suggesting a potential function of RNA-DNA triplexes in transcriptional regulation. Here, we have quantitatively evaluated the potential of different sequences of human and mouse ribosomal RNA genes (rDNA) to form triplexes at different salt and pH conditions. We show by biochemical and biophysical approaches that some of these predicted sequences form triplexes with high affinity, following the canonical rules for triplex formation. We further show that RNA triplex-forming oligos (TFOs) are more stable than their DNA counterpart, and point mutations strongly affect triplex formation. We further show differential sequence requirements of pyrimidine and purine TFO sequences for efficient binding, depending on the G-C content of the TTS. The unexpected sequence specificity, revealing distinct sequence requirements for purine and pyrimidine TFOs, shows that in addition to the Hoogsteen pairing rules, a sequence code and mutations have to be taken into account to predict genomic TTS.

LNA effects on DNA binding and conformation: from single strand to duplex and triplex structures

The anti-gene strategy is based on sequence-specific recognition of double-strand DNA by triplex forming (TFOs) or DNA strand invading oligonucleotides to modulate gene expression. To be efficient, the oligonucleotides (ONs) should target DNA selectively, with high affinity. Here we combined hybridization analysis and electrophoretic mobility shift assay with molecular dynamics (MD) simulations to better understand the underlying structural features of modified ONs in stabilizing duplex- and triplex structures. Particularly, we investigated the role played by the position and number of locked nucleic acid (LNA) substitutions in the ON when targeting a c-MYC or FXN (Frataxin) sequence. We found that LNA-containing single strand TFOs are conformationally pre-organized for major groove binding. Reduced content of LNA at consecutive positions at the 3'-end of a TFO destabilizes the triplex structure, whereas the presence of Twisted Intercalating Nucleic Acid (TINA) at the 3'-end of the TFO increases the rate and extent of triplex formation. A triplex-specific intercalating benzoquinoquinoxaline (BQQ) compound highly stabilizes LNA-containing triplex structures. Moreover, LNA-substitution in the duplex pyrimidine strand alters the double helix structure, affecting x-displacement, slide and twist favoring triplex formation through enhanced TFO major groove accommodation. Collectively, these findings should facilitate the design of potent anti-gene ONs.

Role of Pseudoisocytidine Tautomerization in Triplex-Forming Oligonucleotides: In Silico and in Vitro Studies

Pseudoisocytidine (^ΨC) is a synthetic cytidine analogue that can target DNA duplex to form parallel triplex at neutral pH. Pseudoisocytidine has mainly two tautomers, of which only one is favorable for triplex formation. In this study, we investigated the effect of sequence on ^ΨC tautomerization using λ-dynamics simulation, which takes into account transitions between states. We also performed in vitro binding experiments with sequences containing ^ΨC and furthermore characterized the structure of the formed triplex using molecular dynamics simulation. We found that the neighboring methylated or protonated cytidine promotes the formation of the favorable tautomer, whereas the neighboring thymine or locked nucleic acid has a poor effect, and consecutive ^ΨC has a negative influence. The deleterious effect of consecutive ^ΨC in a triplex formation was confirmed using in vitro binding experiments. Our findings contribute to improving the design of ^ΨC-containing triplex-forming oligonucleotides directed to target G-rich DNA sequences.

Sequence dependency of canonical base pair opening in the DNA double helix

The flipping-out of a DNA base from the double helical structure is a key step of many cellular processes, such as DNA replication, modification and repair. Base pair opening is the first step of base flipping and the exact mechanism is still not well understood. We investigate sequence effects on base pair opening using extensive classical molecular dynamics simulations targeting the opening of 11 different canonical base pairs in two DNA sequences. Two popular biomolecular force fields are applied. To enhance sampling and calculate free energies, we bias the simulation along a simple distance coordinate using a newly developed adaptive sampling algorithm. The simulation is guided back and forth along the coordinate, allowing for multiple opening pathways. We compare the calculated free energies with those from an NMR study and check assumptions of the model used for interpreting the NMR data. Our results further show that the neighboring sequence is an important factor for the opening free energy, but also indicates that other sequence effects may play a role. All base pairs are observed to have a propensity for opening toward the major groove. The preferred opening base is cytosine for GC base pairs, while for AT there is sequence dependent competition between the two bases. For AT opening, we identify two non-canonical base pair interactions contributing to a local minimum in the free energy profile. For both AT and CG we observe long-lived interactions with water and with sodium ions at specific sites on the open base pair.

The DNA double helix, a molecule that stores biological information, has become an iconic image of biomedical research. In order to use or repair the information it carries, the bases that are stacked in the helix need to be chemically exposed. This can happen either by separating the two strands in the helix or by flipping out individual bases. Here, we focus on the latter process. Usually proteins are involved in interactions with bases, but it is still unclear if bases are pulled out actively by proteins or if they act on spontaneously flipped bases. Although experiments can detect base pair opening, it is difficult to detect which base moves in which direction. Here, we present results from molecular dynamics simulations using a recently developed sampling method which improves the statistics in the simulations by enhancing the probability of the base pair opening event. We observe differences in probability, modes and mechanism of opening that depend not only on the types of the bases in the pair, but also strongly on their neighbors. This provides essential information for understanding how DNA functions.

Next-generation bis-locked nucleic acids with stacking linker and 2'-glycylamino-LNA show enhanced DNA invasion into supercoiled duplexes

Targeting and invading double-stranded DNA with synthetic oligonucleotides under physiological conditions remain a challenge. Bis-locked nucleic acids (bisLNAs) are clamp-forming oligonucleotides able to invade into supercoiled DNA via combined Hoogsteen and Watson–Crick binding. To improve the bisLNA design, we investigated its mechanism of binding. Our results suggest that bisLNAs bind via Hoogsteen-arm first, followed by Watson–Crick arm invasion, initiated at the tail. Based on this proposed hybridization mechanism, we designed next-generation bisLNAs with a novel linker able to stack to adjacent nucleobases, a new strategy previously not applied for any type of clamp-constructs. Although the Hoogsteen-arm limits the invasion, upon incorporation of the stacking linker, bisLNA invasion is significantly more efficient than for non-clamp, or nucleotide-linker containing LNA-constructs. Further improvements were obtained by substituting LNA with 2′-glycylamino-LNA, contributing a positive charge. For regular bisLNAs a 14-nt tail significantly enhances invasion. However, when two stacking linkers were incorporated, tail-less bisLNAs were able to efficiently invade. Finally, successful targeting of plasmids inside bacteria clearly demonstrates that strand invasion can take place in a biologically relevant context.

Sharp switches between regular and swinger mitochondrial replication: 16S rDNA systematically exchanging nucleotides A<->T+C<->G in the mitogenome of Kamimuria wangi

Swinger DNAs are sequences whose homology with known sequences is detected only by assuming systematic exchanges between nucleotides. Nine symmetric (X<->Y, i.e. A<->C) and fourteen asymmetric (X->Y->Z, i.e. A->C->G) exchanges exist. All swinger DNA previously detected in GenBank follow the A<->T+C<->G exchange, while mitochondrial swinger RNAs distribute among different swinger types. Here different alignment criteria detect 87 additional swinger mitochondrial DNAs (86 from insects), including the first swinger gene embedded within a complete genome, corresponding to the mitochondrial 16S rDNA of the stonefly Kamimuria wangi. Other Kamimuria mt genome regions are "regular", stressing unanswered questions on (a) swinger polymerization regulation; (b) swinger 16S rDNA functions; and (c) specificity to rDNA, in particular 16S rDNA. Sharp switches between regular and swinger replication, together with previous observations on swinger transcription, suggest that swinger replication might be due to a switch in polymerization mode of regular polymerases and the possibility of swinger-encoded information, predicted in primordial genes such as rDNA.