Strategies for RNA folding and assembly

The FXR1 network acts as signaling scaffold for actomyosin remodeling

It is currently not known whether mRNAs fulfill structural roles in the cytoplasm. Here, we report the FXR1 network, an mRNA-protein (mRNP) network present throughout the cytoplasm, formed by FXR1-mediated packaging of exceptionally long mRNAs. These mRNAs serve as underlying condensate scaffold and concentrate FXR1 molecules. The FXR1 network contains multiple protein binding sites and functions as a signaling scaffold for interacting proteins. We show that it is necessary for RhoA signaling-induced actomyosin reorganization to provide spatial proximity between kinases and their substrates. Point mutations in FXR1, found in its homolog FMR1, where they cause Fragile X syndrome, disrupt the network. FXR1 network disruption prevents actomyosin remodeling-an essential and ubiquitous process for the regulation of cell shape, migration, and synaptic function. These findings uncover a structural role for cytoplasmic mRNA and show how the FXR1 RNA-binding protein as part of the FXR1 network acts as organizer of signaling reactions.

Structure of the SARS-CoV-2 Frameshift Stimulatory Element with an Upstream Multibranch Loop

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) frameshift stimulatory element (FSE) is necessary for programmed -1 ribosomal frameshifting (-1 PRF) and optimized viral efficacy. The FSE has an abundance of context-dependent alternate conformations, but two of the structures most crucial to -1 PRF are an attenuator hairpin and a three-stem H-type pseudoknot structure. A crystal structure of the pseudoknot alone features three RNA stems in a helically stacked linear structure, whereas a 6.9 Å cryo-EM structure including the upstream heptameric slippery site resulted in a bend between two stems. Our previous research alluded to an extended upstream multibranch loop that includes both the attenuator hairpin and the slippery site-a conformation not previously modeled. We aim to provide further context to the SARS-CoV-2 FSE via computational and medium resolution cryo-EM approaches, by presenting a 6.1 Å cryo-EM structure featuring a linear pseudoknot structure and a dynamic upstream multibranch loop.

Determinant of m6A regional preference by transcriptional dynamics

N6-Methyladenosine (m6A) is the most abundant chemical modification occurring on eukaryotic mRNAs, and has been reported to be involved in almost all stages of mRNA metabolism. The distribution of m6A sites is notably asymmetric along mRNAs, with a strong preference toward the 3' terminus of the transcript. How m6A regional preference is shaped remains incompletely understood. In this study, by performing m6A-seq on chromatin-associated RNAs, we found that m6A regional preference arises during transcription. Nucleosome occupancy is remarkedly increased in the region downstream of m6A sites, suggesting an intricate interplay between m6A methylation and nucleosome-mediated transcriptional dynamics. Notably, we found a remarkable slowdown of Pol-II movement around m6A sites. In addition, inhibiting Pol-II movement increases nearby m6A methylation levels. By analyzing massively parallel assays for m6A, we found that RNA secondary structures inhibit m6A methylation. Remarkably, the m6A sites associated with Pol-II pausing tend to be embedded within RNA secondary structures. These results suggest that Pol-II pausing could affect the accessibility of m6A motifs to the methyltransferase complex and subsequent m6A methylation by mediating RNA secondary structure. Overall, our study reveals a crucial role of transcriptional dynamics in the formation of m6A regional preference.

Crystallization as a selection force at the polymerization of nucleotides in a prebiotic context

Accumulation and selection of nucleotides is one of the most challenging problems surrounding the origin of the first RNA molecules on our planet. In the current work we propose that guanosine 3',5' cyclic monophosphate could selectively crystallize upon evaporation of an acidic prebiotic pool containing various other nucleotides. The conditions of the evaporative crystallization are fully compatible with the subsequent acid catalyzed polymerization of this cyclic nucleotide reported in earlier studies and may be relevant in a broad range of possible prebiotic environments. Albeit cytidine 3',5' cyclic monophosphate has the ability to selectively accumulate under the same conditions, its crystal structure is not likely to support polymer formation.

Intrinsically Disordered Proteins: An Overview

Many proteins and protein segments cannot attain a single stable three-dimensional structure under physiological conditions; instead, they adopt multiple interconverting conformational states. Such intrinsically disordered proteins or protein segments are highly abundant across proteomes, and are involved in various effector functions. This review focuses on different aspects of disordered proteins and disordered protein regions, which form the basis of the so-called "Disorder-function paradigm" of proteins. Additionally, various experimental approaches and computational tools used for characterizing disordered regions in proteins are discussed. Finally, the role of disordered proteins in diseases and their utility as potential drug targets are explored.

Alternative RNA Conformations: Companion or Combatant

RNA molecules, in one form or another, are involved in almost all aspects of cell physiology, as well as in disease development. The diversity of the functional roles of RNA comes from its intrinsic ability to adopt complex secondary and tertiary structures, rivaling the diversity of proteins. The RNA molecules form dynamic ensembles of many interconverting conformations at a timescale of seconds, which is a key for understanding how they execute their cellular functions. Given the crucial role of RNAs in various cellular processes, we need to understand the RNA molecules from a structural perspective. Central to this review are studies aimed at revealing the regulatory role of conformational equilibria in RNA in humans to understand genetic diseases such as cancer and neurodegenerative diseases, as well as in pathogens such as bacteria and viruses so as to understand the progression of infectious diseases. Furthermore, we also summarize the prior studies on the use of RNA structures as platforms for the rational design of small molecules for therapeutic applications.

Role of NS2 specific RNA binding and phosphorylation in liquid-liquid phase separation and virus assembly

Liquid-liquid phase separation (LLPS) has assumed a prominent role in biological cell systems, where it underpins the formation of subcellular compartments necessary for cell function. We investigated the underlying mechanism of LLPS in virus infected cells, where virus inclusion bodies are formed by an RNA-binding phosphoprotein (NS2) of Bluetongue virus to serve as sites for subviral particle assembly and virus maturation. We show that NS2 undergoes LLPS that is dependent on protein phosphorylation and RNA-binding and that LLPS occurrence is accompanied by a change in protein secondary structure. Site-directed mutagenesis identified two critical arginine residues in NS2 responsible for specific RNA binding and thus for NS2-RNA complex driven LLPS. Reverse genetics identified the same residues as essential for VIB assembly in infected cells and virus viability. Our findings suggest that a specific arginine-RNA interaction in the context of a phosphorylated state drives LLPS in this, and possibly other, virus infections.

An in vitro cytomimetic of in-cell RNA folding

