Isolation and characterization of a family of stable RNA tetraloops with the motif YNMG that participate in tertiary interactions

Sarbecovirus programmed ribosome frameshift RNA element folding studied by NMR spectroscopy and comparative analyses

The programmed ribosomal frameshift (PRF) region is found in the RNA genome of all coronaviruses and shifts the ribosome reading frame through formation of a three-stem pseudoknot structure, allowing the translation of essential viral proteins. Using NMR spectroscopy, comparative sequence analyses and functional assays we show that, in the absence of the ribosome, a 123-nucleotide sequence encompassing the PRF element of SARS-CoV-2 adopts a well-defined two-stem loop structure that is conserved in all SARS-like coronaviruses. In this conformation, the attenuator hairpin and slippery site nucleotides are exposed in the first stem-loop and two pseudoknot stems are present in the second stem-loop, separated by an 8-nucleotide bulge. Formation of the third pseudoknot stem depends on pairing between bulge nucleotides and base-paired nucleotides of the upstream stem-loop, as shown by a PRF construct where residues of the upstream stem were removed, which formed the pseudoknot structure and had increased frameshifting activity in a dual-luciferase assay. The base-pair switch driving PRF pseudoknot folding was found to be conserved in several human non-SARS coronaviruses. The collective results suggest that the frameshifting pseudoknot structure of these viruses only forms transiently in the presence of the translating ribosome. These findings clarify the frameshifting mechanism in coronaviruses and can have a beneficial impact on antiviral drug discovery.

sincFold: end-to-end learning of short- and long-range interactions in RNA secondary structure

MOTIVATION

Coding and noncoding RNA molecules participate in many important biological processes. Noncoding RNAs fold into well-defined secondary structures to exert their functions. However, the computational prediction of the secondary structure from a raw RNA sequence is a long-standing unsolved problem, which after decades of almost unchanged performance has now re-emerged due to deep learning. Traditional RNA secondary structure prediction algorithms have been mostly based on thermodynamic models and dynamic programming for free energy minimization. More recently deep learning methods have shown competitive performance compared with the classical ones, but there is still a wide margin for improvement.

RESULTS

In this work we present sincFold, an end-to-end deep learning approach, that predicts the nucleotides contact matrix using only the RNA sequence as input. The model is based on 1D and 2D residual neural networks that can learn short- and long-range interaction patterns. We show that structures can be accurately predicted with minimal physical assumptions. Extensive experiments were conducted on several benchmark datasets, considering sequence homology and cross-family validation. sincFold was compared with classical methods and recent deep learning models, showing that it can outperform the state-of-the-art methods.

High-resolution structure of stem-loop 4 from the 5'-UTR of SARS-CoV-2 solved by solution state NMR

We present the high-resolution structure of stem-loop 4 of the 5'-untranslated region (5_SL4) of the severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) genome solved by solution state nuclear magnetic resonance spectroscopy. 5_SL4 adopts an extended rod-like structure with a single flexible looped-out nucleotide and two mixed tandem mismatches, each composed of a G•U wobble base pair and a pyrimidine•pyrimidine mismatch, which are incorporated into the stem-loop structure. Both the tandem mismatches and the looped-out residue destabilize the stem-loop structure locally. Their distribution along the 5_SL4 stem-loop suggests a role of these non-canonical elements in retaining functionally important structural plasticity in particular with regard to the accessibility of the start codon of an upstream open reading frame located in the RNA's apical loop. The apical loop-although mostly flexible-harbors residual structural features suggesting an additional role in molecular recognition processes. 5_SL4 is highly conserved among the different variants of SARS-CoV-2 and can be targeted by small molecule ligands, which it binds with intermediate affinity in the vicinity of the non-canonical elements within the stem-loop structure.

Transcriptome-wide probing reveals RNA thermometers that regulate translation of glycerol permease genes in Bacillus subtilis

RNA structure regulates bacterial gene expression by several distinct mechanisms via environmental and cellular stimuli, one of which is temperature. While some genome-wide studies have focused on heat shock treatments and the subsequent transcriptomic changes, soil bacteria are less likely to experience such rapid and extreme temperature changes. Though RNA thermometers (RNATs) have been found in 5' untranslated leader regions (5' UTRs) of heat shock and virulence-associated genes, this RNA-controlled mechanism could regulate other genes as well. Using Structure-seq2 and the chemical probe dimethyl sulfate (DMS) at four growth temperatures ranging from 23°C to 42°C, we captured a dynamic response of the Bacillus subtilis transcriptome to temperature. Our transcriptome-wide results show RNA structural changes across all four temperatures and reveal nonmonotonic reactivity trends with increasing temperature. Then, focusing on subregions likely to contain regulatory RNAs, we examined 5' UTRs to identify large, local reactivity changes. This approach led to the discovery of RNATs that control the expression of glpF (glycerol permease) and glpT (glycerol-3-phosphate permease); expression of both genes increased with increased temperature. Results with mutant RNATs indicate that both genes are controlled at the translational level. Increased import of glycerols at high temperatures could provide thermoprotection to proteins.

Insights into the secondary and tertiary structure of the Bovine Viral Diarrhea Virus Internal Ribosome Entry Site

The internal ribosome entry site (IRES) RNA of bovine viral diarrhoea virus (BVDV), an economically significant Pestivirus, is required for the cap-independent translation of viral genomic RNA. Thus, it is essential for viral replication and pathogenesis. We applied a combination of high-throughput biochemical RNA structure probing (SHAPE-MaP) and in silico modelling approaches to gain insight into the secondary and tertiary structures of BVDV IRES RNA. Our study demonstrated that BVDV IRES RNA in solution forms a modular architecture composed of three distinct structural domains (I-III). Two regions within domain III are represented in tertiary interactions to form an H-type pseudoknot. Computational modelling of the pseudoknot motif provided a fine-grained picture of the tertiary structure and local arrangement of helices in the BVDV IRES. Furthermore, comparative genomics and consensus structure predictions revealed that the pseudoknot is evolutionarily conserved among many Pestivirus species. These studies provide detailed insight into the structural arrangement of BVDV IRES RNA H-type pseudoknot and encompassing motifs that likely contribute to the optimal functionality of viral cap-independent translation element.

