Structure of an unusually stable RNA hairpin

Structure of Essential RNA Regulatory Elements in the West Nile Virus 3'-Terminal Stem Loop

Flaviviruses, such as West Nile and Dengue Virus, pose a significant and growing threat to global health. Central to the flavivirus life cycle are highly structured 5'- and 3'-untranslated regions (UTRs), which harbor conserved cis-acting RNA elements critical for viral replication and host adaptation. Despite their essential roles, detailed molecular insights into these RNA elements have been limited. By employing nuclear magnetic resonance (NMR) spectroscopy in conjunction with SAXS experiments, we determined the three-dimensional structure of the West Nile Virus (WNV) 3'-terminal stem-loop core, a highly conserved element critical for viral genome cyclization and replication. Single nucleotide mutations at several sites within this RNA abolish the ability of the virus to replicate. These critical sites are located within a short 18-nucleotide hairpin stem, a substructure notable for its conformational flexibility, while the adjoining main stem-loop adopts a well-defined extended helix interrupted by three non-Watson-Crick pairs. This study enhances our understanding of several metastable RNA structures that play key roles in regulating the flavivirus lifecycle, and thereby also opens up potential new avenues for the development of antivirals targeting these conserved RNA structures. In particular, the structure we observe suggests that the plastic junction between the small hairpin and the tail of the longer stem-loop could provide a binding pocket for small molecules, for example potentially stabilizing the RNA in a conformation which hinders the conformational rearrangements critical for viral replication.

Structure of an internal loop motif with three consecutive U•U mismatches from stem-loop 1 in the 3'-UTR of the SARS-CoV-2 genomic RNA

The single-stranded RNA genome of SARS-CoV-2 is highly structured. Numerous helical stem-loop structures interrupted by mismatch motifs are present in the functionally important 5'- and 3'-UTRs. These mismatches modulate local helical geometries and feature unusual arrays of hydrogen bonding donor and acceptor groups. However, their conformational and dynamical properties cannot be directly inferred from chemical probing and are difficult to predict theoretically. A mismatch motif (SL1-motif) consisting of three consecutive U•U base pairs is located in stem-loop 1 of the 3'-UTR. We combined NMR-spectroscopy and MD-simulations to investigate its structure and dynamics. All three U•U base pairs feature two direct hydrogen bonds and are as stable as Watson-Crick A:U base pairs. Plasmodium falciparum 25S rRNA contains a triple U•U mismatch motif (Pf-motif) differing from SL1-motif only with respect to the orientation of the two closing base pairs. Interestingly, while the geometry of the outer two U•U mismatches was identical in both motifs the preferred orientation of the central U•U mismatch was different. MD simulations and potassium ion titrations revealed that the potassium ion-binding mode to the major groove is connected to the different preferred geometries of the central base pair in the two motifs.

RNA adapts its flexibility to efficiently fold and resist unfolding

Recent studies have demonstrated that the mechanisms through which biopolymers like RNA interconvert between multiple folded structures are critical for their cellular functions. A major obstacle to elucidating these mechanisms is the lack of experimental approaches that can resolve these interconversions between functionally relevant biomolecular structures. Here, using a nano-electronic device with microsecond time resolution, we dissect the complete set of structural rearrangements executed by an ultra-stable RNA, the UUCG stem-loop, at the single-molecule level. We show that the stem-loop samples at least four conformations along two folding pathways leading to two distinct folded structures, only one of which has been previously observed. By modulating its flexibility, the stem-loop can adaptively select between these pathways, enabling it to both fold rapidly and resist unfolding. This paradigm of stabilization through compensatory changes in flexibility broadens our understanding of stable RNA structures and is expected to serve as a general strategy employed by all biopolymers.

Identification and characterization of RNA pentaloop sequence families

One of the current methods for predicting RNA tertiary structure is fragment-based homology, which predicts tertiary structure from secondary structure. For a successful prediction, this method requires a library of the tertiary structures of small motifs clipped from previously solved RNA 3D structures. Because of the limited number of available tertiary structures, it is not practical to find structures for all sequences of all motifs. Identifying sequence families for motifs can fill the gaps because all sequences within a family are expected to have similar structural features. Currently, a collection of well-characterized sequence families has been identified for tetraloops. Because of their prevalence and biological functions, pentaloop structures should also be well-characterized. In this study, 10 pentaloop sequence families are identified. For each family, the common and distinguishing structural features are highlighted. These sequence families can be used to predict the tertiary structure of pentaloop sequences for which a solved structure is not available.

Shedding light on the base-pair opening dynamics of nucleic acids in living human cells

Base-pair opening is a fundamental property of nucleic acids that plays important roles in biological functions. However, studying the base-pair opening dynamics inside living cells has remained challenging. Here, to determine the base-pair opening kinetics inside living human cells, the exchange rate constant ([Formula: see text]) of the imino proton with the proton of solvent water involved in hairpin and G-quadruplex (GQ) structures is determined by the in-cell NMR technique. It is deduced on determination of [Formula: see text] values that at least some G-C base pairs of the hairpin structure and all G-G base-pairs of the GQ structure open more frequently in living human cells than in vitro. It is suggested that interactions with endogenous proteins could be responsible for the increase in frequency of base-pair opening. Our studies demonstrate a difference in dynamics of nucleic acids between in-cell and in vitro conditions.

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

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

Structure of the dengue virus RNA promoter

Dengue virus, a single-stranded positive sense RNA virus, is the most prevalent mosquito-borne pathogen in the world. Like all RNA viruses, it uses conserved structural elements within its genome to control essential replicative steps. A 70 nt stem-loop RNA structure (called SLA), found at the 5'-end of the genome of all flaviviruses, functions as the promoter for viral replication. This highly conserved structure interacts with the viral polymerase NS5 to initiate RNA synthesis. Here, we report the NMR structure of a monomeric SLA from dengue virus serotype 1, assembled to high-resolution from independently folded structural elements. The DENV1 SLA has an L-shaped structure, where the top and side helices are coaxially stacked, and the bottom helix is roughly perpendicular to them. Because the sequence is highly conserved among different flavivirus genomes, it is very likely that the three-dimensional fold and local structure of SLA are also conserved among flaviviruses and required for efficient replication. This work provides structural insight into the dengue promoter and provides the foundation for the discovery of new antiviral drugs that target this essential replicative step.

