The sarcin/ricin loop, a modular RNA

Antifungal Activity of Ageritin, a Ribotoxin-like Protein from Cyclocybe aegerita Edible Mushroom, against Phytopathogenic Fungi

Ageritin from poplar mushrooms is a specific endonuclease that hydrolyzes a single phosphodiester bond located in the sarcin-ricin loop (SRL) of the large rRNA, thereby blocking protein synthesis. Considering the possible biotechnological use of this enzyme, here we report its antifungal activity against virulent fungi affecting crops of economic interest. Our results show that ageritin (200 µg/plug; ~13.5 nmole) inhibits the growth of Botrytis cinerea (57%), Colletotrichum truncatum (42%), and Alternaria alternata (57%), when tested on potato dextrose agar plates. At the same time, no effect was observed against Trichoderma harzianum (a fungus promoting beneficial effects in plants). To verify whether the antifungal action of ageritin against B. cinerea and T. harzianum was due to ribosome damage, we tested ageritin in vitro on partially isolated B. cinerea and T. harzianum ribosomes. Interestingly, ageritin was able to release the Endo's fragment from both tested fungal ribosomes. We therefore decided to test the antifungal effect of ageritin on B. cinerea and T. harzianum using a different growth condition (liquid medium). Differently from the result in solid medium, ageritin can inhibit both B. cinerea and T. harzianum fungal growth in liquid medium in a concentration-dependent manner up to 35.7% and 38.7%, respectively, at the highest concentration tested (~200 µg/mL; 12 µM), and the analysis of RNA isolated from ageritin-treated cells revealed the presence of Endo's fragment, highlighting its ability to cross the fungal cell wall and reach the ribosomes. Overall, these data highlight that the efficacy of antifungal treatment to prevent or treat a potential fungal disease may depend not only on the fungal species but also on the conditions of toxin application.

Analysis of the Sequence Preference of Saporin by Deep Sequencing

Ribosome-inactivating proteins (RIPs) are RNA:adenosine glycosidases that inactivate eukaryotic ribosomes by depurinating the sarcin-ricin loop (SRL) in 28S rRNA. The GAGA sequence at the top of the SRL or at the top of a hairpin loop is assumed to be their target motif. Saporin is a RIP widely used to develop immunotoxins for research and medical applications, but its sequence specificity has not been investigated. Here, we combine the conventional aniline cleavage assay for depurinated nucleic acids with high-throughput sequencing to study sequence-specific depurination of oligonucleotides caused by saporin. Our data reveal the sequence preference of saporin for different substrates and show that the GAGA motif is not efficiently targeted by this protein, neither in RNA nor in DNA. Instead, a preference of saporin for certain hairpin DNAs was observed. The observed sequence-specific activity of saporin may be relevant to antiviral or apoptosis-inducing effects of RIPs. The developed method could also be useful for studying the sequence specificity of depurination by other RIPs or enzymes.

Context-sensitivity of isosteric substitutions of non-Watson-Crick basepairs in recurrent RNA 3D motifs

Sequence variation in a widespread, recurrent, structured RNA 3D motif, the Sarcin/Ricin (S/R), was studied to address three related questions: First, how do the stabilities of structured RNA 3D motifs, composed of non-Watson–Crick (non-WC) basepairs, compare to WC-paired helices of similar length and sequence? Second, what are the effects on the stabilities of such motifs of isosteric and non-isosteric base substitutions in the non-WC pairs? And third, is there selection for particular base combinations in non-WC basepairs, depending on the temperature regime to which an organism adapts? A survey of large and small subunit rRNAs from organisms adapted to different temperatures revealed the presence of systematic sequence variations at many non-WC paired sites of S/R motifs. UV melting analysis and enzymatic digestion assays of oligonucleotides containing the motif suggest that more stable motifs tend to be more rigid. We further found that the base substitutions at non-Watson–Crick pairing sites can significantly affect the thermodynamic stabilities of S/R motifs and these effects are highly context specific indicating the importance of base-stacking and base-phosphate interactions on motif stability. This study highlights the significance of non-canonical base pairs and their contributions to modulating the stability and flexibility of RNA molecules.

The Arabidopsis 2'-O-Ribose-Methylation and Pseudouridylation Landscape of rRNA in Comparison to Human and Yeast

Eukaryotic ribosome assembly starts in the nucleolus, where the ribosomal DNA (rDNA) is transcribed into the 35S pre-ribosomal RNA (pre-rRNA). More than two-hundred ribosome biogenesis factors (RBFs) and more than two-hundred small nucleolar RNAs (snoRNA) catalyze the processing, folding and modification of the rRNA in Arabidopsis thaliana. The initial pre-ribosomal 90S complex is formed already during transcription by association of ribosomal proteins (RPs) and RBFs. In addition, small nucleolar ribonucleoprotein particles (snoRNPs) composed of snoRNAs and RBFs catalyze the two major rRNA modification types, 2'-O-ribose-methylation and pseudouridylation. Besides these two modifications, rRNAs can also undergo base methylations and acetylation. However, the latter two modifications have not yet been systematically explored in plants. The snoRNAs of these snoRNPs serve as targeting factors to direct modifications to specific rRNA regions by antisense elements. Today, hundreds of different sites of modifications in the rRNA have been described for eukaryotic ribosomes in general. While our understanding of the general process of ribosome biogenesis has advanced rapidly, the diversities appearing during plant ribosome biogenesis is beginning to emerge. Today, more than two-hundred RBFs were identified by bioinformatics or biochemical approaches, including several plant specific factors. Similarly, more than two hundred snoRNA were predicted based on RNA sequencing experiments. Here, we discuss the predicted and verified rRNA modification sites and the corresponding identified snoRNAs on the example of the model plant Arabidopsis thaliana. Our summary uncovers the plant modification sites in comparison to the human and yeast modification sites.

Characteristic chemical probing patterns of loop motifs improve prediction accuracy of RNA secondary structures

