Exposing the kinetic traps in RNA folding

A Deep Learning-Driven Sampling Technique to Explore the Phase Space of an RNA Stem-Loop

The folding and unfolding of RNA stem-loops are critical biological processes; however, their computational studies are often hampered by the ruggedness of their folding landscape, necessitating long simulation times at the atomistic scale. Here, we adapted DeepDriveMD (DDMD), an advanced deep learning-driven sampling technique originally developed for protein folding, to address the challenges of RNA stem-loop folding. Although tempering- and order parameter-based techniques are commonly used for similar rare-event problems, the computational costs or the need for a priori knowledge about the system often present a challenge in their effective use. DDMD overcomes these challenges by adaptively learning from an ensemble of running MD simulations using generic contact maps as the raw input. DeepDriveMD enables on-the-fly learning of a low-dimensional latent representation and guides the simulation toward the undersampled regions while optimizing the resources to explore the relevant parts of the phase space. We showed that DDMD estimates the free energy landscape of the RNA stem-loop reasonably well at room temperature. Our simulation framework runs at a constant temperature without external biasing potential, hence preserving the information on transition rates, with a computational cost much lower than that of the simulations performed with external biasing potentials. We also introduced a reweighting strategy for obtaining unbiased free energy surfaces and presented a qualitative analysis of the latent space. This analysis showed that the latent space captures the relevant slow degrees of freedom for the RNA folding problem of interest. Finally, throughout the manuscript, we outlined how different parameters are selected and optimized to adapt DDMD for this system. We believe this compendium of decision-making processes will help new users adapt this technique for the rare-event sampling problems of their interest.

Identical sequences, different behaviors: Protein diversity captured at the single-molecule level

The classical "one sequence, one structure, one function" paradigm has shaped much of our intuition of how proteins work inside the cell. Partially due to the insight provided by bulk biochemical assays, individual biomolecules are often assumed to behave as identical entities, and their characterization relies on ensemble averages that flatten any conformational diversity into a unique phenotype. While the emergence of single-molecule techniques opened the gates to interrogating individual molecules, technical shortcomings typically limit the duration of these measurements, which precludes a complete characterization of an individual protein and, hence, capturing the heterogeneity among molecular populations. Here, we introduce an ultrastable magnetic tweezers design, which enables us to measure the folding dynamics of a single protein during several uninterrupted days with high temporal and spatial resolution. Thanks to this instrumental development, we fully characterize the nanomechanics of two proteins with a very distinct force response, the talin R3IVVI domain and protein L. Days-long recordings on the same protein individual accumulate thousands of folding transitions with submicrosecond resolution, allowing us to reconstruct their free energy landscapes and describe how they evolve with force. By mapping the nanomechanical identity of many different protein individuals, we directly capture their molecular diversity as a quantifiable dispersion on their force response and folding kinetics. By significantly expanding the measurable timescales, our instrumental development offers a tool for profiling individual molecules, opening the gates to directly characterizing biomolecular heterogeneity.

Charge transport in individual short base stacked single-stranded RNA molecules

Charge transport in biomolecules is crucial for many biological and technological applications, including biomolecular electronics devices and biosensors. RNA has become the focus of research because of its importance in biomedicine, but its charge transport properties are not well understood. Here, we use the Scanning Tunneling Microscopy-assisted molecular break junction method to measure the electrical conductance of particular 5-base and 10-base single-stranded (ss) RNA sequences capable of base stacking. These ssRNA sequences show single-molecule conductance values around [Formula: see text] ([Formula: see text]), while equivalent-length ssDNAs result in featureless conductance histograms. Circular dichroism (CD) spectra and MD simulations reveal the existence of extended ssRNA conformations versus folded ssDNA conformations, consistent with their different electrical behaviors. Computational molecular modeling and Machine Learning-assisted interpretation of CD data helped us to disentangle the structural and electronic factors underlying CT, thus explaining the observed electrical behavior differences. RNA with a measurable conductance corresponds to sequences with overall extended base-stacking stabilized conformations characterized by lower HOMO energy levels delocalized over a base-stacking mediating CT pathway. In contrast, DNA and a control RNA sequence without significant base-stacking tend to form closed structures and thus are incapable of efficient CT.

Comparative analysis of RNA secondary structure accuracy on predicted RNA 3D models

RNA structure is conformationally dynamic, and accurate all-atom tertiary (3D) structure modeling of RNA remains challenging with the prevailing tools. Secondary structure (2D) information is the standard prerequisite for most RNA 3D modeling. Despite several 2D and 3D structure prediction tools proposed in recent years, one of the challenges is to choose the best combination for accurate RNA 3D structure prediction. Here, we benchmarked seven small RNA PDB structures (40 to 90 nucleotides) with different topologies to understand the effects of different 2D structure predictions on the accuracy of 3D modeling. The current study explores the blind challenge of 2D to 3D conversions and highlights the performances of de novo RNA 3D modeling from their predicted 2D structure constraints. Our results show that conformational sampling-based methods such as SimRNA and IsRNA1 depend less on 2D accuracy, whereas motif-based methods account for 2D evidence. Our observations illustrate the disparities in available 3D and 2D prediction methods and may further offer insights into developing topology-specific or family-specific RNA structure prediction pipelines.

Odd-even disparity in the population of slipped hairpins in RNA repeat sequences with implications for phase separation

Low-complexity nucleotide repeat sequences, which are implicated in several neurological disorders, undergo liquid-liquid phase separation (LLPS) provided the number of repeat units, n, exceeds a critical value. Here, we establish a link between the folding landscapes of the monomers of trinucleotide repeats and their propensity to self-associate. Simulations using a coarse-grained Self-Organized Polymer (SOP) model for (CAG)n repeats in monovalent salt solutions reproduce experimentally measured melting temperatures, which are available only for small n. By extending the simulations to large n, we show that the free-energy gap, ΔGS, between the ground state (GS) and slipped hairpin (SH) states is a predictor of aggregation propensity. The GS for even n is a perfect hairpin (PH), whereas it is a SH when n is odd. The value of ΔGS (zero for odd n) is larger for even n than for odd n. As a result, the rate of dimer formation is slower in (CAG)30 relative to (CAG)31, thus linking ΔGS to RNA-RNA association. The yield of the dimer decreases dramatically, compared to the wild type, in mutant sequences in which the population of the SH decreases substantially. Association between RNA chains is preceded by a transition to the SH even if the GS is a PH. The finding that the excitation spectrum-which depends on the exact sequence, n, and ionic conditions-is a predictor of self-association should also hold for other RNAs (mRNA for example) that undergo LLPS.

Regulation of Gene Expression Through Effector-dependent Conformational Switching by Cobalamin Riboswitches

