Atomic resolution structure of CAG RNA repeats: structural insights and implications for the trinucleotide repeat expansion diseases

CircHTT(2,3,4,5,6) - co-evolving with the HTT CAG-repeat tract - modulates Huntington's disease phenotypes

Circular RNA (circRNA) molecules have critical functions during brain development and in brain-related disorders. Here, we identified and validated a circRNA, circHTT(2,3,4,5,6), stemming from the Huntington's disease (HD) gene locus that is most abundant in the central nervous system (CNS). We uncovered its evolutionary conservation in diverse mammalian species, and a correlation between circHTT(2,3,4,5,6) levels and the length of the CAG-repeat tract in exon-1 of HTT in human and mouse HD model systems. The mouse orthologue, circHtt(2,3,4,5,6), is expressed during embryogenesis, increases during nervous system development, and is aberrantly upregulated in the presence of the expanded CAG tract. While an IRES-like motif was predicted in circH TT (2,3,4,5,6), the circRNA does not appear to be translated in adult mouse brain tissue. Nonetheless, a modest, but consistent fraction of circHtt(2,3,4,5,6) associates with the 40S ribosomal subunit, suggesting a possible role in the regulation of protein translation. Finally, circHtt(2,3,4,5,6) overexpression experiments in HD-relevant STHdh striatal cells revealed its ability to modulate CAG expansion-driven cellular defects in cell-to-substrate adhesion, thus uncovering an unconventional modifier of HD pathology.

Stick-slip unfolding favors self-association of expanded HTT mRNA

In Huntington's Disease (HD) and related disorders, expansion of CAG trinucleotide repeats produces a toxic gain of function in affected neurons. Expanded huntingtin (expHTT) mRNA forms aggregates that sequester essential RNA binding proteins, dysregulating mRNA processing and translation. The physical basis of RNA aggregation has been difficult to disentangle owing to the heterogeneous structure of the CAG repeats. Here, we probe the folding and unfolding pathways of expHTT mRNA using single-molecule force spectroscopy. Whereas normal HTT mRNAs unfold reversibly and cooperatively, expHTT mRNAs with 20 or 40 CAG repeats slip and unravel non-cooperatively at low tension. Slippage of CAG base pairs is punctuated by concerted rearrangement of adjacent CCG trinucleotides, trapping partially folded structures that readily base pair with another RNA strand. We suggest that the conformational entropy of the CAG repeats, combined with stable CCG base pairs, creates a stick-slip behavior that explains the aggregation propensity of expHTT mRNA.

Binding of small molecule inhibitors to RNA polymerase-Spt5 complex impacts RNA and DNA stability

Spt5 is an elongation factor that associates with RNA polymerase II (Pol II) during transcription and has important functions in promoter-proximal pausing and elongation processivity. Spt5 was also recognized for its roles in the transcription of expanded-repeat genes that are related to neurodegenerative diseases. Recently, a set of Spt5-Pol II small molecule inhibitors (SPIs) were reported, which selectively inhibit mutant huntingtin gene transcription. Inhibition mechanisms as well as interaction sites of these SPIs with Pol II and Spt5 are not entirely known. In this study, we predicted the binding sites of three selected SPIs at the Pol II-Spt5 interface by docking and molecular dynamics simulations. Two molecules out of three demonstrated strong binding with Spt5 and Pol II, while the other molecule was more loosely bound and sampled multiple binding sites. Strongly bound SPIs indirectly affected RNA and DNA dynamics at the exit site as DNA became more flexible while RNA was stabilized by increased interactions with Spt5. Our results suggest that the transcription inhibition mechanism induced by SPIs can be related to Spt5-nucleic acid interactions, which were altered to some extent with strong binding of SPIs.

Between Order and Chaos: Understanding the Mechanism and Pathology of RAN Translation

Repeat-associated non-AUG (RAN) translation is a pathogenic mechanism in which repetitive sequences are translated into aggregation-prone proteins from multiple reading frames, even without a canonical AUG start codon. Since its discovery in spinocerebellar ataxia type 8 (SCA8) and myotonic dystrophy type 1 (DM1), RAN translation is now known to occur in the context of 12 disease-linked repeat expansions. This review discusses recent advances in understanding the regulatory mechanisms controlling RAN translation and its contribution to the pathophysiology of repeat expansion diseases. We discuss the key findings in the context of Fragile X Tremor Ataxia Syndrome (FXTAS), a neurodegenerative disorder caused by a CGG repeat expansion in the 5' untranslated region of FMR1.

Uncovering the mechanism for aggregation in repeat expanded RNA reveals a reentrant transition

RNA molecules aggregate under certain conditions. The resulting condensates are implicated in human neurological disorders, and can potentially be designed towards specified bulk properties in vitro. However, the mechanism for aggregation-including how aggregation properties change with sequence and environmental conditions-remains poorly understood. To address this challenge, we introduce an analytical framework based on multimer enumeration. Our approach reveals the driving force for aggregation to be the increased configurational entropy associated with the multiplicity of ways to form bonds in the aggregate. Our model uncovers rich phase behavior, including a sequence-dependent reentrant phase transition, and repeat parity-dependent aggregation. We validate our results by comparison to a complete computational enumeration of the landscape, and to previously published molecular dynamics simulations. Our work unifies and extends published results, both explaining the behavior of CAG-repeat RNA aggregates implicated in Huntington's disease, and enabling the rational design of programmable RNA condensates.

Molecular insights into the interaction of CAG trinucleotide RNA repeats with nucleolin and its implication in polyglutamine diseases

Polyglutamine (polyQ) diseases are a type of inherited neurodegenerative disorders caused by cytosine-adenine-guanine (CAG) trinucleotide expansion within the coding region of the disease-associated genes. We previously demonstrated that a pathogenic interaction between expanded CAG RNA and the nucleolin (NCL) protein triggers the nucleolar stress and neuronal cell death in polyQ diseases. However, mechanisms behind the molecular interaction remain unknown. Here, we report a 1.45 Å crystal structure of the r(CAG)5 oligo that comprises a full A'-form helical turn with widened grooves. Based on this structure, we simulated a model of r(CAG)5 RNA complexed with the RNA recognition motif 2 (RRM2) of NCL and identified NCL residues that are critical for its binding to CAG RNA. Combined with in vitro and in vivo site-directed mutagenesis studies, our model reveals that CAG RNA binds to NCL sites that are not important for other cellular functions like gene expression and rRNA synthesis regulation, indicating that toxic CAG RNA interferes with NCL functions by sequestering it. Accordingly, an NCL mutant that is aberrant in CAG RNA-binding could rescue RNA-induced cytotoxicity effectively. Taken together, our study provides new molecular insights into the pathogenic mechanism of polyQ diseases mediated by NCL-CAG RNA interaction.

A Structural Potential of Rare Trinucleotide Repeat Tracts in RNA

Among types of trinucleotide repeats, there is some disproportion in the frequency of their occurrence in the human exome. This research presents new data describing the folding and thermodynamic stability of short, tandem RNA repeats of 23 types, focusing on the rare, yet poorly analyzed ones. UV-melting experiments included the presence of PEG or potassium and magnesium ions to determine their effect on the stability of RNA repeats structures. Rare repeats predominantly stayed single-stranded but had the potential for base pairing with other partially complementary repeat tracts. A coexistence of suitably complementary repeat types in a single RNA creates opportunities for interaction in the context of the secondary structure of RNA. We searched the human transcriptome for model RNAs in which different, particularly rare trinucleotide repeats coexist and selected the GABRA4 and CHIC1 RNAs to study intramolecular interactions between the repeat tracts that they contain. In vitro secondary structure probing results showed that the UAA and UUG repeat tracts, present in GABRA4 3' UTR, form a double helix, which separates one of its structural domains. For the RNA CHIC1 ORF fragment containing four short AGG repeat tracts and the CGU tract, we proved the formation of quadruplexes that blocked reverse transcription.

