NMR Studies of HAR1 RNA Secondary Structures Reveal Conformational Dynamics in the Human RNA

The expression profile of HAR1A and HAR1B in the peripheral blood cells of multiple sclerosis patients

BACKGROUND

Multiple sclerosis (MS) is a progressive neurodegenerative disease of the central nervous system (CNS) with varying degrees of axonal and neuronal damage. The onset and progression of the disease are influenced by several environmental and genetic variables. Long non-coding RNAs (lncRNAs) have a crucial role in the pathophysiology of MS. Our study aimed to assess the levels of HAR1A and HAR1B lncRNA expression in the blood samples of MS patients and investigate the relationship between these lncRNAs and disease activity.

METHODS AND RESULTS

The blood samples of 100 MS patients, including 82 relapsing-remitting (RR), 8 primary progressive (PP), and 10 secondary progressive (SP) MS cases, and 100 healthy controls were collected. Quantitative real-time PCR was used for the evaluation of gene expression. ROC curve analysis was performed to evaluate the diagnostic potential of lncRNA levels. A significant decrease was detected in HAR1A expressions (P < 0.0001), and a moderate increase was also shown in HAR1B of SPMS patients (P value = 0.0189). HAR1A showed different expression levels in patients over forty (P value = 0.034). The expression levels of HAR1A and HAR1B were positively correlated in MS patients (r = 0.2003, P value = 0.0457). In addition, ROC curve results suggested that HAR1A can be introduced as a novel biomarker for MS diagnosis (AUC = 0.776).

CONCLUSION

The low serum level of HAR1A may be a potential molecular biomarker for MS diagnosis; however, no discernible difference was detected in the expression level of HAR1B in the blood samples of MS patients.

Evolutionary Conservation of RNA Secondary Structure

Noncoding RNAs, ncRNAs, naturally fold into structures, which allow them to perform their functions in the cell. Evolutionarily close species share structures and functions. This occurs because of shared selective pressures, resulting in conserved groups. Previous efforts in finding functional RNAs have been made in detecting conserved structures in genomes or alignments. It may occur that, within a conserved group, species-specific structures arise after species split due to positive selection. Detecting positive selection in ncRNAs is a hard problem in biology as well as bioinformatics. To detect positive selection, one should find species-specific structures within a conserved set. This chapter provides protocols to detect and analyze positive selection in ncRNA structures with the SSS-test and other free software.

RNA Chemical Labeling with Site-Specific, Relative Quantification by Mass Spectrometry for the Structural Study of a Neomycin-Sensing Riboswitch Aptamer Domain

High-resolution mass spectrometry was used for the label-free, direct localization and relative quantification of CMC+ -modifications of a neomycin-sensing riboswitch aptamer domain in the absence and presence of the aminoglycoside ligands neomycin B, ribostamycin, and paromomycin. The chemical probing and MS data for the free riboswitch show high exposure to solvent of the uridine nucleobases U7, U8, U13, U14, U18 as part of the proposed internal and apical loops, but those of U10 and U21 as part of the proposed internal loop were found to be far less exposed than expected. Thus, our data are in better agreement with the proposed secondary structure of the riboswitch in complexes with aminoglycosides than with that of free RNA. For the riboswitch in complexes with neomycin B, ribostamycin, and paromomycin, we found highly similar CMC+ -modification patterns and excellent agreement with previous NMR studies. Differences between the chemical probing and MS data in the absence and presence of the aminoglycoside ligands were quantitative rather than qualitative (i. e., the same nucleobases were labeled, but to different extents) and can be rationalized by stabilization of both the proposed bulge and the apical loop by aminoglycoside binding. Our study shows that chemical probing and mass spectrometry can provide important structural information and complement other techniques such as NMR spectroscopy.

An integrative NMR-SAXS approach for structural determination of large RNAs defines the substrate-free state of a trans-cleaving Neurospora Varkud Satellite ribozyme

