Utilizing the GAAA tetraloop/receptor to facilitate crystal packing and determination of the structure of a CUG RNA helix

Structures of small molecules bound to RNA repeat expansions that cause Huntington's disease-like 2 and myotonic dystrophy type 1

Trinucleotide repeat expansions fold into long, stable hairpins and cause a variety of incurable RNA gain-of-function diseases such as Huntington's disease, the myotonic dystrophies, and spinocerebellar ataxias. One approach for treating these diseases is to bind small molecules to the structured RNAs. Both Huntington's disease-like 2 (HDL2) and myotonic dystrophy type 1 (DM1) are caused by a r(CUG) repeat expansion, or r(CUG)exp. The RNA folds into a hairpin structure with a periodic array of 1 × 1 nucleotide UU loops (5'CUG/3'GUC; where the underlined nucleotides indicate the Us in the internal loop) that sequester various RNA-binding proteins (RBP) and hence the source of its gain-of-function. Here, we report NMR-refined structures of single 5'CUG/3'GUC motifs in complex with three different small molecules, a di-guandinobenzoate (1), a derivative of 1 where the guanidino groups have been exchanged for imidazole (2), and a quinoline with improved drug-like properties (3). These structures were determined using nuclear magnetic resonance (NMR) spectroscopy and simulated annealing with restrained molecular dynamics (MD). Compounds 1, 2, and 3 formed stacking and hydrogen bonding interactions with the 5'CUG/3'GUC motif. Compound 3 also formed van der Waals interactions with the internal loop. The global structure of each RNA-small molecule complexes retains an A-form conformation, while the internal loops are still dynamic but to a lesser extent compared to the unbound form. These results aid our understanding of ligand-RNA interactions and enable structure-based design of small molecules with improved binding affinity for and biological activity against r(CUG)exp. As the first ever reported structures of RNA r(CUG) repeats bound to ligands, these structures can enable virtual screening campaigns combined with machine learning assisted de novo design.

NMR structures of small molecules bound to a model of an RNA CUG repeat expansion

Trinucleotide repeat expansions fold into long, stable hairpins and cause a variety of incurable RNA gain-of-function diseases such as Huntington's disease, the myotonic dystrophies, and spinocerebellar ataxias. One approach for treating these diseases is to bind small molecules to the structured RNAs. Both Huntington's disease-like 2 (HDL2) and myotonic dystrophy type 1 (DM1) are caused by a r(CUG) repeat expansion, or r(CUG)exp. The RNA folds into a hairpin structure with a periodic array of 1×1 nucleotide UU loops (5'CUG/3'GUC; where the underlined nucleotides indicate the Us in the internal loop) that sequester various RNA-binding proteins (RBP) and hence the source of its gain-of-function. Here, we report NMR-refined structures of single 5'CUG/3'GUC motifs in complex with three different small molecules, a di-guandinobenzoate (1), a derivative of 1 where the guanidino groups have been exchanged for imidazole (2), and a quinoline with improved drug-like properties (3). These structures were determined using nuclear magnetic resonance (NMR) spectroscopy and simulated annealing with restrained molecular dynamics (MD). Compounds 1, 2, and 3 formed stacking and hydrogen bonding interactions with the 5'CUG/3'GUC motif. Compound 3 also formed van der Waals interactions with the internal loop. The global structure of each RNA-small molecule complexes retains an A-form conformation, while the internal loops are still dynamic but to a lesser extent compared to the unbound form. These results aid our understanding of ligand-RNA interactions and enable structure-based design of small molecules with improved binding affinity for and biological activity against r(CUG)exp. As the first ever reported structures of RNA r(CUG) repeats bound to ligands, these structures can enable virtual screening campaigns combined with machine learning assisted de novo design.

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.

Structural basis for water modulating RNA duplex formation in the CUG repeats of myotonic dystrophy type 1