Riboswitches are an outstanding example of genetic regulation mediated by RNA conformational switching. In these non-coding RNA elements, the occupancy status of a ligand-binding domain governs the mRNA's decision to form one of two mutually exclusive structures in the downstream expression platform. Temporal constraints upon the function of many riboswitches, requiring folding of complex architectures and conformational switching in a limited co-transcriptional timeframe, make them ideal model systems for studying these processes. In this review, we focus on the mechanism of ligand-directed conformational changes in one of the most widely distributed riboswitches in bacteria: the cobalamin family. We describe the architectural features of cobalamin riboswitches whose structures have been determined by x-ray crystallography, which suggest a direct physical role of cobalamin in effecting the regulatory switch. Next, we discuss a series of experimental approaches applied to several model cobalamin riboswitches that interrogate these structural models. As folding is central to riboswitch function, we consider the differences in folding landscapes experienced by RNAs that are produced in vitro and those that are allowed to fold co-transcriptionally. Finally, we highlight a set of studies that reveal the difficulties of studying cobalamin riboswitches outside the context of transcription and that co-transcriptional approaches are essential for developing a more accurate picture of their structure-function relationships in these switches. This understanding will be essential for future advancements in the use of small-molecule guided RNA switches in a range of applications such as biosensors, RNA imaging tools, and nucleic acid-based therapies.

Cryo-EM reveals an entangled kinetic trap in the folding of a catalytic RNA

Functional RNAs fold through complex pathways that can contain misfolded "kinetic traps." A complete model of RNA folding requires understanding the formation of these misfolded states, but they are difficult to characterize because of their transient and potentially conformationally dynamic nature. We used cryo-electron microscopy (cryo-EM) to visualize a long-lived misfolded state in the folding pathway of the Tetrahymena thermophila group I intron, a paradigmatic RNA structure-function model system. The structure revealed how this state forms native-like secondary structure and tertiary contacts but contains two incorrectly crossed strands, consistent with a previous model. This incorrect topology mispositions a critical catalytic domain and cannot be resolved locally as extensive refolding is required. This work provides a structural framework for interpreting decades of biochemical and functional studies and demonstrates the power of cryo-EM for the exploration of RNA folding pathways.

Nucleic Acid Thermodynamics Derived from Mechanical Unzipping Experiments

Force-spectroscopy techniques have led to significant progress in studying the physicochemical properties of biomolecules that are not accessible in bulk assays. The application of piconewton forces with laser optical tweezers to single nucleic acids has permitted the characterization of molecular thermodynamics and kinetics with unprecedented accuracy. Some examples are the hybridization reaction between complementary strands in DNA and the folding of secondary, tertiary, and other heterogeneous structures, such as intermediate and misfolded states in RNA. Here we review the results obtained in our lab on deriving the nearest-neighbor free energy parameters in DNA and RNA duplexes from mechanical unzipping experiments. Remarkable nonequilibrium effects are also observed, such as the large irreversibility of RNA unzipping and the formation of non-specific secondary structures in single-stranded DNA. These features originate from forming stem-loop structures along the single strands of the nucleic acid. The recently introduced barrier energy landscape model quantifies kinetic trapping effects due to stem-loops being applicable to both RNA and DNA. The barrier energy landscape model contains the essential features to explain the many behaviors observed in heterogeneous nucleic-acid folding.

TPP Riboswitch Populates Holo-Form-like Structure Even in the Absence of Cognate Ligand at High Mg2+ Concentration

Riboswitches are noncoding RNA that regulate gene expression by folding into specific three-dimensional structures (holo-form) upon binding by their cognate ligand in the presence of Mg²⁺. Riboswitch functioning is also hypothesized to be under kinetic control requiring large cognate ligand concentrations. We ask the question under thermodynamic conditions, can the riboswitches populate structures similar to the holo-form only in the presence of Mg²⁺ and absence of cognate ligand binding. We addressed this question using thiamine pyrophosphate (TPP) riboswitch as a model system and computer simulations using a coarse-grained model for RNA. The folding free energy surface (FES) shows that with the initial increase in Mg²⁺ concentration ([Mg²⁺]), the aptamer domain (AD) of TPP riboswitch undergoes a barrierless collapse in its dimensions. On further increase in [Mg²⁺], intermediates separated by barriers appear on the FES, and one of the intermediates has a TPP ligand-binding competent structure. We show that site-specific binding of the Mg²⁺ aids in the formation of tertiary contacts. For [Mg²⁺] greater than physiological concentration, AD folds into a structure similar to the crystal structure of the TPP holo-form even in the absence of the TPP ligand. The folding kinetics shows that TPP AD populates an intermediate due to the misalignment of two arms present in the structure, which acts as a kinetic trap, leading to larger folding timescales. The predictions of the intermediate structures from the simulations are amenable for experimental verification.

A Computational Study of RNA Tetraloop Thermodynamics, Including Misfolded States

An important characteristic of RNA folding is the adoption of alternative configurations of similar stability, often referred to as misfolded configurations. These configurations are considered to compete with correctly folded configurations, although their rigorous thermodynamic and structural characterization remains elusive. Tetraloop motifs found in large ribozymes are ideal systems for an atomistically detailed computational quantification of folding free energy landscapes and the structural characterization of their constituent free energy basins, including nonnative states. In this work, we studied a group of closely related 10-mer tetraloops using a combined parallel tempering and metadynamics technique that allows a reliable sampling of the free energy landscapes, requiring only knowledge that the stem folds into a canonical A-RNA configuration. We isolated and analyzed unfolded, folded, and misfolded populations that correspond to different free energy basins. We identified a distinct misfolded state that has a stability very close to that of the correctly folded state. This misfolded state contains a predominant population that shares the same structural features across all tetraloops studied here and lacks the noncanonical A-G base pair in its loop portion. Further analysis performed with biased trajectories showed that although this competitive misfolded state is not an essential intermediate, it is visited in most of the transitions from unfolded to correctly folded states. Moreover, the tetraloops can transition from this misfolded state to the correctly folded state without requiring extensive unfolding.

Action and function of helicases on RNA G-quadruplexes

•

Methodological progresses and piling evidence prove the rG4 biology in vivo.

•

rG4s step in virtually every aspect of RNA biology.

•

Helicases unwinding of rG4s is a fine regulatory layer to the downstream processes and general cell homeostasis.

•

The current knowledge is however limited to a few cell lines.

•

The regulation of helicases themselves is delineating as a important question.

•

Non-helicase rG4-processing proteins likely play a role.