The divide-and-conquer strategy is commonly used for protein structure determination, but its applications to high-resolution structure determination of RNAs have been limited. Here, we introduce an integrative approach based on the divide-and-conquer strategy that was undertaken to determine the solution structure of an RNA model system, the Neurospora VS ribozyme. NMR and SAXS studies were conducted on a minimal trans VS ribozyme as well as several isolated subdomains. A multi-step procedure was used for structure determination that first involved pairing refined NMR structures with SAXS data to obtain structural subensembles of the various subdomains. These subdomain structures were then assembled to build a large set of structural models of the ribozyme, which was subsequently filtered using SAXS data. The resulting NMR-SAXS structural ensemble shares several similarities with the reported crystal structures of the VS ribozyme. However, a local structural difference is observed that affects the global fold by shifting the relative orientation of the two three-way junctions. Thus, this finding highlights a global conformational change associated with substrate binding in the VS ribozyme that is likely critical for its enzymatic activity. Structural studies of other large RNAs should benefit from similar integrative approaches that allow conformational sampling of assembled fragments.

HAR1: an insight into lncRNA genetic evolution

Long noncoding RNAs (lncRNAs) have a wide range of functions in health and disease, but many remain uncharacterized because of their complex expression patterns and structures. The genetic loci encoding lncRNAs can be subject to accelerated evolutionary changes within the human lineage. HAR1 is a region that has a significantly altered sequence compared to other primates and is a component of two overlapping lncRNA loci, HAR1A and HAR1B. Although the functions of these lncRNAs are unknown, they have been associated with neurological disorders and cancer. Here, we explore the current state of understanding of evolution in human lncRNA genes, using the HAR1 locus as the case study.

NMR Spectroscopy of Large Functional RNAs: From Sample Preparation to Low-Gamma Detection

NMR spectroscopy is a potent method for the structural and biophysical characterization of RNAs. The application of NMR spectroscopy is restricted in RNA size and most often requires isotope-labeled or even selectively labeled RNAs. Additionally, new NMR pulse sequences, such as the heteronuclear-detected NMR experiments, are introduced. We herein provide detailed protocols for the preparation of isotope-labeled RNA for NMR spectroscopy via in vitro transcription. This protocol covers all steps, from the preparation of DNA template to the transcription of milligram RNA quantities. Moreover, we present a protocol for a chemo-enzymatic approach to introduce a single modified nucleotide at any position of any RNA. Regarding NMR methodology, we share protocols for the implementation of a suite of heteronuclear-detected NMR experiments including ¹³ C-detected experiments for ribose assignment and amino groups, the CN-spin filter heteronuclear single quantum coherence (HSQC) for imino groups and the ¹⁵ N-detected band-selective excitation short transient transverse-relaxation-optimized spectroscopy (BEST-TROSY) experiment. © 2020 The Authors. Basic Protocol 1: Preparation of isotope-labeled RNA samples with in vitro transcription using T7 RNAP, DEAE chromatography, and RP-HPLC purification Alternate Protocol 1: Purification of isotope-labeled RNA from in vitro transcription with preparative PAGE Alternate Protocol 2: Purification of isotope-labeled RNA samples from in vitro transcription via centrifugal concentration Support Protocol 1: Preparation of DNA template from plasmid Support Protocol 2: Preparation of PCR DNA as template Support Protocol 3: Preparation of T7 RNA Polymerase (T7 RNAP) Support Protocol 4: Preparation of yeast inorganic pyrophosphatase (YIPP) Basic Protocol 2: Preparation of site-specific labeled RNAs using a chemo-enzymatic synthesis Support Protocol 5: Synthesis of modified nucleoside 3',5'-bisphosphates Support Protocol 6: Preparation of T4 RNA Ligase 2 Support Protocol 7: Setup of NMR spectrometer for heteronuclear-detected NMR experiments Support Protocol 8: IPAP and DIPAP for homonuclear decoupling Basic Protocol 3: ¹³ C-detected 3D (H)CC-TOCSY, (H)CPC, and (H)CPC-CCH-TOCSY experiments for ribose assignment Basic Protocol 4: ¹³ C-detected 2D CN-spin filter HSQC experiment Basic Protocol 5: ¹³ C-detected C(N)H-HDQC experiment for the detection of amino groups Support Protocol 9: ¹³ C-detected CN-HSQC experiment for amino groups Basic Protocol 6: ¹³ C-detected "amino"-NOESY experiment Basic Protocol 7: ¹⁵ N-detected BEST-TROSY experiment.