To discover the cytomimetic that accounts for cytoplasmic crowding and sticking on RNA stability, we conducted a two-dimensional scan of mixtures of artificial crowding and sticking agents, PEG10k and M-PERTM. As our model RNA, we investigate the fourU RNA thermometer motif of Salmonella, a hairpin-structured RNA that regulates translation by unfolding and exposing its RBS in response to temperature perturbations. We found that the addition of artificial crowding and sticking agents leads to a stabilization and destabilization of RNA folding, respectively, through the excluded volume effect and surface interactions. FRET-labels were added to the fourU RNA and Fast Relaxation Imaging (FReI), fluorescence microscopy coupled to temperature-jump spectroscopy, probed differences between folding stability of RNA inside single living cells and in vitro. Our results suggest that the cytoplasmic environment affecting RNA folding is comparable to a combination of 20% v/v M-PERTM and 150 g/L PEG10k.

Prp43/DHX15 exemplify RNA helicase multifunctionality in the gene expression network

Dynamic regulation of RNA folding and structure is critical for the biogenesis and function of RNAs and ribonucleoprotein (RNP) complexes. Through their nucleotide triphosphate-dependent remodelling functions, RNA helicases are key modulators of RNA/RNP structure. While some RNA helicases are dedicated to a specific target RNA, others are multifunctional and engage numerous substrate RNAs in different aspects of RNA metabolism. The discovery of such multitasking RNA helicases raises the intriguing question of how these enzymes can act on diverse RNAs but also maintain specificity for their particular targets within the RNA-dense cellular environment. Furthermore, the identification of RNA helicases that sit at the nexus between different aspects of RNA metabolism raises the possibility that they mediate cross-regulation of different cellular processes. Prominent and extensively characterized multifunctional DEAH/RHA-box RNA helicases are DHX15 and its Saccharomyces cerevisiae (yeast) homologue Prp43. Due to their central roles in key cellular processes, these enzymes have also served as prototypes for mechanistic studies elucidating the mode of action of this type of enzyme. Here, we summarize the current knowledge on the structure, regulation and cellular functions of Prp43/DHX15, and discuss the general concept and implications of RNA helicase multifunctionality.

Orchestrating ribosomal RNA folding during ribosome assembly

Construction of the eukaryotic ribosome is a complex process in which a nascent ribosomal RNA (rRNA) emerging from RNA Polymerase I hierarchically folds into a native three-dimensional structure. Modular assembly of individual RNA domains through interactions with ribosomal proteins and a myriad of assembly factors permit efficient disentanglement of the error-prone RNA folding process. Following these dynamic events, long-range tertiary interactions are orchestrated to compact rRNA. A combination of genetic, biochemical, and structural studies is now providing clues into how a nascent rRNA is transformed into a functional ribosome with high precision. With this essay, we aim to draw attention to the poorly understood process of establishing correct RNA tertiary contacts during ribosome formation.

Diffuse Ions Coordinate Dynamics in a Ribonucleoprotein Assembly

Proper ionic concentrations are required for the functional dynamics of RNA and ribonucleoprotein (RNP) assemblies. While experimental and computational techniques have provided many insights into the properties of chelated ions, less is known about the energetic contributions of diffuse ions to large-scale conformational rearrangements. To address this, we present a model that is designed to quantify the influence of diffuse monovalent and divalent ions on the dynamics of biomolecular assemblies. This model employs all-atom (non-H) resolution and explicit ions, where effective potentials account for hydration effects. We first show that the model accurately predicts the number of excess Mg²⁺ ions for prototypical RNA systems, at a level comparable to modern coarse-grained models. We then apply the model to a complete ribosome and show how the balance between diffuse Mg²⁺ and K⁺ ions can control the dynamics of tRNA molecules during translation. The model predicts differential effects of diffuse ions on the free-energy barrier associated with tRNA entry and the energy of tRNA binding to the ribosome. Together, this analysis reveals the direct impact of diffuse ions on the dynamics of an RNP assembly.

Protein cofactors and substrate influence Mg2+-dependent structural changes in the catalytic RNA of archaeal RNase P

The ribonucleoprotein (RNP) form of archaeal RNase P comprises one catalytic RNA and five protein cofactors. To catalyze Mg²⁺-dependent cleavage of the 5′ leader from pre-tRNAs, the catalytic (C) and specificity (S) domains of the RNase P RNA (RPR) cooperate to recognize different parts of the pre-tRNA. While ∼250–500 mM Mg²⁺ renders the archaeal RPR active without RNase P proteins (RPPs), addition of all RPPs lowers the Mg²⁺ requirement to ∼10–20 mM and improves the rate and fidelity of cleavage. To understand the Mg²⁺- and RPP-dependent structural changes that increase activity, we used pre-tRNA cleavage and ensemble FRET assays to characterize inter-domain interactions in Pyrococcus furiosus (Pfu) RPR, either alone or with RPPs ± pre-tRNA. Following splint ligation to doubly label the RPR (Cy3-RPR_{C domain} and Cy5-RPR_{S domain}), we used native mass spectrometry to verify the final product. We found that FRET correlates closely with activity, the Pfu RPR and RNase P holoenzyme (RPR + 5 RPPs) traverse different Mg²⁺-dependent paths to converge on similar functional states, and binding of the pre-tRNA by the holoenzyme influences Mg²⁺ cooperativity. Our findings highlight how Mg²⁺ and proteins in multi-subunit RNPs together favor RNA conformations in a dynamic ensemble for functional gains.

Formation of frameshift-stimulating RNA pseudoknots is facilitated by remodeling of their folding intermediates

Programmed –1 ribosomal frameshifting is an essential regulation mechanism of translation in viruses and bacteria. It is stimulated by mRNA structures inside the coding region. As the structure is unfolded repeatedly by consecutive translating ribosomes, whether it can refold properly each time is important in performing its function. By using single-molecule approaches and molecular dynamics simulations, we found that a frameshift-stimulating RNA pseudoknot folds sequentially through its upstream stem S1 and downstream stem S2. In this pathway, S2 folds from the downstream side and tends to be trapped in intermediates. By masking the last few nucleotides to mimic their gradual emergence from translating ribosomes, S2 can be directed to fold from the upstream region. The results show that the intermediates are greatly suppressed, suggesting that mRNA refolding may be modulated by ribosomes. Moreover, masking the first few nucleotides of S1 favors the folding from S2 and yields native pseudoknots, which are stable enough to retrieve the masked nucleotides. We hypothesize that translating ribosomes can remodel an intermediate mRNA structure into a stable conformation, which may in turn stimulate backward slippage of the ribosome. This supports an interactive model of ribosomal frameshifting and gives an insightful account addressing previous experimental observations.