RNA multimerization as an organizing force for liquid-liquid phase separation

RNA interactions are exceptionally strong and highly redundant. As such, nearly any two RNAs have the potential to interact with one another over relatively short stretches, especially at high RNA concentrations. This is especially true for pairs of RNAs that do not form strong self-structure. Such phenomena can drive liquid-liquid phase separation, either solely from RNA-RNA interactions in the presence of divalent or organic cations, or in concert with proteins. RNA interactions can drive multimerization of RNA strands via both base-pairing and tertiary interactions. In this article, we explore the tendency of RNA to form stable monomers, dimers, and higher order structures as a function of RNA length and sequence through a focus on the intrinsic thermodynamic, kinetic, and structural properties of RNA. The principles we discuss are independent of any specific type of biomolecular condensate, and thus widely applicable. We also speculate how external conditions experienced by living organisms can influence the formation of nonmembranous compartments, again focusing on the physical and structural properties of RNA. Plants, in particular, are subject to diverse abiotic stresses including extreme temperatures, drought, and salinity. These stresses and the cellular responses to them, including changes in the concentrations of small molecules such as polyamines, salts, and compatible solutes, have the potential to regulate condensate formation by melting or strengthening base-pairing. Reversible condensate formation, perhaps including regulation by circadian rhythms, could impact biological processes in plants, and other organisms.

Importance of Anion-π Interactions in RNA GAAA and GGAG Tetraloops: A Combined MD and QM Study

In this study, we demonstrate that anion−π interactions (an attractive noncovalent force between electron deficient π-systems and anions) are involved in the stabilization of GAAA and GGAG RNA tetraloops. Using the single recognition particle (SRP)–RNA complexes as a case of study, we combined molecular dynamics (MD) and quantum mechanics (QM) calculations to shed light on the structural influence of phosphate–G anion−π interactions and hydrogen bonds (HBs) involving K⁺/Mg²⁺ water clusters. In addition, the RNA assemblies herein were further characterized by means of the “atoms in molecules” (AIM) and noncovalent interactions plot (NCIplot) methodologies. We believe the results derived from this study might be important in the fields of chemical biology (RNA folding and engineering) and supramolecular chemistry (anion−π interactions) as well as to further expand the current knowledge regarding RNA structural motifs.

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

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

Monolithic alkylsilane column: A promising separation medium for oligonucleotides by ion-pair reversed-phase liquid chromatography

In this paper, a monolithic octadecylsilane column and particle-packed octadecylsilane columns were used to investigate the retention behaviors of oligonucleotides by ion-pair reversed-phase liquid chromatography (IP-RPLC). Results showed that, with same base composition, hairpin oligonucleotides always had weaker retention than corresponding random coil oligonucleotides on the monolithic column, but not on the particle-packed columns. In addition, the linear correlation between the retention factor k of oligonucleotides and the reciprocal of temperature (1/T), especially for hairpins, was relatively weaker on the particle-packed columns, as compared to the correlation on the monolithic column. The correlation between k and 1/T became weaker with decreasing particle size of the particle-packed columns. Moreover, results revealed that the overall retention order on the particle-packed column with small particles (3 μm) differed greatly from that on the monolithic column. In contrast, the retention order on the 10 μm particle-packed column was very close to that on the monolithic column. From the above, we inferred that oligonucleotides could keep their primary conformations unchanged when passing through the monolithic column, attributed to the special pore structures of the monolith. However, the conformations of oligonucleotides were suppressed or even destroyed when oligonucleotides passed through the particle-packed columns. This because the narrow and tortuous channels created by the stacked stationary phase particles could lead to more complex and unequable retention behaviors. Therefore, the monolithic column exhibited better retention regularity for oligonucleotides of secondary structure especially for hairpins than the particle-packed columns. It is noteworthy that the monolith-based IP-RPLC opens an intriguing prospect in accurately elucidating the retention behaviors of oligonucleotides.

Interactions of Protonated Guanidine and Guanidine Derivatives with Multiply Deprotonated RNA Probed by Electrospray Ionization and Collisionally Activated Dissociation

Interactions of ribonucleic acid (RNA) with guanidine and guanidine derivatives are important features in RNA-protein and RNA-drug binding. Here we have investigated noncovalently bound complexes of an 8-nucleotide RNA and six different ligands, all of which have a guanidinium moiety, by using electrospray ionization (ESI) and collisionally activated dissociation (CAD) mass spectrometry (MS). The order of complex stability correlated almost linearly with the number of ligand atoms that can potentially be involved in hydrogen-bond or salt-bridge interactions with the RNA, but not with the proton affinity of the ligands. However, ligand dissociation of the complex ions in CAD was generally accompanied by proton transfer from ligand to RNA, which indicated conversion of salt-bridge into hydrogen-bond interactions. The relative stabilities and dissociation pathways of [RNA+m L-n H] n- complexes with different stoichiometries (m=1-5) and net charge (n= 2-5) revealed both specific and unspecific ligand binding to the RNA.

Effects of osmolytes on stable UUCG tetraloops and their preference for a CG closing base pair