Lone Pair…π Contacts and Structure Signatures of r(UNCG) Tetraloops, Z-Turns, and Z-Steps: A WebFR3D Survey

Z-DNA and Z-RNA have long appeared as oddities to nucleic acid scientists. However, their Z-step constituents are recurrently observed in all types of nucleic acid systems including ribosomes. Z-steps are NpN steps that are isostructural to Z-DNA CpG steps. Among their structural features, Z-steps are characterized by the presence of a lone pair…π contact that involves the stacking of the ribose O4′ atom of the first nucleotide with the 3′-face of the second nucleotide. Recently, it has been documented that the CpG step of the ubiquitous r(UNCG) tetraloops is a Z-step. Accordingly, such r(UNCG) conformations were called Z-turns. It has also been recognized that an r(GAAA) tetraloop in appropriate conditions can shapeshift to an unusual Z-turn conformation embedding an ApA Z-step. In this report, we explore the multiplicity of RNA motifs based on Z-steps by using the WebFR3D tool to which we added functionalities to be able to retrieve motifs containing lone pair…π contacts. Many examples that underscore the diversity and universality of these motifs are provided as well as tutorial guidance on using WebFR3D. In addition, this study provides an extensive survey of crystallographic, cryo-EM, NMR, and molecular dynamics studies on r(UNCG) tetraloops with a critical view on how to conduct database searches and exploit their results.

Improved gRNA secondary structures allow editing of target sites resistant to CRISPR-Cas9 cleavage

The first step in CRISPR-Cas9-mediated genome editing is the cleavage of target DNA sequences that are complementary to so-called spacer sequences in CRISPR guide RNAs (gRNAs). However, some DNA sequences are refractory to CRISPR-Cas9 cleavage, which is at least in part due to gRNA misfolding. To overcome this problem, we have engineered gRNAs with highly stable hairpins in their constant parts and further enhanced their stability by chemical modifications. The 'Genome-editing Optimized Locked Design' (GOLD)-gRNA increases genome editing efficiency up to around 1000-fold (from 0.08 to 80.5%) with a mean increase across different other targets of 7.4-fold. We anticipate that this improved gRNA will allow efficient editing regardless of spacer sequence composition and will be especially useful if a desired genomic site is difficult to edit.

Polysomes Bypass a 50-Nucleotide Coding Gap Less Efficiently Than Monosomes Due to Attenuation of a 5' mRNA Stem-Loop and Enhanced Drop-off

Efficient translational bypassing of a 50-nt non-coding gap in a phage T4 topoisomerase subunit gene (gp60) requires several recoding signals. Here we investigate the function of the mRNA stem–loop 5′ of the take-off codon, as well as the importance of ribosome loading density on the mRNA for efficient bypassing. We show that polysomes are less efficient at mediating bypassing than monosomes, both in vitro and in vivo, due to their preventing formation of a stem–loop 5′ of the take-off codon and allowing greater peptidyl-tRNA drop off. A ribosome profiling analysis of phage T4-infected Escherichia coli yielded protected mRNA fragments within the normal size range derived from ribosomes stalled at the take-off codon. However, ribosomes at this position also yielded some 53-nucleotide fragments, 16 longer. These were due to protection of the nucleotides that form the 5′ stem–loop. NMR shows that the 5′ stem–loop is highly dynamic. The importance of different nucleotides in the 5′ stem–loop is revealed by mutagenesis studies. These data highlight the significance of the 5′ stem–loop for the 50-nt bypassing and further enhance appreciation of relevance of the extent of ribosome loading for recoding.

•

Monosomes are more efficient than polysome in mediating 50-nt translational bypassing.

•

A 5′ mRNA stem–loop facilitates translational bypassing by monosomes.

•

Ribosome profiling yields an extra-long, 53-nt, protected fragment of mRNA.

Structure of the RNA Specialized Translation Initiation Element that Recruits eIF3 to the 5'-UTR of c-Jun

Specialized translation initiation is a novel form of regulation of protein synthesis, whereby RNA structures within the 5'-UTR regulate translation rates of specific mRNAs. Similar to internal ribosome entry sites (IRESs), specialized translation initiation requires the recruitment of eukaryotic initiation factor 3 (eIF3), but also requires cap recognition by eIF3d, a new 5'-m7GTP recognizing protein. How these RNA structures mediate eIF3 recruitment to affect translation of specific mRNAs remains unclear. Here, we report the nuclear magnetic resonance (NMR) structure of a stem-loop within the c-JUN 5' UTR recognized by eIF3 and essential for specialized translation initiation of this well-known oncogene. The structure exhibits similarity to eIF3 recognizing motifs found in hepatitis C virus (HCV)-like IRESs, suggesting mechanistic similarities. This work establishes the RNA structural features involved in c-JUN specialized translation initiation and provides a basis to search for small molecule inhibitors of aberrant expression of the proto-oncogenic c-JUN.

Identification and Characterization of New RNA Tetraloop Sequence Families

There is an abundance of RNA sequence information available due to the efforts of sequencing projects. However, current techniques implemented to solve the tertiary structures of RNA, such as NMR and X-ray crystallography, are difficult and time-consuming. Therefore, biophysical techniques are not able to keep pace with the abundance of sequence information available. Because of this, there is a need to develop quick and efficient ways to predict RNA tertiary structure from sequence. One promising approach is to identify structural patterns within previously solved 3D structures and apply these patterns to new sequences. RNA tetraloops are one of the most common naturally occurring secondary structure motifs. Here, we use RNA Characterization of Secondary Structure Motifs (CoSSMos), Dissecting the Spatial Structure of RNA (DSSR), and a bioinformatic approach to search for and characterize tertiary structure patterns among tetraloops. Not surprising, we identified the well-known GNRA and UNCG tetraloops, as well as the previously identified RNYA tetraloop. However, some previously identified characteristics of these families were not observed in this data set, and some new characteristics were identified. In addition, we also identified and characterized three new tetraloop sequence families: YGAR, UGGU, and RMSA. This new structural information sheds light on the tertiary structure of tetraloops and contributes to the efforts of RNA tertiary structure prediction from sequence.

Conformational adaptation of UNCG loops upon crowding