Dynamically probing ATP-dependent RNA helicase A-assisted RNA structure conversion using single molecule fluorescence resonance energy transfer

RNA helicase A (RHA) as a member of DExH-box subgroup of helicase superfamily II, participates in diverse biological processes involved in RNA metabolism in organisms, and these RNA-mediated biological processes rely on RNA structure conversion. However, how RHA regulate the RNA structure conversion was still unknown. In order to unveil the mechanism of RNA structure conversion mediated by RHA, single molecule fluorescence resonance energy transfer was adopted to in our assay, and substrates RNA were from internal ribosome entry site of foot-and-mouth disease virus genome. We first found that the RNA structure conversion by RHA against thermodynamic equilibrium in vitro, and the process of dsRNA YZ converted to dsRNA XY through a tripartite intermediate state. In addition, the rate of the RNA structure conversion and the distribution of dsRNA YZ and XY were affected by ATP concentrations. Our study provides real-time insight into ATP-dependent RHA-assisted RNA structure conversion at the single molecule level, the mechanism displayed by RHA may help in understand how RHA contributes to many biological functions, and the basic mechanistic features illustrated in our work also underlay more complex protein-assisted RNA structure conversions.

Core Hairpin Structure of SpCas9 sgRNA Functions in a Sequence- and Spatial Conformation-Dependent Manner

The clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9 system is a widely used genome-editing tool with great clinical potential. However, its application is limited because of low editing efficiency of some target sequences and off-target effects. As this system contains only the Cas9 protein and a single-guide RNA (sgRNA; engineered from crRNA and tracrRNA), the structure and function of these components should be studied in detail to address the current clinical needs. Consequently, we investigated the structural and sequence features of the core hairpin (the first stem loop of sgRNA) of SpCas9 sgRNA. We showed that the core hairpin structure of sgRNA is essential for SpCas9/sgRNA-mediated DNA cleavage and that the internal loop structure in the core hairpin plays a vital role in target DNA cleavage. We observed that the root stem structure within the core hairpin preferentially forms Watson-Crick base pairs and should be of a specific length to maintain an appropriate spatial conformation for Cas9 binding. However, the length of the leaf stem structure of the core hairpin is flexible, having a variable nucleotide composition. Furthermore, extension of the leaf stem structure enhances the DNA cleavage activity of the Cas9/sgRNA complex, and this could be used to enhance the efficiency of gene editing. These observations provide insight into the sgRNA/Cas9 interaction, indicating that sgRNA modification could be a strategy for improved DNA editing efficiency, and optimized sgRNA can be further used for genome-wide functional screening and clinical application.

RNA-assisted self-assembly of monomeric antigens into virus-like particles as a recombinant vaccine platform

Representing highly ordered repetitive structures of antigen macromolecular assemblies, virus-like particles (VLPs) serve as a high-priority vaccine platform against emerging viral infections, as alternatives to traditional cell culture-based vaccines. RNAs can function as chaperones (Chaperna) and are extremely effective in promoting protein folding. Beyond their canonical function as translational adaptors, tRNAs may moonlight as chaperones for the kinetic control of macromolecular antigen assembly. Capitalizing on genomic RNA co-assembly in infectious virions, we present the first report of a biomimetic assembly of viral capsids that was assisted by non-viral host RNAs into genome-free, non-infectious empty particles. Here, we demonstrate the assembly of bacterially-produced soluble norovirus VP1 forming VLPs (n = 180) in vitro. A tRNA-interacting domain (tRID) was genetically fused with the VP1 capsid protein, as a tRNA docking tag, in the bacterial host to transduce chaperna function for de novo viral antigen folding. tRID/tRNA removal prompted the in vitro assembly of monomeric antigens into highly ordered repetitive structures that elicited robust protective immune responses after immunization. The chaperna-based assembly of monomeric antigens will impact the development and deployment of VLP vaccines for emerging and re-emerging viral infections.

Quality-Control Mechanism for Telomerase RNA Folding in the Cell

Long non-coding RNAs can often fold into different conformations. Telomerase RNA, an essential component of the telomerase ribonucleoprotein (RNP) enzyme, must fold into a defined structure to fulfill its function with the protein catalytic subunit (TERT) and other accessory factors. However, the mechanism by which the correct folding of telomerase RNA is warranted in a cell is still unknown. Here we show that La-related protein Pof8 specifically recognizes the conserved pseudoknot region of telomerase RNA and instructs the binding of the Lsm2-8 complex to its mature 3' end, thus selectively protecting the correctly folded RNA from exonucleolytic degradation. In the absence of Pof8, TERT assembles with misfolded RNA and produces little telomerase activity. Therefore, Pof8 plays a key role in telomerase RNA folding quality control, ensuring that TERT only assembles with functional telomerase RNA to form active telomerase. Our finding reveals a mechanism for non-coding RNA folding quality control.

Key Points to Consider When Studying RNA Remodeling by Proteins

Cellular RNAs depend on proteins for efficient folding to specific functional structures and for transitions between functional structures. This dependence arises from intrinsic properties of RNA structure. Specifically, RNAs possess stable local structure, largely in the form of helices, and there are abundant opportunities for RNAs to form alternative helices and tertiary contacts and therefore to populate alternative structures. Proteins with RNA chaperone activity, either ATP-dependent or ATP-independent, can promote structural transitions by interacting with single-stranded RNA (ssRNA) to compete away partner interactions and then release ssRNA so that it can form new interactions. In this chapter we review the basic properties of RNA and the proteins that function as chaperones and remodelers. We then use these properties as a foundation to explore key points for the design and interpretation of experiments that probe RNA rearrangements and their acceleration by proteins.

Probing Transcriptome-Wide RNA Structural Changes Dependent on the DEAD-box Helicase Dbp2

RNA helicases function in all aspects of RNA biology mainly through remodeling structures of RNA and RNA-protein (RNP) complexes. Among them, DEAD-box proteins form the largest family in eukaryotes and have been shown to remodel RNA/RNP structures and clamping of RNA-binding proteins, both in vitro and in vivo. Nevertheless, for the majority of these enzymes, it is largely unclear what RNAs are targeted and where they modulate RNA/RNP structures to promote RNA metabolism. Several methods have been developed to probe secondary and tertiary structures of specific transcripts or whole transcriptomes in vivo. In this chapter, we describe a protocol for identification of RNA structural changes that are dependent on a Saccharomyces cerevisiae DEAD-box helicase Dbp2. Experiments detailed here can be adapted to the study of other RNA helicases and identification of putative remodeling targets in vivo.