Condensates in RNA repeat sequences are heterogeneously organized and exhibit reptation dynamics

Although it is known that RNA undergoes liquid-liquid phase separation, the interplay between the molecular driving forces and the emergent features of the condensates, such as their morphologies and dynamic properties, is not well understood. We introduce a coarse-grained model to simulate phase separation of trinucleotide repeat RNAs, which are implicated in neurological disorders. After establishing that the simulations reproduce key experimental findings, we show that once recruited inside the liquid droplets, the monomers transition from hairpin-like structures to extended states. Interactions between the monomers in the condensates result in the formation of an intricate and dense intermolecular network, which severely restrains the fluctuations and mobilities of the RNAs inside large droplets. In the largest densely packed high-viscosity droplets, the mobility of RNA chains is best characterized by reptation, reminiscent of the dynamics in polymer melts. Our work provides a microscopic framework for understanding liquid-liquid phase separation in RNA, which is not easily discernible in current experiments.

RNA Toxicity and Perturbation of rRNA Processing in Spinocerebellar Ataxia Type 2

BACKGROUND

Spinocerebellar ataxia type 2 (SCA2) is a neurodegenerative disease caused by expansion of a CAG repeat in Ataxin-2 (ATXN2) gene. The mutant ATXN2 protein with a polyglutamine tract is known to be toxic and contributes to the SCA2 pathogenesis.

OBJECTIVE

Here, we tested the hypothesis that the mutant ATXN2 transcript with an expanded CAG repeat (expATXN2) is also toxic and contributes to SCA2 pathogenesis.

METHODS

The toxic effect of expATXN2 transcripts on SK-N-MC neuroblastoma cells and primary mouse cortical neurons was evaluated by caspase 3/7 activity and nuclear condensation assay, respectively. RNA immunoprecipitation assay was performed to identify RNA binding proteins (RBPs) that bind to expATXN2 RNA. Quantitative PCR was used to examine if ribosomal RNA (rRNA) processing is disrupted in SCA2 and Huntington's disease (HD) human brain tissue.

RESULTS

expATXN2 RNA induces neuronal cell death, and aberrantly interacts with RBPs involved in RNA metabolism. One of the RBPs, transducin β-like protein 3 (TBL3), involved in rRNA processing, binds to both expATXN2 and expanded huntingtin (expHTT) RNA in vitro. rRNA processing is disrupted in both SCA2 and HD human brain tissue.

CONCLUSION

These findings provide the first evidence of a contributory role of expATXN2 transcripts in SCA2 pathogenesis, and further support the role of expHTT transcripts in HD pathogenesis. The disruption of rRNA processing, mediated by aberrant interaction of RBPs with expATXN2 and expHTT transcripts, suggest a point of convergence in the pathogeneses of repeat expansion diseases with potential therapeutic implications. © 2021 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

The MID1 Protein: A Promising Therapeutic Target in Huntington's Disease

Huntington's disease (HD) is caused by an expansion mutation of a CAG repeat in exon 1 of the huntingtin (HTT) gene, that encodes an expanded polyglutamine tract in the HTT protein. HD is characterized by progressive psychiatric and cognitive symptoms associated with a progressive movement disorder. HTT is ubiquitously expressed, but the pathological changes caused by the mutation are most prominent in the central nervous system. Since the mutation was discovered, research has mainly focused on the mutant HTT protein. But what if the polyglutamine protein is not the only cause of the neurotoxicity? Recent studies show that the mutant RNA transcript is also involved in cellular dysfunction. Here we discuss the abnormal interaction of the mutant HTT transcript with a protein complex containing the MID1 protein. MID1 aberrantly binds to CAG repeats and this binding increases with CAG repeat length. Since MID1 is a translation regulator, association of the MID1 complex stimulates translation of mutant HTT mRNA, resulting in an overproduction of polyglutamine protein. Thus, blocking the interaction between MID1 and mutant HTT mRNA is a promising therapeutic approach. Additionally, we show that MID1 expression in the brain of both HD patients and HD mice is aberrantly increased. This finding further supports the concept of blocking the interaction between MID1 and mutant HTT mRNA to counteract mutant HTT translation as a valuable therapeutic strategy. In line, recent studies in which either compounds affecting the assembly of the MID1 complex or molecules targeting HTT RNA, show promising results.

Conformational distortions induced by periodically recurring A…A in d(CAG).d(CAG) provide stereochemical rationale for the trapping of MSH2.MSH3 in polyQ disorders

CAG repeat instability causes a number of neurodegenerative disorders. The unusual hairpin stem structure formed by the CAG repeats in DNA traps the human mismatch repair MSH2.MSH3 (Mutsβ) complex. To understand the mechanism behind the abnormal binding of Mutsβ with the imperfect hairpin stem structure formed by CAG repeats, molecular dynamics simulations have been carried out for Mutsβ-d(CAG)₂(CAG)(CAG)₂.d(CTG)₂(CAG)(CTG)₂ (1 A…A mismatch) and Mutsβ-d(CAG)₅.d(CAG)₅ (5 mismatches, wherein, A…A occurs periodically) complexes. The interaction of MSH3 residue Tyr₂₄₅ at the minor groove side of A…A, an essential interaction responsible for the recognition by Mutsβ, are retained in both the cases. Nevertheless, the periodic unwinding caused by the nonisostericity of A…A with the flanking canonical base pairs in d(CAG)₅.d(CAG)₅ distorts the regular B-form geometry. Such an unwinding exposes one of the A…A mismatches (that interacts with Tyr₂₄₅) at the major groove side and also facilitates the on and off hydrogen bonding interaction with Lys₅₄₆ sidechain (MSH2-domain-IV). In contrast, kinking of the DNA towards the major groove in Mutsβ-d(CAG)₂(CAG)(CAG)₂.d(CTG)₂(CAG)(CTG)₂ doesn’t facilitate such an exposure of the bases at the major groove. Further, the unwinding of the helix in d(CAG)₅.d(CAG)₅ enhances the tighter binding between MSH2-domain-I and d(CAG)₅.d(CAG)₅ at the major groove side as well as between MSH3-domain-I and MSH3-domain-IV. Markedly, such enhanced interactions are absent in Mutsβ-d(CAG)₂(CAG)(CAG)₂.d(CTG)₂(CAG)(CTG)₂ that has a single A…A mismatch. Thus, the above-mentioned enhancement in intra- and inter- molecular interactions in Mutsβ-d(CAG)₅.d(CAG)₅ provide the stereochemical rationale for the trapping of Mutsβ in CAG repeat expansion disorders.

Differential efficacies of Cas nucleases on microsatellites involved in human disorders and associated off-target mutations