The nucleic acid structure called G-quadruplex (G4) is currently discussed to function in nucleic acid-based mechanisms that influence several cellular processes. They can modulate the cellular machinery either positively or negatively, both at the DNA and RNA level. The majority of what we know about G4 biology comes from DNA G4 (dG4) research. RNA G4s (rG4), on the other hand, are gaining interest as researchers become more aware of their role in several aspects of cellular homeostasis. In either case, the correct regulation of G4 structures within cells is essential and demands specialized proteins able to resolve them. Small changes in the formation and unfolding of G4 structures can have severe consequences for the cells that could even stimulate genome instability, apoptosis or proliferation. Helicases are the most relevant negative G4 regulators, which prevent and unfold G4 formation within cells during different pathways. Yet, and despite their importance only a handful of rG4 unwinding helicases have been identified and characterized thus far. This review addresses the current knowledge on rG4s-processing helicases with a focus on methodological approaches. An example of a non-helicase rG4s-unwinding protein is also briefly described.

Effects of chain length of polyethylene glycol molecular crowders on a mutant Tetrahymena group I ribozyme lacking large peripheral module

While current group I ribozymes use several distinct strategies to function under conditions of low Mg²⁺ concentration (≤ 3 mM), a deletion mutant of the Tetrahymena ribozyme (ΔP5 ribozyme) is virtually inactive with 3 mM Mg²⁺ due to removal of the large peripheral module, P5abc, supporting the active conformation of the core module. We investigated the molecular crowding effects of synthetic polyethylene glycols (PEGs) on the activity of the ΔP5 ribozyme. Among PEG molecules with different chain lengths, PEG600 improved the activity of the ΔP5 ribozyme most effectively in the presence of 3 mM Mg²⁺.

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.

The Thermus thermophilus DEAD-box protein Hera is a general RNA binding protein and plays a key role in tRNA metabolism

RNA helicases catalyze the ATP-dependent destabilization of RNA duplexes. DEAD-box helicases share a helicase core that mediates ATP binding and hydrolysis, RNA binding and unwinding. Most members of this family contain domains flanking the core that can confer RNA substrate specificity and guide the helicase to a specific RNA. However, the in vivo RNA substrates of most helicases are currently not defined. The DEAD-box helicase Hera from Thermus thermophilus contains a helicase core, followed by a dimerization domain and an RNA binding domain that folds into an RNA recognition motif (RRM). The RRM mediates high affinity binding to an RNA hairpin, and an adjacent duplex is then unwound by the helicase core. Hera is a cold-shock protein, and has been suggested to act as an RNA chaperone under cold-shock conditions. Using crosslinking immunoprecipitation of Hera/RNA complexes and sequencing, we show that Hera binds to a large fraction of T. thermophilus RNAs under normal-growth and cold-shock conditions without a strong sequence preference, in agreement with a structure-specific recognition of RNAs and a general function in RNA metabolism. Under cold-shock conditions, Hera is recruited to RNAs with high propensities to form stable secondary structures. We show that selected RNAs identified, including a set of tRNAs, bind to Hera in vitro, and activate the Hera helicase core. Gene ontology analysis reveals an enrichment of genes related to translation, including mRNAs of ribosomal proteins, tRNAs, tRNA ligases, and tRNA-modifying enzymes, consistent with a key role of Hera in ribosome and tRNA metabolism.

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.

Kinetic trapping in protein folding

The founding principles of protein folding introduced by Christian Anfinsen, together with the numerous mechanistic investigations that followed, assume that protein folding is a thermodynamically controlled process. On the other hand, this review underscores the fact that thermodynamic control is far from being the norm in protein folding, as long as one considers an extended chemical-potential landscape encompassing aggregates, in addition to native, unfolded and intermediate states. Here, we highlight the key role of kinetic trapping of the protein native state relative to unfolded, intermediate and, most importantly, aggregated states. We propose that kinetic trapping serves an important role in biology by protecting the bioactive states of a large number of proteins from deleterious aggregation. In the event that undesired aggregates were somehow formed, specialized intracellular disaggregation machines have evolved to convert any aberrant populations back to the native state, thus restoring a fully bioactive and aggregation-protected protein cohort.

Iterative annealing mechanism explains the functions of the GroEL and RNA chaperones

Molecular chaperones are ATP-consuming machines, which facilitate the folding of proteins and RNA molecules that are kinetically trapped in misfolded states. Unassisted folding occurs by the kinetic partitioning mechanism according to which folding to the native state, with low probability as well as misfolding to one of the many metastable states, with high probability, occur rapidly. GroEL is an all-purpose stochastic machine that assists misfolded substrate proteins to fold. The RNA chaperones such as CYT-19, which are ATP-consuming enzymes, help the folding of ribozymes that get trapped in metastable states for long times. GroEL does not interact with the folded proteins but CYT-19 disrupts both the folded and misfolded ribozymes. The structures of GroEL and RNA chaperones are strikingly different. Despite these differences, the iterative annealing mechanism (IAM) quantitatively explains all the available experimental data for assisted folding of proteins and ribozymes. Driven by ATP binding and hydrolysis and GroES binding, GroEL undergoes a catalytic cycle during which it samples three allosteric states, T (apo), R (ATP bound), and R^″ (ADP bound). Analyses of the experimental data show that the efficiency of the GroEL-GroES machinery and mutants is determined by the resetting rate k _{R ″  → T} , which is largest for the wild-type (WT) GroEL. Generalized IAM accurately predicts the folding kinetics of Tetrahymena ribozyme and its variants. Chaperones maximize the product of the folding rate and the steady-state native state fold by driving the substrates out of equilibrium. Neither the absolute yield nor the folding rate is optimized.

Real-time assembly of ribonucleoprotein complexes on nascent RNA transcripts

Cellular protein-RNA complexes assemble on nascent transcripts, but methods to observe transcription and protein binding in real time and at physiological concentrations are not available. Here, we report a single-molecule approach based on zero-mode waveguides that simultaneously tracks transcription progress and the binding of ribosomal protein S15 to nascent RNA transcripts during early ribosome biogenesis. We observe stable binding of S15 to single RNAs immediately after transcription for the majority of the transcripts at 35 °C but for less than half at 20 °C. The remaining transcripts exhibit either rapid and transient binding or are unable to bind S15, likely due to RNA misfolding. Our work establishes the foundation for studying transcription and its coupled co-transcriptional processes, including RNA folding, ligand binding, and enzymatic activity such as in coupling of transcription to splicing, ribosome assembly or translation.

The early steps of ribosome assembly occur co-transcriptionally on the nascent ribosomal RNA. Here the authors demonstrate an approach that allows simultaneous monitoring of transcription and ribosomal protein assembly at the single-molecule level in real time.

The Story of RNA Folding, as Told in Epochs

