Minor groove RNA triplex in the crystal structure of a ribosomal frameshifting viral pseudoknot

Antisense RNA C9orf72 hexanucleotide repeat associated with amyotrophic lateral sclerosis and frontotemporal dementia forms a triplex-like structure and binds small synthetic ligand

The abnormal expansion of GGGGCC/GGCCCC hexanucleotide repeats (HR) in C9orf72 is associated with amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Structural polymorphisms of HR result in the multifactorial pathomechanism of ALS/FTD. Consequently, many ongoing studies are focused at developing therapies targeting pathogenic HR RNA. One of them involves small molecules blocking sequestration of important proteins, preventing formation of toxic nuclear foci. However, rational design of potential therapeutics is hindered by limited number of structural studies of RNA-ligand complexes. We determined the crystal structure of antisense HR RNA in complex with ANP77 ligand (1.1 Å resolution) and in the free form (0.92 and 1.5 Å resolution). HR RNA folds into a triplex structure composed of four RNA chains. ANP77 interacted with two neighboring single-stranded cytosines to form pseudo-canonical base pairs by adopting sandwich-like conformation and adjusting the position of its naphthyridine units to the helical twist of the RNA. In the unliganded structure, the cytosines formed a peculiar triplex i-motif, assembled by trans C•C+ pair and a third cytosine located at the Hoogsteen edge of the C•C+ pair. These results extend our knowledge of the structural polymorphisms of HR and can be used for rational design of small molecules targeting disease-related RNAs.

Assessment of Nucleobase Protomeric and Tautomeric States in Nucleic Acid Structures for Interaction Analysis and Structure-Based Ligand Design

With increasing interest in RNA as a therapeutic and a potential target, the role of RNA structures has become more important. Even slight changes in nucleobases, such as modifications or protomeric and tautomeric states, can have a large impact on RNA structure and function, while local environments in turn affect protonation and tautomerization. In this work, the application of empirical tools for pKa and tautomer prediction for RNA modifications was elucidated and compared with ab initio quantum mechanics (QM) methods and expanded toward macromolecular RNA structures, where QM is no longer feasible. In this regard, the Protonate3D functionality within the molecular operating environment (MOE) was expanded for nucleobase protomer and tautomer predictions and applied to reported examples of altered protonation states depending on the local environment. Overall, observations of nonstandard protomers and tautomers were well reproduced, including structural C+G:C(A) and A+GG motifs, several mismatches, and protonation of adenosine or cytidine as the general acid in nucleolytic ribozymes. Special cases, such as cobalt hexamine-soaked complexes or the deprotonation of guanosine as the general base in nucleolytic ribozymes, proved to be challenging. The collected set of examples shall serve as a starting point for the development of further RNA protonation prediction tools, while the presented Protonate3D implementation already delivers reasonable protonation predictions for RNA and DNA macromolecules. For cases where higher accuracy is needed, like following catalytic pathways of ribozymes, incorporation of QM-based methods can build upon the Protonate3D-generated starting structures. Likewise, this protonation prediction can be used for structure-based RNA-ligand design approaches.

Structural and Functional Insights into Viral Programmed Ribosomal Frameshifting

Protein synthesis by the ribosome is the final stage of biological information transfer and represents an irreversible commitment to gene expression. Accurate translation of messenger RNA is therefore essential to all life, and spontaneous errors by the translational machinery are highly infrequent (∼1/100,000 codons). Programmed -1 ribosomal frameshifting (-1PRF) is a mechanism in which the elongating ribosome is induced at high frequency to slip backward by one nucleotide at a defined position and to continue translation in the new reading frame. This is exploited as a translational regulation strategy by hundreds of RNA viruses, which rely on -1PRF during genome translation to control the stoichiometry of viral proteins. While early investigations of -1PRF focused on virological and biochemical aspects, the application of X-ray crystallography and cryo-electron microscopy (cryo-EM), and the advent of deep sequencing and single-molecule approaches have revealed unexpected structural diversity and mechanistic complexity. Molecular players from several model systems have now been characterized in detail, both in isolation and, more recently, in the context of the elongating ribosome. Here we provide a summary of recent advances and discuss to what extent a general model for -1PRF remains a useful way of thinking.

Graded impact of obstacle size on scanning by RNase E

In countless bacterial species, the lifetimes of most mRNAs are controlled by the regulatory endonuclease RNase E, which preferentially degrades RNAs bearing a 5' monophosphate and locates cleavage sites within them by scanning linearly from the 5' terminus along single-stranded regions. Consequently, its rate of cleavage at distal sites is governed by any obstacles that it may encounter along the way, such as bound proteins or ribosomes or base pairing that is coaxial with the path traversed by this enzyme. Here, we report that the protection afforded by such obstacles is dependent on the size and persistence of the structural discontinuities they create, whereas the molecular composition of obstacles to scanning is of comparatively little consequence. Over a broad range of sizes, incrementally larger discontinuities are incrementally more protective, with corresponding effects on mRNA stability. The graded impact of such obstacles suggests possible explanations for why their effect on scanning is not an all-or-none phenomenon dependent simply on whether the size of the resulting discontinuity exceeds the step length of RNase E.

Dynamic effects of the spine of hydrated magnesium on viral RNA pseudoknot structure

In the cellular environment, a viral RNA Pseudoknot (PK) structure is responsive to its prevailing ion atmosphere created by a mixture of monovalent and divalent cations. We investigate the influence of such a mixed-salt environment on RNA-PK structure at an atomic resolution through three sets of 1.5 μs explicit solvent molecular dynamics (MD) simulations and also by building a dynamic counterion-condensation (DCC) model at varying divalent Mg²⁺ concentrations. The DCC model includes explicit interaction of the ligand and adjacent chelated Mg²⁺ by extending the recently developed generalized Manning condensation model. Its implementation within an all-atom structure-based molecular dynamics framework bolsters its opportunity to explore large-length scale and long-timescale phenomena associated with complex RNA systems immersed in its dynamic ion environment. In the present case of RNA-PK, both explicit MD and DCC simulations reveal a spine of hydrated-Mg²⁺ to induce stem-I and stem-II closure where the minor groove between these stems is akin to breathing. Mg²⁺ mediated minor-groove narrowing is coupled with local base-flipping dynamics of a base triple and quadruple, changing the stem structure of such RNA. Contrary to the conversational view of the indispensable need for Mg²⁺ for the tertiary structure of RNA, the study warns about the plausible detrimental effect of specific Mg²⁺-phosphate interactions on the RNA-PK structure beyond a certain concentration of Mg²⁺.

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.