Osmolytes have the potential to affect the stability of secondary structure motifs and alter preferences for conserved nucleic acid sequences in the cell. To contribute to the understanding of the in vivo function of RNA we observed the effects of different classes of osmolytes on the UNCG tetraloop motif. UNCG tetraloops are the most common and stable of the RNA tetraloops and are nucleation sites for RNA folding. They also have a significant thermodynamic preference for a CG closing base pair. The thermal denaturation of model hairpins containing UUCG loops was monitored using UV-Vis spectroscopy in the presence of osmolytes with different chemical properties. Interestingly, all of the osmolytes tested destabilized the hairpins, but all had little effect on the thermodynamic preference for a CG base pair, except for polyethylene glycol (PEG) 200. PEG 200 destabilized the loop with the CG closing base pair relative to the loop with a GC closing base pair. The destabilization was linear with increasing concentrations of PEG 200, and the slope of this relationship was not perturbed by changes in the hairpin stem outside of the closing pair. This result suggests that in the presence of PEG 200, the UUCG loop with a GC closing base pair may retain some preferential interactions with the cosolute that are lost in the presence of the CG closing base pair. These results reveal that relatively small structural changes may influence how osmolytes tune the stability, and thus the function of a secondary structure motif in vivo.

Independent evolution of tetraloop in enterovirus oriL replicative element and its putative binding partners in virus protein 3C

Background

Enteroviruses are small non-enveloped viruses with a (+) ssRNA genome with one open reading frame. Enterovirus protein 3C (or 3CD for some species) binds the replicative element oriL to initiate replication. The replication of enteroviruses features a low-fidelity process, which allows the virus to adapt to the changing environment on the one hand, and requires additional mechanisms to maintain the genome stability on the other. Structural disturbances in the apical region of oriL domain d can be compensated by amino acid substitutions in positions 154 or 156 of 3C (amino acid numeration corresponds to poliovirus 3C), thus suggesting the co-evolution of these interacting sequences in nature. The aim of this work was to understand co-evolution patterns of two interacting replication machinery elements in enteroviruses, the apical region of oriL domain d and its putative binding partners in the 3C protein.

Methods

To evaluate the variability of the domain d loop sequence we retrieved all available full enterovirus sequences (>6, 400 nucleotides), which were present in the NCBI database on February 2017 and analysed the variety and abundance of sequences in domain d of the replicative element oriL and in the protein 3C.

Results

A total of 2,842 full genome sequences was analysed. The majority of domain d apical loops were tetraloops, which belonged to consensus YNHG (Y = U/C, N = any nucleotide, H = A/C/U). The putative RNA-binding tripeptide 154–156 (Enterovirus C 3C protein numeration) was less diverse than the apical domain d loop region and, in contrast to it, was species-specific.

Discussion

Despite the suggestion that the RNA-binding tripeptide interacts with the apical region of domain d, they evolve independently in nature. Together, our data indicate the plastic evolution of both interplayers of 3C-oriL recognition.

All-Atom MD Predicts Magnesium-Induced Hairpin in Chemically Perturbed RNA Analog of F10 Therapeutic

Given their increasingly frequent usage, understanding the chemical and structural properties which allow therapeutic nucleic acids to promote the death of cancer cells is critical for medical advancement. One molecule of interest is a 10-mer of FdUMP (5-fluoro-2'-deoxyuridine-5'-O-monophosphate) also called F10. To investigate causes of structural stability, we have computationally restored the 2' oxygen on each ribose sugar of the phosphodiester backbone, creating FUMP[10]. Microsecond time-scale, all-atom, simulations of FUMP[10] in the presence of 150 mM MgCl₂ predict that the strand has a 45% probability of folding into a stable hairpin-like secondary structure. Analysis of 16 μs of data reveals phosphate interactions as likely contributors to the stability of this folded state. Comparison with polydT and polyU simulations predicts that FUMP[10]'s lowest order structures last for one to 2 orders of magnitude longer than similar nucleic acid strands. Here we provide a brief structural and conformational analysis of the predicted structures of FUMP[10], and suggest insights into its stability via comparison to F10, polydT, and polyU.

Mapping the Universe of RNA Tetraloop Folds

We report a map of RNA tetraloop conformations constructed by calculating pairwise distances among all experimentally determined four-nucleotide hairpin loops. Tetraloops with similar structures are clustered together and, as expected, the two largest clusters are the canonical GNRA and UNCG folds. We identify clusters corresponding to known tetraloop folds such as GGUG, RNYA, AGNN, and CUUG. These clusters are represented in a simple two-dimensional projection that recapitulates the relationship among the different folds. The cluster analysis also identifies 20 novel tetraloop folds that are peculiar to specific positions in ribosomal RNAs and that are stabilized by tertiary interactions. In our RNA tetraloop database we find a significant number of non-GNRA and non-UNCG sequences adopting the canonical GNRA and UNCG folds. Conversely, we find a significant number of GNRA and UNCG sequences adopting non-GNRA and non-UNCG folds. Our analysis demonstrates that there is not a simple one-to-one, but rather a many-to-many mapping between tetraloop sequence and tetraloop fold.

Understanding in-line probing experiments by modeling cleavage of nonreactive RNA nucleotides

Ribonucleic acid (RNA) is involved in many regulatory and catalytic processes in the cell. The function of any RNA molecule is intimately related with its structure. In-line probing experiments provide valuable structural data sets for a variety of RNAs and are used to characterize conformational changes in riboswitches. However, the structural determinants that lead to differential reactivities in unpaired nucleotides have not been investigated yet. In this work, we used a combination of theoretical approaches, i.e., classical molecular dynamics simulations, multiscale quantum mechanical/molecular mechanical calculations, and enhanced sampling techniques in order to compute and interpret the differential reactivity of individual residues in several RNA motifs, including members of the most important GNRA and UNCG tetraloop families. Simulations on the multinanosecond timescale are required to converge the related free-energy landscapes. The results for uGAAAg and cUUCGg tetraloops and double helices are compared with available data from in-line probing experiments and show that the introduced technique is able to distinguish between nucleotides of the uGAAAg tetraloop based on their structural predispositions toward phosphodiester backbone cleavage. For the cUUCGg tetraloop, more advanced ab initio calculations would be required. This study is the first attempt to computationally classify chemical probing experiments and paves the way for an identification of tertiary structures based on the measured reactivity of nonreactive nucleotides.