Synthesis, purification and crystallization of a putative critical bulge of HAR1 RNA

Non-coding RNAs have raised a lot of interest because of their capabilities to perform enzymatic reactions and regulate gene expression in various ways. Human Accelerated Region 1 (HAR1) has been identified during the search for highly conserved regions in mammalian genomes, over one hundred base pairs long, and with high rates of substitution in the human genome. Its potential for coding for a protein is very minimal. However, the HAR1 transcript has been computationally predicted to have a stable secondary structure. Previous structure-probing experiments have suggested that the majority of differences between human and chimp constructs are in helices, designated C and D. For this reason, a 47nt construct consisting of the C and D helices along with two additional C-G pairs was synthesized, purified, and crystallized, and its x-ray structure is reported in this study. The final structure is an artificial dimer, with a bulge that forms different conformations on each monomer. This bulge has been observed in predicted secondary structures, footprinting assays, enzymatic degradation assays, NMR studies, in silico studies, and in this crystalized dimer structure. It is proposed that the HAR1 transcript is a non-coding RNA that interacts with an unknown binding partner responsible for brain development through this inherent structural motif of bulged adenosines.

Pairwise visual comparison of small RNA secondary structures with base pair probabilities

BACKGROUND

Predicted RNA secondary structures are typically visualized using dot-plots for base pair binding probabilities and planar graphs for unique structures, such as the minimum free energy structure. These are however difficult to analyze simultaneously.

RESULTS

This work introduces a compact unified view of the most stable conformation of an RNA secondary structure and its base pair probabilities, which is called the Circular Secondary Structure Base Pairs Probabilities Plot (CS2BP2-Plot). Along with our design we provide access to a web server implementation of our solution that facilitates pairwise comparison of short RNA (and DNA) sequences up to 200 base pairs. The web server first calculates the minimum free energy secondary structure and the base pair probabilities for up to 10 RNA or DNA sequences using RNAfold and then provides a two panel comparative view that includes CS2BP2-Plots along with the traditional graph, planar and circular diagrams obtained with VARNA. The CS2BP2-Plots include highlighting of the nucleotide differences between two selected sequences using ClustalW local alignments. We also provide descriptive statistics, dot-bracket secondary structure representations and ClustalW local alignments for compared sequences.

CONCLUSIONS

Using circular diagrams and colour and weight-coded arcs, we demonstrate how a single image can replace the state-of-the-art dual representations (dot-plots and minimum free energy structures) for base-pair probabilities of RNA secondary structures while allowing efficient exploration and comparison of different RNA conformations via a web server front end. With that, we provide the community, especially the biologically oriented, with an intuitive tool for ncRNA visualization. Web-server: https://nrcmonsrv01.nrc.ca/cs2bp2plot.

Advances that facilitate the study of large RNA structure and dynamics by nuclear magnetic resonance spectroscopy

The characterization of functional yet nonprotein coding (nc) RNAs has expanded the role of RNA in the cell from a passive player in the central dogma of molecular biology to an active regulator of gene expression. The misregulation of ncRNA function has been linked with a variety of diseases and disorders ranging from cancers to neurodegeneration. However, a detailed molecular understanding of how ncRNAs function has been limited; due, in part, to the difficulties associated with obtaining high-resolution structures of large RNAs. Tertiary structure determination of RNA as a whole is hampered by various technical challenges, all of which are exacerbated as the size of the RNA increases. Namely, RNAs tend to be highly flexible and dynamic molecules, which are difficult to crystallize. Biomolecular nuclear magnetic resonance (NMR) spectroscopy offers a viable alternative to determining the structure of large RNA molecules that do not readily crystallize, but is itself hindered by some technical limitations. Recently, a series of advancements have allowed the biomolecular NMR field to overcome, at least in part, some of these limitations. These advances include improvements in sample preparation strategies as well as methodological improvements. Together, these innovations pave the way for the study of ever larger RNA molecules that have important biological function. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