RNA structures play a fundamental role in nearly every aspect of cellular physiology and pathology. Gaining insights into the functions of RNA molecules requires accurate predictions of RNA secondary structures. However, the existing thermodynamic folding models remain less accurate than desired, even when chemical probing data, such as selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) reactivities, are used as restraints. Unlike most SHAPE-directed algorithms that only consider SHAPE restraints for base pairing, we extract two-dimensional structural features encoded in SHAPE data and establish robust relationships between characteristic SHAPE patterns and loop motifs of various types (hairpin, internal, and bulge) and lengths (2-11 nucleotides). Such characteristic SHAPE patterns are closely related to the sugar pucker conformations of loop residues. Based on these patterns, we propose a computational method, SHAPELoop, which refines the predicted results of the existing methods, thereby further improving their prediction accuracy. In addition, SHAPELoop can provide information about local or global structural rearrangements (including pseudoknots) and help researchers to easily test their hypothesized secondary structures.

Antifungal resistance-modifying multiplexing action of Momordica charantia protein and phosphorylated derivatives on the basis of growth-dependent gene coregulation in Candida albicans

Fungal growth-dependent gene coregulation is strongly implicated in alteration of gene-encoding target proteases ruling with an antifungal resistance niche and biology of resistant mutants. On the basis of multi-alterative processes in this platform, the resistance-modifying strategy is designed in ketoconazole resistant Candida albicans and evaluated with less selective Momordica charantia protein and allosterically phosphorylated derivatives at the Thr102, Thr24 and Thr255 sites, respectively. We demonstrate absolutely chemo-sensitizing efficacy regarding stepwise-modifying resistance in sensitivity, by a load of only 26.23-40.00 μg/l agents in Sabouraud's dextrose broth. Five successive modifying-steps realize the decreasing of ketoconazole E-test MIC50 from 11.10 to a lower level than 0.10 mg/l. With the ketoconazole resistance-modifying, colony undergoes a high-frequency morphological switch between high ploidy (opaque) and small budding haploid (white). A cellular event in the first modifying-step associates with relatively slow exponential growth (ie, a 4-h delay)-dependent action, mediated by agents adsorption. Moreover, multiple molecular roles are coupled with intracellularly and extracellularly binding to ATP-dependent RNA helicase dbp6; the 0.08-2.45 fold upregulation of TATA-box-binding protein, rRNA-processing protein and translation initiation factor 5A; and the 7.52-55.33% decrease of cytochrome P450 lanosterol 14α-demethylase, glucan 1, 3-β glucosidase, candidapepsin-1 and 1-acylglycerol-3-phosphate O-acyltransferase. Spatial and temporal gene coregulation, in the transcription and translation initiation stages with rRNA-processing, is a new coprocessing platform enabling target protease attenuations for resistance-impairing. An updated resistance-modifying measure of these agents in the low-dose antifungal strategic design may provide opportunities to a virtually safe therapy that is in high dose-dependency.

LAY SUMMARY

A new platform to modify resistance is fungal growth-dependent gene coregulation. MAP30 and phosphorylated derivatives are candidate resistance-modifying agents. Low-dose stepwise treatment absolutely modifies azole resistance in model fungus.

Mayahuelin, a Type I Ribosome Inactivating Protein: Characterization, Evolution, and Utilization in Phylogenetic Analyses of Agave

Agaves resist extreme heat and drought. In A. tequilana var. azul, the central spike of the rosette -containing the shoot apical meristem and folded leaves in early stages of development- is remarkably heat tolerant. We found that the most abundant protein in this organ is a 27 kDa protein. This protein was named mayahuelin to honor Mayáhuel, the agave goddess in the Aztec pantheon. LC-MS/MS analyses identified mayahuelin as a type I RIP (Ribosome Inactivating Protein). In addition to the spike, mayahuelin was expressed in the peduncle and in seeds, whereas in mature leaves, anthers, filaments, pistils, and tepals was absent. Anti-mayahuelin antibody raised against the A. tequilana var. azul protein revealed strong signals in spike leaves of A. angustifolia, A. bracteosa, A. rhodacantha, and A. vilmoriniana, and moderate signals in A. isthmensis, A. kerchovei, A. striata ssp. falcata, and A. titanota, indicating conservation at the protein level throughout the Agave genus. As in charybdin, a type I RIP characterized in Drimia maritima, mayahuelin from A. tequilana var. azul contains a natural aa substitution (Y76D) in one out of four aa comprising the active site. The RIP gene family in A. tequilana var. azul consists of at least 12 genes and Mayahuelin is the only member encoding active site substitutions. Unlike canonical plant RIPs, expression of Mayahuelin gene in S. cerevisiae did not compromise growth. The inhibitory activity of the purified protein on a wheat germ in vitro translation system was moderate. Mayahuelin orthologs from other Agave species displayed one of six alleles at Y76: (Y/Y, D/D, S/S, Y/D, Y/S, D/S) and proved to be useful markers for phylogenetic analysis. Homozygous alleles were more frequent in wild accessions whereas heterozygous alleles were more frequent in cultivars. Mayahuelin sequences from different wild populations of A. angustifolia and A. rhodacantha allowed the identification of accessions closely related to azul, manso, sigüín, mano larga, and bermejo varieties of A. tequilana and var. espadín of A. angustifolia. Four A. rhodacantha accessions and A. angustifolia var. espadín were closer relatives of A. tequilana var. azul than A. angustifolia wild accessions or other A. tequilana varieties.

How Ricin Damages the Ribosome

Ricin belongs to the group of ribosome-inactivating proteins (RIPs), i.e., toxins that have evolved to provide particular species with an advantage over other competitors in nature. Ricin possesses RNA N-glycosidase activity enabling the toxin to eliminate a single adenine base from the sarcin-ricin RNA loop (SRL), which is a highly conserved structure present on the large ribosomal subunit in all species from the three domains of life. The SRL belongs to the GTPase associated center (GAC), i.e., a ribosomal element involved in conferring unidirectional trajectory for the translational apparatus at the expense of GTP hydrolysis by translational GTPases (trGTPases). The SRL represents a critical element in the GAC, being the main triggering factor of GTP hydrolysis by trGTPases. Enzymatic removal of a single adenine base at the tip of SRL by ricin blocks GTP hydrolysis and, at the same time, impedes functioning of the translational machinery. Here, we discuss the consequences of SRL depurination by ricin for ribosomal performance, with emphasis on the mechanistic model overview of the SRL modus operandi.