Revisiting GNRA and UNCG folds: U-turns versus Z-turns in RNA hairpin loops

When thinking about RNA three-dimensional structures, coming across GNRA and UNCG tetraloops is perceived as a boon since their folds have been extensively described. Nevertheless, analyzing loop conformations within RNA and RNP structures led us to uncover several instances of GNRA and UNCG loops that do not fold as expected. We noticed that when a GNRA does not assume its "natural" fold, it adopts the one we typically associate with a UNCG sequence. The same folding interconversion may occur for loops with UNCG sequences, for instance within tRNA anticodon loops. Hence, we show that some structured tetranucleotide sequences starting with G or U can adopt either of these folds. The underlying structural basis that defines these two fold types is the mutually exclusive stacking of a backbone oxygen on either the first (in GNRA) or the last nucleobase (in UNCG), generating an oxygen-π contact. We thereby propose to refrain from using sequences to distinguish between loop conformations. Instead, we suggest using descriptors such as U-turn (for "GNRA-type" folds) and a newly described Z-turn (for "UNCG-type" folds). Because tetraloops adopt for the largest part only two (inter)convertible turns, we are better able to interpret from a structural perspective loop interchangeability occurring in ribosomes and viral RNA. In this respect, we propose a general view on the inclination for a given sequence to adopt (or not) a specific fold. We also suggest how long-noncoding RNAs may adopt discrete but transient structures, which are therefore hard to predict.

Mutational robustness and resilience of a replicative cis-element of RNA virus: Promiscuity, limitations, relevance

Since replication of RNA-viruses is generally a low-fidelity process, it would be advantageous, if specific interactions of their genomic cis-elements with dedicated ligands are relatively tolerant to mutations. The specificity/promiscuity trade-off of such interactions was addressed here by investigating structural requirements of the oriL (also known as the clover leaf-like element), of poliovirus RNA, a replicative cis-element containing a conserved essential tetraloop functionally interacting with the viral protein 3CD. The sequence of this tetraloop and 2 adjacent base-pairs was randomized in the viral genome, and viable viruses were selected in susceptible cells. Strikingly, each position of this octanucleotide in 62 investigated viable viruses could be occupied by any nucleotide (with the exception of one position, which lacked U), though with certain sequence preferences, confirmed by engineering mutant viral genomes whose phenotypic properties were found to correlate with the strength of the cis-element/ligand interaction. The results were compatible with a hypothesis that functional recognition by 3CD requires that this tetraloop should stably or temporarily adopt a YNMG-like (Y=U/C, N=any nucleotide, M=A/C) fold. The fitness of "weak" viruses could be increased by compensatory mutations "improving" the tetraloops. Otherwise, the recognition of "bad" tetraloops might be facilitated by alterations in the 3CD protein. The virus appeared to tolerate mutations in its cis-element relaying on either robustness (spatial structure degeneracy) or resilience (a combination of dynamic RNA folding, low-fidelity replication modifying the cis-element or its ligand, and negative selection). These mechanisms (especially resilience involving metastable low-fit intermediates) can also contribute to the viral evolvability.

Nuclear magnetic resonance structure of the III-IV-V three-way junction from the Varkud satellite ribozyme and identification of magnesium-binding sites using paramagnetic relaxation enhancement

The VS ribozyme is a catalytic RNA found within some natural isolates of Neurospora that is being used as a model system to improve our understanding of RNA structure, catalysis, and engineering. The catalytic domain contains five helical domains (SLII-SLVI) that are organized by two three-way junctions. The III-IV-V junction is required for high-affinity binding of the substrate domain (SLI) through formation of a kissing loop interaction with SLV. Here, we determine the high-resolution nuclear magnetic resonance (NMR) structure of a 47-nucleotide RNA containing the III-IV-V junction (J345). The J345 RNA adopts a Y-shaped fold typical of the family C three-way junctions, with coaxial stacking between stems III and IV and an acute angle between stems III and V. The NMR structure reveals that the core of the III-IV-V junction contains four stacked base triples, a U-turn motif, a cross-strand stacking interaction, an A-minor interaction, and a ribose zipper. In addition, the NMR structure shows that the cCUUGg tetraloop used to stabilize stem IV adopts a novel RNA tetraloop fold, different from the known gCUUGc tetraloop structure. Using Mn(2+)-induced paramagnetic relaxation enhancement, we identify six Mg(2+)-binding sites within J345, including one associated with the cCUUGg tetraloop and two with the junction core. The NMR structure of J345 likely represents the conformation of the III-IV-V junction in the context of the active VS ribozyme and suggests that this junction functions as a dynamic hinge that contributes to substrate recognition and catalysis. Moreover, this study highlights a new role for family C three-way junctions in long-range tertiary interactions.

Tospovirus ambisense genomic RNA segments use almost complete repertoire of stable tetraloops in the intergenic region

The intergenic regions of the ambisense RNA segments of viruses from the Tospovirus genus form large extended RNA structures that regulate virus replication. Using comparative structure analysis, we show the presence of conserved alternative conformations at the apical parts of these structures. In one conformation, a branched Y-shape, the 5'-proximal hairpin arms are mostly capped by exceptionally stable tetraloop motifs. The tetraloop hairpins are folded in both virus and virus-complementary sense RNAs, and different tetraloops can functionally replace each other. Folding simulations show that the branched Y-shape structures can undergo a conformational transition to alternative extended rod-like conformations. Functional importance of both alternatives is supported by nucleotide covariations. The balanced equilibrium between alternative structures is evidenced by native gel electrophoresis of mutant RNA transcripts with shifted equilibria. The tetraloops play a role in the stability and dynamics of structures but may also be recognized by proteins involved in translation and/or replication.

Counterion and polythymidine loop-length-dependent folding and thermodynamic stability of DNA hairpins reveal the unusual counterion-dependent stability of tetraloop hairpins