Maximizing accuracy of RNA structure in refinement against residual dipolar couplings

Structural information about ribonucleic acid (RNA) is lagging behind that of proteins, in part due to its high charge and conformational variability. Molecular dynamics (MD) has played an important role in describing RNA structure, complementing information from both nuclear magnetic resonance (NMR), or X-ray crystallography. We examine the impact of the choice of the empirical force field for RNA structure refinement using cross-validation against residual dipolar couplings (RDCs) as structural accuracy reporter. Four force fields, representing both the state-of-the art in RNA simulation and the most popular selections in NMR structure determination, are compared for a prototypical A-RNA helix. RNA structural accuracy is also evaluated as a function of both density and nature of input NMR data including RDCs, anisotropic chemical shifts, and distance restraints. Our results show a complex interplay between the experimental restraints and the force fields indicating two best-performing choices: high-fidelity refinement in explicit solvent, and the conformational database-derived potentials. Accuracy of RNA models closely tracks the density of 1-bond C-H RDCs, with other data types having beneficial, but smaller effects. At lower RDC density, or when refining against NOEs only, the two selected force fields are capable of accurately describing RNA helices with little or no experimental RDC data, making them available for the higher order structure assembly or better quantification of the intramolecular dynamics. Unrestrained simulations of simple RNA motifs with state-of-the art MD force fields appear to capture the flexibility inherent in nucleic acids while also maintaining a good agreement with the experimental observables.

SSS-test: a novel test for detecting positive selection on RNA secondary structure

BACKGROUND

Long non-coding RNAs (lncRNAs) play an important role in regulating gene expression and are thus important for determining phenotypes. Most attempts to measure selection in lncRNAs have focused on the primary sequence. The majority of small RNAs and at least some parts of lncRNAs must fold into specific structures to perform their biological function. Comprehensive assessments of selection acting on RNAs therefore must also encompass structure. Selection pressures acting on the structure of non-coding genes can be detected within multiple sequence alignments. Approaches of this type, however, have so far focused on negative selection. Thus, a computational method for identifying ncRNAs under positive selection is needed.

RESULTS

We introduce the SSS-test (test for Selection on Secondary Structure) to identify positive selection and thus adaptive evolution. Benchmarks with biological as well as synthetic controls yield coherent signals for both negative and positive selection, demonstrating the functionality of the test. A survey of a lncRNA collection comprising 15,443 families resulted in 110 candidates that appear to be under positive selection in human. In 26 lncRNAs that have been associated with psychiatric disorders we identified local structures that have signs of positive selection in the human lineage.

CONCLUSIONS

It is feasible to assay positive selection acting on RNA secondary structures on a genome-wide scale. The detection of human-specific positive selection in lncRNAs associated with cognitive disorder provides a set of candidate genes for further experimental testing and may provide insights into the evolution of cognitive abilities in humans.

AVAILABILITY

The SSS-test and related software is available at: https://github.com/waltercostamb/SSS-test . The databases used in this work are available at: http://www.bioinf.uni-leipzig.de/Software/SSS-test/ .

Temporal ordering of substitutions in RNA evolution: Uncovering the structural evolution of the Human Accelerated Region 1

The Human Accelerated Region 1 (HAR1) is the most rapidly evolving region in the human genome. It is part of two overlapping long non-coding RNAs, has a length of only 118 nucleotides and features 18 human specific changes compared to an ancestral sequence that is extremely well conserved across non-human primates. The human HAR1 forms a stable secondary structure that is strikingly different from the one in chimpanzee as well as other closely related species, again emphasizing its human-specific evolutionary history. This suggests that positive selection has acted to stabilize human-specific features in the ensemble of HAR1 secondary structures. To investigate the evolutionary history of the human HAR1 structure, we developed a computational model that evaluates the relative likelihood of evolutionary trajectories as a probabilistic version of a Hamiltonian path problem. The model predicts that the most likely last step in turning the ancestral primate HAR1 into the human HAR1 was exactly the substitution that distinguishes the modern human HAR1 sequence from that of Denisovan, an archaic human, providing independent support for our model. The MutationOrder software is available for download and can be applied to other instances of RNA structure evolution.