If the A-form helix is the major structural motif found in RNA, the loops that cap them constitute the second most important family of motifs. Among those, two are overrepresented, GNRA and UNCG tetraloops. Recent surveys of RNA structures deposited in the PDB show that GNRA and UNCG tetraloops can adopt tertiary folds that are very different from their canonical conformations, characterized by the presence of a U-turn of a Z-turn, respectively. Crystallographic data from both a lariat-capping (LC) ribozyme and a group II intron ribozyme reveal that a given UUCG tetraloop can adopt a distinct fold depending on its structural environment. Specifically, when the crystal packing applies relaxed constraints on the loop, the canonical Z-turn conformation is observed. In contrast, a highly packed environment induces "squashing" of the tetraloop by distorting its sugar-phosphate backbone in a specific way that expels the first and fourth nucleobases out of the loop, and falls in van der Waals distance of the last base pair of the helix, taking the place of the pair formed between the first and fourth residues in Z-turn loops. The biological relevance of our observations is supported by the presence of similarly deformed loops in the highly packed environment of the ribosome and in a complex between a dsRNA and a RNase III. The finding that Z-turn loops change conformation under higher molecular packing suggests that, in addition to their demonstrated role in stabilizing RNA folding, they may contribute to the three-dimensional structure of RNA by mediating tertiary interactions with distal residues.

The yeast C/D box snoRNA U14 adopts a "weak" K-turn like conformation recognized by the Snu13 core protein in solution

Non-coding RNAs associate with proteins to form ribonucleoproteins (RNPs), such as ribosome, box C/D snoRNPs, H/ACA snoRNPs, ribonuclease P, telomerase and spliceosome to ensure cell viability. The assembly of these RNA-protein complexes relies on the ability of the RNA to adopt the correct bound conformation. K-turn motifs represent ubiquitous binding platform for proteins found in several cellular environment. This structural motif has an internal three-nucleotide bulge flanked on its 3' side by a G•A/A•G tandem pairs followed by one or two non-Watson-Crick pairs, and on its 5' side by a classical RNA helix. This peculiar arrangement induces a strong curvature of the phosphodiester backbone, which makes it conducive to multiple tertiary interactions. SNU13/Snu13p (Human/Yeast) binds specifically the U14 C/D box snoRNA K-turn sequence motif. This event is the prerequisite to promote the assembly of the RNP, which contains NOP58/Nop58 and NOP56/Nop56 core proteins and the 2'-O-methyl-transferase, Fibrillarin/Nop1p. The U14 small nucleolar RNA is a conserved non-coding RNA found in yeast and vertebrates required for the pre-rRNA maturation and ribose methylation. Here, we report the solution structure of the free U14 snoRNA K-turn motif (kt-U14) as determined by Nuclear Magnetic Resonance. We demonstrate that a major fraction of free kt-U14 adopts a pre-folded conformation similar to protein bound K-turn, even in the absence of divalent ions. In contrast to the kt-U4 or tyrS RNA, kt-U14 displays a sharp bent in the phosphodiester backbone. The U•U and G•A tandem base pairs are formed through weak hydrogen bonds. Finally, we show that the structure of kt-U14 is stabilized upon Snu13p binding. The structure of the free U14 RNA is the first reference example for the canonical motifs of the C/D box snoRNA family.

High-resolution small RNA structures from exact nuclear Overhauser enhancement measurements without additional restraints

RNA not only translates the genetic code into proteins, but also carries out important cellular functions. Understanding such functions requires knowledge of the structure and dynamics at atomic resolution. Almost half of the published RNA structures have been solved by nuclear magnetic resonance (NMR). However, as a result of severe resonance overlap and low proton density, high-resolution RNA structures are rarely obtained from nuclear Overhauser enhancement (NOE) data alone. Instead, additional semi-empirical restraints and labor-intensive techniques are required for structural averages, while there are only a few experimentally derived ensembles representing dynamics. Here we show that our exact NOE (eNOE) based structure determination protocol is able to define a 14-mer UUCG tetraloop structure at high resolution without other restraints. Additionally, we use eNOEs to calculate a two-state structure, which samples its conformational space. The protocol may open an avenue to obtain high-resolution structures of small RNA of unprecedented accuracy with moderate experimental efforts.

Relative Populations of Some Tautomeric Forms of 2'-Deoxyguanosine-5-Fluorouridine Mismatch

The importance of the 2'-deoxyguanosine-uridine mispair as the most occurring mismatch in transcriptional studies of RNAs from DNAs is multiplied when 5-halo-substituted uridine species cause a serious increase in the probability of its occurrence. Many studies relate this higher probability to the existence of possible tautomeric and ionic forms of its constituent bases. According to these statements, relative populations of mismatches between 5-fluorouridine and both keto and enol forms of 2'-deoxyguanosine are computed by using a conformational search. In order to have a complete scan of all of the highly probable conformers in a moderate computational time, an extensive conformational search methodology is employed here, which benefits from the advantages of both the molecular dynamics simulations and quantum mechanics calculations. The population of an enolic tautomer of normal wobble orientation is about 0.057% of that of its keto tautomer, whereas the population of an enolic tautomer of reverse wobble orientation is about 0.0054% of that of its keto tautomer. Totally, the reverse wobble orientation is about six times more populated than the normal wobble orientation. Calculated populations are in good agreement with experimental populations of closely related compounds. The reliability of the applied methodology is certified, in part, by a good agreement obtained between some experimental data and corresponding Boltzmann-weighted average data of most probable conformers such as NMR parameters. The validation of this methodology is certified with high accuracy by applying it on the substituted diuridine pairs, where experimental populations are available. Not only are the calculated populations and NMR parameters of this test in very good agreement with the experimental data, but also they are free of the ambiguities mentioned by experimentalists.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Emergency Services of Viral RNAs: Repair and Remodeling