Secondary structures formed by expanded CUG RNA are involved in the pathobiology of myotonic dystrophy type 1. Understanding the molecular basis of toxic RNA structures can provide insights into the mechanism of disease pathogenesis and accelerate the drug discovery process. Here, we report the crystal structure of CUG repeat RNA containing three U-U mismatches between C-G and G-C base pairs. The CUG RNA crystallizes as an A-form duplex, with the first and third U-U mismatches adopting a water-mediated asymmetric mirror isoform geometry. We found for the first time that a symmetric, water-bridged U-H₂O-U mismatch is well tolerated within the CUG RNA duplex, which was previously suspected but not observed. The new water-bridged U-U mismatch resulted in high base-pair opening and single-sided cross-strand stacking interactions, which in turn dominate the CUG RNA structure. Furthermore, we performed molecular dynamics (MD) simulations that complemented the structural findings and proposed that the first and third U-U mismatches are interchangeable conformations, while the central water-bridged U-U mismatch represents an intermediate state that modulates the RNA duplex conformation. Collectively, the new structural features provided in this work are important for understanding the recognition of U-U mismatches in CUG repeats by external ligands such as proteins or small molecules.

General Strategies for RNA X-ray Crystallography

An extremely small proportion of the X-ray crystal structures deposited in the Protein Data Bank are of RNA or RNA-protein complexes. This is due to three main obstacles to the successful determination of RNA structure: (1) low yields of pure, properly folded RNA; (2) difficulty creating crystal contacts due to low sequence diversity; and (3) limited methods for phasing. Various approaches have been developed to address these obstacles, such as native RNA purification, engineered crystallization modules, and incorporation of proteins to assist in phasing. In this review, we will discuss these strategies and provide examples of how they are used in practice.

Engineering Crystal Packing in RNA Structures I: Past and Future Strategies for Engineering RNA Packing in Crystals

X-ray crystallography remains a powerful method to gain atomistic insights into the catalytic and regulatory functions of RNA molecules. However, the technique requires the preparation of diffraction-quality crystals. This is often a resource- and time-consuming venture because RNA crystallization is hindered by the conformational heterogeneity of RNA, as well as the limited opportunities for stereospecific intermolecular interactions between RNA molecules. The limited success at crystallization explains in part the smaller number of RNA-only structures in the Protein Data Bank. Several approaches have been developed to aid the formation of well-ordered RNA crystals. The majority of these are construct-engineering techniques that aim to introduce crystal contacts to favor the formation of well-diffracting crystals. A typical example is the insertion of tetraloop-tetraloop receptor pairs into non-essential RNA segments to promote intermolecular association. Other methods of promoting crystallization involve chaperones and crystallization-friendly molecules that increase RNA stability and improve crystal packing. In this review, we discuss the various techniques that have been successfully used to facilitate crystal packing of RNA molecules, recent advances in construct engineering, and directions for future research in this vital aspect of RNA crystallography.

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.

Nop9 recognizes structured and single-stranded RNA elements of preribosomal RNA

Nop9 is an essential factor in the processing of preribosomal RNA. Its absence in yeast is lethal, and defects in the human ortholog are associated with breast cancer, autoimmunity, and learning/language impairment. PUF family RNA-binding proteins are best known for sequence-specific RNA recognition, and most contain eight α-helical repeats that bind to the RNA bases of single-stranded RNA. Nop9 is an unusual member of this family in that it contains eleven repeats and recognizes both RNA structure and sequence. Here we report a crystal structure of Saccharomyces cerevisiae Nop9 in complex with its target RNA within the 20S preribosomal RNA. This structure reveals that Nop9 brings together a carboxy-terminal module recognizing the 5' single-stranded region of the RNA and a bifunctional amino-terminal module recognizing the central double-stranded stem region. We further show that the 3' single-stranded region of the 20S target RNA adds sequence-independent binding energy to the RNA-Nop9 interaction. Both the amino- and carboxy-terminal modules retain the characteristic sequence-specific recognition of PUF proteins, but the amino-terminal module has also evolved a distinct interface, which allows Nop9 to recognize either single-stranded RNA sequences or RNAs with a combination of single-stranded and structured elements.