The past decades have witnessed tremendous developments in our understanding of RNA biology. At the core of these advances have been studies aimed at discerning RNA structure and at understanding the forces that influence the RNA folding process. It is easy to take the present state of understanding for granted, but there is much to be learned by considering the path to our current understanding, which has been tortuous, with the birth and death of models, the adaptation of experimental tools originally developed for characterization of protein structure and catalysis, and the development of novel tools for probing RNA. In this review we tour the stages of RNA folding studies, considering them as "epochs" that can be generalized across scientific disciplines. These epochs span from the discovery of catalytic RNA, through biophysical insights into the putative primordial RNA World, to characterization of structured RNAs, the building and testing of models, and, finally, to the development of models with the potential to yield generalizable predictive and quantitative models for RNA conformational, thermodynamic, and kinetic behavior. We hope that this accounting will aid others as they navigate the many fascinating questions about RNA and its roles in biology, in the past, present, and future.

Effects of flanking regions on HDV cotranscriptional folding kinetics

Hepatitis delta virus (HDV) ribozyme performs the self-cleavage activity through folding to a double pseudoknot structure. The folding of functional RNA structures is often coupled with the transcription process. In this work, we developed a new approach for predicting the cotranscriptional folding kinetics of RNA secondary structures with pseudoknots. We theoretically studied the cotranscriptional folding behavior of the 99-nucleotide (nt) HDV sequence, two upstream flanking sequences, and one downstream flanking sequence. During transcription, the 99-nt HDV can effectively avoid the trap intermediates and quickly fold to the cleavage-active state. It is different from its refolding kinetics, which folds into an intermediate trap state. For all the sequences, the ribozyme regions (from 1 to 73) all fold to the same structure during transcription. However, the existence of the 30-nt upstream flanking sequence can inhibit the ribozyme region folding into the active native state through forming an alternative helix Alt1 with the segments 70-90. The longer upstream flanking sequence of 54 nt itself forms a stable hairpin structure, which sequesters the formation of the Alt1 helix and leads to rapid formation of the cleavage-active structure. Although the 55-nt downstream flanking sequence could invade the already folded active structure during transcription by forming a more stable helix with the ribozyme region, the slow transition rate could keep the structure in the cleavage-active structure to perform the activity.

Single-round isolation of diverse RNA aptamers from a random sequence pool

Aptamers are oligonucleotide ligands with specific binding affinity to target molecules. Generally, RNA aptamers are selected from an RNA pool with random sequences, using the technique termed SELEX, in which the target-binding RNA molecules are repeatedly isolated and exponentially amplified. Despite several advantages, SELEX often produces uncertain results during the iterative amplifications of the rare target-binding RNA molecules. Here, we develop a non-repeated, primer-less and target immobilization-free isolation method for generating RNA aptamers, which is robust to experimental noise. Uniquely, this method focuses on finding and removal of non-aptamer sequences from the RNA pool by RNase digestion leaving target-bound aptamer molecules, and thus is independent of aptamer types. The undigested RNA sequences remaining are so few in number that they must be mixed with a large excess of a known sequence for further manipulations and this sequence is then removed by restriction digestion followed by high-throughput sequencing analysis to identify aptamers. Using this method, we generated multiple RNA aptamers targeting α-thrombin and TGFβ1 proteins, independently. This method potentially generates thousands of sequences as aptamer candidates, which may enable us to predict a common average sequence or structural property of these aptamers that is different from input RNA.

Synergistic SHAPE/Single-Molecule Deconvolution of RNA Conformation under Physiological Conditions

Structural RNA domains are widely involved in the regulation of biological functions, such as gene expression, gene modification, and gene repair. Activity of these dynamic regions depends sensitively on the global fold of the RNA, in particular, on the binding affinity of individual conformations to effector molecules in solution. Consequently, both the 1) structure and 2) conformational dynamics of noncoding RNAs prove to be essential in understanding the coupling that results in biological function. Toward this end, we recently reported observation of three conformational states in the metal-induced folding pathway of the tRNA-like structure domain of Brome Mosaic Virus, via single-molecule fluorescence resonance energy transfer studies. We report herein selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE)-directed structure predictions as a function of metal ion concentrations ([Mn+]) to confirm the three-state folding model, as well as test 2° structure models from the literature. Specifically, SHAPE reactivity data mapped onto literature models agrees well with the secondary structures observed at 0-10 mM [Mg2+], with only minor discrepancies in the E hairpin domain at low [Mg2+]. SHAPE probing and SHAPE-directed structure predictions further confirm the stepwise unfolding pathway previously observed in our single-molecule studies. Of special relevance, this means that reduction in the metal-ion concentration unfolds the 3' pseudoknot interaction before unfolding the long-range stem interaction. This work highlights the synergistic power of combining 1) single-molecule Förster resonance energy transfer and 2) SHAPE-directed structure-probing studies for detailed analysis of multiple RNA conformational states. In particular, single-molecule guided deconvolution of the SHAPE reactivities permits 2° structure predictions of isolated RNA conformations, thereby substantially improving on traditional limitations associated with current structure prediction algorithms.

Marginal stability drives irreversible unfolding of large multi-domain family 3 glycosylhydrolases from thermo-tolerant yeast

Protein folding is an extremely complex and fast, yet perfectly defined process, involving interplay of many intra and inter-molecular forces. In vitro, these molecular interactions are reversible for many proteins e.g., smaller and monomeric, organized into single domains. However, refolding of larger multi-domain/multimeric proteins is much more complicated, proceeds in a hierarchal way and is often irreversible. In a comparative study on two large, multi-domain and multimeric isozymes, β-glucosidase I (BGLI) and β-glucosidase II (BGLII) from Pichia etchellsii, we studied spontaneous and assisted refolding under three denaturing conditions viz. GdnHCl, alkaline pH and heat. During refolding, higher refolding yields were obtained for BGLII in case of pH induced unfolding (13.89%±0.25) than BGLI (6%±0.85) while for GdnHCl induced unfolding, refolding was marginal (BGLI=5%±0.5; BGLII=6%±0.69). Thermal unfolding was irreversible while assisted refolding also showed little structural gain for both proteins. When the apparent free energies of unfolding (ΔG_U^app) were calculated from GdnHCl unfolding data, their values were strikingly found to be lower (BGLI ΔG_U^app=3.02kcal/mol; BGLII ΔG_U^app=2.99kcal/mol) than reported for globular (ΔG_U=5-15kcal/mol)/multimeric proteins (ΔG_U=23-29kcal/mol) indicating marginal stability results in low refolding.

In vivo probing of nascent RNA structures reveals principles of cotranscriptional folding

Defining the in vivo folding pathway of cellular RNAs is essential to understand how they reach their final native conformation. We here introduce a novel method, named Structural Probing of Elongating Transcripts (SPET-seq), that permits single-base resolution analysis of transcription intermediates' secondary structures on a transcriptome-wide scale, enabling base-resolution analysis of the RNA folding events. Our results suggest that cotranscriptional RNA folding in vivo is a mixture of cooperative folding events, in which local RNA secondary structure elements are formed as they get transcribed, and non-cooperative events, in which 5'-halves of long-range helices get sequestered into transient non-native interactions until their 3' counterparts have been transcribed. Together our work provides the first transcriptome-scale overview of RNA cotranscriptional folding in a living organism.

