E-motif formed by extrahelical cytosine bases in DNA homoduplexes of trinucleotide and hexanucleotide repeats

Structure and Dynamics of DNA and RNA Double Helices Formed by d(CTG), d(GTC), r(CUG), and r(GUC) Trinucleotide Repeats and Associated DNA-RNA Hybrids

Myotonic dystrophy type 1 is the most frequent form of muscular dystrophy in adults caused by an abnormal expansion of the CTG trinucleotide. Both the expanded DNA and the expanded CUG RNA transcript can fold into hairpins. Co-transcriptional formation of stable RNA·DNA hybrids can also enhance the instability of repeat tracts. We performed molecular dynamics simulations of homoduplexes associated with the disease, d(CTG)n and r(CUG)n, and their corresponding r(CAG)n:d(CTG)n and r(CUG)n:d(CAG)n hybrids that can form under bidirectional transcription and of non-pathological d(GTC)n and d(GUC)n homoduplexes. We characterized their conformations, stability, and dynamics and found that the U·U and T·T mismatches are dynamic, favoring anti-anti conformations inside the helical core, followed by anti-syn and syn-syn conformations. For DNA, the secondary minima in the non-expanding d(GTC)n helices are deeper, wider, and longer-lived than those in d(CTG)n, which constitutes another biophysical factor further differentiating the expanding and non-expanding sequences. The hybrid helices are closer to A-RNA, with the A-T and A-U pairs forming two stable Watson-Crick hydrogen bonds. The neutralizing ion distribution around the non-canonical pairs is also described.

Frustration Between Preferred States of Complementary Trinucleotide Repeat DNA Hairpins Anticorrelates with Expansion Disease Propensity

DNA trinucleotide repeat (TRs) expansion beyond a threshold often results in human neurodegenerative diseases. The mechanisms causing expansions remain unknown, although the tendency of TR ssDNA to self-associate into hairpins that slip along their length is widely presumed related. Here we apply single molecule FRET (smFRET) experiments and molecular dynamics simulations to determine conformational stabilities and slipping dynamics for CAG, CTG, GAC and GTC hairpins. Tetraloops are favored in CAG (89%), CTG (89%) and GTC (69%) while GAC favors triloops. We also determined that TTG interrupts near the loop in the CTG hairpin stabilize the hairpin against slipping. The different loop stabilities have implications for intermediate structures that may form when TR-containing duplex DNA opens. Opposing hairpins in the (CAG) ∙ (CTG) duplex would have matched stability whereas opposing hairpins in a (GAC) ∙ (GTC) duplex would have unmatched stability, introducing frustration in the (GAC) ∙ (GTC) opposing hairpins that could encourage their resolution to duplex DNA more rapidly than in (CAG) ∙ (CTG) structures. Given that the CAG and CTG TR can undergo large, disease-related expansion whereas the GAC and GTC sequences do not, these stability differences can inform and constrain models of expansion mechanisms of TR regions.

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.

Revisiting recent unusual drug-DNA complex structures: Implications for cancer and neurological disease diagnostics and therapeutics

DNA plays a crucial role in various biological processes such as protein production, replication, recombination etc. by adopting different conformations. Targeting these conformations by small molecules is not only important for disease therapy, but also improves our understanding of the mechanisms of disease development. In this review, we provide an overview of some of the most recent ligand-DNA complexes that have diagnostic and therapeutic applications in neurological diseases caused by abnormal repeat expansions and in cancer associated with mismatches. In addition, we have discussed important implications of ligands targeting higher-order structures, such as four-way junctions, G-quadruplexes and triplexes for drug discovery and DNA nanotechnology. We provide an overview of the results and perspectives of such structural studies on ligand-DNA interactions.

RNA as a Major-Groove Ligand: RNA-RNA and RNA-DNA Triplexes Formed by GAA and UUC or TTC Sequences