Reproduction of RNA viruses is typically error-prone due to the infidelity of their replicative machinery and the usual lack of proofreading mechanisms. The error rates may be close to those that kill the virus. Consequently, populations of RNA viruses are represented by heterogeneous sets of genomes with various levels of fitness. This is especially consequential when viruses encounter various bottlenecks and new infections are initiated by a single or few deviating genomes. Nevertheless, RNA viruses are able to maintain their identity by conservation of major functional elements. This conservatism stems from genetic robustness or mutational tolerance, which is largely due to the functional degeneracy of many protein and RNA elements as well as to negative selection. Another relevant mechanism is the capacity to restore fitness after genetic damages, also based on replicative infidelity. Conversely, error-prone replication is a major tool that ensures viral evolvability. The potential for changes in debilitated genomes is much higher in small populations, because in the absence of stronger competitors low-fit genomes have a choice of various trajectories to wander along fitness landscapes. Thus, low-fit populations are inherently unstable, and it may be said that to run ahead it is useful to stumble. In this report, focusing on picornaviruses and also considering data from other RNA viruses, we review the biological relevance and mechanisms of various alterations of viral RNA genomes as well as pathways and mechanisms of rehabilitation after loss of fitness. The relationships among mutational robustness, resilience, and evolvability of viral RNA genomes are discussed.

The first successful observation of in-cell NMR signals of DNA and RNA in living human cells

In order to understand intracellular biological events, information on the structure, dynamics and interaction of proteins and nucleic acids in living cells is of crucial importance. In-cell NMR is a promising method to obtain this information. Although NMR signals of proteins in human cells have been reported, those of nucleic acids were reported only in Xenopus laevis oocytes, i.e., not in human cells. Here, DNA and RNA were introduced into human cells by means of pore formation by bacterial toxin streptolysin O and subsequent resealing. Then, NMR signals of DNA and RNA were successfully observed for the first time in living human cells. The observed signals directly suggested the formation of DNA and RNA hairpin structures in living human cells.

Effect of GCAA stabilizing loops on three- and four-way intramolecular junctions

Tetraloops are a common way of changing the melting behavior of a DNA or RNA structure without changing the sequence of the stem. Because of the ubiquitous nature of tetraloops, our goal is to understand the effect a GCAA tetraloop, which belongs to the GNRA family of tetraloops, has on the unfolding thermodynamics of intramolecular junctions. Specifically, we have described the melting behavior of intramolecular three-way and four-way junctions where a T₅ loop has been replaced with a GCAA tetraloops in different positions. Their thermodynamic profiles, including Δn_Na⁺ and Δn_W, were analyzed based on the position of the tetraloop. We obtained between -16.7 and -27.5 kcal mol^-1 for all junctions studied. The experimental data indicates the influence of the GCAA tetraloop is primarily dictated by the native unfolding of the junction; if the tetraloop is placed on a stem that unfolds as a single domain when the tetraloop is not present, it will unfold as a single domain when the tetraloop is present but with a higher thermal stability. Conversely, if the tetraloop is placed on a stem which unfolds cooperatively with other stems when the tetraloop is not present, the tetraloop will increase the thermal stability of all the stems in the melting domain. The oligonucleotide structure and not the tetraloop itself affects ion uptake; three-way junctions do not gain an increase in ion uptake, but four-way junctions do. This is not the case for water immobilization, where the position of the tetraloop dictates the amount of water immobilized.

Identification of 15 candidate structured noncoding RNA motifs in fungi by comparative genomics

BACKGROUND

With the development of rapid and inexpensive DNA sequencing, the genome sequences of more than 100 fungal species have been made available. This dataset provides an excellent resource for comparative genomics analyses, which can be used to discover genetic elements, including noncoding RNAs (ncRNAs). Bioinformatics tools similar to those used to uncover novel ncRNAs in bacteria, likewise, should be useful for searching fungal genomic sequences, and the relative ease of genetic experiments with some model fungal species could facilitate experimental validation studies.

RESULTS

We have adapted a bioinformatics pipeline for discovering bacterial ncRNAs to systematically analyze many fungal genomes. This comparative genomics pipeline integrates information on conserved RNA sequence and structural features with alternative splicing information to reveal fungal RNA motifs that are candidate regulatory domains, or that might have other possible functions. A total of 15 prominent classes of structured ncRNA candidates were identified, including variant HDV self-cleaving ribozyme representatives, atypical snoRNA candidates, and possible structured antisense RNA motifs. Candidate regulatory motifs were also found associated with genes for ribosomal proteins, S-adenosylmethionine decarboxylase (SDC), amidase, and HexA protein involved in Woronin body formation. We experimentally confirm that the variant HDV ribozymes undergo rapid self-cleavage, and we demonstrate that the SDC RNA motif reduces the expression of SAM decarboxylase by translational repression. Furthermore, we provide evidence that several other motifs discovered in this study are likely to be functional ncRNA elements.

CONCLUSIONS

Systematic screening of fungal genomes using a computational discovery pipeline has revealed the existence of a variety of novel structured ncRNAs. Genome contexts and similarities to known ncRNA motifs provide strong evidence for the biological and biochemical functions of some newly found ncRNA motifs. Although initial examinations of several motifs provide evidence for their likely functions, other motifs will require more in-depth analysis to reveal their functions.

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.

NMR experiments for the rapid identification of P=O···H-X type hydrogen bonds in nucleic acids

Hydrogen bonds involving the backbone phosphate groups occur with high frequency in functional RNA molecules. They are often found in well-characterized tertiary structural motifs presenting powerful probes for the rapid identification of these motifs for structure elucidation purposes. We have shown recently that stable hydrogen bonds to the phosphate backbone can in principle be detected by relatively simple NMR-experiments, providing the identity of both the donor hydrogen and the acceptor phosphorous within the same experiment (Duchardt-Ferner et al., Angew Chem Int Ed Engl 50:7927-7930, 2011). However, for imino and hydroxyl hydrogen bond donor groups rapidly exchanging with the solvent as well as amino groups broadened by conformational exchange experimental sensitivity is severely hampered by extensive line broadening. Here, we present improved methods for the rapid identification of hydrogen bonds to phosphate groups in nucleic acids by NMR. The introduction of the SOFAST technique into ¹H,³¹P-correlation experiments as well as a BEST-HNP experiment exploiting ^3hJ_N,P rather than ^2hJ_H,P coupling constants enables the rapid and sensitive identification of these hydrogen bonds in RNA. The experiments are applicable for larger RNAs (up to ~ 100-nt), for donor groups influenced by conformational exchange processes such as amino groups and for hydrogen bonds with rather labile hydrogens such as 2'-OH groups as well as for moderate sample concentrations. Interestingly, the size of the through-hydrogen bond scalar coupling constants depends not only on the type of the donor group but also on the structural context. The largest coupling constants were measured for hydrogen bonds involving the imino groups of protonated cytosine nucleotides as donors.