Conditioning and Robustness of RNA Boltzmann Sampling under Thermodynamic Parameter Perturbations

Understanding how RNA secondary structure prediction methods depend on the underlying nearest-neighbor thermodynamic model remains a fundamental challenge in the field. Minimum free energy (MFE) predictions are known to be "ill conditioned" in that small changes to the thermodynamic model can result in significantly different optimal structures. Hence, the best practice is now to sample from the Boltzmann distribution, which generates a set of suboptimal structures. Although the structural signal of this Boltzmann sample is known to be robust to stochastic noise, the conditioning and robustness under thermodynamic perturbations have yet to be addressed. We present here a mathematically rigorous model for conditioning inspired by numerical analysis, and also a biologically inspired definition for robustness under thermodynamic perturbation. We demonstrate the strong correlation between conditioning and robustness and use its tight relationship to define quantitative thresholds for well versus ill conditioning. These resulting thresholds demonstrate that the majority of the sequences are at least sample robust, which verifies the assumption of sampling's improved conditioning over the MFE prediction. Furthermore, because we find no correlation between conditioning and MFE accuracy, the presence of both well- and ill-conditioned sequences indicates the continued need for both thermodynamic model refinements and alternate RNA structure prediction methods beyond the physics-based ones.

Molecular simulations and Markov state modeling reveal the structural diversity and dynamics of a theophylline-binding RNA aptamer in its unbound state

RNA aptamers are oligonucleotides that bind with high specificity and affinity to target ligands. In the absence of bound ligand, secondary structures of RNA aptamers are generally stable, but single-stranded and loop regions, including ligand binding sites, lack defined structures and exist as ensembles of conformations. For example, the well-characterized theophylline-binding aptamer forms a highly stable binding site when bound to theophylline, but the binding site is unstable and disordered when theophylline is absent. Experimental methods have not revealed at atomic resolution the conformations that the theophylline aptamer explores in its unbound state. Consequently, in the present study we applied 21 microseconds of molecular dynamics simulations to structurally characterize the ensemble of conformations that the aptamer adopts in the absence of theophylline. Moreover, we apply Markov state modeling to predict the kinetics of transitions between unbound conformational states. Our simulation results agree with experimental observations that the theophylline binding site is found in many distinct binding-incompetent states and show that these states lack a binding pocket that can accommodate theophylline. The binding-incompetent states interconvert with binding-competent states through structural rearrangement of the binding site on the nanosecond to microsecond timescale. Moreover, we have simulated the complete theophylline binding pathway. Our binding simulations supplement prior experimental observations of slow theophylline binding kinetics by showing that the binding site must undergo a large conformational rearrangement after the aptamer and theophylline form an initial complex, most notably, a major rearrangement of the C27 base from a buried to solvent-exposed orientation. Theophylline appears to bind by a combination of conformational selection and induced fit mechanisms. Finally, our modeling indicates that when Mg²⁺ ions are present the population of binding-competent aptamer states increases more than twofold. This population change, rather than direct interactions between Mg²⁺ and theophylline, accounts for altered theophylline binding kinetics.

The RNA Epistructurome: Uncovering RNA Function by Studying Structure and Post-Transcriptional Modifications

A large fraction of higher metazoan genomes transcribe RNA molecules whose functions extend far beyond carrying instructions for protein synthesis. Although RNA is apparently a simple molecule, the ways in which it performs many of its functions have remained highly elusive for decades. As learned from studying ribosomal and transfer RNAs, two of the key features influencing the function of RNA are its structure and post-transcriptional modifications. A deep understanding of RNA function therefore requires rapid and straightforward approaches to study the complex and intricate landscape of RNA structures and modifications. In this review we summarize and discuss the most recent methods and findings in the field of RNA biology, with an eye toward new frontiers and open questions.

Kinetic and thermodynamic framework for P4-P6 RNA reveals tertiary motif modularity and modulation of the folding preferred pathway

The past decade has seen a wealth of 3D structural information about complex structured RNAs and identification of functional intermediates. Nevertheless, developing a complete and predictive understanding of the folding and function of these RNAs in biology will require connection of individual rate and equilibrium constants to structural changes that occur in individual folding steps and further relating these steps to the properties and behavior of isolated, simplified systems. To accomplish these goals we used the considerable structural knowledge of the folded, unfolded, and intermediate states of P4-P6 RNA. We enumerated structural states and possible folding transitions and determined rate and equilibrium constants for the transitions between these states using single-molecule FRET with a series of mutant P4-P6 variants. Comparisons with simplified constructs containing an isolated tertiary contact suggest that a given tertiary interaction has a stereotyped rate for breaking that may help identify structural transitions within complex RNAs and simplify the prediction of folding kinetics and thermodynamics for structured RNAs from their parts. The preferred folding pathway involves initial formation of the proximal tertiary contact. However, this preference was only ∼10 fold and could be reversed by a single point mutation, indicating that a model akin to a protein-folding contact order model will not suffice to describe RNA folding. Instead, our results suggest a strong analogy with a modified RNA diffusion-collision model in which tertiary elements within preformed secondary structures collide, with the success of these collisions dependent on whether the tertiary elements are in their rare binding-competent conformations.

Visualizing the formation of an RNA folding intermediate through a fast highly modular secondary structure switch

Intermediates play important roles in RNA folding but can be difficult to characterize when short-lived or not significantly populated. By combining (15)N relaxation dispersion NMR with chemical probing, we visualized a fast (kex=k1+k-1≈423 s(-1)) secondary structural switch directed towards a low-populated (∼3%) partially folded intermediate in tertiary folding of the P5abc subdomain of the 'Tetrahymena' group I intron ribozyme. The secondary structure switch changes the base-pairing register across the P5c hairpin, creating a native-like structure, and occurs at rates of more than two orders of magnitude faster than tertiary folding. The switch occurs robustly in the absence of tertiary interactions, Mg(2+) or even when the hairpin is excised from the three-way junction. Fast, highly modular secondary structural switches may be quite common during RNA tertiary folding where they may help smoothen the folding landscape by allowing folding to proceed efficiently via additional pathways.

Polyacrylamide Gel Electrophoresis for Purification of Large Amounts of RNA

Polyacrylamide gel electrophoresis (PAGE) constitutes a powerful technique for the efficient purification of RNA molecules dedicated to applications that require high purity levels. PAGE allows for the fractionation of RNA obtained from cell extracts, chemical or enzymatic synthesis, or modification experiments. Native or denaturing conditions can be chosen for analytical or preparative-scale separations and the nucleotide resolution can be tuned by changing the percentage and reticulation of the gel material. In this protocol, we focus on the preparation of milligram-scale amounts of ~200 nucleotides (nt) RNA molecules that were used in subsequent crystallization experiments.