Friedreich's ataxia is associated with noncanonical nucleic acid structures that emerge when GAA:TTC repeats in the first intron of the FXN gene expand beyond a critical number of repeats. Specifically, the noncanonical repeats are associated with both triplexes and R-loops. Here, we present an in silico investigation of all possible triplexes that form by attaching a third RNA strand to an RNA:RNA or DNA:DNA duplex, complementing previous DNA-based triplex studies. For both new triplexes results are similar. For a pyridimine UUC⁺ third strand, the parallel orientation is stable while its antiparallel counterpart is unstable. For a neutral GAA third strand, the parallel conformation is stable. A protonated GA⁺A third strand is stable in both parallel and antiparallel orientations. We have also investigated Na⁺ and Mg²⁺ ion distributions around the triplexes. The presence of Mg²⁺ ions helps stabilize neutral, antiparallel GAA triplexes. These results (along with previous DNA-based studies) allow for the emergence of a complete picture of the stability and structural characteristics of triplexes based on the GAA and TTC/UUC sequences, thereby contributing to the field of trinucleotide repeats and the associated unusual structures that trigger expansion.

Staggered intercalation of DNA duplexes with base-pair modulation by two distinct drug molecules induces asymmetric backbone twisting and structure polymorphism

The use of multiple drugs simultaneously targeting DNA is a promising strategy in cancer therapy for potentially overcoming single drug resistance. In support of this concept, we report that a combination of actinomycin D (ActD) and echinomycin (Echi), can interact in novel ways with native and mismatched DNA sequences, distinct from the structural effects produced by either drug alone. Changes in the former with GpC and CpG steps separated by a A:G or G:A mismatch or in a native DNA with canonical G:C and C:G base pairs, result in significant asymmetric backbone twists through staggered intercalation and base pair modulations. A wobble or Watson–Crick base pair at the two drug-binding interfaces can result in a single-stranded ‘chair-shaped’ DNA duplex with a straight helical axis. However, a novel sugar-edged hydrogen bonding geometry in the G:A mismatch leads to a ‘curved-shaped’ duplex. Two non-canonical G:C Hoogsteen base pairings produce a sharply kinked duplex in different forms and a four-way junction-like superstructure, respectively. Therefore, single base pair modulations on the two drug-binding interfaces could significantly affect global DNA structure. These structures thus provide a rationale for atypical DNA recognition via multiple DNA intercalators and a structural basis for the drugs’ potential synergetic use.

Novel eGZ-motif formed by regularly extruded guanine bases in a left-handed Z-DNA helix as a major motif behind CGG trinucleotide repeats

The expansion of d(CGG) trinucleotide repeats (TRs) lies behind several important neurodegenerative diseases. Atypical DNA secondary structures have been shown to trigger TR expansion: their characterization is important for a molecular understanding of TR disease. CD spectroscopy experiments in the last decade have unequivocally demonstrated that CGG runs adopt a left-handed Z-DNA conformation, whose features remain uncertain because it entails accommodating GG mismatches. In order to find this missing motif, we have carried out molecular dynamics (MD) simulations to explore all the possible Z-DNA helices that potentially form after the transition from B- to Z-DNA. Such helices combine either CpG or GpC Watson-Crick steps in Z-DNA form with GG-mismatch conformations set as either intrahelical or extrahelical; and participating in BZ or ZZ junctions or in alternately extruded conformations. Characterization of the stability and structural features (especially overall left-handedness, higher-temperature and steered MD simulations) identified two novel Z-DNA helices: the most stable one displays alternately extruded Gs, and is followed by a helix with symmetrically extruded ZZ junctions. The G-extrusion favors a seamless stacking of the Watson-Crick base pairs; extruded Gs favor syn conformations and display hydrogen-bonding and stacking interactions. Such conformations could have the potential to hijack the MMR complex, thus triggering further expansion.

Selectively recognizing extrahelical conformations of DNA trinucleotide repeats by a hydroxylated porphyrin ligand

Trinucleotide repeats (TRs) with abnormal lengths and atypical folding are implicated in various neurodegenerative diseases. The least stable cytosine-cytosine (C-C) mismatches in TRs when structuring into homoduplexes/hairpins have more chance in certain sequence contexts to preferentially adopt an extrahelical (E-motif) conformation with respect to those in polarity-inverted intrahelical counterparts. Herein, we designed a trihydroxyphenyl porphyrin ligand (POH3) to meet the challenge towards resolving the E-motif conformation. POH3 exhibited a specific 2:1 binding with DNAs adopting the E-motif cytosine conformation, independent of the TRs length. The trihydroxyl pattern was very crucial to gain the E-motif selectivity over the polarity-inverted counterparts via the complementary hydrogen bonding that occurred in the minor groove. Our work first elucidates the rationale in designing ligands to selectively resolve the E-motif nucleotides within TRs.