The crystal structure of the 5΄ functional domain of the transcription riboregulator 7SK

In vertebrates, the 7SK RNA forms the scaffold of a complex, which regulates transcription pausing of RNA-polymerase II. By binding to the HEXIM protein, the complex comprising proteins LARP7 and MePCE captures the positive transcription elongation factor P-TEFb and prevents phosphorylation of pausing factors. The HEXIM-binding site embedded in the 5΄-hairpin of 7SK (HP1) encompasses a short signature sequence, a GAUC repeat framed by single-stranded uridines. The present crystal structure of HP1 shows a remarkably straight helical stack involving several unexpected triples formed at a central region. Surprisingly, two uridines of the signature sequence make triple interactions in the major groove of the (GAUC)₂. The third uridine is turned outwards or inward, wedging between the other uridines, thus filling the major groove. A molecular dynamics simulation indicates that these two conformations of the signature sequence represent stable alternatives. Analyses of the interaction with the HEXIM protein confirm the importance of the triple interactions at the signature sequence. Altogether, the present structural analysis of 7SK HP1 highlights an original mechanism of swapping bases, which could represent a possible ‘7SK signature’ and provides new insight into the functional importance of the plasticity of RNA.

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.

Direct Observation of Folding Energy Landscape of RNA Hairpin at Mechanical Loading Rates

By applying a controlled mechanical load using optical tweezers, we measured the diffusive barrier crossing in a 49 nt long P5ab RNA hairpin. We find that in the free-energy landscape the barrier height (G^‡) and transition distance (x^‡) are dependent on the loading rate (r) along the pulling direction, x, as predicted by Bell. The barrier shifted toward the initial state, whereas ΔG^‡ reduced significantly from 50 to 5 kT, as r increased from 0 to 32 pN/s. However, the equilibrium work (ΔG) during strand separation, as estimated by Crook's fluctuation theorem, remained unchanged at different rates. Previously, helix formation and denaturation have been described as two-state (F ↔ U) transitions for P5ab. Herein, we report three intermediate states I₁, I, and I₂ located at 4, 11, and 16 nm respectively, from the folded conformation. The intermediates were observed only when the hairpin was subjected to an optimal r, 7.6 pN/s. The results indicate that the complementary strands in P5ab can zip and unzip through complex routes, whereby mismatches act as checkpoints and often impose barriers. The study highlights the significance of loading rates in force-spectroscopy experiments that are increasingly being used to measure the folding properties of biomolecules.

Revised RNA Dihedral Parameters for the Amber Force Field Improve RNA Molecular Dynamics

The backbone dihedral parameters of the Amber RNA force field were improved by fitting using multiple linear regression to potential energies determined by quantum chemistry calculations. Five backbone and four glycosidic dihedral parameters were fit simultaneously to reproduce the potential energies determined by a high-level density functional theory calculation (B97D3 functional with the AUG-CC-PVTZ basis set). Umbrella sampling was used to determine conformational free energies along the dihedral angles, and these better agree with the population of conformations observed in the protein data bank for the new parameters than for the conventional parameters. Molecular dynamics simulations performed on a set of hairpin loops, duplexes and tetramers with the new parameter set show improved modeling for the structures of tetramers CCCC, CAAU, and GACC, and an RNA internal loop of noncanonical pairs, as compared to the conventional parameters. For the tetramers, the new parameters largely avoid the incorrect intercalated structures that dominate the conformational samples from the conventional parameters. For the internal loop, the major conformation solved by NMR is stable with the new parameters, but not with the conventional parameters. The new force field performs similarly to the conventional parameters for the UUCG and GCAA hairpin loops and the [U(UA)6A]2 duplex.

Solution structure of the 5'-terminal hairpin of the 7SK small nuclear RNA

The small nuclear 7SK RNA regulates RNA polymerase II (RNA Pol II) transcription, by sequestering and inhibiting the positive transcription elongation factor b (P-TEFb). P-TEFb is stored in the 7SK ribonucleoprotein (RNP) that contains the three nuclear proteins Hexim1, LaRP7, and MePCE. P-TEFb interacts with the protein Hexim1 and the 7SK RNA. Once P-TEFb is released from the 7SK RNP, it activates transcription by phosphorylating the C-terminal domain of RNA Pol II. P-TEFb also plays a crucial role in the replication of the human immunodeficiency virus HIV-1, through its recruitment by the viral transactivator Tat. Previous work demonstrated that the protein Tat promotes the release of P-TEFb from the 7SK RNP through direct binding to the 7SK RNA. Hexim1 and Tat proteins both comprise conserved and similar arginine-rich motifs that were identified to bind the 7SK RNA at a repeated GAUC site located at the top of the 5'-terminal hairpin (HPI). Here, we report the solution structure of this region as determined by nuclear magnetic resonance, to identify HPI structural features recognized by Hexim1 and Tat. The HPI solution structure displays an elongated shape featuring four helical segments interrupted by one internal loop and three bulges with distinct folds. In particular, the repeated GAUC motif adopts a pre-organized geometry. Our results suggest that the binding of Hexim1 and Tat to the 7SK RNA could originate from a conformational selection of this motif, highlighting how RNA local structure could lead to an adaptive recognition of their partners.

An intermolecular G-quadruplex as the basis for GTP recognition in the class V-GTP aptamer

Many naturally occurring or artificially created RNAs are capable of binding to guanine or guanine derivatives with high affinity and selectivity. They bind their ligands using very different recognition modes involving a diverse set of hydrogen bonding and stacking interactions. Apparently, the potential structural diversity for guanine, guanosine, and guanine nucleotide binding motifs is far from being fully explored. Szostak and coworkers have derived a large set of different GTP-binding aptamer families differing widely in sequence, secondary structure, and ligand specificity. The so-called class V-GTP aptamer from this set binds GTP with very high affinity and has a complex secondary structure. Here we use solution NMR spectroscopy to demonstrate that the class V aptamer binds GTP through the formation of an intermolecular two-layered G-quadruplex structure that directly incorporates the ligand and folds only upon ligand addition. Ligand binding and G-quadruplex formation depend strongly on the identity of monovalent cations present with a clear preference for potassium ions. GTP binding through direct insertion into an intermolecular G-quadruplex is a previously unobserved structural variation for ligand-binding RNA motifs and rationalizes the previously observed specificity pattern of the class V aptamer for GTP analogs.

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.