A DEAD-box RNA helicase promotes thermodynamic equilibration of kinetically trapped RNA structures in vivo

RNAs must assemble into specific structures in order to carry out their biological functions, but in vitro RNA folding reactions produce multiple misfolded structures that fail to exchange with functional structures on biological time scales. We used carefully designed self-cleaving mRNAs that assemble through well-defined folding pathways to identify factors that differentiate intracellular and in vitro folding reactions. Our previous work showed that simple base-paired RNA helices form and dissociate with the same rate and equilibrium constants in vivo and in vitro. However, exchange between adjacent secondary structures occurs much faster in vivo, enabling RNAs to quickly adopt structures with the lowest free energy. We have now used this approach to probe the effects of an extensively characterized DEAD-box RNA helicase, Mss116p, on a series of well-defined RNA folding steps in yeast. Mss116p overexpression had no detectable effect on helix formation or dissociation kinetics or on the stability of interdomain tertiary interactions, consistent with previous evidence that intracellular factors do not affect these folding parameters. However, Mss116p overexpression did accelerate exchange between adjacent helices. The nonprocessive nature of RNA duplex unwinding by DEAD-box RNA helicases is consistent with a branch migration mechanism in which Mss116p lowers barriers to exchange between otherwise stable helices by the melting and annealing of one or two base pairs at interhelical junctions. These results suggest that the helicase activity of DEAD-box proteins like Mss116p distinguish intracellular RNA folding pathways from nonproductive RNA folding reactions in vitro and allow RNA structures to overcome kinetic barriers to thermodynamic equilibration in vivo.

Genome-wide profiling of mouse RNA secondary structures reveals key features of the mammalian transcriptome

Background

The understanding of RNA structure is a key feature toward the comprehension of RNA functions and mechanisms of action. In particular, non-coding RNAs are thought to exert their functions by specific secondary structures, but an efficient annotation on a large scale of these structures is still missing.

Results

By using a novel high-throughput method, named chemical inference of RNA structures, CIRS-seq, that uses dimethyl sulfate, and N-cyclohexyl- N'-(2-morpholinoethyl)carbodiimide metho-p-toluenesulfonate to modify RNA residues in single-stranded conformation within native deproteinized RNA secondary structures, we investigate the structural features of mouse embryonic stem cell transcripts. Our analysis reveals an unexpected higher structuring of the 5′ and 3′ untranslated regions compared to the coding regions, a reduced structuring at the Kozak sequence and stop codon, and a three-nucleotide periodicity across the coding region of messenger RNAs. We also observe that ncRNAs exhibit a higher degree of structuring with respect to protein coding transcripts. Moreover, we find that the Lin28a binding protein binds selectively to RNA motifs with a strong preference toward a single stranded conformation.

Conclusions

This work defines for the first time the complete RNA structurome of mouse embryonic stem cells, revealing an extremely distinct RNA structural landscape. These results demonstrate that CIRS-seq constitutes an important tool for the identification of native deproteinized RNA structures.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0491-2) contains supplementary material, which is available to authorized users.

Collapse and hybridization of RNA: view from replica technique approach

The replica technique method is applied to investigate the kinetic behavior of the coarse-grained model for the RNA molecule. A non-equilibrium phase transition of second order between the glassy phase and the ensemble of freely fluctuating structures has been observed. The non-equilibrium steady state is investigated as well and the thermodynamic characteristics of the system have been evaluated. The non-equilibrium behavior of the specific heat is discussed. Based on our analysis, we point out the state in the kinetic pathway in which the RNA molecule is most prone to hybridization.

The Globular State of the Single-Stranded RNA: Effect of the Secondary Structure Rearrangements

The mutual influence of the slow rearrangements of secondary structure and fast collapse of the long single-stranded RNA (ssRNA) in approximation of coarse-grained model is studied with analytic calculations. It is assumed that the characteristic time of the secondary structure rearrangement is much longer than that for the formation of the tertiary structure. A nonequilibrium phase transition of the 2nd order has been observed.

MD Simulations of tRNA and Aminoacyl-tRNA Synthetases: Dynamics, Folding, Binding, and Allostery

While tRNA and aminoacyl-tRNA synthetases are classes of biomolecules that have been extensively studied for decades, the finer details of how they carry out their fundamental biological functions in protein synthesis remain a challenge. Recent molecular dynamics (MD) simulations are verifying experimental observations and providing new insight that cannot be addressed from experiments alone. Throughout the review, we briefly discuss important historical events to provide a context for how far the field has progressed over the past few decades. We then review the background of tRNA molecules, aminoacyl-tRNA synthetases, and current state of the art MD simulation techniques for those who may be unfamiliar with any of those fields. Recent MD simulations of tRNA dynamics and folding and of aminoacyl-tRNA synthetase dynamics and mechanistic characterizations are discussed. We highlight the recent successes and discuss how important questions can be addressed using current MD simulations techniques. We also outline several natural next steps for computational studies of AARS:tRNA complexes.

Unzippers, resolvers and sensors: a structural and functional biochemistry tale of RNA helicases

The centrality of RNA within the biological world is an irrefutable fact that currently attracts increasing attention from the scientific community. The panoply of functional RNAs requires the existence of specific biological caretakers, RNA helicases, devoted to maintain the proper folding of those molecules, resolving unstable structures. However, evolution has taken advantage of the specific position and characteristics of RNA helicases to develop new functions for these proteins, which are at the interface of the basic processes for transference of information from DNA to proteins. RNA helicases are involved in many biologically relevant processes, not only as RNA chaperones, but also as signal transducers, scaffolds of molecular complexes, and regulatory elements. Structural biology studies during the last decade, founded in X-ray crystallography, have characterized in detail several RNA-helicases. This comprehensive review summarizes the structural knowledge accumulated in the last two decades within this family of proteins, with special emphasis on the structure-function relationships of the most widely-studied families of RNA helicases: the DEAD-box, RIG-I-like and viral NS3 classes.

Mapping RNA structure in vitro using nucleobase-specific probes

RNAs have to adopt specific three-dimensional structures to fulfill their biological functions. Therefore exploring RNA structure is of interest to understand RNA-dependent processes. Chemical probing in vitro is a very powerful tool to investigate RNA molecules under a variety of conditions. Among the most frequently used chemical reagents are the nucleobase-specific probes dimethyl sulfate (DMS), 1-cyclohexyl-3-(2-morpholinoethyl) carbodiimide metho-p-toluenesulfonate (CMCT) and β-ethoxy-α-ketobutyraldehyde (kethoxal). These chemical reagents modify nucleotides which are not involved in hydrogen bonding or protected by a ligand, such as proteins or metabolites. Upon performing modification reactions with all three chemicals the accessibility of all four nucleobases can be determined. With this fast and inexpensive method local changes in RNA secondary and tertiary structure, as well as the formation of contacts between RNA and its ligands can be detected independent of the RNA's length.