Cooperativity and Interdependency between RNA Structure and RNA-RNA Interactions

Complex RNA–RNA interactions are increasingly known to play key roles in numerous biological processes from gene expression control to ribonucleoprotein granule formation. By contrast, the nature of these interactions and characteristics of their interfaces, especially those that involve partially or wholly structured RNAs, remain elusive. Herein, we discuss different modalities of RNA–RNA interactions with an emphasis on those that depend on secondary, tertiary, or quaternary structure. We dissect recently structurally elucidated RNA–RNA complexes including RNA triplexes, riboswitches, ribozymes, and reverse transcription complexes. These analyses highlight a reciprocal relationship that intimately links RNA structure formation with RNA–RNA interactions. The interactions not only shape and sculpt RNA structures but also are enabled and modulated by the structures they create. Understanding this two-way relationship between RNA structure and interactions provides mechanistic insights into the expanding repertoire of noncoding RNA functions, and may inform the design of novel therapeutics that target RNA structures or interactions.

Graph Theoretical Methods and Workflows for Searching and Annotation of RNA Tertiary Base Motifs and Substructures

The increasing number and complexity of structures containing RNA chains in the Protein Data Bank (PDB) have led to the need for automated structure annotation methods to replace or complement expert visual curation. This is especially true when searching for tertiary base motifs and substructures. Such base arrangements and motifs have diverse roles that range from contributions to structural stability to more direct involvement in the molecule's functions, such as the sites for ligand binding and catalytic activity. We review the utility of computational approaches in annotating RNA tertiary base motifs in a dataset of PDB structures, particularly the use of graph theoretical algorithms that can search for such base motifs and annotate them or find and annotate clusters of hydrogen-bond-connected bases. We also demonstrate how such graph theoretical algorithms can be integrated into a workflow that allows for functional analysis and comparisons of base arrangements and sub-structures, such as those involved in ligand binding. The capacity to carry out such automatic curations has led to the discovery of novel motifs and can give new context to known motifs as well as enable the rapid compilation of RNA 3D motifs into a database.

Tertiary Base Triple Formation in the SRV-1 Frameshifting Pseudoknot Stabilizes Secondary Structure Components

Minor-groove base triples formed between stem 1 and loop 2 of the simian retrovirus type 1 (SRV-1) mRNA frameshifting pseudoknot are essential in stimulating -1 ribosomal frameshifting. How tertiary base triple formation affects the local stabilities of secondary structures (stem 1 and stem 2) and thus ribosomal frameshifting efficiency is not well understood. We made a short peptide nucleic acid (PNA) that is expected to invade stem 1 of the SRV-1 pseudoknot by PNA-RNA duplex formation to mimic the stem 1 unwinding process by a translating ribosome. In addition, we used a PNA for invading stem 2 in the SRV-1 pseudoknot. Our nondenaturing polyacrylamide gel electrophoresis data for the binding of PNA to the SRV-1 pseudoknot and mutants reveal that mutations in loop 2 disrupting base triple formation between loop 2 and stem 1 in the SRV-1 pseudoknot result in enhanced invasion by both PNAs. Our data suggest that tertiary stem 1-loop 2 base triple interactions in the SRV-1 pseudoknot can stabilize both of the secondary structural components, stem 1 and stem 2. Stem 2 stability is thus coupled to the structural stability of stem 1-loop 2 base triples, mediated through a long-range effect. The apparent dissociation constants of both PNAs are positively correlated with the pseudoknot mechanical stabilities and frameshifting efficiencies. The relatively simple PNA local invasion experiment may be used to characterize the energetic contribution of tertiary interactions and ligand binding in many other RNA and DNA structures.

The molecular structure of long non-coding RNAs: emerging patterns and functional implications

Long non-coding RNAs (lncRNAs) are recently-discovered transcripts that regulate vital cellular processes and are crucially connected to diseases. Despite their unprecedented molecular complexity, it is emerging that lncRNAs possess distinct structural motifs. Remarkably, the 3D shape and topology of full-length, native lncRNAs have been visualized for the first time in the last year. These studies reveal that lncRNA structures dictate lncRNA functions. Here, we review experimentally determined lncRNA structures and emphasize that lncRNA structural characterization requires synergistic integration of computational, biochemical and biophysical approaches. Based on these emerging paradigms, we discuss how to overcome the challenges posed by the complex molecular architecture of lncRNAs, with the goal of obtaining a detailed understanding of lncRNA functions and molecular mechanisms in the future.

Unraveling the structure and biological functions of RNA triple helices

It has been nearly 63 years since the first characterization of an RNA triple helix in vitro by Gary Felsenfeld, David Davies, and Alexander Rich. An RNA triple helix consists of three strands: A Watson–Crick RNA double helix whose major‐groove establishes hydrogen bonds with the so‐called “third strand”. In the past 15 years, it has been recognized that these major‐groove RNA triple helices, like single‐stranded and double‐stranded RNA, also mediate prominent biological roles inside cells. Thus far, these triple helices are known to mediate catalysis during telomere synthesis and RNA splicing, bind to ligands and ions so that metabolite‐sensing riboswitches can regulate gene expression, and provide a clever strategy to protect the 3′ end of RNA from degradation. Because RNA triple helices play important roles in biology, there is a renewed interest in better understanding the fundamental properties of RNA triple helices and developing methods for their high‐throughput discovery. This review provides an overview of the fundamental biochemical and structural properties of major‐groove RNA triple helices, summarizes the structure and function of naturally occurring RNA triple helices, and describes prospective strategies to isolate RNA triple helices as a means to establish the “triplexome”.

This article is categorized under:

RNA Structure and Dynamics > RNA Structure and Dynamics

RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry

RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems

Alternate RNA Structures

RNA molecules fold into complex three-dimensional structures that sample alternate conformations ranging from minor differences in tertiary structure dynamics to major differences in secondary structure. This allows them to form entirely different substructures with each population potentially giving rise to a distinct biological outcome. The substructures can be partitioned along an existing energy landscape given a particular static cellular cue or can be shifted in response to dynamic cues such as ligand binding. We review a few key examples of RNA molecules that sample alternate conformations and how these are capitalized on for control of critical regulatory functions.

Competitive Binding of Mg2+ and Na+ Ions to Nucleic Acids: From Helices to Tertiary Structures