Structure of an RNA helix with pyrimidine mismatches and cross-strand stacking

The structure of a 22-base-pair RNA helix with mismatched pyrimidine base pairs is reported. The helix contains two symmetry-related CUG sequences: a triplet-repeat motif implicated in myotonic dystrophy type 1. The CUG repeat contains a U-U mismatch sandwiched between Watson-Crick pairs. Additionally, the center of the helix contains a dimerized UUCG motif with tandem pyrimidine (U-C/C-U) mismatches flanked by U-G wobble pairs. This region of the structure is significantly different from previously observed structures that share the same sequence and neighboring base pairs. The tandem pyrimidine mismatches are unusual and display sheared, cross-strand stacking geometries that locally constrict the helical width, a type of stacking previously associated with purines in internal loops. Thus, pyrimidine-rich regions of RNA have a high degree of structural diversity.

Repeat-associated RNA structure and aberrant splicing

Over 30 hereditary disorders attributed to the expansion of microsatellite repeats have been identified. Despite variant nucleotide content, number of consecutive repeats, and different locations in the genome, many of these diseases have pathogenic RNA gain-of-function mechanisms. The repeat-containing RNAs can form structures in vitro predicted to contribute to the disease through assembly of intracellular RNA aggregates termed foci. The expanded repeat RNAs within these foci sequester RNA binding proteins (RBPs) with important roles in the regulation of RNA metabolism, most notably alternative splicing (AS). These deleterious interactions lead to downstream alterations in transcriptome-wide AS directly linked with disease symptoms. This review summarizes existing knowledge about the association between the repeat RNA structures and RBPs as well as the resulting aberrant AS patterns, specifically in the context of myotonic dystrophy. The connection between toxic, structured RNAs and dysregulation of AS in other repeat expansion diseases is also discussed. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Identification of receptors for UNCG and GNRA Z-turns and their occurrence in rRNA

In contrast to GNRA tetraloop receptors that are common in RNA, receptors for the more thermostable UNCG loops have remained elusive for almost three decades. An analysis of all RNA structures with resolution ≤3.0 Å from the PDB allowed us to identify three previously unnoticed receptors for UNCG and GNRA tetraloops that adopt a common UNCG fold, named ‘Z-turn’ in agreement with our previously published nomenclature. These receptors recognize the solvent accessible second Z-turn nucleotide in different but specific ways. Two receptors participating in a complex network of tertiary interactions are associated with the rRNA UUCG and GAAA Z-turns capping helices H62 and H35a in rRNA large subunits. Structural comparison of fully assembled ribosomes and comparative sequence analysis of >6500 rRNA sequences helped us recognize that these motifs are almost universally conserved in rRNA, where they may contribute to organize the large subunit around the subdomain-IV four-way junction. The third UCCG receptor was identified in a rRNA/protein construct crystallized at acidic pH. These three non-redundant Z-turn receptors are relevant for our understanding of the assembly of rRNA and other long-non-coding RNAs, as well as for the design of novel folding motifs for synthetic biology.

Structural Rationale for the Enhanced Catalysis of Nonenzymatic RNA Primer Extension by a Downstream Oligonucleotide