Construction of DNA/RNA Triplex Helices Based on GAA/TTC Trinucleotide Repeats

Atypical DNA and RNA secondary structures play a crucial role in simple sequence repeat (SSR) diseases, which are associated with a class of neurological and neuromuscular disorders known as "anticipation diseases," where the age of disease onset decreases and the severity of the disease is increased as the intergenerational expansion of the SSR increases. While the mechanisms underlying these diseases are complex and remain elusive, there is a consensus that stable, non-B-DNA atypical secondary structures play an important - if not causative - role. These structures include single-stranded DNA loops and hairpins, G-quartets, Z-DNA, triplex nucleic acid structures, and others. While all of these structures are of interest, structures based on nucleic acid triplexes have recently garnered increased attention as they have been implicated in gene regulation, gene repair, and gene engineering. Our work here focuses on the construction of DNA triplexes and RNA/DNA hybrids formed from GAA/TTC trinucleotide repeats, which underlie Friedreich's ataxia. While there is some software, such as the Discovery Studio Visualizer, that can aid in the initial construction of DNA triple helices, the only option for the triple helix is constrained to be that of an antiparallel pyrimidine for the third strand. In this protocol, we illustrate how to build up more generalized DNA triplexes and DNA/RNA mixed hybrids. We make use of both the Discovery Studio Visualizer and the AMBER simulation package to construct the initial triplexes. Using the steps outlined here, one can - in principle - build up any triple nucleic acid helix with a desired sequence for large-scale molecular dynamics simulation studies.

Molecular conformations and dynamics of nucleotide repeats associated with neurodegenerative diseases: double helices and CAG hairpin loops

Pathogenic DNA secondary structures have been identified as a common and causative factor for expansion in trinucleotide, hexanucleotide, and other simple sequence repeats. These expansions underlie about fifty neurological and neuromuscular disorders known as “anticipation diseases”. Cell toxicity and death have been linked to the pathogenic conformations and functional changes of the RNA transcripts, of DNA itself and, when trinucleotides are present in exons, of the translated proteins. We review some of our results for the conformations and dynamics of pathogenic structures for both RNA and DNA, which include mismatched homoduplexes formed by trinucleotide repeats CAG and GAC; CCG and CGG; CTG(CUG) and GTC(GUC); the dynamics of DNA CAG hairpins; mismatched homoduplexes formed by hexanucleotide repeats (GGGGCC) and (GGCCCC); and G-quadruplexes formed by (GGGGCC) and (GGGCCT). We also discuss the dynamics of strand slippage in DNA hairpins formed by CAG repeats as observed with single-molecule Fluorescence Resonance Energy Transfer. This review focuses on the rich behavior exhibited by the mismatches associated with these simple sequence repeat noncanonical structures.

Atypical structures of GAA/TTC trinucleotide repeats underlying Friedreich's ataxia: DNA triplexes and RNA/DNA hybrids

Expansion of the GAA/TTC repeats in the first intron of the FXN gene causes Friedreich's ataxia. Non-canonical structures are linked to this expansion. DNA triplexes and R-loops are believed to arrest transcription, which results in frataxin deficiency and eventual neurodegeneration. We present a systematic in silico characterization of the possible DNA triplexes that could be assembled with GAA and TTC strands; the two hybrid duplexes [r(GAA):d(TTC) and d(GAA):r(UUC)] in an R-loop; and three hybrid triplexes that could form during bidirectional transcription when the non-template DNA strand bonds with the hybrid duplex (collapsed R-loops, where the two DNA strands remain antiparallel). For both Y·R:Y and R·R:Y DNA triplexes, the parallel third strand orientation is more stable; both parallel and antiparallel protonated d(GA+A)·d(GAA):d(TTC) triplexes are stable. Apparent contradictions in the literature about the R·R:Y triplex stability is probably due to lack of molecular resolution, since shifting the third strand by a single nucleotide alters the stability ranking. In the collapsed R-loops, antiparallel d(TTC+)·d(GAA):r(UUC) is unstable, while parallel d(GAA)·r(GAA):d(TTC) and d(GA+A)·r(GAA):d(TTC) are stable. In addition to providing new structural perspectives for specific therapeutic aims, our results contribute to a systematic structural basis for the emerging field of quantitative R-loop biology.