José M. The Ancient History of Peptidyl Transferase Center Formation as Told by Conservation and Information Analyses

The peptidyl transferase center (PTC) is the catalytic center of the ribosome and forms part of the 23S ribosomal RNA. The PTC has been recognized as the earliest ribosomal part and its origins embodied the First Universal Common Ancestor (FUCA). The PTC is frequently assumed to be highly conserved along all living beings. In this work, we posed the following questions: (i) How many 100% conserved bases can be found in the PTC? (ii) Is it possible to identify clusters of informationally linked nucleotides along its sequence? (iii) Can we propose how the PTC was formed? (iv) How does sequence conservation reflect on the secondary and tertiary structures of the PTC? Aiming to answer these questions, all available complete sequences of 23S ribosomal RNA from Bacteria and Archaea deposited on GenBank database were downloaded. Using a sequence bait of 179 bp from the PTC of Thermus termophilus, we performed an optimum pairwise alignment to retrieve the PTC region from 1424 filtered 23S rRNA sequences. These PTC sequences were multiply aligned, and the conserved regions were assigned and observed along the primary, secondary, and tertiary structures. The PTC structure was observed to be more highly conserved close to the adenine located at the catalytical site. Clusters of interrelated, co-evolving nucleotides reinforce previous assumptions that the PTC was formed by the concatenation of proto-tRNAs and important residues responsible for its assembly were identified. The observed sequence variation does not seem to significantly affect the 3D structure of the PTC ribozyme.

Stress proteins: the biological functions in virus infection, present and challenges for target-based antiviral drug development

Stress proteins (SPs) including heat-shock proteins (HSPs), RNA chaperones, and ER associated stress proteins are molecular chaperones essential for cellular homeostasis. The major functions of HSPs include chaperoning misfolded or unfolded polypeptides, protecting cells from toxic stress, and presenting immune and inflammatory cytokines. Regarded as a double-edged sword, HSPs also cooperate with numerous viruses and cancer cells to promote their survival. RNA chaperones are a group of heterogeneous nuclear ribonucleoproteins (hnRNPs), which are essential factors for manipulating both the functions and metabolisms of pre-mRNAs/hnRNAs transcribed by RNA polymerase II. hnRNPs involve in a large number of cellular processes, including chromatin remodelling, transcription regulation, RNP assembly and stabilization, RNA export, virus replication, histone-like nucleoid structuring, and even intracellular immunity. Dysregulation of stress proteins is associated with many human diseases including human cancer, cardiovascular diseases, neurodegenerative diseases (e.g., Parkinson’s diseases, Alzheimer disease), stroke and infectious diseases. In this review, we summarized the biologic function of stress proteins, and current progress on their mechanisms related to virus reproduction and diseases caused by virus infections. As SPs also attract a great interest as potential antiviral targets (e.g., COVID-19), we also discuss the present progress and challenges in this area of HSP-based drug development, as well as with compounds already under clinical evaluation.

Phosphorylation of the Arginine-Rich C-Terminal Domains of the Hepatitis B Virus (HBV) Core Protein as a Fine Regulator of the Interaction between HBc and Nucleic Acid

The morphogenesis of Hepatitis B Virus (HBV) viral particles is nucleated by the oligomerization of HBc protein molecules, resulting in the formation of an icosahedral capsid shell containing the replication-competent nucleoprotein complex made of the viral polymerase and the pre-genomic RNA (pgRNA). HBc is a phospho-protein containing two distinct domains acting together throughout the viral replication cycle. The N-terminal domain, (residues 1–140), shown to self-assemble, is linked by a short flexible domain to the basic C-terminal domain (residues 150–183) that interacts with nucleic acids (NAs). In addition, the C-terminal domain contains a series of phospho-acceptor residues that undergo partial phosphorylation and de-phosphorylation during virus replication. This highly dynamic process governs the homeostatic charge that is essential for capsid stability, pgRNA packaging and to expose the C-terminal domain at the surface of the particles for cell trafficking. In this review, we discuss the roles of the N-terminal and C-terminal domains of HBc protein during HBV morphogenesis, focusing on how the C-terminal domain phosphorylation dynamics regulate its interaction with nucleic acids throughout the assembly and maturation of HBV particles.

Ribosomal RNA Methyltransferase RsmC Moonlights as an RNA Chaperone

Ribosomes are ribonucleoprotein particles that are essential for protein biosynthesis in all forms of life. During ribosome biogenesis, transcription, folding, modification, and processing of rRNA are coupled to the assembly of proteins. Various assembly factors are required to synchronize all different processes that occur during ribosome biogenesis. Herein, the RNA chaperone and RNA strand annealing activity of rRNA modification enzyme ribosome small subunit methyltransferase C (RsmC), which modifies guanine to 2-methylguanosine (m² G) at position 1207 of 16S rRNA (Escherichia coli nucleotide numbering) located at helix 34 (h34), are reported. A 25-fold increase in the h34 RNA strand annealing rates is observed in the presence of RsmC. Single-molecule FRET experiments confirmed the ability of protein RsmC to denature a non-native structure formed by one of the two h34 strands and to form a native-like duplex. This observed RNA chaperone activity of protein RsmC might play a vital role in the rapid generation of functional ribosomes.

Against Expectations: Unassisted RNA Adsorption onto Negatively Charged Lipid Bilayers

The composition and physicochemical properties of biological membranes can be altered by diverse membrane integral and peripheral proteins as well as by small molecules, natural and synthetic. Diverse oligonucleotides have been shown to electrostatically interact with cationic and bivalent ion loaded zwitterionic liposomes, leading to the formation of oligonucleotide-liposome aggregates. However, interaction of RNAs with other membrane surfaces remains ill understood. We used the nonnatural RNA10 to investigate RNA binding to anionic and net-uncharged membrane surfaces. RNA10 had initially been selected in a screen for nonnatural RNA motives that bind to phosphatidylcholine liposomes in the presence of Mg²⁺. Here we show that interaction of defined RNA molecules with membrane surfaces crucially depends on electrostatic surface properties. Furthermore, RNA10 electrostatically binds to anionic lipid bilayers in the absence of Mg²⁺ or other bivalent cations, and this interaction leads to measurably changed physicochemical properties of the bilayer and the oligonucleotide. Thus, the structure of polyanionic RNA can be modulated via contact with negatively charged membrane surfaces and vice versa.

A Modified In Vitro Transcription Approach to Improve RNA Synthesis and Ribozyme Cleavage Efficiency