Nucleic acids generally reside in cellular aqueous solutions with mixed divalent/monovalent ions, and the competitive binding of divalent and monovalent ions is critical to the structures of nucleic acids because of their polyanionic nature. In this work, we first proposed a general and effective method for simulating a nucleic acid in mixed divalent/monovalent ion solutions with desired bulk ion concentrations via molecular dynamics (MD) simulations and investigated the competitive binding of Mg²⁺/Na⁺ ions to various nucleic acids by all-atom MD simulations. The extensive MD-based examinations show that single MD simulations conducted using the proposed method can yield desired bulk divalent/monovalent ion concentrations for various nucleic acids, including RNA tertiary structures. Our comprehensive analyses show that the global binding of Mg²⁺/Na⁺ to a nucleic acid is mainly dependent on its structure compactness, as well as Mg²⁺/Na⁺ concentrations, rather than the specific structure of the nucleic acid. Specifically, the relative global binding of Mg²⁺ over Na⁺ is stronger for a nucleic acid with higher effective surface charge density and higher relative Mg²⁺/Na⁺ concentrations. Furthermore, the local binding of Mg²⁺/Na⁺ to a phosphate of a nucleic acid mainly depends on the local phosphate density in addition to Mg²⁺/Na⁺ concentrations.

jViz.RNA 4.0-Visualizing pseudoknots and RNA editing employing compressed tree graphs

Previously, we have introduced an improved version of jViz.RNA which enabled faster and more stable RNA visualization by employing compressed tree graphs. However, the new RNA representation and visualization method required a sophisticated mechanism of pseudoknot visualization. In this work, we present our novel pseudoknot classification and implementation of pseudoknot visualization in the context of the new RNA graph model. We then compare our approach with other RNA visualization software, and demonstrate jViz.RNA 4.0's benefits compared to other software. Additionally, we introduce interactive editing functionality into jViz.RNA and demonstrate its benefits in exploring and building RNA structures. The results presented highlight the new high degree of utility jViz.RNA 4.0 now offers. Users are now able to visualize pseudoknotted RNA, manipulate the resulting automatic layouts to suit their individual needs, and change both positioning and connectivity of the RNA molecules examined. Care was taken to limit overlap between structural elements, particularly in the case of pseudoknots to ensure an intuitive and informative layout of the final RNA structure. Availability: The software is freely available at: https://jviz.cs.sfu.ca/.

Hexahydrated Mg2+ Binding and Outer-Shell Dehydration on RNA Surface

The interaction between metal ions, especially Mg²⁺ ions, and RNA plays a critical role in RNA folding. Upon binding to RNA, a metal ion that is fully hydrated in bulk solvent can become dehydrated. Here we use molecular dynamics simulation to investigate the dehydration of bound hexahydrated Mg²⁺ ions. We find that a hydrated Mg²⁺ ion in the RNA groove region can involve significant dehydration in the outer hydration shell. The first or innermost hydration shell of the Mg²⁺ ion, however, is retained during the simulation because of the strong ion-water electrostatic attraction. As a result, water-mediated hydrogen bonding remains an important form for Mg²⁺-RNA interaction. Analysis for ions at different binding sites shows that the most pronounced water deficiency relative to the fully hydrated state occurs at a radial distance of around 11 Å from the center of the ion. Based on the independent 200 ns molecular dynamics simulations for three different RNA structures (Protein Data Bank: 1TRA, 2TPK, and 437D), we find that Mg²⁺ ions overwhelmingly dominate over monovalent ions such as Na⁺ and K⁺ in ion-RNA binding. Furthermore, application of the free energy perturbation method leads to a quantitative relationship between the Mg²⁺ dehydration free energy and the local structural environment. We find that ΔΔG_hyd, the change of the Mg²⁺ hydration free energy upon binding to RNA, varies linearly with the inverse distance between the Mg²⁺ ion and the nearby nonbridging oxygen atoms of the phosphate groups, and ΔΔG_hyd can reach -2.0 kcal/mol and -3.0 kcal/mol for an Mg²⁺ ion bound to the surface and to the groove interior, respectively. In addition, the computation results in an analytical formula for the hydration ratio as a function of the average inverse Mg²⁺-O distance. The results here might be useful for further quantitative investigations of ion-RNA interactions in RNA folding.

Frictional Effects on RNA Folding: Speed Limit and Kramers Turnover

We investigated frictional effects on the folding rates of a human telomerase hairpin (hTR HP) and H-type pseudoknot from the Beet Western Yellow Virus (BWYV PK) using simulations of the Three Interaction Site (TIS) model for RNA. The heat capacity from TIS model simulations, calculated using temperature replica exchange simulations, reproduces nearly quantitatively the available experimental data for the hTR HP. The corresponding results for BWYV PK serve as predictions. We calculated the folding rates ( k_F) from more than 100 folding trajectories for each value of the solvent viscosity (η) at a fixed salt concentration of 200 mM. By using the theoretical estimate (∝ √N where N is the number of nucleotides) for folding free energy barrier, k_F data for both the RNAs are quantitatively fit using one-dimensional Kramers's theory with two parameters specifying the curvatures in the unfolded basin and the barrier top. In the high-friction regime (η ≳ 10^-5 Pa·s), for both HP and PK, k_F values decrease as 1/η, whereas in the low friction regime, k_F values increase as η increases, leading to a maximum folding rate at a moderate viscosity (∼10^-6 Pa·s), which is the Kramers turnover. From the fits, we find that the speed limit to RNA folding at water viscosity is between 1 and 4 μs, which is in accord with our previous theoretical prediction as well as results from several single molecule experiments. Both the RNA constructs fold by parallel pathways. Surprisingly, we find that the flux through the pathways could be altered by changing solvent viscosity, a prediction that is more easily testable in RNA than in proteins.

Predicting 3D structure and stability of RNA pseudoknots in monovalent and divalent ion solutions

RNA pseudoknots are a kind of minimal RNA tertiary structural motifs, and their three-dimensional (3D) structures and stability play essential roles in a variety of biological functions. Therefore, to predict 3D structures and stability of RNA pseudoknots is essential for understanding their functions. In the work, we employed our previously developed coarse-grained model with implicit salt to make extensive predictions and comprehensive analyses on the 3D structures and stability for RNA pseudoknots in monovalent/divalent ion solutions. The comparisons with available experimental data show that our model can successfully predict the 3D structures of RNA pseudoknots from their sequences, and can also make reliable predictions for the stability of RNA pseudoknots with different lengths and sequences over a wide range of monovalent/divalent ion concentrations. Furthermore, we made comprehensive analyses on the unfolding pathway for various RNA pseudoknots in ion solutions. Our analyses for extensive pseudokonts and the wide range of monovalent/divalent ion concentrations verify that the unfolding pathway of RNA pseudoknots is mainly dependent on the relative stability of unfolded intermediate states, and show that the unfolding pathway of RNA pseudoknots can be significantly modulated by their sequences and solution ion conditions.