Microsatellite expansions are the cause of >20 neurological or developmental human disorders. Shortening expanded repeats using specific DNA endonucleases may be envisioned as a gene editing approach. Here, we measured the efficacy of several CRISPR-Cas nucleases to induce recombination within disease-related microsatellites, in Saccharomyces cerevisiae. Broad variations in nuclease performances were detected on all repeat tracts. Wild-type Streptococcus pyogenes Cas9 (SpCas9) was more efficient than Staphylococcus aureus Cas9 on all repeats tested, except (CAG)33. Cas12a (Cpf1) was the most efficient on GAA trinucleotide repeats, whereas GC-rich repeats were more efficiently cut by SpCas9. The main genetic factor underlying Cas efficacy was the propensity of the recognition part of the sgRNA to form a stable secondary structure, independently of its structural part. This suggests that such structures form in vivo and interfere with sgRNA metabolism. The yeast genome contains 221 natural CAG/CTG and GAA/CTT trinucleotide repeats. Deep sequencing after nuclease induction identified three of them as carrying statistically significant low frequency mutations, corresponding to SpCas9 off-target double-strand breaks.

Role and Perspective of Molecular Simulation-Based Investigation of RNA-Ligand Interaction: From Small Molecules and Peptides to Photoswitchable RNA Binding

Aberrant RNA–protein complexes are formed in a variety of diseases. Identifying the ligands that interfere with their formation is a valuable therapeutic strategy. Molecular simulation, validated against experimental data, has recently emerged as a powerful tool to predict both the pose and energetics of such ligands. Thus, the use of molecular simulation may provide insight into aberrant molecular interactions in diseases and, from a drug design perspective, may allow for the employment of less wet lab resources than traditional in vitro compound screening approaches. With regard to basic research questions, molecular simulation can support the understanding of the exact molecular interaction and binding mode. Here, we focus on examples targeting RNA–protein complexes in neurodegenerative diseases and viral infections. These examples illustrate that the strategy is rather general and could be applied to different pharmacologically relevant approaches. We close this study by outlining one of these approaches, namely the light-controllable association of small molecules with RNA, as an emerging approach in RNA-targeting therapy.

Spontaneous and frequent conformational dynamics induced by A…A mismatch in d(CAA)·d(TAG) duplex

Base pair mismatches in DNA can erroneously be incorporated during replication, recombination, etc. Here, the influence of A…A mismatch in the context of 5′CAA·5′TAG sequence is explored using molecular dynamics (MD) simulation, umbrella sampling MD, circular dichroism (CD), microscale thermophoresis (MST) and NMR techniques. MD simulations reveal that the A…A mismatch experiences several transient events such as base flipping, base extrusion, etc. facilitating B–Z junction formation. A…A mismatch may assume such conformational transitions to circumvent the effect of nonisostericity with the flanking canonical base pairs so as to get accommodated in the DNA. CD and 1D proton NMR experiments further reveal that the extent of B–Z junction increases when the number of A…A mismatch in d(CAA)·d(T(A/T)G) increases (1–5). CD titration studies of d(CAA)·d(TAG)_n=5 with the hZα_ADAR1 show the passive binding between the two, wherein, the binding of protein commences with B–Z junction recognition. Umbrella sampling simulation indicates that the mismatch samples anti…+ syn/+ syn…anti, anti…anti & + syn…+ syn glycosyl conformations. The concomitant spontaneous transitions are: a variety of hydrogen bonding patterns, stacking and minor or major groove extrahelical movements (with and without the engagement of hydrogen bonds) involving the mismatch adenines. These transitions frequently happen in anti…anti conformational region compared with the other three regions as revealed from the lifetime of these states. Further, 2D-NOESY experiments indicate that the number of cross-peaks diminishes with the increasing number of A…A mismatches implicating its dynamic nature. The spontaneous extrahelical movement seen in A…A mismatch may be a key pre-trapping event in the mismatch repair due to the accessibility of the base(s) to the sophisticated mismatch repair machinery.

An Overview of Alternative Splicing Defects Implicated in Myotonic Dystrophy Type I

Myotonic dystrophy type I (DM1) is the most common form of adult muscular dystrophy, caused by expansion of a CTG triplet repeat in the 3′ untranslated region (3′UTR) of the myotonic dystrophy protein kinase (DMPK) gene. The pathological CTG repeats result in protein trapping by expanded transcripts, a decreased DMPK translation and the disruption of the chromatin structure, affecting neighboring genes expression. The muscleblind-like (MBNL) and CUG-BP and ETR-3-like factors (CELF) are two families of tissue-specific regulators of developmentally programmed alternative splicing that act as antagonist regulators of several pre-mRNA targets, including troponin 2 (TNNT2), insulin receptor (INSR), chloride channel 1 (CLCN1) and MBNL2. Sequestration of MBNL proteins and up-regulation of CELF1 are key to DM1 pathology, inducing a spliceopathy that leads to a developmental remodelling of the transcriptome due to an adult-to-foetal splicing switch, which results in the loss of cell function and viability. Moreover, recent studies indicate that additional pathogenic mechanisms may also contribute to disease pathology, including a misregulation of cellular mRNA translation, localization and stability. This review focuses on the cause and effects of MBNL and CELF1 deregulation in DM1, describing the molecular mechanisms underlying alternative splicing misregulation for a deeper understanding of DM1 complexity. To contribute to this analysis, we have prepared a comprehensive list of transcript alterations involved in DM1 pathogenesis, as well as other deregulated mRNA processing pathways implications.

Dynamics of strand slippage in DNA hairpins formed by CAG repeats: roles of sequence parity and trinucleotide interrupts

DNA trinucleotide repeats (TRs) can exhibit dynamic expansions by integer numbers of trinucleotides that lead to neurodegenerative disorders. Strand slipped hairpins during DNA replication, repair and/or recombination may contribute to TR expansion. Here, we combine single-molecule FRET experiments and molecular dynamics studies to elucidate slipping dynamics and conformations of (CAG)_n TR hairpins. We directly resolve slipping by predominantly two CAG units. The slipping kinetics depends on the even/odd repeat parity. The populated states suggest greater stability for 5′-AGCA-3′ tetraloops, compared with alternative 5′-CAG-3′ triloops. To accommodate the tetraloop, even(odd)-numbered repeats have an even(odd) number of hanging bases in the hairpin stem. In particular, a paired-end tetraloop (no hanging TR) is stable in (CAG)_{n = even}, but such situation cannot occur in (CAG)_{n = odd}, where the hairpin is “frustrated’’ and slips back and forth between states with one TR hanging at the 5′ or 3′ end. Trinucleotide interrupts in the repeating CAG pattern associated with altered disease phenotypes select for specific conformers with favorable loop sequences. Molecular dynamics provide atomic-level insight into the loop configurations. Reducing strand slipping in TR hairpins by sequence interruptions at the loop suggests disease-associated variations impact expansion mechanisms at the level of slipped hairpins.

Nuclear Magnetic Resonance of Single-Stranded RNAs and DNAs of CAAU and UCAAUC as Benchmarks for Molecular Dynamics Simulations

RNA and DNA are rapidly emerging as targets for therapeutics and as potential frameworks for nanotechnology. Accurate methods for predicting and designing structures and dynamics of nucleic acids would accelerate progress in these and other applications. Suitable approximations for modeling nucleic acids are being developed but require validation against disparate experimental observations. Here, nuclear magnetic resonance spectra for RNA and DNA single strands, CAAU and UCAAUC, are used as benchmarks to test molecular dynamics simulations with AMBER force fields OL3 and ROC-RNA for RNA and BSC1 for DNA. A detailed scheme for making comparisons is also presented. The results reflect recent progress in approximations and reveal remaining challenges.

Combination of coarse-grained molecular dynamics simulations and small-angle X-ray scattering experiments