RNA elements such as catalytic RNA, riboswitch, microRNA, and long non-coding RNA perform a major role in cellular processes. A complete understanding of cellular processes is impossible without knowing the structure-function relationship of participating RNA molecules that ultimately requires large quantities of pure RNAs. Thus, structural/functional analyses of emerging RNAs necessitate revised protocols for improved RNA quantity and quality. Here we present a modified in vitro transcription protocol to enhance ribozyme cleaving efficiency and RNA yield by working on two variables, i.e., incubation temperature and limiting GTPs. Following an improved RNA synthesis, the target RNA is purified from transcription mixture components through denaturing size-exclusion chromatography. The protocol confirms that cyclic elevated incubation temperatures during transcription and increased concentrations of GTPs improve the production rate of RNA. Our modified in vitro transcription method improves the ribozyme cleaving efficiency and targets RNA yield by four- to fivefold that can benefit almost any RNA-related study from protein-RNA interaction analysis to crystallography.

About the Crystallization of Abiotic Coded Matter

The crystallization of digitally encoded polyurethanes was studied by electron diffraction. A series of oligomers with different primary structures was analyzed in this work. They all form hydrogen-bonding-directed lamellar single crystals with a base-centered orthorhombic unit cell. Although crystal morphology was the same in all cases, the digital coding of the oligomers has a small influence on the intersheet distance in the crystals. The crystal lattices allow calculation of the volume occupied by one basic information unit, which is in the range 148-188 ų. Interestingly, this volume is about 3× smaller than that occupied by a coded nucleotide in a DNA double helix. Furthermore, crystallization of blends of oligourethanes with different coded primary structures was investigated. Oligomers with drastically different monomer compositions form structures that are not cocrystals but more probably segregated crystals containing distinct domains of different composition.

Computational prediction of functions of intrinsically disordered regions

Intrinsically disorder regions (IDRs) are abundant in nature, particularly among Eukaryotes. While they facilitate a wide spectrum of cellular functions including signaling, molecular assembly and recognition, translation, transcription and regulation, only several hundred IDRs are annotated functionally. This annotation gap motivates the development of fast and accurate computational methods that predict IDR functions directly from protein sequences. We introduce and describe a comprehensive collection of 25 methods that provide accurate predictions of IDRs that interact with proteins and nucleic acids, that function as flexible linkers and that moonlight multiple functions. Virtually all of these predictors can be accessed online and many were developed in the last few years. They utilize a wide range of predictive architectures and take advantage of modern machine learning algorithms. Our empirical analysis shows that predictors that are available as webservers enjoy high rates of citations, attesting to their practical value and popularity. The most cited methods include DISOPRED3, ANCHOR, alpha-MoRFpred, MoRFpred, fMoRFpred and MoRFCHiBi. We present two case studies to demonstrate that predictions produced by these computational tools are relatively easy to interpret and that they deliver valuable functional clues. However, the current computational tools cover a relatively narrow range of disorder functions. Further development efforts that would cover a broader range of functions should be pursued. We demonstrate that a sufficient amount of functionally annotated IDRs that are associated with several other disorder functions is already available and can be used to design and validate novel predictors.

Proteins That Chaperone RNA Regulation

RNA-binding proteins chaperone the biological functions of noncoding RNA by reducing RNA misfolding, improving matchmaking between regulatory RNA and targets, and exerting quality control over RNP biogenesis. Recent studies of Escherichia coli CspA, HIV NCp, and E. coli Hfq are beginning to show how RNA-binding proteins remodel RNA structures. These different protein families use common strategies for disrupting or annealing RNA double helices, which can be used to understand the mechanisms by which proteins chaperone RNA-dependent regulation in bacteria.

Discovery of a Branched Peptide That Recognizes the Rev Response Element (RRE) RNA and Blocks HIV-1 Replication

We synthesized and screened a unique 46 656-member library composed of unnatural amino acids that revealed several hits against RRE IIB RNA. Among the hit peptides identified, peptide 4A5 was found to be selective against competitor RNAs and inhibited HIV-1 Rev-RRE RNA interaction in cell culture in a p24 ELISA assay. Biophysical characterization in a ribonuclease protection assay suggested that 4A5 bound to the stem-loop region in RRE IIB while SHAPE MaP probing with 234 nt RRE RNA indicated additional interaction with secondary Rev binding sites. Taken together, our investigation suggests that HIV replication is inhibited by 4A5 blocking binding of Rev and subsequent multimerization.

How the Conformations of an Internal Junction Contribute to Fold an RNA Domain

Like proteins, some RNAs fold to compact structures. We can model functional RNAs as a series of short, rigid, base-paired elements, connected by non-base-paired nucleotides that serve as junctions. These connecting regions bend and twist, facilitating the formation of tertiary contacts that stabilize compact states. Here, we explore the roles of salt and junction sequence in determining the structures of a ubiquitous connector: an asymmetric internal loop. We focus on the J5/5a junction from the widely studied P4-P6 domain of the Tetrahymena ribozyme. Following the addition of magnesium ions to fold P4-P6, this junction bends dramatically, bringing the two halves of the RNA domain together for tertiary contact engagement. Using single-molecule fluorescence resonance energy transfer (smFRET), we examine the role of sequence and salt on model RNA constructs that contain these junction regions. We explore the wild-type J5/5a junction as well as two sequence variants. These junctions display distinct, salt-dependent conformations. Small-angle X-ray scattering (SAXS) measurements verify that these effects persist in the full-length P4-P6 domain. These measurements underscore the importance of junction sequence and interactions with ions in facilitating RNA folding.

RNA as a Complex Polymer with Coupled Dynamics of Ions and Water in the Outer Solvation Sphere

We unravel the internal and collective modes of a widely studied 58-nucleotide rRNA fragment in solvent using atomically detailed molecular dynamics simulations. The variation of lifetimes for water hydrogen bonds with nucleotide groups indicates heterogeneity of water dynamics on the RNA surface. The time scales of interactions of the discrete water molecules with RNA nucleotides extend from several hundred picoseconds to a few nanoseconds. We determine all of the association sites and the spatial distribution of residence times for Mg²⁺, K⁺, and water molecules in those sites. We provide insights into the population of Mg²⁺ and K⁺ ions and water molecules in the outer sphere and how their fluctuations are intricately linked with the kinetics of the 58-mer. We find that many of the long-lived Mg²⁺ sites identified from the simulations agree with the locations of ions in the X-ray structure. We determine the excess ion atmosphere around the rRNA and compare it with experimental data. We investigate the collective behavior of RNA, ions, and water, by performing a joint principle component analysis for the Cartesian coordinates of the RNA phosphorus atoms and for the occupation counts of the association sites. Our results indicate that the 58-mer system is a complex polymer, composed of RNA that is encased by a fluctuating network of associated counterions, co-ions, and water.