Know Your Enemy: Successful Bioinformatic Approaches to Predict Functional RNA Structures in Viral RNAs

Structured RNA elements may control virus replication, transcription and translation, and their distinct features are being exploited by novel antiviral strategies. Viral RNA elements continue to be discovered using combinations of experimental and computational analyses. However, the wealth of sequence data, notably from deep viral RNA sequencing, viromes, and metagenomes, necessitates computational approaches being used as an essential discovery tool. In this review, we describe practical approaches being used to discover functional RNA elements in viral genomes. In addition to success stories in new and emerging viruses, these approaches have revealed some surprising new features of well-studied viruses e.g., human immunodeficiency virus, hepatitis C virus, influenza, and dengue viruses. Some notable discoveries were facilitated by new comparative analyses of diverse viral genome alignments. Importantly, comparative approaches for finding RNA elements embedded in coding and non-coding regions differ. With the exponential growth of computer power we have progressed from stem-loop prediction on single sequences to cutting edge 3D prediction, and from command line to user friendly web interfaces. Despite these advances, many powerful, user friendly prediction tools and resources are underutilized by the virology community.

Ensemble simulations: folding, unfolding and misfolding of a high-efficiency frameshifting RNA pseudoknot

Massive all-atom molecular dynamics simulations were conducted across a distributed computing network to study the folding, unfolding, misfolding and conformational plasticity of the high-efficiency frameshifting double mutant of the 26 nt potato leaf roll virus RNA pseudoknot. Our robust sampling, which included over 40 starting structures spanning the spectrum from the extended unfolded state to the native fold, yielded nearly 120 μs of cumulative sampling time. Conformational microstate transitions on the 1.0 ns to 10.0 μs timescales were observed, with post-equilibration sampling providing detailed representations of the conformational free energy landscape and the complex folding mechanism inherent to the pseudoknot motif. Herein, we identify and characterize two alternative native structures, three intermediate states, and numerous misfolded states, the latter of which have not previously been characterized via atomistic simulation techniques. While in line with previous thermodynamics-based models of a general RNA folding mechanism, our observations indicate that stem-strand-sequence-separation may serve as an alternative predictor of the order of stem formation during pseudoknot folding. Our results contradict a model of frameshifting based on structural rigidity and resistance to mechanical unfolding, and instead strongly support more recent studies in which conformational plasticity is identified as a determining factor in frameshifting efficiency.

Coordination among tertiary base pairs results in an efficient frameshift-stimulating RNA pseudoknot

Frameshifting is an essential process that regulates protein synthesis in many viruses. The ribosome may slip backward when encountering a frameshift motif on the messenger RNA, which usually contains a pseudoknot structure involving tertiary base pair interactions. Due to the lack of detailed molecular explanations, previous studies investigating which features of the pseudoknot are important to stimulate frameshifting have presented diverse conclusions. Here we constructed a bimolecular pseudoknot to dissect the interior tertiary base pairs and used single-molecule approaches to assess the structure targeted by ribosomes. We found that the first ribosome target stem was resistant to unwinding when the neighboring loop was confined along the stem; such constrained conformation was dependent on the presence of consecutive adenosines in this loop. Mutations that disrupted the distal base triples abolished all remaining tertiary base pairs. Changes in frameshifting efficiency correlated with the stem unwinding resistance. Our results demonstrate that various tertiary base pairs are coordinated inside a highly efficient frameshift-stimulating RNA pseudoknot and suggest a mechanism by which mechanical resistance of the pseudoknot may persistently act on translocating ribosomes.

Detection and discrimination of maintenance and de novo CpG methylation events using MethylBreak

Understanding the principles governing the establishment and maintenance activities of DNA methyltransferases (DNMTs) can help in the development of predictive biomarkers associated with genetic disorders and diseases. A detection system was developed that distinguishes and quantifies methylation events using methylation-sensitive endonucleases and molecular beacon technology. MethylBreak (MB) is a 22-mer oligonucleotide with one hemimethylated and two unmethylated CpG sites, which are also recognition sites for Sau96I and SacII, and is attached to a fluorophore and a quencher. Maintenance methylation was quantified by fluorescence emission due to the digestion of SacII when the hemimethylated CpG site is methylated, which inhibits Sau96I cleavage. The signal difference between SacII digestion of both MB substrate and maintenance methylated MB corresponds to de novo methylation event. Our technology successfully discriminated and measured both methylation activities at different concentrations of MB and achieved a high correlation coefficient of R²=0.997. Additionally, MB was effectively applied to normal and cancer cell lines and in the analysis of enzymatic kinetics and RNA inhibition of recombinant human DNMT1.

Structure and folding of the Tetrahymena telomerase RNA pseudoknot

Telomerase maintains telomere length at the ends of linear chromosomes using an integral telomerase RNA (TER) and telomerase reverse transcriptase (TERT). An essential part of TER is the template/pseudoknot domain (t/PK) which includes the template, for adding telomeric repeats, template boundary element (TBE), and pseudoknot, enclosed in a circle by stem 1. The Tetrahymena telomerase holoenzyme catalytic core (p65-TER-TERT) was recently modeled in our 9 Å resolution cryo-electron microscopy map by fitting protein and TER domains, including a solution NMR structure of the Tetrahymena pseudoknot. Here, we describe in detail the structure and folding of the isolated pseudoknot, which forms a compact structure with major groove U•A-U and novel C•G-A⁺ base triples. Base substitutions that disrupt the base triples reduce telomerase activity in vitro. NMR studies also reveal that the pseudoknot does not form in the context of full-length TER in the absence of TERT, due to formation of a competing structure that sequesters pseudoknot residues. The residues around the TBE remain unpaired, potentially providing access by TERT to this high affinity binding site during an early step in TERT-TER assembly. A model for the assembly pathway of the catalytic core is proposed.

Ribosomal frameshifting and transcriptional slippage: From genetic steganography and cryptography to adventitious use