Functional analysis of RIP toxins from the Drosophila endosymbiont Spiroplasma poulsonii

Background

Insects frequently live in close relationship with symbiotic bacteria that carry out beneficial functions for their host, like protection against parasites and viruses. However, in some cases, the mutualistic nature of such associations is put into question because of detrimental phenotypes caused by the symbiont. One example is the association between the vertically transmitted facultative endosymbiont Spiroplasma poulsonii and its natural host Drosophila melanogaster. Whereas S. poulsonii protects its host against parasitoid wasps and nematodes by the action of toxins from the family of Ribosome Inactivating Proteins (RIPs), the presence of S. poulsonii has been reported to reduce host’s life span and to kill male embryos by a toxin called Spaid. In this work, we investigate the harmful effects of Spiroplasma RIPs on Drosophila in the absence of parasite infection.

Results

We show that only two Spiroplasma RIPs (SpRIP1 and SpRIP2) among the five RIP genes encoded in the S. poulsonii genome are significantly expressed during the whole Drosophila life cycle. Heterologous expression of SpRIP1 and 2 in uninfected flies confirms their toxicity, as indicated by a reduction of Drosophila lifespan and hemocyte number. We also show that RIPs can cause the death of some embryos, including females.

Conclusion

Our results indicate that RIPs released by S. poulsonii contribute to the reduction of host lifespan and embryo mortality. This suggests that SpRIPs may impact the insect-symbiont homeostasis beyond their protective function against parasites.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1410-1) contains supplementary material, which is available to authorized users.

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.

A pseudouridylation switch in rRNA is implicated in ribosome function during the life cycle of Trypanosoma brucei

The protozoan parasite Trypanosoma brucei, which causes devastating diseases in humans and animals in sub-Saharan Africa, undergoes a complex life cycle between the mammalian host and the blood-feeding tsetse fly vector. However, little is known about how the parasite performs most molecular functions in such different environments. Here, we provide evidence for the intriguing possibility that pseudouridylation of rRNA plays an important role in the capacity of the parasite to transit between the insect midgut and the mammalian bloodstream. Briefly, we mapped pseudouridines (Ψ) on rRNA by Ψ-seq in procyclic form (PCF) and bloodstream form (BSF) trypanosomes. We detected 68 Ψs on rRNA, which are guided by H/ACA small nucleolar RNAs (snoRNA). The small RNome of both life cycle stages was determined by HiSeq and 83 H/ACAs were identified. We observed an elevation of 21 Ψs modifications in BSF as a result of increased levels of the guiding snoRNAs. Overexpression of snoRNAs guiding modification on H69 provided a slight growth advantage to PCF parasites at 30 °C. Interestingly, these modifications are predicted to significantly alter the secondary structure of the large subunit (LSU) rRNA suggesting that hypermodified positions may contribute to the adaption of ribosome function during cycling between the two hosts.

RAG-3D: a search tool for RNA 3D substructures

To address many challenges in RNA structure/function prediction, the characterization of RNA's modular architectural units is required. Using the RNA-As-Graphs (RAG) database, we have previously explored the existence of secondary structure (2D) submotifs within larger RNA structures. Here we present RAG-3D—a dataset of RNA tertiary (3D) structures and substructures plus a web-based search tool—designed to exploit graph representations of RNAs for the goal of searching for similar 3D structural fragments. The objects in RAG-3D consist of 3D structures translated into 3D graphs, cataloged based on the connectivity between their secondary structure elements. Each graph is additionally described in terms of its subgraph building blocks. The RAG-3D search tool then compares a query RNA 3D structure to those in the database to obtain structurally similar structures and substructures. This comparison reveals conserved 3D RNA features and thus may suggest functional connections. Though RNA search programs based on similarity in sequence, 2D, and/or 3D structural elements are available, our graph-based search tool may be advantageous for illuminating similarities that are not obvious; using motifs rather than sequence space also reduces search times considerably. Ultimately, such substructuring could be useful for RNA 3D structure prediction, structure/function inference and inverse folding.

RNAMotifScanX: a graph alignment approach for RNA structural motif identification

RNA structural motifs are recurrent three-dimensional (3D) components found in the RNA architecture. These RNA structural motifs play important structural or functional roles and usually exhibit highly conserved 3D geometries and base-interaction patterns. Analysis of the RNA 3D structures and elucidation of their molecular functions heavily rely on efficient and accurate identification of these motifs. However, efficient RNA structural motif search tools are lacking due to the high complexity of these motifs. In this work, we present RNAMotifScanX, a motif search tool based on a base-interaction graph alignment algorithm. This novel algorithm enables automatic identification of both partially and fully matched motif instances. RNAMotifScanX considers noncanonical base-pairing interactions, base-stacking interactions, and sequence conservation of the motifs, which leads to significantly improved sensitivity and specificity as compared with other state-of-the-art search tools. RNAMotifScanX also adopts a carefully designed branch-and-bound technique, which enables ultra-fast search of large kink-turn motifs against a 23S rRNA. The software package RNAMotifScanX is implemented using GNU C++, and is freely available from http://genome.ucf.edu/RNAMotifScanX.

Rfam 12.0: updates to the RNA families database