The combination of molecular dynamics (MD) simulations and small-angle X-ray scattering (SAXS), called the MD-SAXS method, is efficient for investigating protein dynamics. To overcome the time-scale limitation of all-atom MD simulations, coarse-grained (CG) representations are often utilized for biomolecular simulations. In this study, we propose a method to combine CG MD simulations with SAXS, termed the CG-MD-SAXS method. In the CG-MD-SAXS method, the scattering factors of CG particles for proteins and nucleic acids are evaluated using high-resolution structural data in the Protein Data Bank, and the excluded volume and the hydration shell are modeled using two adjustable parameters to incorporate solvent effects. To avoid overfitting, only the two parameters are adjusted for an entire structure ensemble. To verify the developed method, theoretical SAXS profiles for various proteins, DNA/RNA, and a protein-RNA complex are compared with both experimental profiles and theoretical profiles obtained by the all-atom representation. In the present study, we applied the CG-MD-SAXS method to the Swi5-Sfr1 complex and three types of nucleosomes to obtain reliable ensemble models consistent with the experimental SAXS data.

Structural insights into synthetic ligands targeting A-A pairs in disease-related CAG RNA repeats

The trinucleotide repeat expansion disorders (TREDs) constitute of a group of >40 hereditary neurodegenerative human diseases associated with abnormal expansion of repeated sequences, such as CAG repeats. The pathogenic factor is a transcribed RNA or protein whose function in the cell is compromised. The disorders are progressive and incurable. Consequently, many ongoing studies are oriented at developing therapies. We have analyzed crystal structures of RNA containing CAG repeats in complex with synthetic cyclic mismatch-binding ligands (CMBLs). The models show well-defined interactions between the molecules in which the CMBLs mimic nucleobases as they form pseudo-canonical base pairs with adenosine residues and engage in extensive stacking interactions with neighboring nucleotides. The binding of ligands is associated with major structural changes of the CAG repeats, which is consistent with results of biochemical studies. The results constitute an early characterization of the first lead compounds in the search for therapy against TREDs. The crystallographic data indicate how the compounds could be further refined in future biomedical studies.

Reducing Mutant Huntingtin Protein Expression in Living Cells by a Newly Identified RNA CAG Binder

Expanded CAG trinucleotide repeats in Huntington's disease (HD) are causative for neurotoxicity. The mutant CAG repeat RNA encodes neurotoxic polyglutamine proteins and can lead to a toxic gain of function by aberrantly recruiting RNA-binding proteins. One of these is the MID1 protein, which induces aberrant Huntingtin (HTT) protein translation upon binding. Here we have identified a set of CAG repeat binder candidates by in silico methods. One of those, furamidine, reduces the level of binding of HTT mRNA to MID1 and other target proteins in vitro. Metadynamics calculations, fairly consistent with experimental data measured here, provide hints about the binding mode of the ligand. Importantly, furamidine also decreases the protein level of HTT in a HD cell line model. This shows that small molecules masking RNA-MID1 interactions may be active against mutant HTT protein in living cells.

Design of Bivalent Nucleic Acid Ligands for Recognition of RNA-Repeated Expansion Associated with Huntington's Disease

We report the development of a new class of nucleic acid ligands that is comprised of Janus bases and the MPγPNA backbone and is capable of binding rCAG repeats in a sequence-specific and selective manner via, inference, bivalent H-bonding interactions. Individually, the interactions between ligands and RNA are weak and transient. However, upon the installation of a C-terminal thioester and an N-terminal cystine and the reduction of disulfide bond, they undergo template-directed native chemical ligation to form concatenated oligomeric products that bind tightly to the RNA template. In the absence of an RNA target, they self-deactivate by undergoing an intramolecular reaction to form cyclic products, rendering them inactive for further binding. The work has implications for the design of ultrashort nucleic acid ligands for targeting rCAG-repeat expansion associated with Huntington's disease and a number of other related neuromuscular and neurodegenerative disorders.

A B-Z junction induced by an A … A mismatch in GAC repeats in the gene for cartilage oligomeric matrix protein promotes binding with the hZαADAR1 protein

GAC repeat expansion from five to seven in the exonic region of the gene for cartilage oligomeric matrix protein (COMP) leads to pseudoachondroplasia, a skeletal abnormality. However, the molecular mechanism by which GAC expansions in the COMP gene lead to skeletal dysplasias is poorly understood. Here we used molecular dynamics simulations, which indicate that an A … A mismatch in a d(GAC)6·d(GAC)6 duplex induces negative supercoiling, leading to a local B-to-Z DNA transition. This transition facilitates the binding of d(GAC)7·d(GAC)7 with the Zα-binding domain of human adenosine deaminase acting on RNA 1 (ADAR1, hZαADAR1), as confirmed by CD, NMR, and microscale thermophoresis studies. The CD results indicated that hZαADAR1 recognizes the zigzag backbone of d(GAC)7·d(GAC)7 at the B-Z junction and subsequently converts it into Z-DNA via the so-called passive mechanism. Molecular dynamics simulations carried out for the modeled hZαADAR1-d(GAC)6d(GAC)6 complex confirmed the retention of previously reported important interactions between the two molecules. These findings suggest that hZαADAR1 binding with the GAC hairpin stem in COMP can lead to a non-genetic, RNA editing-mediated substitution in COMP that may then play a crucial role in the development of pseudoachondroplasia.

Influence of mismatched and bulged nucleotides on SNP-preferential RNase H cleavage of RNA-antisense gapmer heteroduplexes

This study focused on determining design rules for gapmer-type antisense oligonucleotides (ASOs), that can differentiate cleavability of two SNP variants of RNA in the presence of ribonuclease H based on the mismatch type and position in the heteroduplex. We describe the influence of structural motifs formed by several arrangements of multiple mismatches (various types of mismatches and their position within the ASO/target RNA duplex) on RNase H cleavage selectivity of five different SNP types. The targets were mRNA fragments of APP, SCA3, SNCA and SOD1 genes, carrying C-to-G, G-to-C, G-to-A, A-to-G and C-to-U substitutions. The results show that certain arrangements of mismatches enhance discrimination between wild type and mutant SNP alleles of RNA in vitro as well as in HeLa cells. Among the over 120 gapmers tested, we found two gapmers that caused preferential degradation of the mutant allele APP 692 G and one that led to preferential cleavage of the mutant SNCA 53 A allele, both in vitro and in cells. However, several gapmers promoted selective cleavage of mRNA mutant alleles in in vitro experiments only.

Structure and Dynamics of DNA and RNA Double Helices of CAG and GAC Trinucleotide Repeats

CAG trinucleotide repeats are known to cause 10 late-onset progressive neurodegenerative disorders as the repeats expand beyond a threshold, whereas GAC repeats are associated with skeletal dysplasias and expand from the normal five to a maximum of seven repeats. The TR secondary structure is believed to play a role in CAG expansions. We have carried out free energy and molecular dynamics studies to determine the preferred conformations of the A-A noncanonical pairs in (CAG)_n and (GAC)_n trinucleotide repeats (n = 1, 4) and the consequent changes in the overall structure of the RNA and DNA duplexes. We find that the global free energy minimum corresponds to A-A pairs stacked inside the core of the helix with anti-anti conformations in RNA and (high-anti)-(high-anti) conformations in DNA. The next minimum corresponds to anti-syn conformations, whereas syn-syn conformations are higher in energy. Transition rates of the A-A conformations are higher for RNA than DNA. Mechanisms for these various transitions are identified. Additional structural and dynamical aspects of the helical conformations are explored, with a focus on contrasting CAG and GAC duplexes. The neutralizing ion distribution around the noncanonical pairs is described.