Nonenzymatic RNA primer extension by activated mononucleotides has long served as a model for the study of prebiotic RNA copying. We have recently shown that the rate of primer extension is greatly enhanced by the formation of an imidazolium-bridged dinucleotide between the incoming monomer and a second, downstream activated monomer. However, the rate of primer extension is further enhanced if the downstream monomer is replaced by an activated oligonucleotide. Even an unactivated downstream oligonucleotide provides a modest enhancement in the rate of reaction of a primer with a single activated monomer. Here we study the mechanism of these effects through crystallographic studies of RNA complexes with the recently synthesized nonhydrolyzable substrate analog, guanosine 5'-(4-methylimidazolyl)-phosphonate (ICG). ICG mimics 2-methylimidazole activated guanosine-5'-phosphate (2-MeImpG), a commonly used substrate in nonenzymatic primer extension experiments. We present crystal structures of primer-template complexes with either one or two ICG residues bound downstream of a primer. In both cases, the aryl-phosphonate moiety of the ICG adjacent to the primer is disordered. To investigate the effect of a downstream oligonucleotide, we transcribed a short RNA oligonucleotide with either a 5'-ICG residue, a 5'-phosphate or a 5'-hydroxyl. We then determined crystal structures of primer-template complexes with a bound ICG monomer sandwiched between the primer and each of the three downstream oligonucleotides. Surprisingly, all three oligonucleotides rigidify the ICG monomer conformation and position it for attack by the primer 3'-hydroxyl. Furthermore, when GpppG, an analog of the imidazolium-bridged intermediate, is sandwiched between an upstream primer and a downstream helper oligonucleotide, or covalently linked to the 5'-end of the downstream oligonucleotide, the complex is better preorganized for primer extension than in the absence of a downstream oligonucleotide. Our results suggest that a downstream helper oligonucleotide contributes to the catalysis of primer extension by favoring a reactive conformation of the primer-template-intermediate complex.

Effects of osmolytes and macromolecular crowders on stable GAAA tetraloops and their preference for a CG closing base pair

Osmolytes and macromolecular crowders have the potential to influence the stability of secondary structure motifs and alter preferences for conserved nucleic acid sequences in vivo. To further understand the cellular function of RNA we observed the effects of a model osmolyte, polyethylene glycol (PEG) 200, and a model macromolecular crowding agent, PEG 8000, on the GAAA tetraloop motif. GAAA tetraloops are conserved, stable tetraloops, and are critical participants in RNA tertiary structure. They also have a thermodynamic preference for a CG closing base pair. The thermal denaturation of model hairpins containing GAAA loops was monitored using UV-Vis spectroscopy in the presence and absence of PEG 200 or PEG 8000. Both of the cosolutes tested influenced the thermodynamic preference for a CG base pair by destabilizing the loop with a CG closing base pair relative to the loop with a GC closing base pair. This result also extended to a related DNA triloop, which provides further evidence that the interactions between the loop and closing base pair are identical for the d(GCA) triloop and the GAAA tetraloop. Our results suggest that in the presence of model PEG molecules, loops with a GC closing base pair may retain some preferential interactions with the cosolutes that are lost in the presence of the CG closing base pair. These results reveal that relatively small structural changes could influence how neutral cosolutes tune the stability and function of secondary structure motifs in vivo.

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.

Martini Coarse-Grained Force Field: Extension to RNA

RNA has an important role not only as the messenger of genetic information but also as a regulator of gene expression. Given its central role in cell biology, there is significant interest in studying the structural and dynamic behavior of RNA in relation to other biomolecules. Coarse-grain molecular dynamics simulations are a key tool to that end. Here, we have extended the coarse-grain Martini force field to include RNA after our recent extension to DNA. In the same way DNA was modeled, the tertiary structure of RNA is constrained using an elastic network. This model, therefore, is not designed for applications involving RNA folding but rather offers a stable RNA structure for studying RNA interactions with other (bio)molecules. The RNA model is compatible with all other Martini models and opens the way to large-scale explicit-solvent molecular dynamics simulations of complex systems involving RNA.

In silico discovery of substituted pyrido[2,3-d]pyrimidines and pentamidine-like compounds with biological activity in myotonic dystrophy models