Genetic decoding is not ‘frozen’ as was earlier thought, but dynamic. One facet of this is frameshifting that often results in synthesis of a C-terminal region encoded by a new frame. Ribosomal frameshifting is utilized for the synthesis of additional products, for regulatory purposes and for translational ‘correction’ of problem or ‘savior’ indels. Utilization for synthesis of additional products occurs prominently in the decoding of mobile chromosomal element and viral genomes. One class of regulatory frameshifting of stable chromosomal genes governs cellular polyamine levels from yeasts to humans. In many cases of productively utilized frameshifting, the proportion of ribosomes that frameshift at a shift-prone site is enhanced by specific nascent peptide or mRNA context features. Such mRNA signals, which can be 5′ or 3′ of the shift site or both, can act by pairing with ribosomal RNA or as stem loops or pseudoknots even with one component being 4 kb 3′ from the shift site. Transcriptional realignment at slippage-prone sequences also generates productively utilized products encoded trans-frame with respect to the genomic sequence. This too can be enhanced by nucleic acid structure. Together with dynamic codon redefinition, frameshifting is one of the forms of recoding that enriches gene expression.

Salt Effects on the Thermodynamics of a Frameshifting RNA Pseudoknot under Tension

Because of the potential link between -1 programmed ribosomal frameshifting and response of a pseudoknot (PK) RNA to force, a number of single-molecule pulling experiments have been performed on PKs to decipher the mechanism of programmed ribosomal frameshifting. Motivated in part by these experiments, we performed simulations using a coarse-grained model of RNA to describe the response of a PK over a range of mechanical forces (fs) and monovalent salt concentrations (Cs). The coarse-grained simulations quantitatively reproduce the multistep thermal melting observed in experiments, thus validating our model. The free energy changes obtained in simulations are in excellent agreement with experiments. By varying f and C, we calculated the phase diagram that shows a sequence of structural transitions, populating distinct intermediate states. As f and C are changed, the stem-loop tertiary interactions rupture first, followed by unfolding of the 3'-end hairpin (I⇌F). Finally, the 5'-end hairpin unravels, producing an extended state (E⇌I). A theoretical analysis of the phase boundaries shows that the critical force for rupture scales as (logCm)(α) with α=1(0.5) for E⇌I (I⇌F) transition. This relation is used to obtain the preferential ion-RNA interaction coefficient, which can be quantitatively measured in single-molecule experiments, as done previously for DNA hairpins. A by-product of our work is the suggestion that the frameshift efficiency is likely determined by the stability of the 5'-end hairpin that the ribosome first encounters during translation.

Cell cycle control (and more) by programmed -1 ribosomal frameshifting: implications for disease and therapeutics

Like most basic molecular mechanisms, programmed –1 ribosomal frameshifting (−1 PRF) was first identified in viruses. Early observations that global dysregulation of −1 PRF had deleterious effects on yeast cell growth suggested that −1 PRF may be used to control cellular gene expression, and the cell cycle in particular. Collection of sufficient numbers of viral −1 PRF signals coupled with advances in computer sciences enabled 2 complementary computational approaches to identify −1 PRF signals in free living organisms. The unexpected observation that almost all −1 PRF events on eukaryotic mRNAs direct ribosomes to premature termination codons engendered the hypothesis that −1 PRF signals post-transcriptionally regulate gene expression by functioning as mRNA destabilizing elements. Emerging research suggests that some human diseases are associated with global defects in −1 PRF. The recent discovery of −1 PRF signal-specific trans-acting regulators may provide insight into novel therapeutic strategies aimed at treating diseases caused by changes in gene expression patterns.

Generalized Manning Condensation Model Captures the RNA Ion Atmosphere

RNA is highly sensitive to the ionic environment and typically requires Mg(2+) to form compact structures. There is a need for models capable of describing the ion atmosphere surrounding RNA with quantitative accuracy. We present a model of RNA electrostatics and apply it within coarse-grained molecular dynamics simulation. The model treats Mg(2+) ions explicitly to account for ion-ion correlations neglected by mean-field theories. Since mean-field theories capture KCl well, it is treated implicitly by a generalized Manning counterion condensation model. The model extends Manning condensation to deal with arbitrary RNA conformations, nonlimiting KCl concentrations, and the ion inaccessible volume of RNA. The model is tested against experimental measurements of the excess Mg(2+) associated with the RNA, Γ(2+), because Γ(2+) is directly related to the Mg(2+)-RNA interaction free energy. The excellent agreement with experiment demonstrates that the model captures the ionic dependence of the RNA free energy landscape.

Progress and outlook in structural biology of large viral RNAs

The field of viral molecular biology has reached a precipice for which pioneering studies on the structure of viral RNAs are beginning to bridge the gap. It has become clear that viral genomic RNAs are not simply carriers of hereditary information, but rather are active players in many critical stages during replication. Indeed, functions such as cap-independent translation initiation mechanisms are, in some cases, primarily driven by RNA structural determinants. Other stages including reverse transcription initiation in retroviruses, nuclear export and viral packaging are specifically dependent on the proper 3-dimensional folding of multiple RNA domains to recruit necessary viral and host factors required for activity. Furthermore, a large-scale conformational change within the 5'-untranslated region of HIV-1 has been proposed to regulate the temporal switch between viral protein synthesis and packaging. These RNA-dependent functions are necessary for replication of many human disease-causing viruses such as severe acute respiratory syndrome (SARS)-associated coronavirus, West Nile virus, and HIV-1. The potential for antiviral development is currently hindered by a poor understanding of RNA-driven molecular mechanisms, resulting from a lack of structural information on large RNAs and ribonucleoprotein complexes. Herein, we describe the recent progress that has been made on characterizing these large RNAs and provide brief descriptions of the techniques that will be at the forefront of future advances. Ongoing and future work will contribute to a more complete understanding of the lifecycles of retroviruses and RNA viruses and potentially lead to novel antiviral strategies.

Local structural and environmental factors define the efficiency of an RNA pseudoknot involved in programmed ribosomal frameshift process

In programmed -1 ribosomal frameshift, an RNA pseudoknot stalls the ribosome at specific sequence and restarts translation in a new reading frame. A precise understanding of structural characteristics of these pseudoknots and their PRF inducing ability has not been clear to date. To investigate this phenomenon, we have studied various structural aspects of a -1 PRF inducing RNA pseudoknot from BWYV using extensive molecular dynamics simulations. A set of functional and poorly functional forms, for which previous mutational data were available, were chosen for analysis. These structures differ from each other by either single base substitutions or base-pair replacements from the native structure. We have rationalized how certain mutations in RNA pseudoknot affect its function; e.g., a specific base substitution in loop 2 stabilizes the junction geometry by forming multiple noncanonical hydrogen bonds, leading to a highly rigid structure that could effectively resist ribosome-induced unfolding, thereby increasing efficiency. While, a CG to AU pair substitution in stem 1 leads to loss of noncanonical hydrogen bonds between stems and loop, resulting in a less stable structure and reduced PRF inducing ability, inversion of a pair in stem 2 alters specific base-pair geometry that might be required in ribosomal recognition of nucleobase groups, negatively affecting pseudoknot functioning. These observations illustrate that the ability of an RNA pseudoknot to induce -1 PRF with an optimal rate depends on several independent factors that contribute to either the local conformational variability or geometry.