Stem-loop DNA hairpins containing a 5-base-pair (bp) stem and single-stranded polythymidine loop were investigated using thermodynamic melting analysis and stopped-flow kinetics. These studies revealed the thermodynamic stability and folding kinetics as a function of loop length and counterion concentration. Our results show the unusually high thermodynamic stability for tetraloop or 4 poly(dT) loop hairpin as compared with longer loop length hairpins. Furthermore, this exceptional stability is highly counterion-dependent. For example, in the higher counterion concentration regime of 50 mM NaCl and above, the tetraloop hairpin displays enhanced stability as compared with longer loop length hairpins. However, at lower counterion concentration of 25 mM NaCl and below, the thermal stability of tetraloop hairpin is consistent with the longer loop hairpins. The enhanced stability of tetraloop hairpins at higher counterion concentration can be explained on the basis of the combined entropic effect of loop closure as well as base stacking in the loop regions. The stability of longer loop length hairpins at all counterion concentrations as well as tetraloop hairpin at lower counterion concentration can be explained on the basis of entropic effect of loop closure alone. The thermodynamic parameters at lower and higher counterion concentrations were determined to quantify the enhanced stability of base-stacking effects occurring at higher counterion concentrations. For example, for 100 mM NaCl, excess Gibbs energy and enthalpy due to base stacking within the tetraloops were measured to be -1.2 ± 0.14 and -3.28 ± 0.32 kcal/mol, respectively, whereas, no excess of Gibbs energy and enthalpy was observed for 0, 5, 10, and 25 mM NaCl. These findings suggest significant base-stacking interactions occurring in the loop region of the tetraloop hairpins at higher counterion concentration and less significant base-stacking interactions in the lower counterion concentration regime. We suggest that at higher counterion concentrations, hydrophobic collapse of the nucleotides in the loop may be enhanced due to the increased polarity of the solvent, thereby enhancing base-stacking interactions that contribute to unusually high stability.

An RNA folding motif: GNRA tetraloop-receptor interactions

Nearly two decades after Westhof and Michel first proposed that RNA tetraloops may interact with distal helices, tetraloop–receptor interactions have been recognized as ubiquitous elements of RNA tertiary structure. The unique architecture of GNRA tetraloops (N=any nucleotide, R=purine) enables interaction with a variety of receptors, e.g., helical minor grooves and asymmetric internal loops. The most common example of the latter is the GAAA tetraloop–11 nt tetraloop receptor motif. Biophysical characterization of this motif provided evidence for the modularity of RNA structure, with applications spanning improved crystallization methods to RNA tectonics. In this review, we identify and compare types of GNRA tetraloop–receptor interactions. Then we explore the abundance of structural, kinetic, and thermodynamic information on the frequently occurring and most widely studied GAAA tetraloop–11 nt receptor motif. Studies of this interaction have revealed powerful paradigms for structural assembly of RNA, as well as providing new insights into the roles of cations, transition states and protein chaperones in RNA folding pathways. However, further research will clearly be necessary to characterize other tetraloop–receptor and long-range tertiary binding interactions in detail – an important milestone in the quantitative prediction of free energy landscapes for RNA folding.

Evaluation of end-capped DNA modules for pRNA capture and functionalization

The ability of packaging RNA (pRNA) from the phi29 DNA packaging motor to form nanoassemblies and nanostructures has been exploited for the development of the nascent field of RNA nanotechnology and subsequent applications in nanomedicine. For applications in nanomedicine, it is necessary to modify the pRNA structure for the conjugation of active molecules. We have investigated end-capped double-stranded DNA segments as reversible capture reagents for pRNA. These capture agents can be designed to allow the conjugation of any desired molecule for pRNA functionalization. The results of model studies presented in this report show that 5- to 7-nucleotide overhangs on a target RNA can provide efficient handles for the high-affinity association to capped double-stranded DNA.

Computational approaches for RNA energy parameter estimation

Methods for efficient and accurate prediction of RNA structure are increasingly valuable, given the current rapid advances in understanding the diverse functions of RNA molecules in the cell. To enhance the accuracy of secondary structure predictions, we developed and refined optimization techniques for the estimation of energy parameters. We build on two previous approaches to RNA free-energy parameter estimation: (1) the Constraint Generation (CG) method, which iteratively generates constraints that enforce known structures to have energies lower than other structures for the same molecule; and (2) the Boltzmann Likelihood (BL) method, which infers a set of RNA free-energy parameters that maximize the conditional likelihood of a set of reference RNA structures. Here, we extend these approaches in two main ways: We propose (1) a max-margin extension of CG, and (2) a novel linear Gaussian Bayesian network that models feature relationships, which effectively makes use of sparse data by sharing statistical strength between parameters. We obtain significant improvements in the accuracy of RNA minimum free-energy pseudoknot-free secondary structure prediction when measured on a comprehensive set of 2518 RNA molecules with reference structures. Our parameters can be used in conjunction with software that predicts RNA secondary structures, RNA hybridization, or ensembles of structures. Our data, software, results, and parameter sets in various formats are freely available at http://www.cs.ubc.ca/labs/beta/Projects/RNA-Params.

Thermodynamic characterization of RNA triloops

Relatively few thermodynamic parameters are available for RNA triloops. Therefore, 24 stem-loop sequences containing naturally occurring triloops were optically melted, and the thermodynamic parameters ΔH°, ΔS°, ΔG°(37), and T(M) for each stem-loop were determined. These new experimental values, on average, are 0.5 kcal/mol different from the values predicted for these triloops using the model proposed by Mathews et al. [Mathews, D. H., Disney, M. D., Childs, J. L., Schroeder, S. J., Zuker, M., and Turner, D. H. (2004) Proc. Natl. Acad. Sci. U.S.A. 101, 7287-7292]. The data for the 24 triloops reported here were then combined with the data for five triloops that were published previously. A new model was derived to predict the free energy contribution of previously unmeasured triloops. The average absolute difference between the measured values and the values predicted using this proposed model is 0.3 kcal/mol. These new experimental data and updated predictive model allow for more accurate calculations of the free energy of RNA stem-loops containing triloops and, furthermore, should allow for improved prediction of secondary structure from sequence.