Towards Long-Range RNA Structure Prediction in Eukaryotic Genes

The ability to form an intramolecular structure plays a fundamental role in eukaryotic RNA biogenesis. Proximate regions in the primary transcripts fold into a local secondary structure, which is then hierarchically assembled into a tertiary structure that is stabilized by RNA-binding proteins and long-range intramolecular base pairings. While the local RNA structure can be predicted reasonably well for short sequences, long-range structure at the scale of eukaryotic genes remains problematic from the computational standpoint. The aim of this review is to list functional examples of long-range RNA structures, to summarize current comparative methods of structure prediction, and to highlight their advances and limitations in the context of long-range RNA structures. Most comparative methods implement the “first-align-then-fold” principle, i.e., they operate on multiple sequence alignments, while functional RNA structures often reside in non-conserved parts of the primary transcripts. The opposite “first-fold-then-align” approach is currently explored to a much lesser extent. Developing novel methods in both directions will improve the performance of comparative RNA structure analysis and help discover novel long-range structures, their higher-order organization, and RNA⁻RNA interactions across the transcriptome.

An Activity Switch in Human Telomerase Based on RNA Conformation and Shaped by TCAB1

Ribonucleoprotein enzymes require dynamic conformations of their RNA constituents for regulated catalysis. Human telomerase employs a non-coding RNA (hTR) with a bipartite arrangement of domains-a template-containing core and a distal three-way junction (CR4/5) that stimulates catalysis through unknown means. Here, we show that telomerase activity unexpectedly depends upon the holoenzyme protein TCAB1, which in turn controls conformation of CR4/5. Cells lacking TCAB1 exhibit a marked reduction in telomerase catalysis without affecting enzyme assembly. Instead, TCAB1 inactivation causes unfolding of CR4/5 helices that are required for catalysis and for association with the telomerase reverse-transcriptase (TERT). CR4/5 mutations derived from patients with telomere biology disorders provoke defects in catalysis and TERT binding similar to TCAB1 inactivation. These findings reveal a conformational "activity switch" in human telomerase RNA controlling catalysis and TERT engagement. The identification of two discrete catalytic states for telomerase suggests an intramolecular means for controlling telomerase in cancers and progenitor cells.

Quantitative Characterization of Translational Riboregulators Using an in Vitro Transcription-Translation System

Riboregulators are short RNA sequences that, upon binding to a ligand, change their secondary structure and influence the expression rate of a downstream gene. They constitute an attractive alternative to transcription factors for building synthetic gene regulatory networks because they can be engineered de novo. However, riboregulators are generally designed in silico and tested in vivo, which provides little quantitative information about their performances, thus hindering the improvement of design algorithms. Here we show that a cell-free transcription-translation (TX-TL) system provides valuable information about the performances of in silico designed riboregulators. We first propose a simple model that provides a quantitative definition of the dynamic range of a riboregulator. We further characterize two types of translational riboregulators composed of a cis-repressed (cr) and a trans-activating (ta) strand. At the DNA level we demonstrate that high concentrations of taDNA poisoned the activator until total shut off, in agreement with our model, and that relative dynamic ranges of riboregulators determined in vitro are in agreement with published in vivo data. At the RNA level, we show that this approach provides a fast and simple way to measure dissociation constants of functional riboregulators, in contrast to standard mobility-shift assays. Our method opens the route for using cell-free TX-TL systems for the quantitative characterization of functional riboregulators in order to improve their design in silico.

Allele-specific SHAPE-MaP assessment of the effects of somatic variation and protein binding on mRNA structure

The impact of inherited and somatic mutations on messenger RNA (mRNA) structure remains poorly understood. Recent technological advances that leverage next-generation sequencing to obtain experimental structure data, such as SHAPE-MaP, can reveal structural effects of mutations, especially when these data are incorporated into structure modeling. Here, we analyze the ability of SHAPE-MaP to detect the relatively subtle structural changes caused by single-nucleotide mutations. We find that allele-specific sorting greatly improved our detection ability. Thus, we used SHAPE-MaP with a novel combination of clone-free robotic mutagenesis and allele-specific sorting to perform a rapid, comprehensive survey of noncoding somatic and inherited riboSNitches in two cancer-associated mRNAs, TPT1 and LCP1 Using rigorous thermodynamic modeling of the Boltzmann suboptimal ensemble, we identified a subset of mutations that change TPT1 and LCP1 RNA structure, with approximately 14% of all variants identified as riboSNitches. To confirm that these in vitro structures were biologically relevant, we tested how dependent TPT1 and LCP1 mRNA structures were on their environments. We performed SHAPE-MaP on TPT1 and LCP1 mRNAs in the presence or absence of cellular proteins and found that both mRNAs have similar overall folds in all conditions. RiboSNitches identified within these mRNAs in vitro likely exist under biological conditions. Overall, these data reveal a robust mRNA structural landscape where differences in environmental conditions and most sequence variants do not significantly alter RNA structural ensembles. Finally, predicting riboSNitches in mRNAs from sequence alone remains particularly challenging; these data will provide the community with benchmarks for further algorithmic development.

Structural Changes of RNA in Complex with Proteins in the SRP

The structural flexibility of RNA allows it to exist in several shapes and sizes. Thus, RNA is functionally diverse and is known to be involved in processes such as catalysis, ligand binding, and most importantly, protein recognition. RNA can adopt different structures, which can often dictate its functionality. When RNA binds onto protein to form a ribonucleoprotein complex (RNP), multiple interactions and conformational changes occur with the RNA and protein. However, there is the question of whether there is a specific pattern for these changes to occur upon recognition. In particular when RNP complexity increases with the addition of multiple proteins/RNA, it becomes difficult to structurally characterize the overall changes using the current structural determination techniques. Hence, there is a need to use a combination of biochemical, structural and computational modeling to achieve a better understanding of the processes that RNPs are involved. Nevertheless, there are well-characterized systems that are evolutionarily conserved [such as the signal recognition particle (SRP)] that give us important information on the structural changes of RNA and protein upon complex formation.

Single-Molecule Fluorescence Reveals Commonalities and Distinctions among Natural and in Vitro-Selected RNA Tertiary Motifs in a Multistep Folding Pathway