Many-body effect in ion binding to RNA

Ion-mediated electrostatic interactions play an important role in RNA folding stability. For a RNA in a solution with higher Mg(2+) ion concentration, more counterions in the solution can bind to the RNA, causing a strong many-body coupling between the bound ions. The many-body effect can change the effective potential of mean force between the tightly bound ions. This effect tends to dampen ion binding and lower RNA folding stability. Neglecting the many-body effect leads to a systematic error (over-estimation) of RNA folding stability at high Mg(2+) ion concentrations. Using the tightly bound ion model combined with a conformational ensemble model, we investigate the influence of the many-body effect on the ion-dependent RNA folding stability. Comparisons with the experimental data indicate that including the many-body effect led to much improved predictions for RNA folding stability at high Mg(2+) ion concentrations. The results suggest that the many-body effect can be important for RNA folding in high concentrations of multivalent ions. Further investigation showed that the many-body effect can influence the spatial distribution of the tightly bound ions and the effect is more pronounced for compact RNA structures and structures prone to the formation of local clustering of ions.

A genome-wide analysis of RNA pseudoknots that stimulate efficient -1 ribosomal frameshifting or readthrough in animal viruses

Programmed −1 ribosomal frameshifting (PRF) and stop codon readthrough are two translational recoding mechanisms utilized by some RNA viruses to express their structural and enzymatic proteins at a defined ratio. Efficient recoding usually requires an RNA pseudoknot located several nucleotides downstream from the recoding site. To assess the strategic importance of the recoding pseudoknots, we have carried out a large scale genome-wide analysis in which we used an in-house developed program to detect all possible H-type pseudoknots within the genomic mRNAs of 81 animal viruses. Pseudoknots are detected downstream from ~85% of the recoding sites, including many previously unknown pseudoknots. ~78% of the recoding pseudoknots are the most stable pseudoknot within the viral genomes. However, they are not as strong as some designed pseudoknots that exhibit roadblocking effect on the translating ribosome. Strong roadblocking pseudoknots are not detected within the viral genomes. These results indicate that the decoding pseudoknots have evolved to possess optimal stability for efficient recoding. We also found that the sequence at the gag-pol frameshift junction of HIV1 harbors potential elaborated pseudoknots encompassing the frameshift site. A novel mechanism is proposed for possible involvement of the elaborated pseudoknots in the HIV1 PRF event.

pKa shifting in double-stranded RNA is highly dependent upon nearest neighbors and bulge positioning

Shifting of pKa's in RNA is important for many biological processes; however, the driving forces responsible for shifting are not well understood. Herein, we determine how structural environments surrounding protonated bases affect pKa shifting in double-stranded RNA (dsRNA). Using (31)P NMR, we determined the pKa of the adenine in an A(+)·C base pair in various sequence and structural environments. We found a significant dependence of pKa on the base pairing strength of nearest neighbors and the location of a nearby bulge. Increasing nearest neighbor base pairing strength shifted the pKa of the adenine in an A(+)·C base pair higher by an additional 1.6 pKa units, from 6.5 to 8.1, which is well above neutrality. The addition of a bulge two base pairs away from a protonated A(+)·C base pair shifted the pKa by only ~0.5 units less than a perfectly base paired hairpin; however, positioning the bulge just one base pair away from the A(+)·C base pair prohibited formation of the protonated base pair as well as several flanking base pairs. Comparison of data collected at 25 °C and 100 mM KCl to biological temperature and Mg(2+) concentration revealed only slight pKa changes, suggesting that similar sequence contexts in biological systems have the potential to be protonated at biological pH. We present a general model to aid in the determination of the roles protonated bases may play in various dsRNA-mediated processes including ADAR editing, miRNA processing, programmed ribosomal frameshifting, and general acid-base catalysis in ribozymes.

The emerging role of triple helices in RNA biology

The ability of RNA to form sophisticated secondary and tertiary structures enables it to perform a wide variety of cellular functions. One tertiary structure, the RNA triple helix, was first observed in vitro over 50 years ago, but biological activities for triple helices are only beginning to be appreciated. The recent determination of several RNA structures has implicated triple helices in distinct biological functions. For example, the SAM-II riboswitch forms a triple helix that creates a highly specific binding pocket for S-adenosylmethionine. In addition, a triple helix in the conserved pseudoknot domain of the telomerase-associated RNA TER is essential for telomerase activity. A viral RNA cis-acting RNA element called the ENE contributes to the nuclear stability of a viral noncoding RNA by forming a triple helix with the poly(A) tail. Finally, a cellular noncoding RNA, MALAT1, includes a triple helix at its 3'-end that contributes to RNA stability, but surprisingly also supports translation. These examples highlight the diverse roles that RNA triple helices play in biology. Moreover, the dissection of triple helix mechanisms has the potential to uncover fundamental pathways in cell biology.

The complete nucleotide sequence of the genome of Barley yellow dwarf virus-RMV reveals it to be a new Polerovirus distantly related to other yellow dwarf viruses

The yellow dwarf viruses (YDVs) of the Luteoviridae family represent the most widespread group of cereal viruses worldwide. They include the Barley yellow dwarf viruses (BYDVs) of genus Luteovirus, the Cereal yellow dwarf viruses (CYDVs) and Wheat yellow dwarf virus (WYDV) of genus Polerovirus. All of these viruses are obligately aphid transmitted and phloem-limited. The first described YDVs (initially all called BYDV) were classified by their most efficient vector. One of these viruses, BYDV-RMV, is transmitted most efficiently by the corn leaf aphid, Rhopalosiphum maidis. Here we report the complete 5612 nucleotide sequence of the genomic RNA of a Montana isolate of BYDV-RMV (isolate RMV MTFE87, Genbank accession no. KC921392). The sequence revealed that BYDV-RMV is a polerovirus, but it is quite distantly related to the CYDVs or WYDV, which are very closely related to each other. Nor is BYDV-RMV closely related to any other particular polerovirus. Depending on the gene that is compared, different poleroviruses (none of them a YDV) share the most sequence similarity to BYDV-RMV. Because of its distant relationship to other YDVs, and because it commonly infects maize via its vector, R. maidis, we propose that BYDV-RMV be renamed Maize yellow dwarf virus-RMV (MYDV-RMV).