Structure and Dynamics of RNA Repeat Expansions That Cause Huntington's Disease and Myotonic Dystrophy Type 1

RNA repeat expansions cause a host of incurable, genetically defined diseases. The most common class of RNA repeats consists of trinucleotide repeats. These long, repeating transcripts fold into hairpins containing 1 × 1 internal loops that can mediate disease via a variety of mechanism(s) in which RNA is the central player. Two of these disorders are Huntington's disease and myotonic dystrophy type 1, which are caused by r(CAG) and r(CUG) repeats, respectively. We report the structures of two RNA constructs containing three copies of a r(CAG) [r(3×CAG)] or r(CUG) [r(3×CUG)] motif that were modeled with nuclear magnetic resonance spectroscopy and simulated annealing with restrained molecular dynamics. The 1 × 1 internal loops of r(3×CAG) are stabilized by one-hydrogen bond (cis Watson-Crick/Watson-Crick) AA pairs, while those of r(3×CUG) prefer one- or two-hydrogen bond (cis Watson-Crick/Watson-Crick) UU pairs. Assigned chemical shifts for the residues depended on the identity of neighbors or next nearest neighbors. Additional insights into the dynamics of these RNA constructs were gained by molecular dynamics simulations and a discrete path sampling method. Results indicate that the global structures of the RNA are A-form and that the loop regions are dynamic. The results will be useful for understanding the dynamic trajectory of these RNA repeats but also may aid in the development of therapeutics.

Conformational flexibility in the RNA stem-loop structures formed by CAG repeats

The expansion of CAG repeats has been found to be associated with at least nine human genetic disorders. In these disorders, the full-length expanded CAG RNA transcripts are cleaved into small CAG-repeated RNAs which are cytotoxic and known to be capable of forming hairpins. To better understand the RNA pathogenic mechanism, in this study we have performed high-resolution nuclear magnetic resonance structural investigations on the RNA hairpins formed by CAG repeats. Our results show the formation of a type III AGCA tetraloop and reveal the effect of stem rigidity on the loop conformational flexibility.

Structural Characteristics of Simple RNA Repeats Associated with Disease and their Deleterious Protein Interactions

Short Tandem Repeats (STRs) are frequent entities in many transcripts, however, in some cases, pathological events occur when a critical repeat length is reached. This phenomenon is observed in various neurological disorders, such as myotonic dystrophy type 1 (DM1), fragile X-associated tremor/ataxia syndrome, C9orf72-related amyotrophic lateral sclerosis and frontotemporal dementia (C9ALS/FTD), and polyglutamine diseases, such as Huntington's disease (HD) and spinocerebellar ataxias (SCA). The pathological effects of these repeats are triggered by mutant RNA transcripts and/or encoded mutant proteins, which depend on the localization of the expanded repeats in non-coding or coding regions. A growing body of recent evidence revealed that the RNA structures formed by these mutant RNA repeat tracts exhibit toxic effects on cells. Therefore, in this review article, we present existing knowledge on the structural aspects of different RNA repeat tracts as revealed mainly using well-established biochemical and biophysical methods. Furthermore, in several cases, it was shown that these expanded RNA structures are potent traps for a variety of RNA-binding proteins and that the sequestration of these proteins from their normal intracellular environment causes alternative splicing aberration, inhibition of nuclear transport and export, or alteration of a microRNA biogenesis pathway. Therefore, in this review article, we also present the most studied examples of abnormal interactions that occur between mutant RNAs and their associated proteins.

Structure and Dynamics of DNA and RNA Double Helices Obtained from the GGGGCC and CCCCGG Hexanucleotide Repeats That Are the Hallmark of C9FTD/ALS Diseases

A (GGGGCC) hexanucleotide repeat (HR) expansion in the C9ORF72 gene, and its associated antisense (CCCCGG) expansion, are considered the major cause behind frontotemporal dementia and amyotrophic lateral sclerosis. We have performed molecular dynamics simulations to characterize the conformation and dynamics of the 12 duplexes that result from the three different reading frames in sense and antisense HRs for both DNA and RNA. These duplexes display atypical structures relevant not only for a molecular level understanding of these diseases but also for enlarging the repertoire of nucleic-acid structural motifs. G-rich helices share common features. The inner G-G mismatches stay inside the helix in Gsyn-Ganti conformations and form two hydrogen bonds (HBs) between the Watson-Crick edge of Ganti and the Hoogsteen edge of Gsyn. In addition, Gsyn in RNA forms a base-phosphate HB. Inner G-G mismatches cause local unwinding of the helix. G-rich double helices are more stable than C-rich helices due to better stacking and HBs of G-G mismatches. C-rich helix conformations vary wildly. C mismatches flip out of the helix in DNA but not in RNA. Least (most) stable C-rich RNA and DNA helices have single (double) mismatches separated by two (four) Watson-Crick basepairs. The most stable DNA structure displays an "e-motif" where mismatched bases flip toward the minor groove and point in the 5' direction. There are two RNA conformations, where the orientation and HB pattern of the mismatches is coupled to bending of the helix.

Structures of RNA repeats associated with neurological diseases

All RNA molecules possess a 'propensity' to fold into complex secondary and tertiary structures. Although they are composed of only four types of nucleotides, they show an enormous structural richness which reflects their diverse functions in the cell. However, in some cases the folding of RNA can have deleterious consequences. Aberrantly expanded, repeated RNA sequences can exhibit gain-of-function abnormalities and become pathogenic, giving rise to many incurable neurological diseases. Most RNA repeats form long hairpin structures whose stem consists of noncanonical base pairs interspersed among Watson-Crick pairs. The expanded hairpins have an ability to sequester important proteins and form insoluble nuclear foci. The RNA pathology, common to many repeat disorders, has drawn attention to the structures of the RNA repeats. In this review, we summarize secondary structure probing and crystallographic studies of disease-related RNA repeat sequences. We discuss the unique structural features which can contribute to the pathogenic properties of the repeated runs. In addition, we present the newest reports concerning structural data linked to therapeutic approaches. WIREs RNA 2017, 8:e1412. doi: 10.1002/wrna.1412 For further resources related to this article, please visit the WIREs website.

Regulation of mRNA Translation by MID1: A Common Mechanism of Expanded CAG Repeat RNAs

Expansion of CAG repeats, which code for the disease-causing polyglutamine protein, is a common feature in polyglutamine diseases. RNA-mediated mechanisms that contribute to neuropathology in polyglutamine diseases are important. RNA-toxicity describes a phenomenon by which the mutant CAG repeat RNA recruits RNA-binding proteins, thereby leading to aberrant function. For example the MID1 protein binds to mutant huntingtin (HTT) RNA, which is linked to Huntington's disease (HD), at its CAG repeat region and induces protein synthesis of mutant protein. But is this mechanism specific to HD or is it a common mechanism in CAG repeat expansion disorders? To answer this question, we have analyzed the interaction between MID1 and three other CAG repeat mRNAs, Ataxin2 (ATXN2), Ataxin3 (ATXN3), and Ataxin7 (ATXN7), that all differ in the sequence flanking the CAG repeat. We show that ATXN2, ATXN3, and ATXN7 bind to MID1 in a CAG repeat length-dependent manner. Furthermore, we show that functionally, in line with what we have previously observed for HTT, the binding of MID1 to ATXN2, ATXN3, and ATXN7 mRNA induces protein synthesis in a repeat length-dependent manner. Our data suggest that regulation of protein translation by the MID1 complex is a common mechanism for CAG repeat containing mRNAs.