Myotonic dystrophy type 1 (DM1) is a rare multisystemic disorder associated with an expansion of CUG repeats in mutant DMPK (dystrophia myotonica protein kinase) transcripts; the main effect of these expansions is the induction of pre-mRNA splicing defects by sequestering muscleblind-like family proteins (e.g. MBNL1). Disruption of the CUG repeats and the MBNL1 protein complex has been established as the best therapeutic approach for DM1, hence two main strategies have been proposed: targeted degradation of mutant DMPK transcripts and the development of CUG-binding molecules that prevent MBNL1 sequestration. Herein, suitable CUG-binding small molecules were selected using in silico approaches such as scaffold analysis, similarity searching, and druggability analysis. We used polarization assays to confirm the CUG repeat binding in vitro for a number of candidate compounds, and went on to evaluate the biological activity of the two with the strongest affinity for CUG repeats (which we refer to as compounds 1–2 and 2–5) in DM1 mutant cells and Drosophila DM1 models with an impaired locomotion phenotype. In particular, 1–2 and 2–5 enhanced the levels of free MBNL1 in patient-derived myoblasts in vitro and greatly improved DM1 fly locomotion in climbing assays. This work provides new computational approaches for rational large-scale virtual screens of molecules that selectively recognize CUG structures. Moreover, it contributes valuable knowledge regarding two compounds with desirable biological activity in DM1 models.

Structural basis for high-affinity fluorophore binding and activation by RNA Mango

Genetically encoded fluorescent protein tags revolutionized proteome studies, while the lack of intrinsically fluorescent RNAs has hindered transcriptome exploration. Among several RNA-fluorophore complexes that potentially address this problem, RNA Mango has an exceptionally high affinity for its thiazole orange (TO)-derived fluorophore, TO1-Biotin (K_d ~3 nM), and in complex with related ligands, is one of the most red-shifted fluorescent macromolecular tags known. To elucidate how this small aptamer exhibits such properties, which make it well suited for studying low-copy cellular RNAs, we determined its 1.7 Å resolution co-crystal structure. Unexpectedly, the entire ligand, including TO, biotin, and the linker connecting them, abuts one of the near-planar faces of the three-tiered G-quadruplex. The two heterocycles of TO are held in place by two loop adenines and make a 45° angle with respect to each other. Minimizing this angle would increase quantum yield and further improve this tool for in vivo RNA visualization.

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.

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.

Controlled dehydration improves the diffraction quality of two RNA crystals

BACKGROUND

Post-crystallization dehydration methods, applying either vapor diffusion or humidity control devices, have been widely used to improve the diffraction quality of protein crystals. Despite the fact that RNA crystals tend to diffract poorly, there is a dearth of reports on the application of dehydration methods to improve the diffraction quality of RNA crystals.

RESULTS

We use dehydration techniques with a Free Mounting System (FMS, a humidity control device) to recover the poor diffraction quality of RNA crystals. These approaches were applied to RNA constructs that model various RNA-mediated repeat expansion disorders.

CONCLUSION

The method we describe herein could serve as a general tool to improve diffraction quality of RNA crystals to facilitate structure determinations.

Dose-Dependent Regulation of Alternative Splicing by MBNL Proteins Reveals Biomarkers for Myotonic Dystrophy

Alternative splicing is a regulated process that results in expression of specific mRNA and protein isoforms. Alternative splicing factors determine the relative abundance of each isoform. Here we focus on MBNL1, a splicing factor misregulated in the disease myotonic dystrophy. By altering the concentration of MBNL1 in cells across a broad dynamic range, we show that different splicing events require different amounts of MBNL1 for half-maximal response, and respond more or less steeply to MBNL1. Motifs around MBNL1 exon 5 were studied to assess how cis-elements mediate the MBNL1 dose-dependent splicing response. A framework was developed to estimate MBNL concentration using splicing responses alone, validated in the cell-based model, and applied to myotonic dystrophy patient muscle. Using this framework, we evaluated the ability of individual and combinations of splicing events to predict functional MBNL concentration in human biopsies, as well as their performance as biomarkers to assay mild, moderate, and severe cases of DM.