Folding pathways of the Tetrahymena ribozyme

Like many structured RNAs, the Tetrahymena group I intron ribozyme folds through multiple pathways and intermediates. Under standard conditions in vitro, a small fraction reaches the native state (N) with kobs ≈ 0.6 min(-1), while the remainder forms a long-lived misfolded conformation (M) thought to differ in topology. These alternative outcomes reflect a pathway that branches late in folding, after disruption of a trapped intermediate (Itrap). Here we use catalytic activity to probe the folding transitions from Itrap to the native and misfolded states. We show that mutations predicted to weaken the core helix P3 do not increase the rate of folding from Itrap but they increase the fraction that reaches the native state rather than forming the misfolded state. Thus, P3 is disrupted during folding to the native state but not to the misfolded state, and P3 disruption occurs after the rate-limiting step. Interestingly, P3-strengthening mutants also increase native folding. Additional experiments show that these mutants are rapidly committed to folding to the native state, although they reach the native state with approximately the same rate constant as the wild-type ribozyme (~1 min(-1)). Thus, the P3-strengthening mutants populate a distinct pathway that includes at least one intermediate but avoids the M state, most likely because P3 and the correct topology are formed early. Our results highlight multiple pathways in RNA folding and illustrate how kinetic competitions between rapid events can have long-lasting effects because the "choice" is enforced by energy barriers that grow larger as folding progresses.

Single-molecule studies of riboswitch folding

The folding dynamics of riboswitches are central to their ability to modulate gene expression in response to environmental cues. In most cases, a structural competition between the formation of a ligand-binding aptamer and an expression platform (or some other competing off-state) determines the regulatory outcome. Here, we review single-molecule studies of riboswitch folding and function, predominantly carried out using single-molecule FRET or optical trapping approaches. Recent results have supplied new insights into riboswitch folding energy landscapes, the mechanisms of ligand binding, the roles played by divalent ions, the applicability of hierarchical folding models, and kinetic vs. thermodynamic control schemes. We anticipate that future work, based on improved data sets and potentially combining multiple experimental techniques, will enable the development of more complete models for complex RNA folding processes. This article is part of a Special Issue entitled: Riboswitches.

RNA catalytic activity as a probe of chaperone-mediated RNA folding

For structured RNAs that possess catalytic activity, this activity provides a powerful probe for measuring the progress of folding and the effects of RNA chaperone proteins on the folding rate. The crux of this approach is that only the natively folded RNA is able to perform the catalytic reaction. This method can provide a quantitative measure of the fraction of native RNA over time, and it can readily distinguish the native state from all misfolded conformations. Here we describe an activity-based method measuring native folding of ribozymes derived from self-splicing group I introns, and we show how the assay can be used to monitor acceleration of native folding by DEAD-box RNA helicase proteins that function as general RNA chaperones. By measuring the amount of substrate that is converted to product in a rapid first turnover, we describe how to determine the fraction of the ribozyme population that is present in the native state. Further, we describe how to perform a two-stage or discontinuous assay in which folding proceeds in stage one and then solution conditions are changed in stage two to permit catalytic activity and block further folding. This protocol allows folding to be followed under a broad range of solution conditions, including those that do not support catalytic activity, and facilitates studies of chaperone proteins.

A chemogenetic approach to study the structural basis of protein-facilitated RNA folding

Large RNA molecules play important roles in all aspects of cellular metabolism ranging from mRNA splicing and protein biosynthesis to regulation of gene expression. In order to correctly perform its function in the cell, an RNA molecule must fold into a complex tertiary structure. Folding of many large RNAs is slow either due to formation of stable misfolded intermediates or due to high contact order or instability of obligate folding intermediates. Therefore many RNAs use protein cofactors to facilitate their folding in vivo. Folding of the yeast mitochondrial group II intron ai5γ to the native state under physiological conditions is facilitated by the protein cofactor Mss116. This chapter describes the use of Nucleotide Analog Interference Mapping (NAIM) to identify specific substructures within the intron molecule that are directly affected by the protein.

Chemical probing of RNA in living cells

RNAs need to adopt a specific architecture to exert their task in cells. While significant progress has been made in describing RNA folding landscapes in vitro, understanding intracellular RNA structure formation is still in its infancy. This is in part due to the complex nature of the cellular environment but also to the limited availability of suitable methodologies. To assess the intracellular structure of large RNAs, we recently applied a chemical probing technique and a metal-induced cleavage assay in vivo. These methods are based on the fact that small molecules, like dimethyl sulfate (DMS), or metal ions, such as Pb(2+), penetrate and spread throughout the cell very fast. Hence, these chemicals are able to modify accessible RNA residues or to induce cleavage of the RNA strand in the vicinity of a metal ion in living cells. Mapping of these incidents allows inferring information on the intracellular conformation, metal ion binding sites or ligand-induced structural changes of the respective RNA molecule. Importantly, in vivo chemical probing can be easily adapted to study RNAs in different cell types.

The kinetics of ribozyme cleavage: a tool to analyze RNA folding as a function of catalysis

As catalytically active RNAs, ribozymes can be characterized by kinetic measurements similar to classical enzyme kinetics. However, in contrast to standard protein enzymes, for which reactions can usually be started by mixing the enzyme with its substrate, ribozymes are typically self-cleaving. The reaction has to be initiated by folding the RNA into its active conformation. Thus, ribozyme kinetics are influenced by both folding and catalytic components and often enable indirect observation of RNA folding. Here, I describe how to obtain quantitative ribozyme cleavage data via denaturing polyacrylamide gel electrophoresis (PAGE) of radioactively labeled in vitro transcripts and discuss general considerations for subsequent kinetic analysis.

Sequence-dependent base-stacking stabilities guide tRNA folding energy landscapes

The folding of bacterial tRNAs with disparate sequences has been observed to proceed in distinct folding mechanisms despite their structural similarity. To explore the folding landscapes of tRNA, we performed ion concentration-dependent coarse-grained TIS model MD simulations of several E. coli tRNAs to compare their thermodynamic melting profiles to the classical absorbance spectra of Crothers and co-workers. To independently validate our findings, we also performed atomistic empirical force field MD simulations of tRNAs, and we compared the base-to-base distances from coarse-grained and atomistic MD simulations to empirical base-stacking free energies. We then projected the free energies to the secondary structural elements of tRNA, and we observe distinct, parallel folding mechanisms whose differences can be inferred on the basis of their sequence-dependent base-stacking stabilities. In some cases, a premature, nonproductive folding intermediate corresponding to the Ψ hairpin loop must backtrack to the unfolded state before proceeding to the folded state. This observation suggests a possible explanation for the fast and slow phases observed in tRNA folding kinetics.