Watson-Crick-like pairs in CCUG repeats: evidence for tautomeric shifts or protonation

RNA transcripts that include expanded CCUG repeats are associated with myotonic dystrophy type 2. Crystal structures of two CCUG-containing oligomers show that the RNA strands associate into slipped duplexes that contain noncanonical C-U pairs that have apparently undergone tautomeric transition or protonation resulting in an unusual Watson-Crick-like pairing. The overhanging ends of the duplexes interact forming U-U pairs, which also show tautomerism. Duplexes consisting of CCUG repeats are thermodynamically less stable than the trinucleotide repeats involved in the TRED genetic disorders, but introducing LNA residues increases their stability and raises the melting temperature of the studied oligomers by ∼10°C, allowing detailed crystallographic studies. Quantum mechanical calculations were performed to test the possibility of the tautomeric transitions or protonation within the noncanonical pairs. The results indicate that tautomeric or ionic shifts of nucleobases can manifest themselves in biological systems, supplementing the canonical "rules of engagement."

The novel double-folded structure of d(GCATGCATGC): a possible model for triplet-repeat sequences

The structure of the decadeoxyribonucleotide d(GCATGCATGC) is presented at a resolution of 1.8 Å. The decamer adopts a novel double-folded structure in which the direction of progression of the backbone changes at the two thymine residues. Intra-strand stacking interactions (including an interaction between the endocylic O atom of a ribose moiety and the adjacent purine base), hydrogen bonds and cobalt-ion interactions stabilize the double-folded structure of the single strand. Two such double-folded strands come together in the crystal to form a dimer. Inter-strand Watson-Crick hydrogen bonds form four base pairs. This portion of the decamer structure is similar to that observed in other previously reported oligonucleotide structures and has been dubbed a `bi-loop'. Both the double-folded single-strand structure, as well as the dimeric bi-loop structure, serve as starting points to construct models for triplet-repeat DNA sequences, which have been implicated in many human diseases.

Molecular view of ligands specificity for CAG repeats in anti-Huntington therapy

Huntington's disease is a fatal and devastating neurodegenerative genetic disorder for which there is currently no cure. It is characterized by Huntingtin protein's mRNA transcripts with 36 or more CAG repeats. Inhibiting the formation of pathological complexes between these expanded transcripts and target proteins may be a valuable strategy against the disease. Yet, the rational design of molecules specifically targeting the expanded CAG repeats is limited by the lack of structural information. Here, we use well-tempered metadynamics-based free energy calculations to investigate pose and affinity of two ligands targeting CAG repeats for which affinities have been previously measured. The first consists of two 4-guanidinophenyl rings linked by an ester group. It is the most potent ligand identified so far, with Kd = 60(30) nM. The second consists of a 4-phenyl dihydroimidazole and 4-1H-indole dihydroimidazole connected by a C-C bond (Kd = 700(80) nM). Our calculations reproduce the experimental affinities and uncover the recognition pattern between ligands' and their RNA target. They also provide a molecular basis for the markedly different affinity of the two ligands for CAG repeats as observed experimentally. These findings may pave the way for a structure-based hit-to-lead optimization to further improve ligand selectivity toward CAG repeat-containing mRNAs.

Structural Insights Reveal the Dynamics of the Repeating r(CAG) Transcript Found in Huntington's Disease (HD) and Spinocerebellar Ataxias (SCAs)

In humans, neurodegenerative disorders such as Huntington's disease (HD) and many spinocerebellar ataxias (SCAs) have been found to be associated with CAG trinucleotide repeat expansion. An important RNA-mediated mechanism that causes these diseases involves the binding of the splicing regulator protein MBNL1 (Muscleblind-like 1 protein) to expanded r(CAG) repeats. Moreover, mutant huntingtin protein translated from expanded r(CAG) also yields toxic effects. To discern the role of mutant RNA in these diseases, it is essential to gather information about its structure. Detailed insight into the different structures and conformations adopted by these mutant transcripts is vital for developing therapeutics targeting them. Here, we report the crystal structure of an RNA model with a r(CAG) motif, which is complemented by an NMR-based solution structure obtained from restrained Molecular Dynamics (rMD) simulation studies. Crystal structure data of the RNA model resolved at 2.3 Å reveals non-canonical pairing of adenine in 5´-CAG/3´-GAC motif samples in different syn and anti conformations. The overall RNA structure has helical parameters intermediate to the A- and B-forms of nucleic acids due to the global widening of major grooves and base-pair preferences near internal AA loops. The comprehension of structural behaviour by studying the spectral features and the dynamics also supports the flexible nature of the r(CAG) motif.

Twisting right to left: A…A mismatch in a CAG trinucleotide repeat overexpansion provokes left-handed Z-DNA conformation

Conformational polymorphism of DNA is a major causative factor behind several incurable trinucleotide repeat expansion disorders that arise from overexpansion of trinucleotide repeats located in coding/non-coding regions of specific genes. Hairpin DNA structures that are formed due to overexpansion of CAG repeat lead to Huntington’s disorder and spinocerebellar ataxias. Nonetheless, DNA hairpin stem structure that generally embraces B-form with canonical base pairs is poorly understood in the context of periodic noncanonical A…A mismatch as found in CAG repeat overexpansion. Molecular dynamics simulations on DNA hairpin stems containing A…A mismatches in a CAG repeat overexpansion show that A…A dictates local Z-form irrespective of starting glycosyl conformation, in sharp contrast to canonical DNA duplex. Transition from B-to-Z is due to the mechanistic effect that originates from its pronounced nonisostericity with flanking canonical base pairs facilitated by base extrusion, backbone and/or base flipping. Based on these structural insights we envisage that such an unusual DNA structure of the CAG hairpin stem may have a role in disease pathogenesis. As this is the first study that delineates the influence of a single A…A mismatch in reversing DNA helicity, it would further have an impact on understanding DNA mismatch repair.

When a set of 3 nucleotides in a DNA sequence repeats beyond a certain number, it leads to incurable neurological or neuromuscular disorders. Such DNA sequences tend to form unusual DNA structures comprising of base pairing schemes different from the canonical A…T/G…C base pairs. Influence of such unusual base pairing on the overall 3-dimensional structure of DNA and its impact on the pathogenesis of disorder is not well understood. CAG repeat overexpansion that leads to Huntington’s disorder and several spinocerebellar ataxias forms noncanonical A…A base pair in between canonical C…G and G…C base pairs. However, no detailed structural information is available on the influence of an A…A mismatch on a DNA structure under any sequence context. Here, we have shown for the first time that A…A base pairing in a CAG repeat provokes the formation of left-handed Z-DNA due to the pronounced structural dissimilarity of A…A base pair with G…C base pair, leading to periodic B-Z junction. Thus, these results suggest that formation of periodic B-Z junction may be one of the molecular bases for CAG repeat instability.

Shortening trinucleotide repeats using highly specific endonucleases: a possible approach to gene therapy?

Trinucleotide repeat expansions are involved in more than two dozen neurological and developmental disorders. Conventional therapeutic approaches aimed at regulating the expression level of affected genes, which rely on drugs, oligonucleotides, and/or transgenes, have met with only limited success so far. An alternative approach is to shorten repeats to non-pathological lengths using highly specific nucleases. Here, I review early experiments using meganucleases, zinc-finger nucleases (ZFN), and transcription-activator like effector nucleases (TALENs) to contract trinucleotide repeats, and discuss the possibility of using CRISPR-Cas nucleases to the same end. Although this is a nascent field, I explore the possibility of designing nucleases and effectively delivering them in the context of gene therapy.

RNA-mediated pathogenic mechanisms in polyglutamine diseases and amyotrophic lateral sclerosis

Gene transcription produces a wide variety of ribonucleic acid (RNA) species in eukaryotes. Individual types of RNA, such as messenger, structural and regulatory RNA, are known to play distinct roles in the cell. Recently, researchers have identified a large number of RNA-mediated toxicity pathways that play significant pathogenic roles in numerous human disorders. In this article, we describe various common RNA toxicity pathways, namely epigenetic gene silencing, nucleolar stress, nucleocytoplasmic transport, bi-directional gene transcription, repeat-associated non-ATG translation, RNA foci formation and cellular protein sequestration. We emphasize RNA toxicity mechanisms that involve nucleotide repeat expansion, such as those related to polyglutamine (polyQ) disorders and frontotemporal lobar degeneration-amyotrophic lateral sclerosis.

RAN translation and frameshifting as translational challenges at simple repeats of human neurodegenerative disorders

Repeat-associated disorders caused by expansions of short sequences have been classified as coding and noncoding and are thought to be caused by protein gain-of-function and RNA gain-of-function mechanisms, respectively. The boundary between such classifications has recently been blurred by the discovery of repeat-associated non-AUG (RAN) translation reported in spinocerebellar ataxia type 8, myotonic dystrophy type 1, fragile X tremor/ataxia syndrome and C9ORF72 amyotrophic lateral sclerosis and frontotemporal dementia. This noncanonical translation requires no AUG start codon and can initiate in multiple frames of CAG, CGG and GGGGCC repeats of the sense and antisense strands of disease-relevant transcripts. RNA structures formed by the repeats have been suggested as possible triggers; however, the precise mechanism of the translation initiation remains elusive. Templates containing expansions of microsatellites have also been shown to challenge translation elongation, as frameshifting has been recognized across CAG repeats in spinocerebellar ataxia type 3 and Huntington's disease. Determining the critical requirements for RAN translation and frameshifting is essential to decipher the mechanisms that govern these processes. The contribution of unusual translation products to pathogenesis needs to be better understood. In this review, we present current knowledge regarding RAN translation and frameshifting and discuss the proposed mechanisms of translational challenges imposed by simple repeat expansions.

MBNL proteins and their target RNAs, interaction and splicing regulation

Muscleblind-like (MBNL) proteins are key regulators of precursor and mature mRNA metabolism in mammals. Based on published and novel data, we explore models of tissue-specific MBNL interaction with RNA. We portray MBNL domains critical for RNA binding and splicing regulation, and the structure of MBNL's normal and pathogenic RNA targets, particularly in the context of myotonic dystrophy (DM), in which expanded CUG or CCUG repeat transcripts sequester several nuclear proteins including MBNLs. We also review the properties of MBNL/RNA complex, including recent data obtained from UV cross-linking and immunoprecipitation (CLIP-Seq), and discuss how this interaction shapes normal MBNL-dependent alternative splicing regulation. Finally, we review how this acquired knowledge about the pathogenic RNA structure and nature of MBNL sequestration can be translated into the design of therapeutic strategies against DM.

Structural studies of CNG repeats

CNG repeats (where N denotes one of the four natural nucleotides) are abundant in the human genome. Their tendency to undergo expansion can lead to hereditary diseases known as TREDs (trinucleotide repeat expansion disorders). The toxic factor can be protein, if the abnormal gene is expressed, or the gene transcript, or both. The gene transcripts have attracted much attention in the biomedical community, but their molecular structures have only recently been investigated. Model RNA molecules comprising CNG repeats fold into long hairpins whose stems generally conform to an A-type helix, in which the non-canonical N-N pairs are flanked by C-G and G-C pairs. Each homobasic pair is accommodated in the helical context in a unique manner, with consequences for the local helical parameters, solvent structure, electrostatic potential and potential to interact with ligands. The detailed three-dimensional profiles of RNA CNG repeats can be used in screening of compound libraries for potential therapeutics and in structure-based drug design. Here is a brief survey of the CNG structures published to date.

Oligonucleotide-based strategies to combat polyglutamine diseases

Considerable advances have been recently made in understanding the molecular aspects of pathogenesis and in developing therapeutic approaches for polyglutamine (polyQ) diseases. Studies on pathogenic mechanisms have extended our knowledge of mutant protein toxicity, confirmed the toxicity of mutant transcript and identified other toxic RNA and protein entities. One very promising therapeutic strategy is targeting the causative gene expression with oligonucleotide (ON) based tools. This straightforward approach aimed at halting the early steps in the cascade of pathogenic events has been widely tested for Huntington's disease and spinocerebellar ataxia type 3. In this review, we gather information on the use of antisense oligonucleotides and RNA interference triggers for the experimental treatment of polyQ diseases in cellular and animal models. We present studies testing non-allele-selective and allele-selective gene silencing strategies. The latter include targeting SNP variants associated with mutations or targeting the pathologically expanded CAG repeat directly. We compare gene silencing effectors of various types in a number of aspects, including their design, efficiency in cell culture experiments and pre-clinical testing. We discuss advantages, current limitations and perspectives of various ON-based strategies used to treat polyQ diseases.

RNA-protein interactions in unstable microsatellite diseases

A novel RNA-mediated disease mechanism has emerged from studies on dominantly inherited neurological disorders caused by unstable microsatellite expansions in non-coding regions of the genome. These non-coding tandem repeat expansions trigger the production of unusual RNAs that gain a toxic function, which involves the formation of RNA repeat structures that interact with, and alter the activities of, various factors required for normal RNA processing as well as additional cellular functions. In this review, we explore the deleterious effects of toxic RNA expression and discuss the various model systems currently available for studying RNA gain-of-function in neurologic diseases. Common themes, including bidirectional transcription and repeat-associated non-ATG (RAN) translation, have recently emerged from expansion disease studies. These and other discoveries have highlighted the need for further investigations designed to provide the additional mechanistic insights essential for future therapeutic development.

Endonuclease V cleaves at inosines in RNA

Endonuclease V orthologues are highly conserved proteins found in all kingdoms of life. While the prokaryotic enzymes are DNA repair proteins for removal of deaminated adenosine (inosine) from the genome, no clear role for the eukaryotic counterparts has hitherto been described. Here we report that human endonuclease V (ENDOV) and also Escherichia coli endonuclease V are highly active ribonucleases specific for inosine in RNA. Inosines are normal residues in certain RNAs introduced by specific deaminases. Adenosine-to-inosine editing is essential for proper function of these transcripts and defects are linked to various human disease. Here we show that human ENDOV cleaves an RNA substrate containing inosine in a position corresponding to a biologically important site for deamination in the Gabra-3 transcript of the GABA(A) neurotransmitter. Further, human ENDOV specifically incises transfer RNAs with inosine in the wobble position. This previously unknown RNA incision activity may suggest a role for endonuclease V in normal RNA metabolism.

Methods to enable the design of bioactive small molecules targeting RNA

RNA is an immensely important target for small molecule therapeutics or chemical probes of function. However, methods that identify, annotate, and optimize RNA-small molecule interactions that could enable the design of compounds that modulate RNA function are in their infancies. This review describes recent approaches that have been developed to understand and optimize RNA motif-small molecule interactions, including structure-activity relationships through sequencing (StARTS), quantitative structure-activity relationships (QSAR), chemical similarity searching, structure-based design and docking, and molecular dynamics (MD) simulations. Case studies described include the design of small molecules targeting RNA expansions, the bacterial A-site, viral RNAs, and telomerase RNA. These approaches can be combined to afford a synergistic method to exploit the myriad of RNA targets in the transcriptome.

Optimization of an AMBER force field for the artificial nucleic acid, LNA, and benchmarking with NMR of L(CAAU)

Locked Nucleic Acids (LNAs) are RNA analogues with an O2'-C4' methylene bridge which locks the sugar into a C3'-endo conformation. This enhances hybridization to DNA and RNA, making LNAs useful in microarrays and potential therapeutics. Here, the LNA, L(CAAU), provides a simplified benchmark for testing the ability of molecular dynamics (MD) to approximate nucleic acid properties. LNA χ torsions and partial charges were parametrized to create AMBER parm99_LNA. The revisions were tested by comparing MD predictions with AMBER parm99 and parm99_LNA against a 200 ms NOESY NMR spectrum of L(CAAU). NMR indicates an A-Form equilibrium ensemble. In 3000 ns simulations starting with an A-form structure, parm99_LNA and parm99 provide 66% and 35% agreement, respectively, with NMR NOE volumes and (3)J-couplings. In simulations of L(CAAU) starting with all χ torsions in a syn conformation, only parm99_LNA is able to repair the structure. This implies methods for parametrizing force fields for nucleic acid mimics can reasonably approximate key interactions and that parm99_LNA will improve reliability of MD studies for systems with LNA. A method for approximating χ population distribution on the basis of base to sugar NOEs is also introduced.

Is the DPT tautomerization of the long A·G Watson-Crick DNA base mispair a source of the adenine and guanine mutagenic tautomers? A QM and QTAIM response to the biologically important question

Herein, we first address the question posed in the title by establishing the tautomerization trajectory via the double proton transfer of the adenine·guanine (A·G) DNA base mispair formed by the canonical tautomers of the A and G bases into the A*·G* DNA base mispair, involving mutagenic tautomers, with the use of the quantum-mechanical calculations and quantum theory of atoms in molecules (QTAIM). It was detected that the A·G ↔ A*·G* tautomerization proceeds through the asynchronous concerted mechanism. It was revealed that the A·G base mispair is stabilized by the N6H···O6 (5.68) and N1H···N1 (6.51) hydrogen bonds (H-bonds) and the N2H···HC2 dihydrogen bond (DH-bond) (0.68 kcal·mol(-1) ), whereas the A*·G* base mispair-by the O6H···N6 (10.88), N1H···N1 (7.01) and C2H···N2 H-bonds (0.42 kcal·mol(-1) ). The N2H···HC2 DH-bond smoothly and without bifurcation transforms into the C2H···N2 H-bond at the IRC = -10.07 Bohr in the course of the A·G ↔ A*·G* tautomerization. Using the sweeps of the energies of the intermolecular H-bonds, it was observed that the N6H···O6 H-bond is anticooperative to the two others-N1H···N1 and N2H···HC2 in the A·G base mispair, while the latters are significantly cooperative, mutually strengthening each other. In opposite, all three O6H···N6, N1H···N1, and C2H···N2 H-bonds are cooperative in the A*·G* base mispair. All in all, we established the dynamical instability of the А*·G* base mispair with a short lifetime (4.83·10(-14) s), enabling it not to be deemed feasible source of the A* and G* mutagenic tautomers of the DNA bases. The small lifetime of the А*·G* base mispair is predetermined by the negative value of the Gibbs free energy for the A*·G* → A·G transition. Moreover, all of the six low-frequency intermolecular vibrations cannot develop during this lifetime that additionally confirms the aforementioned results. Thus, the A*·G* base mispair cannot be considered as a source of the mutagenic tautomers of the DNA bases, as the A·G base mispair dissociates during DNA replication exceptionally into the A and G monomers in the canonical tautomeric form.

Role of context in RNA structure: flanking sequences reconfigure CAG motif folding in huntingtin exon 1 transcripts

The length of the CAG-repeat region in the huntingtin mRNA is predictive of Huntington's disease. Structural studies of CAG-repeat-containing RNAs suggest that these sequences form simple hairpin structures; however, in the context of the full-length huntingtin mRNA, CAG repeats may form complex structures that could be targeted for therapeutic intervention. We examined the structures of transcripts spanning the first exon of the huntingtin mRNA with both healthy and disease-prone repeat lengths. In transcripts with 17-70 repeats, the CAG sequences base paired extensively with nucleotides in the 5' UTR and with conserved downstream sequences including a CCG-repeat region. In huntingtin transcripts with healthy numbers of repeats, the previously observed CAG hairpin was either absent or short. In contrast, in transcripts with disease-associated numbers of repeats, a CAG hairpin was present and extended from a three-helix junction. Our findings demonstrate the profound importance of sequence context in RNA folding and identify specific structural differences between healthy and disease-inducing huntingtin alleles that may be targets for therapeutic intervention.

Infrared multiple photon dissociation action spectroscopy of proton-bound dimers of cytosine and modified cytosines: effects of modifications on gas-phase conformations

The gas-phase structures of proton-bound dimers of cytosine and modified cytosines and their d6-analogues generated by electrospray ionization are probed via infrared multiple photon dissociation (IRMPD) action spectroscopy and theoretical electronic structure calculations. The modified cytosines examined include the 5-methyl-, 5-fluoro-, and 5-bromo-substituted species. IRMPD action spectra of seven proton-bound dimers exhibit both similar and distinctive spectral features over the range of ∼2600-3700 cm(-1). The IRMPD spectra of all of these proton-bound dimers are relatively simple, but exhibit obvious shifts in the positions of several bands that correlate with the properties of the substituent. The measured IRMPD spectra are compared to linear IR spectra calculated for the stable low-energy tautomeric conformations, determined at the B3LYP/6-31G* level of theory, to identify the conformations accessed in the experiments. Comparison of the measured IRMPD and calculated IR spectra indicates that only a single conformation, the ground-state structure, is accessed for all proton-bound homodimers, whereas the ground-state and a small population of the first-excited tautomeric conformations are accessed for all proton-bound heterodimers. In all cases, three hydrogen-bonding interactions in which the nucleobases are aligned in an antiparallel fashion analogous to that of the DNA i-motif are responsible for stabilizing the base pairing. Thus, base modifications such as 5-methyl- and 5-halo-substitution of cytosine should not alter the structure of the DNA i-motif.

An innate twist between Crick's wobble and Watson-Crick base pairs

Non-Watson-Crick pairs like the G·U wobble are frequent in RNA duplexes. Their geometric dissimilarity (nonisostericity) with the Watson-Crick base pairs and among themselves imparts structural variations decisive for biological functions. Through a novel circular representation of base pairs, a simple and general metric scheme for quantification of base-pair nonisostericity, in terms of residual twist and radial difference that can also envisage its mechanistic effect, is proposed. The scheme is exemplified by G·U and U·G wobble pairs, and their predicable local effects on helical twist angle are validated by MD simulations. New insights into a possible rationale for contextual occurrence of G·U and other non-WC pairs, as well as the influence of a G·U pair on other non-Watson-Crick pair neighborhood and RNA-protein interactions are obtained from analysis of crystal structure data. A few instances of RNA-protein interactions along the major groove are documented in addition to the well-recognized interaction of the G·U pair along the minor groove. The nonisostericity-mediated influence of wobble pairs for facilitating helical packing through long-range interactions in ribosomal RNAs is also reviewed.