Our studies provide insight into the mechanisms of myotonic dystrophy, the most common adult form of muscular dystrophy. In this disease, a family of RNA binding proteins is sequestered by toxic RNA, which leads to mis-regulation and disease symptoms. We have created a cellular model with one of these family members to study how these RNA binding proteins function in the absence of the toxic RNA. In parallel, we analyzed transcriptomic data from over 50 individuals (44 affected by myotonic dystrophy) with a range of disease severity. The results from the transcriptomic data provide a rational approach to select biomarkers for clinical research and therapeutic trials.

On the Applicability of Elastic Network Models for the Study of RNA CUG Trinucleotide Repeat Overexpansion

Non-coding RNAs play a pivotal role in a number of diseases promoting an aberrant sequestration of nuclear RNA-binding proteins. In the particular case of myotonic dystrophy type 1 (DM1), a multisystemic autosomal dominant disease, the formation of large non-coding CUG repeats set up long-tract hairpins able to bind muscleblind-like proteins (MBNL), which trigger the deregulation of several splicing events such as cardiac troponin T (cTNT) and insulin receptor’s, among others. Evidence suggests that conformational changes in RNA are determinant for the recognition and binding of splicing proteins, molecular modeling simulations can attempt to shed light on the structural diversity of CUG repeats and to understand their pathogenic mechanisms. Molecular dynamics (MD) are widely used to obtain accurate results at atomistic level, despite being very time consuming, and they contrast with fast but simplified coarse-grained methods such as Elastic Network Model (ENM). In this paper, we assess the application of ENM (traditionally applied on proteins) for studying the conformational space of CUG repeats and compare it to conventional and accelerated MD conformational sampling. Overall, the results provided here reveal that ANM can provide useful insights into dynamic rCUG structures at a global level, and that their dynamics depend on both backbone and nucleobase fluctuations. On the other hand, ANM fail to describe local U-U dynamics of the rCUG system, which require more computationally expensive methods such as MD. Given that several limitations are inherent to both methods, we discuss here the usefulness of the current theoretical approaches for studying highly dynamic RNA systems such as CUG trinucleotide repeat overexpansions.

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."

Crystal structure of group II intron domain 1 reveals a template for RNA assembly

Although the importance of large noncoding RNAs is increasingly appreciated, our understanding of their structures and architectural dynamics remains limited. In particular, we know little about RNA folding intermediates and how they facilitate the productive assembly of RNA tertiary structures. Here, we report the crystal structure of an obligate intermediate that is required during the earliest stages of group II intron folding. Composed of domain 1 from the Oceanobacillus iheyensis group II intron (266 nucleotides), this intermediate retains native-like features but adopts a compact conformation in which the active site cleft is closed. Transition between this closed and the open (native) conformation is achieved through discrete rotations of hinge motifs in two regions of the molecule. The open state is then stabilized by sequential docking of downstream intron domains, suggesting a 'first come, first folded' strategy that may represent a generalizable pathway for assembly of large RNA and ribonucleoprotein structures.

Computational investigation of RNA CUG repeats responsible for myotonic dystrophy 1

Despite the importance of the knowledge of molecular hydration entropy (ΔShyd) in chemical and biological processes, the exact calculation of ΔShyd is very difficult, because of the complexity in solute–water interactions. Although free-energy perturbation (FEP) methods have been employed quite widely in the literature, the poor convergent behavior of the van der Waals interaction term in the potential function limited the accuracy and robustness. In this study, we propose a new method for estimating ΔShyd by means of combining the FEP approach and the scaled particle theory (or information theory) to separately calculate the electrostatic solute–water interaction term (ΔSelec) and the hydrophobic contribution approximated by the cavity formation entropy (ΔScav), respectively. Decomposition of ΔShyd into ΔScav and ΔSelec terms is found to be very effective with a substantial accuracy enhancement in ΔShyd estimation, when compared to the conventional full FEP calculations. ΔScav appears to dominate over ΔSelec in magnitude, even in the case of polar solutes, implying that the major contribution to the entropic cost for hydration comes from the formation of a solvent-excluded volume. Our hybrid scaled particle theory and FEP method is thus found to enhance the accuracy of ΔShyd prediction by effectively complementing the conventional full FEP method.