Ligand-modulated folding of the full-length adenine riboswitch probed by NMR and single-molecule FRET spectroscopy

The full-length translation-regulating add adenine riboswitch (Asw) from Vibrio vulnificus has a more complex conformational space than its isolated aptamer domain. In addition to the predicted apo (apoA) and holo conformation that feature the conserved three-way junctional purine riboswitch aptamer, it adopts a second apo (apoB) conformation with a fundamentally different secondary structure. Here, we characterized the ligand-dependent conformational dynamics of the full-length add Asw by NMR and by single-molecule FRET (smFRET) spectroscopy. Both methods revealed an adenine-induced secondary structure switch from the apoB-form to the apoA-form that involves no tertiary structural interactions between aptamer and expression platform. This strongly suggests that the add Asw triggers translation by capturing the apoA-form secondary structure in the holo state. Intriguingly, NMR indicated a homogenous, docked aptamer kissing loop fold for apoA and holo, while smFRET showed persistent aptamer kissing loop docking dynamics between comparably stable, undocked and docked substates of the apoA and the holo conformation. Unraveling the folding of large junctional riboswitches thus requires the integration of complementary solution structural techniques such as NMR and smFRET.

Evaluation of 15N-detected H-N correlation experiments on increasingly large RNAs

Recently, ¹⁵N-detected multidimensional NMR experiments have been introduced for the investigation of proteins. Utilization of the slow transverse relaxation of nitrogen nuclei in a ¹⁵N-TROSY experiment allowed recording of high quality spectra for high molecular weight proteins, even in the absence of deuteration. Here, we demonstrate the applicability of three ¹⁵N-detected H-N correlation experiments (TROSY, BEST-TROSY and HSQC) to RNA. With the newly established ¹⁵N-detected BEST-TROSY experiment, which proves to be the most sensitive ¹⁵N-detected H-N correlation experiment, spectra for five RNA molecules ranging in size from 5 to 100 kDa were recorded. These spectra yielded high resolution in the ¹⁵N-dimension even for larger RNAs since the increase in line width with molecular weight is more pronounced in the ¹H- than in the ¹⁵N-dimension. Further, we could experimentally validate the difference in relaxation behavior of imino groups in AU and GC base pairs. Additionally, we showed that ¹⁵N-detected experiments theoretically should benefit from sensitivity and resolution advantages at higher static fields but that the latter is obscured by exchange dynamics within the RNAs.

Nuclear Magnetic Resonance-Assisted Prediction of Secondary Structure for RNA: Incorporation of Direction-Dependent Chemical Shift Constraints

Knowledge of RNA structure is necessary to determine structure–function relationships and to facilitate design of potential therapeutics. RNA secondary structure prediction can be improved by applying constraints from nuclear magnetic resonance (NMR) experiments to a dynamic programming algorithm. Imino proton walks from NOESY spectra reveal double-stranded regions. Chemical shifts of protons in GH1, UH3, and UH5 of GU pairs, UH3, UH5, and AH2 of AU pairs, and GH1 of GC pairs were analyzed to identify constraints for the 5′ to 3′ directionality of base pairs in helices. The 5′ to 3′ directionality constraints were incorporated into an NMR-assisted prediction of secondary structure (NAPSS-CS) program. When it was tested on 18 structures, including nine pseudoknots, the sensitivity and positive predictive value were improved relative to those of three unrestrained programs. The prediction accuracy for the pseudoknots improved the most. The program also facilitates assignment of chemical shifts to individual nucleotides, a necessary step for determining three-dimensional structure.

Signatures of accelerated somatic evolution in gene promoters in multiple cancer types