Exploiting preQ(1) riboswitches to regulate ribosomal frameshifting

Knowing the molecular details of the interaction between riboswitch aptamers and their corresponding metabolites is important to understand gene expression. Here we report on a novel in vitro assay to study preQ(1) riboswitch aptamers upon binding of 7-aminomethyl-7-deazaguanine (preQ(1)). The assay is based on the ability of the preQ(1) aptamer to fold, upon ligand binding, into a pseudoknotted structure that is capable of stimulating -1 ribosomal frameshifting (-1 FS). Aptamers from three different species were found to induce between 7% and 20% of -1 FS in response to increasing preQ(1) levels, whereas preQ(1) analogues were 100-1000-fold less efficient. In depth mutational analysis of the Fusobacterium nucleatum aptamer recapitulates most of the structural details previously identified for preQ(1) aptamers from other bacteria by crystallography and/or NMR spectroscopy. In addition to providing insight into the role of individual nucleotides of the preQ(1) riboswitch aptamer in ligand binding, the presented system provides a valuable tool to screen small molecules against bacterial riboswitches in a eukaryotic background.

Novel cis-acting element within the capsid-coding region enhances flavivirus viral-RNA replication by regulating genome cyclization

cis-Acting elements in the viral genome RNA (vRNA) are essential for the translation, replication, and/or encapsidation of RNA viruses. In this study, a novel conserved cis-acting element was identified in the capsid-coding region of mosquito-borne flavivirus. The downstream of 5' cyclization sequence (5'CS) pseudoknot (DCS-PK) element has a three-stem pseudoknot structure, as demonstrated by structure prediction and biochemical analysis. Using dengue virus as a model, we show that DCS-PK enhances vRNA replication and that its function depends on its secondary structure and specific primary sequence. Mutagenesis revealed that the highly conserved stem 1 and loop 2, which are involved in potential loop-helix interactions, are crucial for DCS-PK function. A predicted loop 1-stem 3 base triple interaction is important for the structural stability and function of DCS-PK. Moreover, the function of DCS-PK depends on its position relative to the 5'CS, and the presence of DCS-PK facilitates the formation of 5'-3' RNA complexes. Taken together, our results reveal that the cis-acting element DCS-PK enhances vRNA replication by regulating genome cyclization, and DCS-PK might interplay with other cis-acting elements to form a functional vRNA cyclization domain, thus playing critical roles during the flavivirus life cycle and evolution.

Molecular mechanism of preQ1 riboswitch action: a molecular dynamics study

Riboswitches often occur in the 5'-untranslated regions of bacterial mRNA where they regulate gene expression. The preQ(1) riboswitch controls the biosynthesis of a hypermodified nucleoside queuosine in response to binding the queuosine metabolic intermediate. Structures of the ligand-bound and ligand-free states of the preQ(1) riboswitch from Thermoanaerobacter tengcongensis were determined recently by X-ray crystallography. We used multiple, microsecond-long molecular dynamics simulations (29 μs in total) to characterize the structural dynamics of preQ(1) riboswitches in both states. We observed different stabilities of the stem in the bound and free states, resulting in different accessibilities of the ribosome-binding site. These differences are related to different stacking interactions between nucleotides of the stem and the associated loop, which itself adopts different conformations in the bound and free states. We suggest that the loop not only serves to bind preQ(1) but also transmits information about ligand binding from the ligand-binding pocket to the stem, which has implications for mRNA accessibility to the ribosome. We explain functional results obscured by a high salt crystallization medium and help to refine regions of disordered electron density, which demonstrates the predictive power of our approach. Besides investigating the functional dynamics of the riboswitch, we have also utilized this unique small folded RNA system for analysis of performance of the RNA force field on the μs time scale. The latest AMBER parmbsc0χ(OL3) RNA force field is capable of providing stable trajectories of the folded molecule on the μs time scale. On the other hand, force fields that are not properly balanced lead to significant structural perturbations on the sub-μs time scale, which could easily lead to inappropriate interpretation of the simulation data.

Programmed -1 frameshifting efficiency correlates with RNA pseudoknot conformational plasticity, not resistance to mechanical unfolding

Programmed -1 frameshifting, whereby the reading frame of a ribosome on messenger RNA is shifted in order to generate an alternate gene product, is often triggered by a pseudoknot structure in the mRNA in combination with an upstream slippery sequence. The efficiency of frameshifting varies widely for different sites, but the factors that determine frameshifting efficiency are not yet fully understood. Previous work has suggested that frameshifting efficiency is related to the resistance of the pseudoknot against mechanical unfolding. We tested this hypothesis by studying the mechanical properties of a panel of pseudoknots with frameshifting efficiencies ranging from 2% to 30%: four pseudoknots from retroviruses, two from luteoviruses, one from a coronavirus, and a nonframeshifting bacteriophage pseudoknot. Using optical tweezers to apply tension across the RNA, we measured the distribution of forces required to unfold each pseudoknot. We found that neither the average unfolding force, nor the unfolding kinetics, nor the parameters describing the energy landscape for mechanical unfolding of the pseudoknot (energy barrier height and distance to the transition state) could be correlated to frameshifting efficiency. These results indicate that the resistance of pseudoknots to mechanical unfolding is not a primary determinant of frameshifting efficiency. However, increased frameshifting efficiency was correlated with an increased tendency to form alternate, incompletely folded structures, suggesting a more complex picture of the role of the pseudoknot involving the conformational dynamics.

Predicting ion-nucleic acid interactions by energy landscape-guided sampling

The recently developed Tightly Bound Ion (TBI) model offers improved predictions for ion effect in nucleic acid systems by accounting for ion correlation and fluctuation effects. However, further application of the model to larger systems is limited by the low computational efficiency of the model. Here, we develop a new computational efficient TBI model using free energy landscape-guided sampling method. The method leads to drastic reduction in the computer time by a factor of 50 for RNAs of 50-100 nucleotides long. The improvement in the computational efficiency would be more significant for larger structures. To test the new method, we apply the model to predict the free energies and the number of bound ions for a series of RNA folding systems. The validity of this new model is supported by the nearly exact agreement with the results from the original TBI model and the agreement with the experimental data. The method may pave the way for further applications of the TBI model to treat a broad range of biologically significant systems such as tetraloop-receptor and riboswitches.

Statistical mechanical modeling of RNA folding: from free energy landscape to tertiary structural prediction