Thermodynamic characterization of naturally occurring RNA tetraloops

Although tetraloops are one of the most frequently occurring secondary structure motifs in RNA, less than one-third of the 30 most frequently occurring RNA tetraloops have been thermodynamically characterized. Therefore, 24 stem-loop sequences containing common tetraloops were optically melted, and the thermodynamic parameters DeltaH degrees , DeltaS degrees , DeltaG degrees (37,) and T(M) for each stem-loop were determined. These new experimental values, on average, are 0.7 kcal/mol different from the values predicted for these tetraloops using the model proposed by Vecenie CJ, Morrow CV, Zyra A, Serra MJ. 2006. Biochemistry 45: 1400-1407. The data for the 24 tetraloops reported here were then combined with the data for 28 tetraloops that were published previously. A new model, independent of terminal mismatch data, was derived to predict the free energy contribution of previously unmeasured tetraloops. The average absolute difference between the measured values and the values predicted using this proposed model is 0.4 kcal/mol. This new experimental data and updated predictive model allow for more accurate calculations of the free energy of RNA stem-loops containing tetraloops and, furthermore, should allow for improved prediction of secondary structure from sequence. It was also shown that tetraloops within the sequence 5'-GCCNNNNGGC-3' are, on average, 0.6 kcal/mol more stable than the same tetraloop within the sequence 5'-GGCNNNNGCC-3'. More systemic studies are required to determine the full extent of non-nearest-neighbor effects on tetraloop stability.

Mouse hepatitis virus stem-loop 2 adopts a uYNMG(U)a-like tetraloop structure that is highly functionally tolerant of base substitutions

Stem-loop 2 (SL2) of the 5'-untranslated region of the mouse hepatitis virus (MHV) contains a highly conserved pentaloop (C47-U48-U49-G50-U51) stacked on a 5-bp stem. Solution nuclear magnetic resonance experiments are consistent with a 5'-uYNMG(U)a or uCUYG(U)a tetraloop conformation characterized by an anti-C47-syn-G50 base-pairing interaction, with U51 flipped out into solution and G50 stacked on A52. Previous studies showed that U48C and U48A substitutions in MHV SL2 were lethal, while a U48G substitution was viable. Here, we characterize viruses harboring all remaining single-nucleotide substitutions in the pentaloop of MHV SL2 and also investigate the degree to which the sequence context of key pentaloop point mutations influences the MHV replication phenotype. U49 or U51 substitution mutants all are viable; C47 substitution mutants also are viable but produce slightly smaller plaques than wild-type virus. In contrast, G50A and G50C viruses are severely crippled and form much smaller plaques. Virus could not be recovered from G50U-containing mutants; rather, only true wild-type revertants or a virus, G50U/C47A, containing a second site mutation were recovered. These functional data suggest that the Watson-Crick edges of C47 and G50 (or A47 and U50 in the G50U/C47A mutant) are in close enough proximity to a hydrogen bond with U51 flipped out of the hairpin. Remarkably, increasing the helical stem stability rescues the previously lethal mutants U48C and G50U. These studies suggest that SL2 functions as an important, but rather plastic, structural element in stimulating subgenomic RNA synthesis in coronaviruses.

Contribution of the closing base pair to exceptional stability in RNA tetraloops: roles for molecular mimicry and electrostatic factors

Hairpins are common RNA secondary structures that play multiple roles in nature. Tetraloops are the most frequent RNA hairpin loops and are often phylogenetically conserved. For both the UNCG and GNRA families, CG closing base pairs (cbps) confer exceptional thermodynamic stability but the molecular basis for this has remained unclear. We propose that, despite having very different overall folds, these two tetraloop families achieve stability by presenting the same functionalities to the major groove edge of the CG cbp. Thermodynamic contributions of this molecular mimicry were investigated using substitutions at the nucleobase and functional group levels. By either interrupting or deleting loop-cbp electrostatic interactions, which were identified by solving the nonlinear Poisson-Boltzmann (NLPB) equation, stability changed in a manner consistent with molecular mimicry. We also observed a linear relationship between DeltaG(o)(37) and log[Na(+)] for both families, and loops with a CG cbp had a decreased dependence of stability on salt. NLPB calculations revealed that, for both UUCG and GAAA tetraloops, the GC cbp form has a higher surface charge density, although it arises from changes in loop compaction for UUCG and changes in loop configuration for GAAA. Higher surface charge density leads to stronger interactions of GC cbp loops with solvent and salt, which explains the correlation between experimental and calculated trends of free energy with salt. Molecular mimicry as evidenced in these two stable but otherwise unrelated tetraloops may underlie common functional roles in other RNA and DNA motifs.

Thinking inside the box: designing, implementing, and interpreting thermodynamic cycles to dissect cooperativity in RNA and DNA folding

Double and triple mutant thermodynamic cycles provide a means to dissect the cooperativity of RNA and DNA folding at both the secondary and tertiary structural levels through use of the thermodynamic box or cube. In this article, we describe three steps for applying thermodynamic cycles to nucleic acid folding, with considerations of both conceptual and experimental features. The first step is design of an appropriate system and development of hypotheses regarding which residues might interact. Next is implementing this design in terms of a tractable experimental strategy, with an emphasis on UV melting. The final step, and the one we emphasize the most, is interpreting mutant cycles in terms of coupling between specific residues in the RNA or DNA. Coupling free energy in the absence and presence of changes elsewhere in the molecule is discussed in terms of specific folding models, including stepwise folding and concerted changes. Last, we provide a practical section on the use of commercially available software (KaleidaGraph) to fit melting data, along with a consideration of error propagation. Along the way, specific examples are chosen from the literature to illustrate the methods. This article is intended to be accessible to the biochemist or biologist without extensive thermodynamics background.