RNA Structures as Mediators of Neurological Diseases and as Drug Targets

RNAs adopt diverse folded structures that are essential for function and thus play critical roles in cellular biology. A striking example of this is the ribosome, a complex, three-dimensionally folded macromolecular machine that orchestrates protein synthesis. Advances in RNA biochemistry, structural and molecular biology, and bioinformatics have revealed other non-coding RNAs whose functions are dictated by their structure. It is not surprising that aberrantly folded RNA structures contribute to disease. In this Review, we provide a brief introduction into RNA structural biology and then describe how RNA structures function in cells and cause or contribute to neurological disease. Finally, we highlight successful applications of rational design principles to provide chemical probes and lead compounds targeting structured RNAs. Based on several examples of well-characterized RNA-driven neurological disorders, we demonstrate how designed small molecules can facilitate the study of RNA dysfunction, elucidating previously unknown roles for RNA in disease, and provide lead therapeutics.

Crystallographic and Computational Analyses of AUUCU Repeating RNA That Causes Spinocerebellar Ataxia Type 10 (SCA10)

Spinocerebellar ataxia type 10 (SCA10) is caused by a pentanucleotide repeat expansion of r(AUUCU) within intron 9 of the ATXN10 pre-mRNA. The RNA causes disease by a gain-of-function mechanism in which it inactivates proteins involved in RNA biogenesis. Spectroscopic studies showed that r(AUUCU) repeats form a hairpin structure; however, there were no high-resolution structural models prior to this work. Herein, we report the first crystal structure of model r(AUUCU) repeats refined to 2.8 Å and analysis of the structure via molecular dynamics simulations. The r(AUUCU) tracts adopt an overall A-form geometry in which 3 × 3 nucleotide (5')UCU(3')/(3')UCU(5') internal loops are closed by AU pairs. Helical parameters of the refined structure as well as the corresponding electron density map on the crystallographic model reflect dynamic features of the internal loop. The computational analyses captured dynamic motion of the loop closing pairs, which can form single-stranded conformations with relatively low energies. Overall, the results presented here suggest the possibility for r(AUUCU) repeats to form metastable A-from structures, which can rearrange into single-stranded conformations and attract proteins such as heterogeneous nuclear ribonucleoprotein K (hnRNP K). The information presented here may aid in the rational design of therapeutics targeting this RNA.

Modifications to toxic CUG RNAs induce structural stability, rescue mis-splicing in a myotonic dystrophy cell model and reduce toxicity in a myotonic dystrophy zebrafish model

CUG repeat expansions in the 3′ UTR of dystrophia myotonica protein kinase (DMPK) cause myotonic dystrophy type 1 (DM1). As RNA, these repeats elicit toxicity by sequestering splicing proteins, such as MBNL1, into protein–RNA aggregates. Structural studies demonstrate that CUG repeats can form A-form helices, suggesting that repeat secondary structure could be important in pathogenicity. To evaluate this hypothesis, we utilized structure-stabilizing RNA modifications pseudouridine (Ψ) and 2′-O-methylation to determine if stabilization of CUG helical conformations affected toxicity. CUG repeats modified with Ψ or 2′-O-methyl groups exhibited enhanced structural stability and reduced affinity for MBNL1. Molecular dynamics and X-ray crystallography suggest a potential water-bridging mechanism for Ψ-mediated CUG repeat stabilization. Ψ modification of CUG repeats rescued mis-splicing in a DM1 cell model and prevented CUG repeat toxicity in zebrafish embryos. This study indicates that the structure of toxic RNAs has a significant role in controlling the onset of neuromuscular diseases.

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.

Targeting toxic RNAs that cause myotonic dystrophy type 1 (DM1) with a bisamidinium inhibitor