Insights into Stability and Folding of GNRA and UNCG Tetraloops Revealed by Microsecond Molecular Dynamics and Well-Tempered Metadynamics

RNA hairpins capped by 5'-GNRA-3' or 5'-UNCG-3' tetraloops (TLs) are prominent RNA structural motifs. Despite their small size, a wealth of experimental data, and recent progress in theoretical simulations of their structural dynamics and folding, our understanding of the folding and unfolding processes of these small RNA elements is still limited. Theoretical description of the folding and unfolding processes requires robust sampling, which can be achieved by either an exhaustive time scale in standard molecular dynamics simulations or sophisticated enhanced sampling methods, using temperature acceleration or biasing potentials. Here, we study structural dynamics of 5'-GNRA-3' and 5'-UNCG-3' TLs by 15-μs-long standard simulations and a series of well-tempered metadynamics, attempting to accelerate sampling by bias in a few chosen collective variables (CVs). Both methods provide useful insights. The unfolding and refolding mechanisms of the GNRA TL observed by well-tempered metadynamics agree with the (reverse) folding mechanism suggested by recent replica exchange molecular dynamics simulations. The orientation of the glycosidic bond of the GL4 nucleobase is critical for the UUCG TL folding pathway, and our data strongly support the hypothesis that GL4-anti forms a kinetic trap along the folding pathway. Along with giving useful insight, our study also demonstrates that using only a few CVs apparently does not capture the full folding landscape of the RNA TLs. Despite using several sophisticated selections of the CVs, formation of the loop appears to remain a hidden variable, preventing a full convergence of the metadynamics. Finally, our data suggest that the unfolded state might be overstabilized by the force fields used.

Free-energy landscape of a hyperstable RNA tetraloop

We report the characterization of the energy landscape and the folding/unfolding thermodynamics of a hyperstable RNA tetraloop obtained through high-performance molecular dynamics simulations at microsecond timescales. Sampling of the configurational landscape is conducted using temperature replica exchange molecular dynamics over three isochores at high, ambient, and negative pressures to determine the thermodynamic stability and the free-energy landscape of the tetraloop. The simulations reveal reversible folding/unfolding transitions of the tetraloop into the canonical A-RNA conformation and the presence of two alternative configurations, including a left-handed Z-RNA conformation and a compact purine Triplet. Increasing hydrostatic pressure shows a stabilizing effect on the A-RNA conformation and a destabilization of the left-handed Z-RNA. Our results provide a comprehensive description of the folded free-energy landscape of a hyperstable RNA tetraloop and highlight the significant advances of all-atom molecular dynamics in describing the unbiased folding of a simple RNA secondary structure motif.

Spectral density mapping at multiple magnetic fields suitable for (13)C NMR relaxation studies

Standard spectral density mapping protocols, well suited for the analysis of (15)N relaxation rates, introduce significant systematic errors when applied to (13)C relaxation data, especially if the dynamics is dominated by motions with short correlation times (small molecules, dynamic residues of macromolecules). A possibility to improve the accuracy by employing cross-correlated relaxation rates and on measurements taken at several magnetic fields has been examined. A suite of protocols for analyzing such data has been developed and their performance tested. Applicability of the proposed protocols is documented in two case studies, spectral density mapping of a uniformly labeled RNA hairpin and of a selectively labeled disaccharide exhibiting highly anisotropic tumbling. Combination of auto- and cross-correlated relaxation data acquired at three magnetic fields was applied in the former case in order to separate effects of fast motions and conformational or chemical exchange. An approach using auto-correlated relaxation rates acquired at five magnetic fields, applicable to anisotropically moving molecules, was used in the latter case. The results were compared with a more advanced analysis of data obtained by interpolation of auto-correlated relaxation rates measured at seven magnetic fields, and with the spectral density mapping of cross-correlated relaxation rates. The results showed that sufficiently accurate values of auto- and cross-correlated spectral density functions at zero and (13)C frequencies can be obtained from data acquired at three magnetic fields for uniformly (13)C-labeled molecules with a moderate anisotropy of the rotational diffusion tensor. Analysis of auto-correlated relaxation rates at five magnetic fields represents an alternative for molecules undergoing highly anisotropic motions.

The RNA 3D Motif Atlas: Computational methods for extraction, organization and evaluation of RNA motifs

RNA 3D motifs occupy places in structured RNA molecules that correspond to the hairpin, internal and multi-helix junction "loops" of their secondary structure representations. As many as 40% of the nucleotides of an RNA molecule can belong to these structural elements, which are distinct from the regular double helical regions formed by contiguous AU, GC, and GU Watson-Crick basepairs. With the large number of atomic- or near atomic-resolution 3D structures appearing in a steady stream in the PDB/NDB structure databases, the automated identification, extraction, comparison, clustering and visualization of these structural elements presents an opportunity to enhance RNA science. Three broad applications are: (1) identification of modular, autonomous structural units for RNA nanotechnology, nanobiology and synthetic biology applications; (2) bioinformatic analysis to improve RNA 3D structure prediction from sequence; and (3) creation of searchable databases for exploring the binding specificities, structural flexibility, and dynamics of these RNA elements. In this contribution, we review methods developed for computational extraction of hairpin and internal loop motifs from a non-redundant set of high-quality RNA 3D structures. We provide a statistical summary of the extracted hairpin and internal loop motifs in the most recent version of the RNA 3D Motif Atlas. We also explore the reliability and accuracy of the extraction process by examining its performance in clustering recurrent motifs from homologous ribosomal RNA (rRNA) structures. We conclude with a summary of remaining challenges, especially with regard to extraction of multi-helix junction motifs.

Synthesis of specifically deuterated nucleotides for NMR studies on RNA