Structures, kinetics, thermodynamics, and biological functions of RNA hairpins

Most RNA comprises one strand and therefore can fold back on itself to form complex structures. At the heart of these structures is the hairpin, which is composed of a stem having Watson-Crick base pairing and a loop wherein the backbone changes directionality. First, we review the structure of hairpins including diversity in the stem, loop, and closing base pair. The function of RNA hairpins in biology is discussed next, including roles for isolated hairpins, as well as hairpins in the context of complex tertiary structures. We describe the kinetics and thermodynamics of hairpin folding including models for hairpin folding, folding transition states, and the cooperativity of folding. Lastly, we discuss some ways in which hairpins can influence the folding and function of tertiary structures, both directly and indirectly. RNA hairpins provide a simple means of controlling gene expression that can be understood in the language of physical chemistry.

Nature and magnitude of aromatic stacking of nucleic acid bases

This review summarises recent advances in quantum chemical calculations of base-stacking forces in nucleic acids. We explain in detail the very complex relationship between the gas-phase base-stacking energies, as revealed by quantum chemical (QM) calculations, and the highly variable roles of these interactions in nucleic acids. This issue is rarely discussed in quantum chemical and physical chemistry literature. We further extensively discuss methods that are available for base-stacking studies, complexity of comparison of stacking calculations with gas phase experiments, balance of forces in stacked complexes of nucleic acid bases, and the relation between QM and force field descriptions. We also review all recent calculations on base-stacking systems, including details analysis of the B-DNA stacking. Specific attention is paid to the highest accuracy QM calculations, to the decomposition of the interactions, and development of dispersion-balanced DFT methods. Future prospects of computational studies of base stacking are discussed.

A simple ligation-based method to increase the information density in sequencing reactions used to deconvolute nucleic acid selections

Herein, a method is described to increase the information density of sequencing experiments used to deconvolute nucleic acid selections. The method is facile and should be applicable to any selection experiment. A critical feature of this method is the use of biotinylated primers to amplify and encode a BamHI restriction site on both ends of a PCR product. After amplification, the PCR reaction is captured onto streptavidin resin, washed, and digested directly on the resin. Resin-based digestion affords clean product that is devoid of partially digested products and unincorporated PCR primers. The product's complementary ends are annealed and ligated together with T4 DNA ligase. Analysis of ligation products shows formation of concatemers of different length and little detectable monomer. Sequencing results produced data that routinely contained three to four copies of the library. This method allows for more efficient formulation of structure-activity relationships since multiple active sequences are identified from a single clone.

Using selection to identify and chemical microarray to study the RNA internal loops recognized by 6'-N-acylated kanamycin A

Herein, we describe our initial steps towards identifying the RNA secondary structure motifs that are recognized by small molecules. We selected members of an RNA 3x3 internal loop motif library that bind kanamycin A, an RNA-binding aminoglycoside antibiotic, by using only one round of selection. A small internal-loop library was chosen because members are likely to be present in other larger, biologically relevant RNAs. We have identified several internal loops of various size and base composition that kanamycin A prefers to bind. The highest affinity structures are two 5'-UU/3'-CU 2x2 internal loops closed by AU pairs. Binding is specific for the selected internal loops with the highest affinities, since binding to the RNA cassette used to display the library or to DNA is >150-fold weaker. Enzymatic mapping experiments also confirm binding of kanamycin A to the internal loops. This method lays the foundation for finding RNA secondary structure elements that bind small molecules and for interrogating factors affecting RNA-ligand interactions. Information from these and subsequent studies will: 1) facilitate the rational and modular design of drugs or probes that bind target RNAs with high affinity, provided the secondary structure of the target is known and 2) give insight into the potential bystander RNAs that aminoglycosides bind.

Distant Segments of Saccharomyces cerevisiae scR1 RNA Promote Assembly and Function of the Signal Recognition Particle

The conserved signal recognition particle targets ribosomes synthesizing presecretory proteins to the endoplasmic reticulum membrane. Key to the activity of SRP is its ability to bind the ribosome at distant locations, the signal sequence exit and elongation factor-binding sites. These contacts are made by the S and Alu domains of SRP, respectively. We tested earlier secondary structure predictions of the Saccharomyces cerevisiae SRP RNA, scR1, and provide and test a consensus structure. The structure contains four non-conserved insertions, helices 9-12, into the core SRP RNA fold, and an extended helix 7. Using a series of scR1 mutants lacking part or all of these structural elements, we find that they are important for the RNA in both function and assembly of the RNP. About 20% of the RNA, corresponding to the outer regions of these helices, is dispensable for function. Further, we examined the role of several features within the S-domain section of the core, helix 5, and find that its length and flexibility are important for proper SRP function and become essential in the absence of helix 10, 11 and/or 7 regions. Overall, the genetic data indicate that regions of scR1 distant in both primary sequence and secondary structure have interrelated roles in the function of the complex, and possibly mediate communication between Alu and S domains during targeting.

A set of nearest neighbor parameters for predicting the enthalpy change of RNA secondary structure formation

A complete set of nearest neighbor parameters to predict the enthalpy change of RNA secondary structure formation was derived. These parameters can be used with available free energy nearest neighbor parameters to extend the secondary structure prediction of RNA sequences to temperatures other than 37°C. The parameters were tested by predicting the secondary structures of sequences with known secondary structure that are from organisms with known optimal growth temperatures. Compared with the previous set of enthalpy nearest neighbor parameters, the sensitivity of base pair prediction improved from 65.2 to 68.9% at optimal growth temperatures ranging from 10 to 60°C. Base pair probabilities were predicted with a partition function and the positive predictive value of structure prediction is 90.4% when considering the base pairs in the lowest free energy structure with pairing probability of 0.99 or above. Moreover, a strong correlation is found between the predicted melting temperatures of RNA sequences and the optimal growth temperatures of the host organism. This indicates that organisms that live at higher temperatures have evolved RNA sequences with higher melting temperatures.