In spite of the success of computational methods for predicting RNA secondary structure, the problem of predicting RNA tertiary structure folding remains. Low-resolution structural models show promise as they allow for rigorous statistical mechanical computation for the conformational entropies, free energies, and the coarse-grained structures of tertiary folds. Molecular dynamics refinement of coarse-grained structures leads to all-atom 3D structures. Modeling based on statistical mechanics principles also has the unique advantage of predicting the full free energy landscape, including local minima and the global free energy minimum. The energy landscapes combined with the 3D structures form the basis for quantitative predictions of RNA functions. In this chapter, we present an overview of statistical mechanical models for RNA folding and then focus on a recently developed RNA statistical mechanical model -- the Vfold model. The main emphasis is placed on the physics underpinning the models, the computational strategies, and the connections to RNA biology.

RNA and its ionic cloud: solution scattering experiments and atomically detailed simulations

RNA molecules play critical roles in many cellular processes. Traditionally viewed as genetic messengers, RNA molecules were recently discovered to have diverse functions related to gene regulation and expression. RNA also has great potential as a therapeutic and a tool for further investigation of gene regulation. Metal ions are an integral part of RNA structure and should be considered in any experimental or theoretical study of RNA. Here, we report a multidisciplinary approach that combines anomalous small-angle x-ray scattering and molecular-dynamics (MD) simulations with explicit solvent and ions around RNA. From experiment and simulation results, we find excellent agreement in the number and distribution of excess monovalent and divalent ions around a short RNA duplex. Although similar agreement can be obtained from a continuum description of the solvent and mobile ions (by solving the Poisson-Boltzmann equation and accounting for finite ion size), the use of MD is easily extended to flexible RNA systems with thermal fluctuations. Therefore, we also model a short RNA pseudoknot and find good agreement between the MD results and the experimentally derived solution structures. Surprisingly, both deviate from crystal structure predictions. These favorable comparisons of experiment and simulations encourage work on RNA in all-atom dynamic models.

The synthesis of 2'-methylseleno adenosine and guanosine 5'-triphosphates

Modified nucleoside triphosphates (NTPs) represent powerful building blocks to generate nucleic acids with novel properties by enzymatic synthesis. We have recently demonstrated the access to 2′-SeCH₃-uridine and 2′-SeCH₃-cytidine derivatized RNAs for applications in RNA crystallography, using the corresponding nucleoside triphosphates and distinct mutants of T7 RNA polymerase. In the present note, we introduce the chemical synthesis of the novel 2′-methylseleno-2′-deoxyadenosine and -guanosine 5′-triphosphates (2′-SeCH₃-ATP and 2′-SeCH₃-GTP) that represent further candidates for the enzymatic RNA synthesis with engineered RNA polymerases.

A domain-based model for predicting large and complex pseudoknotted structures

Pseudoknotted structures play important structural and functional roles in RNA cellular functions at the level of transcription, splicing and translation. However, the problem of computational prediction for large pseudoknotted folds remains. Here we develop a domain-based method for predicting complex and large pseudoknotted structures from RNA sequences. The model is based on the observation that large RNAs can be separated into different structural domains. The basic idea is to first identify the domains and then predict the structures for each domain. Assembly of the domain structures gives the full structure. The use of the domain-based approach leads to a reduction of computational time by a factor of about ~N ( 2) for an N-nt sequence. As applications of the model, we predict structures for a variety of RNA systems, such as regions in human telomerase RNA (hTR), internal ribosome entry site (IRES) and HIV genome. The lengths of these sequences range from 200-nt to 400-nt. The results show good agreements with the experiments.

Novel base triples in RNA structures revealed by graph theoretical searching methods

Background

Highly hydrogen bonded base interactions play a major part in stabilizing the tertiary structures of complex RNA molecules, such as transfer-RNAs, ribozymes and ribosomal RNAs.

Results

We describe the graph theoretical identification and searching of highly hydrogen bonded base triples, where each base is involved in at least two hydrogen bonds with the other bases. Our approach correlates theoretically predicted base triples with literature-based compilations and other actual occurrences in crystal structures. The use of ‘fuzzy’ search tolerances has enabled us to discover a number of triple interaction types that have not been previously recorded nor predicted theoretically.

Conclusions

Comparative analyses of different ribosomal RNA structures reveal several conserved base triple motifs in 50S rRNA structures, indicating an important role in structural stabilization and ultimately RNA function.

Ribosomal frameshifting into an overlapping gene in the 2B-encoding region of the cardiovirus genome

The genus Cardiovirus (family Picornaviridae) currently comprises the species Encephalomyocarditis virus (EMCV) and Theilovirus. Cardioviruses have a positive-sense, single-stranded RNA genome that encodes a large polyprotein (L-1ABCD-2ABC-3ABCD) that is cleaved to produce approximately 12 mature proteins. We report on a conserved ORF that overlaps the 2B-encoding sequence of EMCV in the +2 reading frame. The ORF is translated as a 128-129 amino acid transframe fusion (2B*) with the N-terminal 11-12 amino acids of 2B, via ribosomal frameshifting at a conserved GGUUUUY motif. Mutations that knock out expression of 2B* result in a small-plaque phenotype. Curiously, although theilovirus sequences lack a long ORF in the +2 frame at this genomic location, they maintain a conserved GGUUUUU motif just downstream of the 2A-2B junction, and a highly localized peak in conservation at polyprotein-frame synonymous sites suggests that theiloviruses also utilize frameshifting here, albeit into a very short +2-frame ORF. Unlike previous cases of programmed -1 frameshifting, here frameshifting is modulated by virus infection, thus suggesting a novel regulatory role for frameshifting in these viruses.

Salt contribution to RNA tertiary structure folding stability

Accurate quantification of the ionic contribution to RNA folding stability could greatly enhance our ability to understand and predict RNA functions. Recently, motivated by the potential importance of ion correlation and fluctuation in RNA folding, we developed the tightly bound ion (TBI) model. Extensive experimental tests showed that the TBI model can lead to better treatment of multivalent ions than the Poisson-Boltzmann equation. In this study, we use the model to quantify the contribution of salt (Na(+) and Mg(2+)) to the RNA tertiary structure folding free energy. Folding of the RNA tertiary structure often involves intermediates. We focus on the folding transition from an intermediate state to the native state, and compute the electrostatic folding free energy of the RNA. Based on systematic calculations for a variety of RNA molecules, we derive a set of formulas for the electrostatic free energy for tertiary structural folding as a function of the sequence length and compactness of the RNA and the Na(+) and Mg(2+) concentrations. Extensive comparisons with experimental data suggest that our model and the extracted empirical formulas are quite reliable.