Structural Basis for Targeting T:T Mismatch with Triaminotriazine-Acridine Conjugate Induces a U-Shaped Head-to-Head Four-Way Junction in CTG Repeat DNA

The potent DNA-binding compound triaminotriazine-acridine conjugate (Z1) functions by targeting T:T mismatches in CTG trinucleotide repeats that are responsible for causing neurological diseases such as myotonic dystrophy type 1, but its binding mechanism remains unclear. We solved a crystal structure of Z1 in a complex with DNA containing three consecutive CTG repeats with three T:T mismatches. Crystallographic studies revealed that direct intercalation of two Z1 molecules at both ends of the CTG repeat induces thymine base flipping and DNA backbone deformation to form a four-way junction. The core of the complex unexpectedly adopts a U-shaped head-to-head topology to form a crossover of each chain at the junction site. The crossover junction is held together by two stacked G:C pairs at the central core that rotate with respect to each other in an X-shape to form two nonplanar minor-groove-aligned G·C·G·C tetrads. Two stacked G:C pairs on both sides of the center core are involved in the formation of pseudo-continuous duplex DNA. Four metal-mediated base pairs are observed between the N7 atoms of G and Co^II, an interaction that strongly preserves the central junction site. Beyond revealing a new type of ligand-induced, four-way junction, these observations enhance our understanding of the specific supramolecular chemistry of Z1 that is essential for the formation of a noncanonical DNA superstructure. The structural features described here serve as a foundation for the design of new sequence-specific ligands targeting mismatches in the repeat-associated structures.

New Roles for Canonical Transcription Factors in Repeat Expansion Diseases

The presence of microsatellite repeat expansions within genes is associated with >30 neurological diseases. Of interest, (GGGGCC)_>30-repeats within C9orf72 are associated with amyotrophic lateral sclerosis and frontotemporal dementia (ALS/FTD). These expansions can be 100s to 1000s of units long. Thus, it is perplexing how RNA-polymerase II (RNAPII) can successfully transcribe them. Recent investigations focusing on GGGGCC-transcription have identified specific, canonical complexes that may promote RNAPII-transcription at these GC-rich microsatellites: the DSIF complex and PAF1C. These complexes may be important for resolving the unique secondary structures formed by GGGGCC-DNA during transcription. Importantly, this process can produce potentially toxic repeat-containing RNA that can encode potentially toxic peptides, impacting neuron function and health. Understanding how transcription of these repeats occurs has implications for therapeutics in multiple diseases.

Targeted Oligonucleotides for Treating Neurodegenerative Tandem Repeat Diseases

Nucleotide repeat disorders encompass more than 30 diseases, most of which show dominant inheritance, such as Huntington's disease, spinocerebellar ataxias, and myotonic dystrophies. Yet others, including Friedreich's ataxia, are recessively inherited. A common feature is the presence of a DNA tandem repeat in the disease-associated gene and the propensity of the repeats to expand in germ and in somatic cells, with ensuing neurological and frequently also neuromuscular defects. Repeat expansion is the most frequent event in these diseases; however, sequence contractions, deletions, and mutations have also been reported. Nucleotide repeat sequences are predisposed to adopt non-B-DNA conformations, such as hairpins, cruciform, and intramolecular triple-helix structures (triplexes), also known as H-DNA. For gain-of-function disorders, oligonucleotides can be used to target either transcripts or duplex DNA and in diseases with recessive inheritance oligonucleotides may be used to alter repressive DNA or RNA conformations. Most current treatment strategies are aimed at altering transcript levels, but therapies directed against DNA are also emerging, and novel strategies targeting DNA, instead of RNA, are described. Different mechanisms using modified oligonucleotides are discussed along with the structural aspects of repeat sequences, which can influence binding modes and efficiencies.

CoII(Chromomycin)₂ Complex Induces a Conformational Change of CCG Repeats from i-Motif to Base-Extruded DNA Duplex