A working hypothesis for the pathogenesis of myotonic dystrophy type 1 (DM1) involves the aberrant sequestration of an alternative splicing regulator, MBNL1, by expanded CUG repeats, r(CUG)(exp). It has been suggested that a reversal of the myotonia and potentially other symptoms of the DM1 disease can be achieved by inhibiting the toxic MBNL1-r(CUG)(exp) interaction. Using rational design, we discovered an RNA-groove binding inhibitor (ligand 3) that contains two triaminotriazine units connected by a bisamidinium linker. Ligand 3 binds r(CUG)12 with a low micromolar affinity (K(d) = 8 ± 2 μM) and disrupts the MBNL1-r(CUG)12 interaction in vitro (K(i) = 8 ± 2 μM). In addition, ligand 3 is cell and nucleus permeable, exhibits negligible toxicity to mammalian cells, dissolves MBNL1-r(CUG)(exp) ribonuclear foci, and restores misregulated splicing of IR and cTNT in a DM1 cell culture model. Importantly, suppression of r(CUG)(exp) RNA-induced toxicity in a DM1 Drosophila model was observed after treatment with ligand 3. These results suggest ligand 3 as a lead for the treatment of DM1.

Structure of the myotonic dystrophy type 2 RNA and designed small molecules that reduce toxicity

Myotonic dystrophy type 2 (DM2) is an incurable neuromuscular disorder caused by a r(CCUG) expansion (r(CCUG)(exp)) that folds into an extended hairpin with periodically repeating 2×2 nucleotide internal loops (5'CCUG/3'GUCC). We designed multivalent compounds that improve DM2-associated defects using information about RNA-small molecule interactions. We also report the first crystal structure of r(CCUG) repeats refined to 2.35 Å. Structural analysis of the three 5'CCUG/3'GUCC repeat internal loops (L) reveals that the CU pairs in L1 are each stabilized by one hydrogen bond and a water-mediated hydrogen bond, while CU pairs in L2 and L3 are stabilized by two hydrogen bonds. Molecular dynamics (MD) simulations reveal that the CU pairs are dynamic and stabilized by Na(+) and water molecules. MD simulations of the binding of the small molecule to r(CCUG) repeats reveal that the lowest free energy binding mode occurs via the major groove, in which one C residue is unstacked and the cross-strand nucleotides are displaced. Moreover, we modeled the binding of our dimeric compound to two 5'CCUG/3'GUCC motifs, which shows that the scaffold on which the RNA-binding modules are displayed provides an optimal distance to span two adjacent loops.

RNA synthesis and purification for structural studies

RNAs play pivotal roles in the cell, ranging from catalysis (e.g., RNase P), acting as adaptor molecule (tRNA) to regulation (e.g., riboswitches). Precise understanding of its three-dimensional structures has given unprecedented insight into the molecular basis for all of these processes. Nevertheless, structural studies on RNA are still limited by the very special nature of this polymer. The most common methods for the determination of 3D RNA structures are NMR and X-ray crystallography. Both methods have their own set of requirements and give different amounts of information about the target RNA. For structural studies, the major bottleneck is usually obtaining large amounts of highly pure and homogeneously folded RNA. Especially for X-ray crystallography it can be necessary to screen a large number of variants to obtain well-ordered single crystals. In this mini-review we give an overview about strategies for the design, in vitro production, and purification of RNA for structural studies.

RNA and RNA–Protein Complex Crystallography and its Challenges

RNA biology has changed completely in the past decade with the discovery of non-coding RNAs. Unfortunately, obtaining mechanistic information about these RNAs alone or in cellular complexes with proteins has been a major problem. X-ray crystallography of RNA and RNA–protein complexes has suffered from the major problems encountered in preparing and purifying them in large quantity. Here, we review the available techniques and methods in vitro and in vivo used to prepare and purify RNA and RNA–protein complex for crystallographic studies. We also discuss the future directions necessary to explore the vast number of RNA species waiting for their atomic-resolution structure to be determined.