The DEAD-box protein Dbp2 functions with the RNA-binding protein Yra1 to promote mRNP assembly

Eukaryotic gene expression involves numerous biochemical steps that are dependent on RNA structure and ribonucleoprotein (RNP) complex formation. The DEAD-box class of RNA helicases plays fundamental roles in formation of RNA and RNP structure in every aspect of RNA metabolism. In an effort to explore the diversity of biological roles for DEAD-box proteins, our laboratory recently demonstrated that the DEAD-box protein Dbp2 associates with actively transcribing genes and is required for normal gene expression in Saccharomyces cerevisiae. We now provide evidence that Dbp2 interacts genetically and physically with the mRNA export factor Yra1. In addition, we find that Dbp2 is required for in vivo assembly of mRNA-binding proteins Yra1, Nab2, and Mex67 onto poly(A)+ RNA. Strikingly, we also show that Dbp2 is an efficient RNA helicase in vitro and that Yra1 decreases the efficiency of ATP-dependent duplex unwinding. We provide a model whereby messenger ribonucleoprotein (mRNP) assembly requires Dbp2 unwinding activity and once the mRNP is properly assembled, inhibition by Yra1 prevents further rearrangements. Both Yra1 and Dbp2 are conserved in multicellular eukaryotes, suggesting that this constitutes a broadly conserved mechanism for stepwise assembly of mature mRNPs in the nucleus.

The long-range P3 helix of the Tetrahymena ribozyme is disrupted during folding between the native and misfolded conformations

RNAs are prone to misfolding, but how misfolded structures are formed and resolved remains incompletely understood. The Tetrahymena group I intron ribozyme folds in vitro to a long-lived misfolded conformation (M) that includes extensive native structure but is proposed to differ in topology from the native state (N). A leading model predicts that exchange of the topologies requires unwinding of the long-range, core helix P3, despite the presence of P3 in both conformations. To test this model, we constructed 16 mutations to strengthen or weaken P3. Catalytic activity and in-line probing showed that nearly all of the mutants form the M state before folding to N. The P3-weakening mutations accelerated refolding from M (3- to 30-fold) and the P3-strengthening mutations slowed refolding (6- to 1400-fold), suggesting that P3 indeed unwinds transiently. Upon depletion of Mg(2+), the mutations had analogous effects on unfolding from N to intermediates that subsequently fold to M. The magnitudes for the P3-weakening mutations were larger than in refolding from M, and small-angle X-ray scattering showed that the ribozyme expands rapidly to intermediates from which P3 is disrupted subsequently. These results are consistent with previous results indicating unfolding of native peripheral structure during refolding from M, which probably permits rearrangement of the core. Together, our results demonstrate that exchange of the native and misfolded conformations requires loss of a core helix in addition to peripheral structure. Further, the results strongly suggest that misfolding arises from a topological error within the ribozyme core, and a specific topology is proposed.

Coarse-grained model for predicting RNA folding thermodynamics

We present a thermodynamically robust coarse-grained model to simulate folding of RNA in monovalent salt solutions. The model includes stacking, hydrogen bond, and electrostatic interactions as fundamental components in describing the stability of RNA structures. The stacking interactions are parametrized using a set of nucleotide-specific parameters, which were calibrated against the thermodynamic measurements for single-base stacks and base-pair stacks. All hydrogen bonds are assumed to have the same strength, regardless of their context in the RNA structure. The ionic buffer is modeled implicitly, using the concept of counterion condensation and the Debye-Hückel theory. The three adjustable parameters in the model were determined by fitting the experimental data for two RNA hairpins and a pseudoknot. A single set of parameters provides good agreement with thermodynamic data for the three RNA molecules over a wide range of temperatures and salt concentrations. In the process of calibrating the model, we establish the extent of counterion condensation onto the single-stranded RNA backbone. The reduced backbone charge is independent of the ionic strength and is 60% of the RNA bare charge at 37 °C. Our model can be used to predict the folding thermodynamics for any RNA molecule in the presence of monovalent ions.

Molecular dynamics simulation of RNA pseudoknot unfolding pathway

Many biological functions of RNA molecules are related to their pseudoknot structures. It is significant for predicting the structure and function of RNA that learning about the stability and the process of RNA pseudoknot folding and unfolding. The structural features of mouse mammary tumor virus (MMTV) RNA pseudoknot in different ion concentration, the unfolding process of the RNA pseudoknot, and the two hairpin helices that constitute the RNA pseudoknot were studied with all atom molecule dynamics simulation method in this paper. We found that the higher cation concentration can cause structure of the RNA molecules more stable, and ions played an indispensable role in keeping the structure of RNA molecules stable; the unfolding process of hairpin structure was corresponding to the antiprocess of its folding process. The main pathway of pseudoknot unfolding was that the inner base pair opened first, and then, the two helices, which formed the RNA pseudoknot opened decussately, while the folding pathway of the RNA pseudoknot was a helix folding after formation of the other helix. Therefore, the unfolding process of RNA pseudoknot is different from the antiprocess of its folding process, and the unfolding process of each helix in the RNA pseudoknot is similar to the hairpin structure’s unfolding process, which means that both are the unzipping process.

Shape and secondary structure prediction for ncRNAs including pseudoknots based on linear SVM

Background

Accurate secondary structure prediction provides important information to undefirstafinding the tertiary structures and thus the functions of ncRNAs. However, the accuracy of the native structure derivation of ncRNAs is still not satisfactory, especially on sequences containing pseudoknots. It is recently shown that using the abstract shapes, which retain adjacency and nesting of structural features but disregard the length details of helix and loop regions, can improve the performance of structure prediction. In this work, we use SVM-based feature selection to derive the consensus abstract shape of homologous ncRNAs and apply the predicted shape to structure prediction including pseudoknots.

Results

Our approach was applied to predict shapes and secondary structures on hundreds of ncRNA data sets with and without psuedoknots. The experimental results show that we can achieve 18% higher accuracy in shape prediction than the state-of-the-art consensus shape prediction tools. Using predicted shapes in structure prediction allows us to achieve approximate 29% higher sensitivity and 10% higher positive predictive value than other pseudoknot prediction tools.

Conclusions

Extensive analysis of RNA properties based on SVM allows us to identify important properties of sequences and structures related to their shapes. The combination of mass data analysis and SVM-based feature selection makes our approach a promising method for shape and structure prediction. The implemented tools, Knot Shape and Knot Structure are open source software and can be downloaded at: http://www.cse.msu.edu/~achawana/KnotShape.