A GCUA tetranucleotide loop found in the poliovirus oriL by in vivo SELEX (un)expectedly forms a YNMG-like structure: Extending the YNMG family with GYYA

The cloverleaf structure in the 5'-untranslated region of enterovirus RNA that regulates viral RNA replication contains an evolutionarily conserved YNMG tetraloop closed by a Y-G base pair. This loop is believed to interact specifically with the viral protease 3C. To further characterize the specificity of this interaction, the tetraloop and two flanking base pairs of the poliovirus RNA were randomized, and viable viral clones were obtained using in vivo SELEX. Among many different mutants with the canonical YNMG sequences to be described elsewhere, a large-plaque-forming clone contained a deviating uGCUAg sequence. The NMR structure of a small hairpin capped with uGCUAg that we present here shows that the GCUA tetraloop adopts a novel fold, which is highly similar to that of the YNMG tetraloop with common stacking properties and hydrogen-bond interactions including an unusual syn conformation of the adenosine. Thermodynamic studies show moderate stabilities of hairpins with canonical YNMG and the novel GCUA loops, which, together with the similarity of spatial structures, illustrates that the tetraloop structure itself is crucial for the RNA-protein interaction required for the viral replication. A re-evaluation of the ribosomal secondary structure database reveals a hairpin containing a GCUA loop, which covaries with YNMG and is involved in a tertiary interaction, and in the 50S ribosomal subunit from Haloarcula marismortui the structurally comparable apex of stem-loop 35a is a recognition site for protein L2. These observations show a more general occurrence and importance of the so-far unrecognized GYYA hairpin loops.

Structural parameters affecting the kinetics of RNA hairpin formation

There is little experimental knowledge on the sequence dependent rate of hairpin formation in RNA. We have therefore designed RNA sequences that can fold into either of two mutually exclusive hairpins and have determined the ratio of folding of the two conformations, using structure probing. This folding ratio reflects their respective folding rates. Changing one of the two loop sequences from a purine- to a pyrimidine-rich loop did increase its folding rate, which corresponds well with similar observations in DNA hairpins. However, neither changing one of the loops from a regular non-GNRA tetra-loop into a stable GNRA tetra-loop, nor increasing the loop size from 4 to 6 nt did affect the folding rate. The folding kinetics of these RNAs have also been simulated with the program ‘Kinfold’. These simulations were in agreement with the experimental results if the additional stabilization energies for stable tetra-loops were not taken into account. Despite the high stability of the stable tetra-loops, they apparently do not affect folding kinetics of these RNA hairpins. These results show that it is possible to experimentally determine relative folding rates of hairpins and to use these data to improve the computer-assisted simulation of the folding kinetics of stem–loop structures.

Exploring the energy landscape of a small RNA hairpin

The energy landscape of a small RNA tetraloop hairpin is explored by temperature jump kinetics and base-substitution. The folding kinetics are single-exponential near the folding transition midpoint T(m). An additional fast phase appears below the midpoint, and an additional slow phase appears above the midpoint. Stem mutation affects the high-temperature phase, while loop mutation affects the low-temperature phase. An adjusted 2-D lattice model reproduces the temperature-dependent phases, although it oversimplifies the structural interpretation. A four-state free energy landscape model is generated based on the lattice model. This model explains the thermodynamics and multiphase kinetics over the full temperature range of the experiments. An analysis of three variants shows that one of the intermediate RNA structures is a stacking-related trap affected by stem but not loop modification, while the other is an early intermediate that forms some stem and loop structure. Even a very fast-folding 8-mer RNA with an ideal tetraloop sequence has a rugged energy landscape, ideal for testing analytical and computational models.

Transcription of potato spindle tuber viroid by RNA polymerase II starts in the left terminal loop

Viroids are single-stranded, circular RNAs of 250 to 400 bases, that replicate autonomously in their host plants but do not code for a protein. Viroids of the family Pospiviroidae, of which potato spindle tuber viroid (PSTVd) is the type strain, are replicated by the host's DNA-dependent RNA polymerase II in the nucleus. To analyze the initiation site of transcription from the (+)-stranded circles into (-)-stranded replication intermediates, we used a nuclear extract from a non-infected cell culture of the host plant S. tuberosum. The (-)-strands, which were de novo-synthesized in the extract upon addition of circular (+)-PSTVd, were purified by affinity chromatography. This purification avoided contamination by host nucleic acids that had resulted in a misassignment of the start site in an earlier study. Primer-extension analysis of the de novo-synthesized (-)-strands revealed a single start site located in the hairpin loop of the left terminal region in circular PSTVd's secondary structure. This start site is supported further by analysis of the infectivity and replication behavior of site-directed mutants in planta.

In vitro selection of halo-thermophilic RNA reveals two families of resistant RNA

The "RNA world" hypothesis proposes that early in the evolution of life, RNA was responsible both for the storage and transfer of genetic information and for the catalysis of biochemical reactions. One of the problems of the hypothesis is that RNA is known to be temperature sensitive. Nevertheless, different types of sequences with a thermostable phenotype may exist. In order to test this possibility, we applied an in vitro evolution method (SELEX) to isolate RNA molecules that are resistant at high temperatures (80 degrees C for 65 h) and high salt concentrations (2 M NaCl). The sequences of the resulting cloned halo-thermophilic RNAs can be grouped in two families (I and II) possessing very different thermal and chemical stabilities and very different secondary structures. The selected RNA molecules illustrate two different possibilities leading to thermal resistance which may be related to primitive conditions. We propose that members of family I constitute a good means of storing sequence information while members of family II are less efficient but replicate faster in early steps of the SELEX. These selected RNA behaviors may be related to primitive conditions and could allow to define limits for survival, and demonstrate that what is at stake for RNA molecules, as for living organisms, is survival and reproduction.