The Rfam database (available at http://rfam.xfam.org) is a collection of non-coding RNA families represented by manually curated sequence alignments, consensus secondary structures and annotation gathered from corresponding Wikipedia, taxonomy and ontology resources. In this article, we detail updates and improvements to the Rfam data and website for the Rfam 12.0 release. We describe the upgrade of our search pipeline to use Infernal 1.1 and demonstrate its improved homology detection ability by comparison with the previous version. The new pipeline is easier for users to apply to their own data sets, and we illustrate its ability to annotate RNAs in genomic and metagenomic data sets of various sizes. Rfam has been expanded to include 260 new families, including the well-studied large subunit ribosomal RNA family, and for the first time includes information on short sequence- and structure-based RNA motifs present within families.

Chemical probing of RNA with the hydroxyl radical at single-atom resolution

While hydroxyl radical cleavage is widely used to map RNA tertiary structure, lack of mechanistic understanding of strand break formation limits the degree of structural insight that can be obtained from this experiment. Here, we determine how individual ribose hydrogens of sarcin/ricin loop RNA participate in strand cleavage. We find that substituting deuterium for hydrogen at a ribose 5'-carbon produces a kinetic isotope effect on cleavage; the major cleavage product is an RNA strand terminated by a 5'-aldehyde. We conclude that hydroxyl radical abstracts a 5'-hydrogen atom, leading to RNA strand cleavage. We used this approach to obtain structural information for a GUA base triple, a common tertiary structural feature of RNA. Cleavage at U exhibits a large 5' deuterium kinetic isotope effect, a potential signature of a base triple. Others had noted a ribose-phosphate hydrogen bond involving the G 2'-OH and the U phosphate of the GUA triple, and suggested that this hydrogen bond contributes to backbone rigidity. Substituting deoxyguanosine for G, to eliminate this hydrogen bond, results in a substantial decrease in cleavage at G and U of the triple. We conclude that this hydrogen bond is a linchpin of backbone structure around the triple.

Computational identification of RNA functional determinants by three-dimensional quantitative structure-activity relationships

Anti-infection drugs target vital functions of infectious agents, including their ribosome and other essential non-coding RNAs. One of the reasons infectious agents become resistant to drugs is due to mutations that eliminate drug-binding affinity while maintaining vital elements. Identifying these elements is based on the determination of viable and lethal mutants and associated structures. However, determining the structure of enough mutants at high resolution is not always possible. Here, we introduce a new computational method, MC-3DQSAR, to determine the vital elements of target RNA structure from mutagenesis and available high-resolution data. We applied the method to further characterize the structural determinants of the bacterial 23S ribosomal RNA sarcin–ricin loop (SRL), as well as those of the lead-activated and hammerhead ribozymes. The method was accurate in confirming experimentally determined essential structural elements and predicting the viability of new SRL variants, which were either observed in bacteria or validated in bacterial growth assays. Our results indicate that MC-3DQSAR could be used systematically to evaluate the drug-target potentials of any RNA sites using current high-resolution structural data.

RNA structure and dynamics: a base pairing perspective

RNA is now known to possess various structural, regulatory and enzymatic functions for survival of cellular organisms. Functional RNA structures are generally created by three-dimensional organization of small structural motifs, formed by base pairing between self-complementary sequences from different parts of the RNA chain. In addition to the canonical Watson-Crick or wobble base pairs, several non-canonical base pairs are found to be crucial to the structural organization of RNA molecules. They appear within different structural motifs and are found to stabilize the molecule through long-range intra-molecular interactions between basic structural motifs like double helices and loops. These base pairs also impart functional variation to the minor groove of A-form RNA helices, thus forming anchoring site for metabolites and ligands. Non-canonical base pairs are formed by edge-to-edge hydrogen bonding interactions between the bases. A large number of theoretical studies have been done to detect and analyze these non-canonical base pairs within crystal or NMR derived structures of different functional RNA. Theoretical studies of these isolated base pairs using ab initio quantum chemical methods as well as molecular dynamics simulations of larger fragments have also established that many of these non-canonical base pairs are as stable as the canonical Watson-Crick base pairs. This review focuses on the various structural aspects of non-canonical base pairs in the organization of RNA molecules and the possible applications of these base pairs in predicting RNA structures with more accuracy.

The nuclear magnetic resonance of CCCC RNA reveals a right-handed helix, and revised parameters for AMBER force field torsions improve structural predictions from molecular dynamics

The sequence dependence of RNA energetics is important for predicting RNA structure. Hairpins with C(n) loops are consistently less stable than hairpins with other loops, which suggests the structure of C(n) regions could be unusual in the "unfolded" state. For example, previous nuclear magnetic resonance (NMR) evidence suggested that polycytidylic acid forms a left-handed helix. In this study, UV melting experiments show that the hairpin formed by r(5'GGACCCCCGUCC) is less stable than r(5'GGACUUUUGUCC). NMR spectra for single-stranded C(4) oligonucleotide, mimicking the unfolded hairpin loop, are consistent with a right-handed A-form-like helix. Comparisons between NMR spectra and molecular dynamics (MD) simulations suggest that recent reparametrizations, parm99χ_YIL and parm99TOR, of the AMBER parm99 force field improve the agreement between structural features for C(4) determined by NMR and predicted by MD. Evidently, the force field revisions to parm99 improve the modeling of RNA energetics and therefore structure.

Novel conformation of an RNA structural switch

The RNA duplex, (5′GACGAGUGUCA)₂, has two conformations in equilibrium. The nuclear magnetic resonance solution structure reveals that the major conformation of the loop, 5′GAGU/3′UGAG, is novel and contains two unusual Watson–Crick/Hoogsteen GG pairs with G residues in the syn conformation, two A residues stacked on each other in the center of the helix with inverted sugars, and two bulged-out U residues. The structure provides a benchmark for testing approaches for predicting local RNA structure and a sequence that allows the design of a unique arrangement of functional groups and/or a conformational switch into nucleic acids.

A procedure to validate and correct the 13C chemical shift calibration of RNA datasets

Chemical shifts reflect the structural environment of a certain nucleus and can be used to extract structural and dynamic information. Proper calibration is indispensable to extract such information from chemical shifts. Whereas a variety of procedures exist to verify the chemical shift calibration for proteins, no such procedure is available for RNAs to date. We present here a procedure to analyze and correct the calibration of (13)C NMR data of RNAs. Our procedure uses five (13)C chemical shifts as a reference, each of them found in a narrow shift range in most datasets deposited in the Biological Magnetic Resonance Bank. In 49 datasets we could evaluate the (13)C calibration and detect errors or inconsistencies in RNA (13)C chemical shifts based on these chemical shift reference values. More than half of the datasets (27 out of those 49) were found to be improperly referenced or contained inconsistencies. This large inconsistency rate possibly explains that no clear structure-(13)C chemical shift relationship has emerged for RNA so far. We were able to recalibrate or correct 17 datasets resulting in 39 usable (13)C datasets. 6 new datasets from our lab were used to verify our method increasing the database to 45 usable datasets. We can now search for structure-chemical shift relationships with this improved list of (13)C chemical shift data. This is demonstrated by a clear relationship between ribose (13)C shifts and the sugar pucker, which can be used to predict a C2'- or C3'-endo conformation of the ribose with high accuracy. The improved quality of the chemical shift data allows statistical analysis with the potential to facilitate assignment procedures, and the extraction of restraints for structure calculations of RNA.

Identification of a minimal region of the HIV-1 5'-leader required for RNA dimerization, NC binding, and packaging

Assembly of human immunodeficiency virus type 1 (HIV-1) particles is initiated in the cytoplasm by the formation of a ribonucleoprotein complex comprising the dimeric RNA genome and a small number of viral Gag polyproteins. Genomes are recognized by the nucleocapsid (NC) domains of Gag, which interact with packaging elements believed to be located primarily within the 5'-leader (5'-L) of the viral RNA. Recent studies revealed that the native 5'-L exists as an equilibrium of two conformers, one in which dimer-promoting residues and NC binding sites are sequestered and packaging is attenuated, and one in which these sites are exposed and packaging is promoted. To identify the elements within the dimeric 5'-L that are important for packaging, we generated HIV-1 5'-L RNAs containing mutations and deletions designed to eliminate substructures without perturbing the overall structure of the leader and examined effects of the mutations on RNA dimerization, NC binding, and packaging. Our findings identify a 159-residue RNA packaging signal that possesses dimerization and NC binding properties similar to those of the intact 5'-L and contains elements required for efficient RNA packaging.