Cancer-associated somatic mutations outside protein-coding regions remain largely unexplored. Analyses of the TERT locus have indicated that non-coding regulatory mutations can be more frequent than previously suspected and play important roles in oncogenesis. Using a computational method called SASE-hunter, developed here, we identified a novel signature of accelerated somatic evolution (SASE) marked by a significant excess of somatic mutations localized in a genomic locus, and prioritized those loci that carried the signature in multiple cancer patients. Interestingly, even when an affected locus carried the signature in multiple individuals, the mutations contributing to SASE themselves were rarely recurrent at the base-pair resolution. In a pan-cancer analysis of 906 samples from 12 tumor types, we detected SASE in the promoters of several genes, including known cancer genes such as MYC, BCL2, RBM5 and WWOX. Nucleotide substitution patterns consistent with oxidative DNA damage and local somatic hypermutation appeared to contribute to this signature in selected gene promoters (e.g. MYC). SASEs in selected cancer gene promoters were associated with over-expression, and also correlated with the age of onset of cancer, aggressiveness of the disease and survival. Taken together, our work detects a hitherto under-appreciated and clinically important class of regulatory changes in cancer genomes.

The role of gene regulatory factors in the evolutionary history of humans

Deciphering the molecular basis of how modern human phenotypes have evolved is one of the most fascinating challenges in biology. Here, we will focus on the roles of gene regulatory factors (GRFs), in particular transcription factors (TFs) and long non-coding RNAs (lncRNAs) during human evolution. We will present examples of TFs and lncRNAs that have changed or show signs of positive selection in humans compared to chimpanzees, in modern humans compared to archaic humans, or within modern human populations. On the basis of current knowledge about the functions of these GRF genes, we speculate that they have been involved in speciation as well as in shaping phenotypes such as brain functions, skeletal morphology, and metabolic processes.

Structure determination of noncanonical RNA motifs guided by ¹H NMR chemical shifts

Structured noncoding RNAs underlie fundamental cellular processes, but determining their three-dimensional structures remains challenging. We demonstrate that integrating ¹H NMR chemical shift data with Rosetta de novo modeling can be used to consistently determine high-resolution RNA structures. On a benchmark set of 23 noncanonical RNA motifs, including 11 'blind' targets, chemical-shift Rosetta for RNA (CS-Rosetta-RNA) recovered experimental structures with high accuracy (0.6-2.0 Å all-heavy-atom r.m.s. deviation) in 18 cases.

Tackling structures of long noncoding RNAs

RNAs are important catalytic machines and regulators at every level of gene expression. A new class of RNAs has emerged called long non-coding RNAs, providing new insights into evolution, development and disease. Long non-coding RNAs (lncRNAs) predominantly found in higher eukaryotes, have been implicated in the regulation of transcription factors, chromatin-remodeling, hormone receptors and many other processes. The structural versatility of RNA allows it to perform various functions, ranging from precise protein recognition to catalysis and metabolite sensing. While major housekeeping RNA molecules have long been the focus of structural studies, lncRNAs remain the least characterized class, both structurally and functionally. Here, we review common methodologies used to tackle RNA structure, emphasizing their potential application to lncRNAs. When considering the complexity of lncRNAs and lack of knowledge of their structure, chemical probing appears to be an indispensable tool, with few restrictions in terms of size, quantity and heterogeneity of the RNA molecule. Probing is not constrained to in vitro analysis and can be adapted to high-throughput sequencing platforms. Significant efforts have been applied to develop new in vivo chemical probing reagents, new library construction protocols for sequencing platforms and improved RNA prediction software based on the experimental evidence.

Long noncoding RNAs in development and disease of the central nervous system

The central nervous system (CNS) is a complex biological system composed of numerous cell types working in concert. The intricate development and functioning of this highly ordered structure depends upon exquisite spatial and temporal control of gene expression in the cells comprising the CNS. Thus, gene regulatory networks that control cell fates and functions play critical roles in the CNS. Failure to develop and maintain intricate regulatory networks properly leads to impaired development or neural dysfunction, which might manifest as neurological disorders. Long noncoding RNAs (lncRNAs) are emerging as important components of gene regulatory networks, working in concert with transcription factors and epigenetic regulators of gene expression. Interestingly, many lncRNAs are highly expressed in the adult and developing brain, often showing precise temporal and spatial patterns of expression. This specificity of expression and growing awareness of the importance of lncRNAs suggest that they play key roles in CNS development and function. In this review, we highlight the growing evidence for the importance of lncRNAs in the CNS and the indications that their dysregulation underlies some neurological disorders.