Abstract We propose a strategy for NMR structure determination of RNA based on deuteration and use of specific labeling patterns. This strategy involves the use of NTPs that are deuterated in the ribose ring except for specific positions, e.g. H2', and that are either unlabeled or uniformly labeled in (13)C and (15)N in either the ribose or the base or both. Incorporation of these NTPs into an RNA sequence reduces both resonance line-width and spectral overlap. A limited number of combinations of these differently labeled NTPs in an RNA sequence suffices to obtain all relevant proton resonance assignments and structure parameters necessary for structure determination of larger systems (≫ 50 nucleotides). We describe the in vitro synthesis of the deuterated and/or (13)C/(15)N-labeled NTPs from glucose via separate enzymatic reactions. First, enzymes from the pentose-phosphate pathway efficiently convert glucose into ribose and enzymes from nucleotide biosynthesis and salvage pathways subsequently convert the ribose into nucleosides triphosphates (NTPs). The enzymes from the pentosephosphate pathway are all commercially available; the remaining enzymes have been purified from over-expressing strains. Separate enzymatic reactions were used to convert (2)H(7)- (13)C(6)-glucose into [1',3',4',5',5″-(2)H(5)-1',2',3',4',5',2,4,5,6-(13)C(9)-1,3-(15)N(2)]UTP and (2)H(7)-glucose into [1',3',4',5',5″-(2)H(5)]ATP, [1',3',4',5',5″-(2)H(5)]GTP, and [1',3',4',5',5″-(2)H(5)] CTP. The synthesis yields up to 1 gram of NTPs from 1 gram of glucose, which is about 5 to 10 times as efficient extraction for E. Coli grown on glucose. The synthesis presented here, is a modification of the method described by Tolbert & Williamson (1,2). (1)D and (2)D NMR spectra were acquired to demonstrate the utility of the new labeling patterns. The enzymatically synthesized NTPs were used in the synthesis of a 31-nucleotide RNA derived from the primer binding site of Hepatitis B virus genomic RNA to asses their efficiency in transcription.

An introduction to recurrent nucleotide interactions in RNA

RNA secondary structure diagrams familiar to molecular biologists summarize at a glance the folding of RNA chains to form Watson–Crick paired double helices. However, they can be misleading: First of all, they imply that the nucleotides in loops and linker segments, which can amount to 35% to 50% of a structured RNA, do not significantly interact with other nucleotides. Secondly, they give the impression that RNA molecules are loosely organized in three-dimensional (3D) space. In fact, structured RNAs are compactly folded as a result of numerous long-range, sequence-specific interactions, many of which involve loop or linker nucleotides. Here, we provide an introduction for students and researchers of RNA on the types, prevalence, and sequence variations of inter-nucleotide interactions that structure and stabilize RNA 3D motifs and architectures, using Escherichia coli (E. coli) 16S ribosomal RNA as a concrete example. The picture that emerges is that almost all nucleotides in structured RNA molecules, including those in nominally single-stranded loop or linker regions, form specific interactions that stabilize functional structures or mediate interactions with other molecules. The small number of noninteracting, ‘looped-out’ nucleotides make it possible for the RNA chain to form sharp turns. Base-pairing is the most specific interaction in RNA as it involves edge-to-edge hydrogen bonding (H-bonding) of the bases. Non-Watson–Crick base pairs are a significant fraction (30% or more) of base pairs in structured RNAs.

Formation of Tertiary Interactions during rRNA GTPase Center Folding

The 60-nt GTPase center (GAC) of 23S rRNA has a phylogenetically conserved secondary structure with two hairpin loops and a 3-way junction. It folds into an intricate tertiary structure upon addition of Mg(2+) ions, which is stabilized by the L11 protein in cocrystal structures. Here, we monitor the kinetics of its tertiary folding and Mg(2+)-dependent intermediate states by observing selected nucleobases that contribute specific interactions to the GAC tertiary structure in the cocrystals. The fluorescent nucleobase 2-aminopurine replaced three individual adenines, two of which make long-range stacking interactions and one that also forms hydrogen bonds. Each site reveals a unique response to Mg(2+) addition and temperature, reflecting its environmental change from secondary to tertiary structure. Stopped-flow fluorescence experiments revealed that kinetics of tertiary structure formation upon addition of MgCl2 are also site specific, with local conformational changes occurring from 5 ms to 4s and with global folding from 1 to 5s. Site-specific substitution with (15)N-nucleobases allowed observation of stable hydrogen bond formation by NMR experiments. Equilibrium titration experiments indicate that a stable folding intermediate is present at stoichiometric concentrations of Mg(2+) and suggest that there are two initial sites of Mg(2+) ion association.

Structural fidelity and NMR relaxation analysis in a prototype RNA hairpin

RNA hairpins are widespread and very stable motifs that contribute decisively to RNA folding and biological function. The GTP1G2C3A4C5U6U7C8G9G10U11G12C13C14 construct (with a central UUCG tetraloop) has been extensively studied by solution NMR, and offers and excellent opportunity to evaluate the structure and dynamical description afforded by molecular dynamics (MD) simulations. Here, we compare average structural parameters and NMR relaxation rates estimated from a series of multiple independent explicit solvent MD simulations using the two most recent RNA AMBER force fields (ff99 and ff10). Predicted overall tumbling times are ∼20% faster than those inferred from analysis of NMR data and follow the same trend when temperature and ionic strength is varied. The Watson-Crick stem and the "canonical" UUCG loop structure are maintained in most simulations including the characteristic syn conformation along the glycosidic bond of G9, although some key hydrogen bonds in the loop are partially disrupted. Our analysis pinpoints G9-G10 backbone conformations as a locus of discrepancies between experiment and simulation. In general the results for the more recent force-field parameters (ff10) are closer to experiment than those for the older ones (ff99). This work provides a comprehensive and detailed comparison of state of the art MD simulations against a wide variety of solution NMR measurements.

Specific nucleotides at the 3'-terminal promoter of viral hemorrhagic septicemia virus are important for virulence in vitro and in vivo

Viral hemorrhagic septicemia virus (VHSV), a member of the Novirhabdovirus genus, contains an 11-nucleotide conserved sequence at the terminal 3'- and 5'-untranslated regions (UTRs) that are complementary. To study the importance of nucleotides in the 3'-UTR of VHSV for replication of novirhabdoviruses, we performed site-directed mutagenesis of selected residues at the 3'-terminus and generated mutant viruses using a reverse genetics approach. Assessment of growth kinetics and in vitro real-time cytopathogenicity studies showed that the order of two nucleotides (A4G5) of the 3'-terminus of VHSV directly affects growth kinetics in vitro. The mutant A4G-G5A virus has reduced total positive-strand RNA synthesis efficiency (51% of wild-type) at 48h post-transfection and 70h delay in causing complete cytopathic effect in susceptible fish cells, as compared to the WT-VHSV. Furthermore, when the A4G-G5A virus was used to challenge zebrafish, it exhibited reduced pathogenicity (54% lower end-point mortality) compared to the WT-VHSV. From these studies, we infer that specific residues in the 3'-UTR of VHSV have a promoter function and are essential to modulate the virulence in cells and pathogenicity in fish.