Interaction of ricin and Shiga toxins with ribosomes

Ricin and Shiga toxins designated as ribosome inactivating proteins (RIPs) are RNA N-glycosidases that depurinate a specific adenine (A₄₃₂₄ in rat 28S rRNA) in the conserved α-sarcin/ricin loop of the large rRNA, inhibiting protein synthesis. Evidence obtained from a number of studies suggests that interaction with ribosomal proteins plays an important role in the catalytic activity and ribosome specificity of RIPs. This review summarizes the recent developments in identification of the ribosomal proteins that interact with ricin and Shiga toxins and the principles governing these interactions.

Comprehensive survey and geometric classification of base triples in RNA structures

Base triples are recurrent clusters of three RNA nucleobases interacting edge-to-edge by hydrogen bonding. We find that the central base in almost all triples forms base pairs with the other two bases of the triple, providing a natural way to geometrically classify base triples. Given 12 geometric base pair families defined by the Leontis-Westhof nomenclature, combinatoric enumeration predicts 108 potential geometric base triple families. We searched representative atomic-resolution RNA 3D structures and found instances of 68 of the 108 predicted base triple families. Model building suggests that some of the remaining 40 families may be unlikely to form for steric reasons. We developed an on-line resource that provides exemplars of all base triples observed in the structure database and models for unobserved, predicted triples, grouped by triple family, as well as by three-base combination (http://rna.bgsu.edu/Triples). The classification helps to identify recurrent triple motifs that can substitute for each other while conserving RNA 3D structure, with applications in RNA 3D structure prediction and analysis of RNA sequence evolution.

Benchmarking AMBER force fields for RNA: comparisons to NMR spectra for single-stranded r(GACC) are improved by revised χ torsions

Accurately modeling unpaired regions of RNA is important for predicting structure, dynamics, and thermodynamics of folded RNA. Comparisons between NMR data and molecular dynamics simulations provide a test of force fields used for modeling. Here, NMR spectroscopy, including NOESY, (1)H-(31)P HETCOR, DQF-COSY, and TOCSY, was used to determine conformational preferences for single-stranded GACC RNA. The spectra are consistent with a conformational ensemble containing major and minor A-form-like structures. In a series of 50 ns molecular dynamics (MD) simulations with the AMBER99 force field in explicit solvent, initial A-form-like structures rapidly evolve to disordered conformations. A set of 50 ns simulations with revised χ torsions (AMBER99χ force field) gives two primary conformations, consistent with the NMR spectra. A single 1.9 μs MD simulation with the AMBER99χ force field showed that the major and minor conformations are retained for almost 68% of the time in the first 700 ns, with multiple transformations from A-form to non-A-form conformations. For the rest of the simulation, random-coil structures and a stable non-A-form conformation inconsistent with NMR spectra were seen. Evidently, the AMBER99χ force field improves structural predictions for single-stranded GACC RNA compared to the AMBER99 force field, but further force field improvements are needed.

The ribotoxin restrictocin recognizes its RNA substrate by selective engagement of active site residues

Restrictocin and related fungal endoribonucleases from the α-sarcin family site-specifically cleave the sarcin/ricin loop (SRL) on the ribosome to inhibit translation and ultimately trigger cell death. Previous studies showed that the SRL folds into a bulged-G motif and tetraloop, with restrictocin achieving a specificity of ∼1000-fold by recognizing both motifs only after the initial binding step. Here, we identify contacts within the protein−RNA interface and determine the extent to which each one contributes to enzyme specificity by examining the effect of protein mutations on the cleavage of the SRL substrate compared to a variety of other RNA substrates. As with other biomolecular interfaces, only a subset of contacts contributes to specificity. One contact of this subset is critical, with the H49A mutation resulting in quantitative loss of specificity. Maximum catalytic activity occurs when both motifs of the SRL are present, with the major contribution involving the bulged-G motif recognized by three lysine residues located adjacent to the active site: K110, K111, and K113. Our findings support a kinetic proofreading mechanism in which the active site residues H49 and, to a lesser extent, Y47 make greater catalytic contributions to SRL cleavage than to suboptimal substrates. This systematic and quantitative analysis begins to elucidate the principles governing RNA recognition by a site-specific endonuclease and may thus serve as a mechanistic model for investigating other RNA modifying enzymes.

Fluorescence competition and optical melting measurements of RNA three-way multibranch loops provide a revised model for thermodynamic parameters

Three-way multibranch loops (junctions) are common in RNA secondary structures. Computer algorithms such as RNAstructure and MFOLD do not consider the identity of unpaired nucleotides in multibranch loops when predicting secondary structure. There is limited experimental data, however, to parametrize this aspect of these algorithms. In this study, UV optical melting and a fluorescence competition assay are used to measure stabilities of multibranch loops containing up to five unpaired adenosines or uridines or a loop E motif. These results provide a test of our understanding of the factors affecting multibranch loop stability and provide revised parameters for predicting stability. The results should help to improve predictions of RNA secondary structure.

NMR structure and dynamics of the Specifier Loop domain from the Bacillus subtilis tyrS T box leader RNA

Gram-positive bacteria utilize a tRNA-responsive transcription antitermination mechanism, designated the T box system, to regulate expression of many amino acid biosynthetic and aminoacyl-tRNA synthetase genes. The RNA transcripts of genes controlled by this mechanism contain 5′ untranslated regions, or leader RNAs, that specifically bind cognate tRNA molecules through pairing of nucleotides in the tRNA anticodon loop with nucleotides in the Specifier Loop domain of the leader RNA. We have determined the solution structure of the Specifier Loop domain of the tyrS leader RNA from Bacillus subtilis. Fifty percent of the nucleotides in the Specifier Loop domain adopt a loop E motif. The Specifier Sequence nucleotides, which pair with the tRNA anticodon, stack with their Watson–Crick edges rotated toward the minor groove and exhibit only modest flexibility. We also show that a Specifier Loop domain mutation that impairs the function of the B. subtilis glyQS T box RNA disrupts the tyrS loop E motif. Our results suggest a mechanism for tRNA–Specifier Loop binding in which the phosphate backbone kink created by the loop E motif causes the Specifier Sequence bases to rotate toward the minor groove, which increases accessibility for pairing with bases in the anticodon loop of tRNA.

Comparative analysis of depurination catalyzed by ricin A-chain on synthetic 32mer and 25mer oligoribonucleotides mimicking the sarcin/ricin domain of the rat 28S rRNA and E. coli 23S rRNA

Ricin A-chain can inactivate eukaryotic ribosomes, but exhibits no N-glycosidase activity on intact E. coli ribosomes. In the present research, in order to avoid using radiolabeled oligoribonucleotides, two kinds of synthetic 5'-FAM fluorescence-labeled oligoribonucleotide substrates were used to mimic the sarcin/ricin domains of rat 28S rRNA and E. coli 23S rRNA (32mer and 25mer, named as Rat FAM-SRD and E. coli FAM-SRD, respectively). Ricin A-chain was able to specifically release adenine from the first adenosine of the GAGA tetraloop and exhibited specific N-glycosidase activity under neutral and weak acidic conditions with both substrates. However, under more acidic conditions, ricin A-chain was able to release purines from other sites on eukaryotic substrates, but it retained specific depurination activity on prokaryotic substrates. At pH 5.0, the Michaelis constant (K(m)) for the reaction with Rat FAM-SRD (4.57+/-0.28microM) corresponded to that with E. coli FAM-SRD (4.64+/-0.26microM). However, the maximum velocity (V(max)) for ricin A-chain with Rat FAM-SRD was 0.5+/-0.024microM/min, which is higher than that with E. coli FAM-SRD (0.32+/-0.011microM/min).

A putative loop E motif and an H-H kissing loop interaction are conserved and functional features in a group C enterovirus RNA that inhibits ribonuclease L

A phylogenetically conserved RNA structure within the open reading frame of poliovirus and other group C enteroviruses functions as a competitive inhibitor of the antiviral endoribonuclease RNase L. Hence, we call this viral RNA the RNase L competitive inhibitor RNA (RNase L ciRNA). In this investigation we used phylogenetic information, RNA structure prediction software, site-directed mutagenesis, and RNase L activity assays to identify functionally important sequences and structures of the RNase L ciRNA. A putative loop E motif is phylogenetically conserved in the RNA structure and mutations of nucleotides within the putative loop E motif destroyed the ability of the RNA molecule to inhibit RNase L. A putative H-H kissing loop interaction is phylogenetically conserved in the RNA structure and covariant polymorphisms that maintain the Watson-Crick complementarity required for the kissing interaction provide evidence of its importance. Compensatory mutations that disrupted and then restored the putative kissing interaction confirm that it contributes to the ability of the viral RNA to inhibit RNase L. RNase L was activated late during the course of poliovirus replication in HeLa cells, as virus replication and assembly neared completion. We conclude that a putative loop E motif and an H-H kissing loop interaction are key features of the group C enterovirus RNA associated with the inhibition of RNase L.

Electrostatic interactions guide the active site face of a structure-specific ribonuclease to its RNA substrate

Restrictocin, a member of the α-sarcin family of site-specific endoribonucleases, uses electrostatic interactions to bind to the ribosome and to RNA oligonucleotides, including the minimal specific substrate, the sarcin/ricin loop (SRL) of 23S−28S rRNA. Restrictocin binds to the SRL by forming a ground-state E:S complex that is stabilized predominantly by Coulomb interactions and depends on neither the sequence nor structure of the RNA, suggesting a nonspecific complex. The 22 cationic residues of restrictocin are dispersed throughout this protein surface, complicating a priori identification of a Coulomb interacting surface. Structural studies have identified an enzyme−substrate interface, which is expected to overlap with the electrostatic E:S interface. Here, we identified restrictocin residues that contribute to binding in the E:S complex by determining the salt dependence [∂ log(k₂/K_1/2)/∂ log[KCl]] of cleavage of the minimal SRL substrate for eight point mutants within the protein designed to disrupt contacts in the crystallographically defined interface. Relative to the wild-type salt dependence of −4.1, a subset of the mutants clustering near the active site shows significant changes in salt dependence, with differences of magnitude being ≥0.4. This same subset was identified using calculated salt dependencies for each mutant derived from solutions to the nonlinear Poisson−Boltzmann equation. Our findings support a mechanism in which specific residues on the active site face of restrictocin (primarily K110, K111, and K113) contribute to formation of the E:S complex, thereby positioning the SRL substrate for site-specific cleavage. The same restrictocin residues are expected to facilitate targeting of the SRL on the surface of the ribosome.

Introduction to special issue on RNA

In this introduction to the special issue on RNA, we provide a brief overview of some of the novel and exciting biological discoveries concerning diverse roles played by RNA, and subsequently we give a rapid summary of some algorithmic aspects of RNA structure and alignment. Each of the contributions to this special issue is briefly described.

Accommodating the bacterial decoding release factor as an alien protein among the RNAs at the active site of the ribosome

The decoding release factor (RF) triggers termination of protein synthesis by functionally mimicking a tRNA to span the decoding centre and the peptidyl transferase centre (PTC) of the ribosome. Structurally, it must fit into a site crafted for a tRNA and surrounded by five other RNAs, namely the adjacent peptidyl tRNA carrying the completed polypeptide, the mRNA and the three rRNAs. This is achieved by extending a structural domain from the body of the protein that results in a critical conformational change allowing it to contact the PTC. A structural model of the bacterial termination complex with the accommodated RF shows that it makes close contact with the first, second and third bases of the stop codon in the mRNA with two separate loops of structure: the anticodon loop and the loop at the tip of helix alpha5. The anticodon loop also makes contact with the base following the stop codon that is known to strongly influence termination efficiency. It confirms the close contact of domain 3 of the protein with the key RNA structures of the PTC. The mRNA signal for termination includes sequences upstream as well as downstream of the stop codon, and this may reflect structural restrictions for specific combinations of tRNA and RF to be bound onto the ribosome together. An unbiased SELEX approach has been investigated as a tool to identify potential rRNA-binding contacts of the bacterial RF in its different binding conformations within the active centre of the ribosome.

A comparative analysis of the triloops in all high-resolution RNA structures reveals sequence structure relationships

Despite an increasing number of experimentally determined RNA structures, the gap between the number of structures and that of RNA families is still growing. To overcome this limitation, efficient and reliable RNA modeling methodologies must be developed. In order to reach this goal, here, we show how triloop sequence-structure relationships have been inferred through a systematic analysis of all triloops found in available high-resolution structures. The structural annotation of all triloops allowed us to define discrete states of the triloop's conformational space, and therefore an explicit sequence-to-structure relation. The sequence-structure relationships inferred from this explicit relation are presented in a convenient modeling table that provides a limited set of possible three-dimensional structures given any triloop sequence. The table is indexed by the two nucleotides that form the triloop's flanking base pair, since they are shown to provide the most information about the triloop three-dimensional structures. We also report the observations in the X-ray crystallographic structures of important conformational variations, which we believe might be the result of RNA dynamic.

Modeling RNA tertiary structure motifs by graph-grammars

A new approach, graph-grammars, to encode RNA tertiary structure patterns is introduced and exemplified with the classical sarcin–ricin motif. The sarcin–ricin motif is found in the stem of the crucial ribosomal loop E (also referred to as the sarcin–ricin loop), which is sensitive to the α-sarcin and ricin toxins. Here, we generate a graph-grammar for the sarcin-ricin motif and apply it to derive putative sequences that would fold in this motif. The biological relevance of the derived sequences is confirmed by a comparison with those found in known sarcin–ricin sites in an alignment of over 800 bacterial 23S ribosomal RNAs. The comparison raised alternative alignments in few sarcin–ricin sites, which were assessed using tertiary structure predictions and 3D modeling. The sarcin–ricin motif graph-grammar was built with indivisible nucleotide interaction cycles that were recently observed in structured RNAs. A comparison of the sequences and 3D structures of each cycle that constitute the sarcin–ricin motif gave us additional insights about RNA sequence–structure relationships. In particular, this analysis revealed the sequence space of an RNA motif depends on a structural context that goes beyond the single base pairing and base-stacking interactions.

Fungal ribotoxins: molecular dissection of a family of natural killers

RNase T1 is the best known representative of a large family of ribonucleolytic proteins secreted by fungi, mostly Aspergillus and Penicillium species. Ribotoxins stand out among them by their cytotoxic character. They exert their toxic action by first entering the cells and then cleaving a single phosphodiester bond located within a universally conserved sequence of the large rRNA gene, known as the sarcin-ricin loop. This cleavage leads to inhibition of protein biosynthesis, followed by cellular death by apoptosis. Although no protein receptor has been found for ribotoxins, they preferentially kill cells showing altered membrane permeability, such as those that are infected with virus or transformed. Many steps of the cytotoxic process have been elucidated at the molecular level by means of a variety of methodological approaches and the construction and purification of different mutant versions of these ribotoxins. Ribotoxins have been used for the construction of immunotoxins, because of their cytotoxicity. Besides this activity, Aspf1, a ribotoxin produced by Aspergillus fumigatus, has been shown to be one of the major allergens involved in allergic aspergillosis-related pathologies. Protein engineering and peptide synthesis have been used in order to understand the basis of these pathogenic mechanisms as well as to produce hypoallergenic proteins with potential diagnostic and immunotherapeutic applications.

Structural energetics and base-pair opening dynamics in sarcin-ricin domain RNA

The sarcin-ricin domain is a universal element of the RNA from the large ribosomal subunit. The domain is part of the binding site for elongation factors and is specifically cleaved by the toxins alpha-sarcin and ricin. In this work, we have mapped the energetics and dynamics of individual structural motifs in a 29-mer RNA oligomer containing the sarcin-ricin domain. The stability of individual base pairs in the structure was characterized from measurements of the exchange rates of imino protons using nuclear magnetic resonance spectroscopy at 10 degrees C. The measurements also provided the rates of opening and closing for selected base pairs. The results reveal that the structural stabilization free energies in the sarcin-ricin domain are broadly distributed between 2.9 and 10.6 kcal/mol. One of the least stable sites in the structure is the noncanonical G-A base pair located next to the phosphodiester bond that is cleaved by alpha-sarcin. The low stability of this base pair supports the proposal that cleavage by alpha-sarcin occurs by a base flipping mechanism. The opening dynamics of other base pairs is affected by elements of the structure such as the bulged-G motif and its cross-strand stacking. Participation in these motifs increases the lifetimes of the bases in an open, solvent-accessible conformation.

Tertiary structural and functional analyses of a viroid RNA motif by isostericity matrix and mutagenesis reveal its essential role in replication

RNA-templated RNA replication is essential for viral or viroid infection, as well as for regulation of cellular gene expression. Specific RNA motifs likely regulate various aspects of this replication. Viroids of the Pospiviroidae family, as represented by the Potato spindle tuber viroid (PSTVd), replicate in the nucleus by utilizing DNA-dependent RNA polymerase II. We investigated the role of the loop E (sarcin/ricin) motif of the PSTVd genomic RNA in replication. A tertiary-structural model of this motif, inferred by comparative sequence analysis and comparison with nuclear magnetic resonance and X-ray crystal structures of loop E motifs in other RNAs, is presented in which core non-Watson-Crick base pairs are precisely specified. Isostericity matrix analysis of these base pairs showed that the model accounts for the reported natural sequence variations and viable experimental mutations in loop E motifs of PSTVd and other viroids. Furthermore, isostericity matrix analysis allowed us to design disruptive, as well as compensatory, mutations of PSTVd loop E. Functional analyses of such mutants by in vitro and in vivo experiments demonstrated that loop E structural integrity is crucial for replication, specifically during transcription. Our results suggest that the PSTVd loop E motif exists and functions in vivo and provide loss-of-function genetic evidence for the essential role of a viroid RNA three-dimensional motif in rolling-circle replication. The use of isostericity matrix analysis of non-Watson-Crick base pairing to rationalize mutagenesis of tertiary motifs and systematic in vitro and in vivo functional assays of mutants offers a novel, comprehensive approach to elucidate the tertiary-structure-function relationships for RNA motifs of general biological significance.

Molecular dynamics simulations of sarcin-ricin rRNA motif

Explicit solvent molecular dynamics (MD) simulations were carried out for sarcin-ricin domain (SRD) motifs from 23S (Escherichia coli) and 28S (rat) rRNAs. The SRD motif consists of GAGA tetraloop, G-bulged cross-strand A-stack, flexible region and duplex part. Detailed analysis of the overall dynamics, base pairing, hydration, cation binding and other SRD features is presented. The SRD is surprisingly static in multiple 25 ns long simulations and lacks any non-local motions, with root mean square deviation (r.m.s.d.) values between averaged MD and high-resolution X-ray structures of 1-1.4 A. Modest dynamics is observed in the tetraloop, namely, rotation of adenine in its apex and subtle reversible shift of the tetraloop with respect to the adjacent base pair. The deformed flexible region in low-resolution rat X-ray structure is repaired by simulations. The simulations reveal few backbone flips, which do not affect positions of bases and do not indicate a force field imbalance. Non-Watson-Crick base pairs are rigid and mediated by long-residency water molecules while there are several modest cation-binding sites around SRD. In summary, SRD is an unusually stiff rRNA building block. Its intrinsic structural and dynamical signatures seen in simulations are strikingly distinct from other rRNA motifs such as Loop E and Kink-turns.

Photochemical generation of ribose abasic sites in RNA oligonucleotides

[reaction: see text] A general method for the photochemical generation of ribose abasic sites within RNA oligonucleotides is reported. Photochemically caged nucleoside phosphoramidite analogues were prepared and incorporated into RNA oligonucleotides by automated RNA synthesis. Irradiation of the modified RNA at 350 nm efficiently produced ribose abasic sites at specific sites within RNA sequences. The current approach offers a chemical route to RNA abasic lesions for RNA biochemical studies.

Induced fit and "lock and key" recognition of 5S RNA by zinc fingers of transcription factor IIIA

Transcription factor IIIA (TFIIIA) is a Cys2His2 zinc finger protein that regulates expression of the 5 S ribosomal RNA gene by binding specifically to the internal control element. TFIIIA also functions in transport and storage of 5 S RNA by binding directly to the RNA transcript. To obtain insights into the mechanism by which TFIIIA recognizes 5 S RNA, we determined the solution structure of the middle three zinc fingers bound to the central core of 5 S RNA. Finger 4 utilizes "lock and key" recognition to bind in the widened major groove of the pre-structured RNA loop E motif. This interaction is mediated by direct hydrogen bonding interactions with bases. In contrast, recognition of loop A, a flexible junction of three helices, occurs by an induced fit mechanism that involves reorganization of the conserved CAUA motif and structuring of the finger 5-finger 6 interface to form a complementary RNA binding surface.

Modeling the highly specific ribotoxin recognition of ribosomes

The three-dimensional structures of the alpha-sarcin ribotoxin and its delta(7-22) deletion mutant, both complexed with a 20-mer oligonucleotide mimicking the sarcin/ricin loop (SRL) of the ribosome, have been docked into the structure of the Halobacterium marismortui ribosome by fitting the nucleotide atomic coordinates into those of the ribosomal SRL. This study has revealed that two regions of the ribotoxin, residues 11-16 and 84-85, contact the ribosomal proteins L14 (residues 99-105) and L6 (residues 88-92), respectively. The first of these two ribotoxin regions appears to be crucial for its specific ribosome recognition.

Solution probing of metal ion binding by helix 27 from Escherichia coli 16S rRNA

Helix (H)27 from Escherichia coli 16S ribosomal (r)RNA is centrally located within the small (30S) ribosomal subunit, immediately adjacent to the decoding center. Bacterial 30S subunit crystal structures depicting Mg(2+) binding sites resolve two magnesium ions within the vicinity of H27: one in the major groove of the G886-U911 wobble pair, and one within the GCAA tetraloop. Binding of such metal cations is generally thought to be crucial for RNA folding and function. To ask how metal ion-RNA interactions in crystals compare with those in solution, we have characterized, using solution NMR spectroscopy, Tb(3+) footprinting and time-resolved fluorescence resonance energy transfer (tr-FRET), location, and modes of metal ion binding in an isolated H27. NMR and Tb(3+) footprinting data indicate that solution secondary structure and Mg(2+) binding are generally consistent with the ribosomal crystal structures. However, our analyses also suggest that H27 is dynamic in solution and that metal ions localize within the narrow major groove formed by the juxtaposition of the loop E motif with the tandem G894-U905 and G895-U904 wobble pairs. In addition, tr-FRET studies provide evidence that Mg(2+) uptake by the H27 construct results in a global lengthening of the helix. We propose that only a subset of H27-metal ion interactions has been captured in the crystal structures of the 30S ribosomal subunit, and that small-scale structural dynamics afforded by solution conditions may contribute to these differences. Our studies thus highlight an example for differences between RNA-metal ion interactions observed in solution and in crystals.