Decades of study of the RNA folding problem have revealed that diverse and complex structured RNAs are built from a common set of recurring structural motifs, leading to the perspective that a generalizable model of RNA folding may be developed from understanding of the folding properties of individual structural motifs. We used single-molecule fluorescence to dissect the kinetic and thermodynamic properties of a set of variants of a common tertiary structural motif, the tetraloop/tetraloop-receptor (TL/TLR). Our results revealed a multistep TL/TLR folding pathway in which preorganization of the ubiquitous AA-platform submotif precedes the formation of the docking transition state and tertiary A-minor hydrogen bond interactions form after the docking transition state. Differences in ion dependences between TL/TLR variants indicated the occurrence of sequence-dependent conformational rearrangements prior to and after the formation of the docking transition state. Nevertheless, varying the junction connecting the TL/TLR produced a common kinetic and ionic effect for all variants, suggesting that the global conformational search and compaction electrostatics are energetically independent from the formation of the tertiary motif contacts. We also found that in vitro-selected variants, despite their similar stability at high Mg²⁺ concentrations, are considerably less stable than natural variants under near-physiological ionic conditions, and the occurrence of the TL/TLR sequence variants in Nature correlates with their thermodynamic stability in isolation. Overall, our findings are consistent with modular but complex energetic properties of RNA structural motifs and will aid in the eventual quantitative description of RNA folding from its secondary and tertiary structural elements.

RNA Chaperone Function of a Universal Stress Protein in Arabidopsis Confers Enhanced Cold Stress Tolerance in Plants

The physiological function of Arabidopsis thaliana universal stress protein (AtUSP) in plant has remained unclear. Thus, we report here the functional role of the Arabidopsis universal stress protein, AtUSP (At3g53990). To determine how AtUSP affects physiological responses towards cold stress, AtUSP overexpression (AtUSP OE) and T-DNA insertion knock-out (atusp, SALK_146059) mutant lines were used. The results indicated that AtUSP OE enhanced plant tolerance to cold stress, whereas atusp did not. AtUSP is localized in the nucleus and cytoplasm, and cold stress significantly affects RNA metabolism such as by misfolding and secondary structure changes of RNA. Therefore, we investigated the relationship of AtUSP with RNA metabolism. We found that AtUSP can bind nucleic acids, including single- and double-stranded DNA and luciferase mRNA. AtUSP also displayed strong nucleic acid-melting activity. We expressed AtUSP in RL211 Escherichia coli, which contains a hairpin-loop RNA structure upstream of chloramphenicol acetyltransferase (CAT), and observed that AtUSP exhibited anti-termination activity that enabled CAT gene expression. AtUSP expression in the cold-sensitive Escherichia coli (E. coli) mutant BX04 complemented the cold sensitivity of the mutant cells. As these properties are typical characteristics of RNA chaperones, we conclude that AtUSP functions as a RNA chaperone under cold-shock conditions. Thus, the enhanced tolerance of AtUSP OE lines to cold stress is mediated by the RNA chaperone function of AtUSP.

RNA Regulations and Functions Decoded by Transcriptome-wide RNA Structure Probing

RNA folds into intricate structures that are crucial for its functions and regulations. To date, a multitude of approaches for probing structures of the whole transcriptome, i.e., RNA structuromes, have been developed. Applications of these approaches to different cell lines and tissues have generated a rich resource for the study of RNA structure–function relationships at a systems biology level. In this review, we first introduce the designs of these methods and their applications to study different RNA structuromes. We emphasize their technological differences especially their unique advantages and caveats. We then summarize the structural insights in RNA functions and regulations obtained from the studies of RNA structuromes. And finally, we propose potential directions for future improvements and studies.

Double Hydrogen Bonding between Side Chain Carboxyl Groups in Aqueous Solutions of Poly (β-L-Malic Acid): Implication for the Evolutionary Origin of Nucleic Acids

The RNA world hypothesis holds that in the evolutionary events that led to the emergence of life RNA preceded proteins and DNA and is supported by the ability of RNA to act as both a genetic polymer and a catalyst. On the other hand, biosynthesis of nucleic acids requires a large number of enzymes and chemical synthesis of RNA under presumed prebiotic conditions is complicated and requires many sequential steps. These observations suggest that biosynthesis of RNA is the end product of a long evolutionary process. If so, what was the original polymer from which RNA and DNA evolved? In most syntheses of simpler RNA or DNA analogs, the D-ribose phosphate polymer backbone is altered and the purine and pyrimidine bases are retained for hydrogen bonding between complementary base pairs. However, the bases are themselves products of complex biosynthetic pathways and hence they too may have evolved from simpler polymer side chains that had the ability to form hydrogen bonds. We hypothesize that the earliest evolutionary predecessor of nucleic acids was the simple linear polyester, poly (β-D-malic acid), for which the carboxyl side chains could form double hydrogen bonds. In this study, we show that in accord with this hypothesis a closely related polyester, poly (β-L-malic acid), uses carboxyl side chains to form robust intramolecular double hydrogen bonds in moderately acidic solution.

Bioinformatics comparisons of RNA-binding proteins of pathogenic and non-pathogenic Escherichia coli strains reveal novel virulence factors

Background

Pathogenic bacteria have evolved various strategies to counteract host defences. They are also exposed to environments that are undergoing constant changes. Hence, in order to survive, bacteria must adapt themselves to the changing environmental conditions by performing regulations at the transcriptional and/or post-transcriptional levels. Roles of RNA-binding proteins (RBPs) as virulence factors have been very well studied. Here, we have used a sequence search-based method to compare and contrast the proteomes of 16 pathogenic and three non-pathogenic E. coli strains as well as to obtain a global picture of the RBP landscape (RBPome) in E. coli.

Results

Our results show that there are no significant differences in the percentage of RBPs encoded by the pathogenic and the non-pathogenic E. coli strains. The differences in the types of Pfam domains as well as Pfam RNA-binding domains, encoded by these two classes of E. coli strains, are also insignificant. The complete and distinct RBPome of E. coli has been established by studying all known E. coli strains till date. We have also identified RBPs that are exclusive to pathogenic strains, and most of them can be exploited as drug targets since they appear to be non-homologous to their human host proteins. Many of these pathogen-specific proteins were uncharacterised and their identities could be resolved on the basis of sequence homology searches with known proteins. Detailed structural modelling, molecular dynamics simulations and sequence comparisons have been pursued for selected examples to understand differences in stability and RNA-binding.

Conclusions

The approach used in this paper to cross-compare proteomes of pathogenic and non-pathogenic strains may also be extended to other bacterial or even eukaryotic proteomes to understand interesting differences in their RBPomes. The pathogen-specific RBPs reported in this study, may also be taken up further for clinical trials and/or experimental validations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4045-3) contains supplementary material, which is available to authorized users.

Coronavirus nucleocapsid proteins assemble constitutively in high molecular oligomers