We have reported the propensity of a DNA sequence containing CCG repeats to form a stable i-motif tetraplex structure in the absence of ligands. Here we show that an i-motif DNA sequence may transition to a base-extruded duplex structure with a GGCC tetranucleotide tract when bound to the (Co^II)-mediated dimer of chromomycin A3, Co^II(Chro)₂. Biophysical experiments reveal that CCG trinucleotide repeats provide favorable binding sites for Co^II(Chro)₂. In addition, water hydration and divalent metal ion (Co^II) interactions also play a crucial role in the stabilization of CCG trinucleotide repeats (TNRs). Our data furnish useful structural information for the design of novel therapeutic strategies to treat neurological diseases caused by repeat expansions.

Structural and Dynamical Characterization of DNA and RNA Quadruplexes Obtained from the GGGGCC and GGGCCT Hexanucleotide Repeats Associated with C9FTD/ALS and SCA36 Diseases

A (GGGGCC) hexanucleotide repeat (HR) expansion in the C9ORF72 gene has been considered the major cause behind both frontotemporal dementia and amyotrophic lateral sclerosis, while a (GGGCCT) is associated with spinocerebellar ataxia 36. Recent experiments involving NMR, CD, optical melting and 1D 1H NMR spectroscopy, suggest that the r(GGGGCC) HR can adopt a hairpin structure with G-G mismatches in equilibrium with a G-quadruplex structure. G-Quadruplexes have also been identified for d(GGGGCC). As these experiments lack molecular resolution, we have used molecular dynamics microsecond simulations to obtain a structural characterization of the G-quadruplexes associated with both HRs. All DNA G-quadruplexes, parallel or antiparallel, with or without loops are stable, while only parallel and one antiparallel (stabilized by diagonal loops) RNA G-quadruplexes are stable. It is known that antiparallel G-quadruplexes require alternating guanines to be in a syn conformation that is hindered by the C3'-endo pucker preferred by RNA. Initial RNA antiparallel quadruplexes built with C2'-endo sugars evolve such that the transition (C2'-endo)-to-(C3'-endo) triggers unwinding and buckling of the flat G-tetrads, resulting in the unfolding of the RNA antiparallel quadruplex. Finally, a parallel G-quadruplex stabilizes an adjacent C-tetrad in both DNA and RNA (thus effectively becoming a mixed quadruplex of 5 layers). The C-tetrad is stabilized by the stacking interactions with the preceding G-tetrad, by cyclical hydrogen bonds C(N4)-(O2), and by an ion between the G-tetrad and the C-tetrad. In addition, antiparallel DNA G-quadruplexes also stabilize flat C-layers at the ends of the quadruplexes.

Structure and Dynamics of DNA and RNA Double Helices Obtained from the CCG and GGC Trinucleotide Repeats

Expansions of both GGC and CCG sequences lead to a number of expandable, trinucleotide repeat (TR) neurodegenerative diseases. Understanding of these diseases involves, among other things, the structural characterization of the atypical DNA and RNA secondary structures. We have performed molecular dynamics simulations of (GCC) _n and (GGC) _n homoduplexes in order to characterize their conformations, stability, and dynamics. Each TR has two reading frames, which results in eight nonequivalent RNA/DNA homoduplexes, characterized by CpG or GpC steps between the Watson-Crick base pairs. Free energy maps for the eight homoduplexes indicate that the C-mismatches prefer anti-anti conformations, while G-mismatches prefer anti-syn conformations. Comparison between three modifications of the DNA AMBER force field shows good agreement for the mismatch free energy maps. The mismatches in DNA-GCC (but not CCG) are extrahelical, forming an extended e-motif. The mismatched duplexes exhibit characteristic sequence-dependent step twist, with strong variations in the G-rich sequences and the e-motif. The distribution of Na⁺ is highly localized around the mismatches, especially G-mismatches. In the e-motif, there is strong Na⁺ binding by two G(N7) atoms belonging to the pseudo GpC step created when cytosines are extruded and by extrahelical cytosines. Finally, we used a novel technique based on fast melting by means of an infrared laser pulse to classify the relative stability of the different DNA-CCG and -GGC homoduplexes.