A geometrical parametrization of C1'-C5' RNA ribose chemical shifts calculated by density functional theory

It has been recently shown that NMR chemical shifts can be used to determine the structures of proteins. In order to begin to extend this type of approach to nucleic acids, we present an equation that relates the structural parameters and the (13)C chemical shifts of the ribose group. The parameters in the equation were determined by maximizing the agreement between the DFT-derived chemical shifts and those predicted through the equation for a database of ribose structures. Our results indicate that this type of approach represents a promising way of establishing quantitative and computationally efficient analytical relationships between chemical shifts and structural parameters in nucleic acids.

Structure modulation of helix 69 from Escherichia coli 23S ribosomal RNA by pseudouridylations

Helix 69 (H69) is a 19-nt stem-loop region from the large subunit ribosomal RNA. Three pseudouridine (Ψ) modifications clustered in H69 are conserved across phylogeny and known to affect ribosome function. To explore the effects of Ψ on the conformations of Escherichia coli H69 in solution, nuclear magnetic resonance spectroscopy was used to reveal the structural differences between H69 with (ΨΨΨ) and without (UUU) Ψ modifications. Comparison of the two structures shows that H69 ΨΨΨ has the following unique features: (i) the loop region is closed by a Watson-Crick base pair between Ψ1911 and A1919, which is potentially reinforced by interactions involving Ψ1911N1H and (ii) Ψ modifications at loop residues 1915 and 1917 promote base stacking from Ψ1915 to A1918. In contrast, the H69 UUU loop region, which lacks Ψ modifications, is less organized. Structure modulation by Ψ leads to alteration in conformational behavior of the 5' half of the H69 loop region, observed as broadening of C1914 non-exchangeable base proton resonances in the H69 ΨΨΨ nuclear magnetic resonance spectra, and plays an important biological role in establishing the ribosomal intersubunit bridge B2a and mediating translational fidelity.

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.

Recognition modes of RNA tetraloops and tetraloop-like motifs by RNA-binding proteins

RNA hairpins are the most commonly occurring secondary structural elements in RNAs and serve as nucleation sites for RNA folding, RNA-RNA, and RNA-protein interactions. RNA hairpins are frequently capped by tetraloops, and based on sequence similarity, three broad classes of RNA tetraloops have been defined: GNRA, UNCG, and CUYG. Other classes such as the UYUN tetraloop in histone mRNAs, the UGAA in 16S rRNA, the AUUA tetraloop from the MS2 bacteriophage, and the AGNN tetraloop that binds RNase III have also been characterized. The tetraloop structure is compact and is usually characterized by a paired interaction between the first and fourth nucleotides. The two unpaired nucleotides in the loop are usually involved in base-stacking or base-phosphate hydrogen bonding interactions. Several structures of RNA tetraloops, free and complexed to other RNAs or proteins, are now available and these studies have increased our understanding of the diverse mechanisms by which this motif is recognized. RNA tetraloops can mediate RNA-RNA contacts via the tetraloop-receptor motif, kissing hairpin loops, A-minor interactions, and pseudoknots. While these RNA-RNA interactions are fairly well understood, how RNA-binding proteins recognize RNA tetraloops and tetraloop-like motifs remains unclear. In this review, we summarize the structures of RNA tetraloop-protein complexes and the general themes that have emerged on sequence- and structure-specific recognition of RNA tetraloops. We highlight how proteins achieve molecular recognition of this nucleic acid motif, the structural adaptations observed in the tetraloop to accommodate the protein-binding partner, and the role of dynamics in recognition.

Backbone-base interactions critical to quantum stabilization of transfer RNA anticodon structure

Transfer RNA (tRNA) anticodons adopt a highly ordered 3'-stack without significant base overlap. Density functional theory at the M06-2X/6-31+G(d,p) level in combination with natural bond orbital analysis was utilized to calculate the intramolecular interactions within the tRNA anticodon that are responsible for stabilizing the stair-stepped conformation. Ten tRNA X-ray crystal structures were obtained from the PDB databank and were trimmed to include only the anticodon bases. Hydrogenic positions were added and optimized for the structures in the stair-stepped conformation. The sugar-phosphate backbone has been retained for these calculations, revealing the role it plays in RNA structural stability. It was found that electrostatic interactions between the sugar-phosphate backbone and the base provide the most stability, rather than the traditionally studied interbase stacking. Base-stacking interactions, though present, were weak and inconsistent. Aqueous solvation was found to have little effect on the intramolecular interactions.

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.

Structure of the K-turn U4 RNA: a combined NMR and SANS study

K-turn motifs are universal RNA structural elements providing a binding platform for proteins in several cellular contexts. Their characteristic is a sharp kink in the phosphate backbone that puts the two helical stems of the protein-bound RNA at an angle of 60°. However, to date no high-resolution structure of a naked K-turn motif is available. Here, we present the first structural investigation at atomic resolution of an unbound K-turn RNA (the spliceosomal U4-Kt RNA) by a combination of NMR and small-angle neutron scattering data. With this study, we wish to address the question whether the K-turn structural motif assumes the sharply kinked conformation in the absence of protein binders and divalent cations. Previous studies have addressed this question by fluorescence resonance energy transfer, biochemical assays and molecular dynamics simulations, suggesting that the K-turn RNAs exist in equilibrium between a kinked conformation, which is competent for protein binding, and a more extended conformation, with the population distribution depending on the concentration of divalent cations. Our data shows that the U4-Kt RNA predominantly assumes the more extended conformation in the absence of proteins and divalent cations. The internal loop region is well structured but adopts a different conformation from the one observed in complex with proteins. Our data suggests that the K-turn consensus sequence does not per se code for the kinked conformation; instead the sharp backbone kink requires to be stabilized by protein binders.