Coronaviruses (CoV) are enveloped viruses and rely on their nucleocapsid N protein to incorporate the positive-stranded genomic RNA into the virions. CoV N proteins form oligomers but the mechanism and relevance underlying their multimerization remain to be fully understood. Using in vitro pull-down experiments and density glycerol gradients, we found that at least 3 regions distributed over its entire length mediate the self-interaction of mouse hepatitis virus (MHV) and severe acute respiratory syndrome coronavirus (SARS-CoV) N protein. The fact that these regions can bind reciprocally between themselves provides a possible molecular basis for N protein oligomerization. Interestingly, cytoplasmic N molecules of MHV-infected cells constitutively assemble into oligomers through a process that does not require binding to genomic RNA. Based on our data, we propose a model where constitutive N protein oligomerization allows the optimal loading of the genomic viral RNA into a ribonucleoprotein complex via the presentation of multiple viral RNA binding motifs.

The multiple roles of the nucleocapsid in retroviral RNA conversion into proviral DNA by reverse transcriptase

Retroviruses are enveloped plus-strand RNA viruses that can cause cancer, immunodeficiency and neurological disorder in human and animals. Retroviruses have several unique properties, such as a genomic RNA in a dimeric form found in the virus, and a replication strategy called 'copy-and-paste' during which the plus-strand genomic RNA is converted into a double-stranded DNA, subsequently integrated into the cellular genome. Two essential viral enzymes, reverse transcriptase (RT) and integrase (IN), direct this 'copy-and-paste' replication. RT copies the genomic RNA generating the double-stranded proviral DNA, while IN catalyzes proviral DNA integration into the cellular DNA, then called the provirus. In that context, a major component of the virion core, the nucleocapsid protein (NC), was found to be a potent nucleic-acid chaperone that assists RT during the conversion of the genomic RNA into proviral DNA. Here we briefly review the interplay of NC with viral nucleic-acids, which enables rapid and faithful folding and hybridization of complementary sequences, and with active RT thus providing assistance to the synthesis of the complete proviral DNA. Because of its multiple roles in retrovirus replication, NC could be viewed as a two-faced Janus-chaperone acting on viral nucleic-acids and enzymes.

RNA splicing in human disease and in the clinic

Defects at the level of the pre-mRNA splicing process represent a major cause of human disease. Approximately 15-50% of all human disease mutations have been shown to alter functioning of basic and auxiliary splicing elements. These elements are required to ensure proper processing of pre-mRNA splicing molecules, with their disruption leading to misprocessing of the pre-mRNA molecule and disease. The splicing process is a complex process, with much still to be uncovered before we are able to accurately predict whether a reported genomic sequence variant (GV) represents a splicing-associated disease mutation or a harmless polymorphism. Furthermore, even when a mutation is correctly identified as affecting the splicing process, there still remains the difficulty of providing an exact evaluation of the potential impact on disease onset, severity and duration. In this review, we provide a brief overview of splicing diagnostic methodologies, from in silico bioinformatics approaches to wet lab in vitro and in vivo systems to evaluate splicing efficiencies. In particular, we provide an overview of how the latest developments in high-throughput sequencing can be applied to the clinic, and are already changing clinical approaches.

Collagen Membranes with Ribonuclease Inhibitors for Long-Term Stability of Electrochemical Aptamer-Based Sensors Employing RNA

Electrochemical aptamer-based (E-AB) sensors offer advantageous analytical detection abilities due to their rapid response time (seconds to minutes), specificity to a target, and selectivity to function in complex media. Ribonucleic acid (RNA) aptamers employed in this class of sensor offer favorable binding characteristics resulting from the ability of RNA to form stable tertiary folds aided by long-range intermolecular interactions. As a result, RNA aptamers can fold into three-dimensional structures more complex than those of their DNA counterparts and consequently exhibit better binding ability to target analytes. Unfortunately, RNA aptamers are susceptible to degradation by nucleases, and for this reason, RNA-based sensors are scarce or require significant sample pretreatment before use in clinically relevant media. Here, we combine the usefulness of a collagen I hydrogel membrane with entrapped ribonuclease inhibitors (RI) to protect small molecule RNA E-AB sensors from endogenous nucleases in complex media. More specifically, the biocompatibility of the naturally polymerized hydrogel with encapsulated RI promotes the protection of an aminoglycoside-binding RNA E-AB sensor up to 6 h, enabling full sensor function in nuclease-rich environments (undiluted serum) without the need for prior sample preparation or oligonucleotide modification. The use of collagen as a biocompatible membrane represents a general approach to compatibly interface E-AB sensors with complex biological samples.

Capturing RNA Folding Free Energy with Coarse-Grained Molecular Dynamics Simulations

We introduce a coarse-grained RNA model for molecular dynamics simulations, RACER (RnA CoarsE-gRained). RACER achieves accurate native structure prediction for a number of RNAs (average RMSD of 2.93 Å) and the sequence-specific variation of free energy is in excellent agreement with experimentally measured stabilities (R2 = 0.93). Using RACER, we identified hydrogen-bonding (or base pairing), base stacking, and electrostatic interactions as essential driving forces for RNA folding. Also, we found that separating pairing vs. stacking interactions allowed RACER to distinguish folded vs. unfolded states. In RACER, base pairing and stacking interactions each provide an approximate stability of 3-4 kcal/mol for an A-form helix. RACER was developed based on PDB structural statistics and experimental thermodynamic data. In contrast with previous work, RACER implements a novel effective vdW potential energy function, which led us to re-parameterize hydrogen bond and electrostatic potential energy functions. Further, RACER is validated and optimized using a simulated annealing protocol to generate potential energy vs. RMSD landscapes. Finally, RACER is tested using extensive equilibrium pulling simulations (0.86 ms total) on eleven RNA sequences (hairpins and duplexes).

Defective Ribonucleoproteins, Mistakes in RNA Processing, and Diseases

Ribonucleoproteins (RNPs) are vital to many cellular events. To this end, many neurodegenerative diseases and cancers have been linked to RNP malfunction, particularly as this relates to defective processing of cellular RNA. The connection of RNPs and diseases has also propagated a shift of focus onto RNA targeting from traditional protein targeting treatments. However, therapeutic development in this area has been limited by incomplete molecular insight into the specific contributions of RNPs to disease. This review outlines the role of several RNPs in diseases, focusing on molecular defects in processes that affect proper RNA handling in the cell. This work also evaluates the contributions of recently developed methods to understanding RNP association and function. We review progress in this area by focusing on molecular malfunctions of RNPs associated with the onset and progression of several neurodegenerative diseases and cancer and conclude with a brief discussion of RNA-based therapeutic efforts.