RNA folds: insights from recent crystal structures

Accurately Modeling RNA Stem-Loops in an Implicit Solvent Environment

Ribonucleic acid (RNA) molecules can adopt a variety of secondary and tertiary structures in solution, with stem-loops being one of the more common motifs. Here, we present a systematic analysis of 15 RNA stem-loop sequences simulated with molecular dynamics simulations in an implicit solvent environment. Analysis of RNA cluster ensembles showed that the stem-loop structures can generally adopt the A-form RNA in the stem region. Loop structures are more sensitive, and experimental structures could only be reproduced with modification of CH···O interactions in the force field, combined with an implicit solvent nonpolar correction to better model base stacking interactions. Accurately modeling RNA with current atomistic physics-based models remains challenging, but the RNA systems studied herein may provide a useful benchmark set for testing other RNA modeling methods in the future.

Intricate 3D architecture of a DNA mimic of GFP

Numerous studies have shown how RNA molecules can adopt elaborate three-dimensional (3D) architectures1-3. By contrast, whether DNA can self-assemble into complex 3D folds capable of sophisticated biochemistry, independent of protein or RNA partners, has remained mysterious. Lettuce is an in vitro-evolved DNA molecule that binds and activates4 conditional fluorophores derived from GFP. To extend previous structural studies5,6 of fluorogenic RNAs, GFP and other fluorescent proteins7 to DNA, we characterize Lettuce-fluorophore complexes by X-ray crystallography and cryogenic electron microscopy. The results reveal that the 53-nucleotide DNA adopts a four-way junction (4WJ) fold. Instead of the canonical L-shaped or H-shaped structures commonly seen8 in 4WJ RNAs, the four stems of Lettuce form two coaxial stacks that pack co-linearly to form a central G-quadruplex in which the fluorophore binds. This fold is stabilized by stacking, extensive nucleobase hydrogen bonding-including through unusual diagonally stacked bases that bridge successive tiers of the main coaxial stacks of the DNA-and coordination of monovalent and divalent cations. Overall, the structure is more compact than many RNAs of comparable size. Lettuce demonstrates how DNA can form elaborate 3D structures without using RNA-like tertiary interactions and suggests that new principles of nucleic acid organization will be forthcoming from the analysis of complex DNAs.

Gene editing with CRISPR-Cas12a guides possessing ribose-modified pseudoknot handles

CRISPR-Cas12a is a leading technology for development of model organisms, therapeutics, and diagnostics. These applications could benefit from chemical modifications that stabilize or tune enzyme properties. Here we chemically modify ribonucleotides of the AsCas12a CRISPR RNA 5' handle, a pseudoknot structure that mediates binding to Cas12a. Gene editing in human cells required retention of several native RNA residues corresponding to predicted 2'-hydroxyl contacts. Replacing these RNA residues with a variety of ribose-modified nucleotides revealed 2'-hydroxyl sensitivity. Modified 5' pseudoknots with as little as six out of nineteen RNA residues, with phosphorothioate linkages at remaining RNA positions, yielded heavily modified pseudoknots with robust cell-based editing. High trans activity was usually preserved with cis activity. We show that the 5' pseudoknot can tolerate near complete modification when design is guided by structural and chemical compatibility. Rules for modification of the 5' pseudoknot should accelerate therapeutic development and be valuable for CRISPR-Cas12a diagnostics.

Photoswitches for controllable RNA binding: a future approach in the RNA-targeting therapy

RNA is an emerging drug target that opens new perspectives in the treatment of viral and bacterial infections, cancer and a range of so far incurable genetic diseases. Among the various strategies towards the design and development of selective and efficient ligands for targeting and detection of therapeutically relevant RNA, photoswitchable RNA binders represent a very promising approach due to the possibility to control the ligand-RNA and protein-RNA interactions by light with high spatiotemporal resolution. However, the field of photoswitchable RNA binders still remains underexplored due to challenging design of lead structures that should combine high RNA binding selectivity with efficient photochemical performance. The aim of this highlight article is to describe the development of photoswitchable noncovalent RNA binders and to outline the current situation and perspectives of this emerging interdisciplinary field.

Structural analyses of an RNA stability element interacting with poly(A)

Cis-acting RNA elements are crucial for the regulation of polyadenylated RNA stability. The element for nuclear expression (ENE) contains a U-rich internal loop flanked by short helices. An ENE stabilizes RNA by sequestering the poly(A) tail via formation of a triplex structure that inhibits a rapid deadenylation-dependent decay pathway. Structure-based bioinformatic studies identified numerous ENE-like elements in evolutionarily diverse genomes, including a subclass containing two ENE motifs separated by a short double-helical region (double ENEs [dENEs]). Here, the structure of a dENE derived from a rice transposable element (TWIFB1) before and after poly(A) binding (∼24 kDa and ∼33 kDa, respectively) is investigated. We combine biochemical structure probing, small angle X-ray scattering (SAXS), and cryo-electron microscopy (cryo-EM) to investigate the dENE structure and its local and global structural changes upon poly(A) binding. Our data reveal 1) the directionality of poly(A) binding to the dENE, and 2) that the dENE-poly(A) interaction involves a motif that protects the 3'-most seven adenylates of the poly(A). Furthermore, we demonstrate that the dENE does not undergo a dramatic global conformational change upon poly(A) binding. These findings are consistent with the recently solved crystal structure of a dENE+poly(A) complex [S.-F. Torabi et al., Science 371, eabe6523 (2021)]. Identification of additional modes of poly(A)-RNA interaction opens new venues for better understanding of poly(A) tail biology.

Large Phenotypic Enhancement of Structured Random RNA Pools

Laboratory evolution of functional RNAs has applications in many areas of chemical and synthetic biology. In vitro selections critically depend on the presence of functional molecules, such as aptamers and ribozymes, in the starting sequence pools. For selection of novel functions the pools are typically transcribed from random-sequence DNA templates, yielding a highly diverse set of RNAs that contain a multitude of folds and biochemical activities. The phenotypic potential, the frequency of functional RNAs, is very low, requiring large complexity of starting pools, surpassing 10¹⁵ different sequences, to identify highly active isolates. Furthermore, the majority of random sequences is not structured and has a high propensity for aggregation; the in vitro selection process thus involves not just enrichment of functional RNAs, but also their purification from aggregation-prone "free-riders". We reasoned that purification of the nonaggregating, monomeric subpopulation of a random-sequence RNA pool will yield pools of folded, functional RNAs. We performed six rounds of selection for monomeric sequences and show that the enriched population is compactly folded. In vitro selections originating from various mixtures of the compact pool and a fully random pool showed that sequences from the compact pool always dominate the population once a biochemical activity is detectable. A head-to-head competition of the two pools starting from a low (5 × 10¹²) sequence diversity revealed that the phenotypic potential of the compact pool is about 1000-times higher than the fully random pool. A selection for folded and monomeric RNA pools thus greatly increases the frequency of functional RNAs from that seen in random-sequence pools, providing a facile experimental approach to isolation of highly active functional RNAs from low-diversity populations.

Effects of Noncanonical Base Pairing on RNA Folding: Structural Context and Spatial Arrangements of G·A Pairs

Noncanonical base pairs play important roles in assembling the three-dimensional structures critical to the diverse functions of RNA. These associations contribute to the looped segments that intersperse the canonical double-helical elements within folded, globular RNA molecules. They stitch together various structural elements, serve as recognition elements for other molecules, and act as sites of intrinsic stiffness or deformability. This work takes advantage of new software (DSSR) designed to streamline the analysis and annotation of RNA three-dimensional structures. The multiscale structural information gathered for individual molecules, combined with the growing number of unique, well-resolved RNA structures, makes it possible to examine the collective features deeply and to uncover previously unrecognized patterns of chain organization. Here we focus on a subset of noncanonical base pairs involving guanine and adenine and the links between their modes of association, secondary structural context, and contributions to tertiary folding. The rigorous descriptions of base-pair geometry that we employ facilitate characterization of recurrent geometric motifs and the structural settings in which these arrangements occur. Moreover, the numerical parameters hint at the natural motions of the interacting bases and the pathways likely to connect different spatial forms. We draw attention to higher-order multiplexes involving two or more G·A pairs and the roles these associations appear to play in bridging different secondary structural units. The collective data reveal pairing propensities in base organization, secondary structural context, and deformability and serve as a starting point for further multiscale investigations and/or simulations of RNA folding.

Surprising Sequence Effects on GU Closure of Symmetric 2 × 2 Nucleotide RNA Internal Loops

GU base pairs are important RNA structural motifs and often close loops. Accurate prediction of RNA structures relies upon understanding the interactions determining structure. The thermodynamics of some 2 × 2 nucleotide internal loops closed by GU pairs are not well understood. Here, several self-complementary oligonucleotide sequences expected to form duplexes with 2 × 2 nucleotide internal loops closed by GU pairs were investigated. Surprisingly, nuclear magnetic resonance revealed that many of the sequences exist in equilibrium between hairpin and duplex conformations. This equilibrium is not observed with loops closed by Watson-Crick pairs. To measure the thermodynamics of some 2 × 2 nucleotide internal loops closed by GU pairs, non-self-complementary sequences that preclude formation of hairpins were designed. The measured thermodynamics indicate that some internal loops closed by GU pairs are unusually unstable. This instability accounts for the observed equilibria between duplex and hairpin conformations. Moreover, it suggests that future three-dimensional structures of loops closed by GU pairs may reveal interactions that unexpectedly destabilize folding.

In vitro evolution of coenzyme-independent variants from the glmS ribozyme structural scaffold

Uniquely among known natural ribozymes that cleave RNA sequence-specifically, the glmS ribozyme-riboswitch employs a small molecule, glucosamine-6-phosphate (GlcN6P) as a catalytic cofactor. In vitro selection was employed to search for coenzyme-independent variants of this ribozyme. In addition to shedding light on the catalytic mechanism of the ribozyme, such variants could resemble the evolutionary ancestors of the modern, GlcN6P-regulated ribozyme-riboswitch. A mutant pool was constructed such that the secondary structure elements, which define the triply-pseudoknotted global fold of the ribozyme, was preserved. A stringent selection scheme that relies on thiol-mercury affinity chromatography for separating active and inactive sequences ultimately yielded a triple mutant with a cleavage rate exceeding 3min(-1) that only requires divalent cations for activity. Mutational analysis demonstrated that a point reversion of the variant toward the wild-type sequence was sufficient to partially restore GlcN6P-dependence, suggesting that coenzyme dependence can be readily be acquired by RNAs that adopt the glmS ribozyme fold. The methods employed to perform this selection experiment are described in detail in this review.

Progress and Current Challenges in Modeling Large RNAs

Recent breakthroughs in next-generation sequencing technologies have led to the discovery of several classes of non-coding RNAs (ncRNAs). It is now apparent that RNA molecules are not only just carriers of genetic information but also key players in many cellular processes. While there has been a rapid increase in the number of ncRNA sequences deposited in various databases over the past decade, the biological functions of these ncRNAs are largely not well understood. Similar to proteins, RNA molecules carry out a function by forming specific three-dimensional structures. Understanding the function of a particular RNA therefore requires a detailed knowledge of its structure. However, determining experimental structures of RNA is extremely challenging. In fact, RNA-only structures represent just 1% of the total structures deposited in the PDB. Thus, computational methods that predict three-dimensional RNA structures are in high demand. Computational models can provide valuable insights into structure-function relationships in ncRNAs and can aid in the development of functional hypotheses and experimental designs. In recent years, a set of diverse RNA structure prediction tools have become available, which differ in computational time, input data and accuracy. This review discusses the recent progress and challenges in RNA structure prediction methods.

RNA 3D Modules in Genome-Wide Predictions of RNA 2D Structure

Recent experimental and computational progress has revealed a large potential for RNA structure in the genome. This has been driven by computational strategies that exploit multiple genomes of related organisms to identify common sequences and secondary structures. However, these computational approaches have two main challenges: they are computationally expensive and they have a relatively high false discovery rate (FDR). Simultaneously, RNA 3D structure analysis has revealed modules composed of non-canonical base pairs which occur in non-homologous positions, apparently by independent evolution. These modules can, for example, occur inside structural elements which in RNA 2D predictions appear as internal loops. Hence one question is if the use of such RNA 3D information can improve the prediction accuracy of RNA secondary structure at a genome-wide level. Here, we use RNAz in combination with 3D module prediction tools and apply them on a 13-way vertebrate sequence-based alignment. We find that RNA 3D modules predicted by metaRNAmodules and JAR3D are significantly enriched in the screened windows compared to their shuffled counterparts. The initially estimated FDR of 47.0% is lowered to below 25% when certain 3D module predictions are present in the window of the 2D prediction. We discuss the implications and prospects for further development of computational strategies for detection of RNA 2D structure in genomic sequence.

Recognition of the bacterial alarmone ZMP through long-distance association of two RNA subdomains

The bacterial alarmone 5-aminoimidazole-4-carboxamide riboside 5'-triphosphate (ZTP), derived from the monophosphorylated purine precursor ZMP, accumulates during folate starvation. ZTP regulates genes involved in purine and folate metabolism through a cognate riboswitch. The linker connecting this riboswitch’s two sub-domains varies in length by over 100 nucleotides. We report the co-crystal structure of the Fusobacterium ulcerans riboswitch bound to ZMP, which spans the two sub-domains whose interface also comprises a pseudoknot and ribose zipper. The riboswitch recognizes the carboxamide oxygen of ZMP through an unprecedented inner-sphere coordination with a Mg²⁺ ion. We demonstrate that the affinity of the riboswitch for ZMP is modulated by the linker length. Notably, ZMP can bind to the two sub-domains together even when synthesized as separate RNAs. The ZTP riboswitch demonstrates how specific small-molecule binding can drive association of distant non-coding RNA domains to regulate gene expression.

Structure and mechanism of the T-box riboswitches

In most Gram-positive bacteria, including many clinically devastating pathogens from genera such as Bacillus, Clostridium, Listeria, and Staphylococcus, T-box riboswitches sense and regulate intracellular availability of amino acids through a multipartite messenger RNA (mRNA)-transfer RNA (tRNA) interaction. The T-box mRNA leaders respond to nutrient starvation by specifically binding cognate tRNAs and sensing whether the bound tRNA is aminoacylated, as a proxy for amino acid availability. Based on this readout, T-boxes direct a transcriptional or translational switch to control the expression of downstream genes involved in various aspects of amino acid metabolism: biosynthesis, transport, aminoacylation, transamidation, and so forth. Two decades after its discovery, the structural and mechanistic underpinnings of the T-box riboswitch were recently elucidated, producing a wealth of insights into how two structured RNAs can recognize each other with robust affinity and exquisite selectivity. The T-box paradigm exemplifies how natural noncoding RNAs can interact not just through sequence complementarity but can add molecular specificity by precisely juxtaposing RNA structural motifs, exploiting inherently flexible elements and the biophysical properties of post-transcriptional modifications, ultimately achieving a high degree of shape complementarity through mutually induced fit. The T-box also provides a proof-of-principle that compact RNA domains can recognize minute chemical changes (such as tRNA aminoacylation) on another RNA. The unveiling of the structure and mechanism of the T-box system thus expands our appreciation of the range of capabilities and modes of action of structured noncoding RNAs, and hints at the existence of networks of noncoding RNAs that communicate through both, structural and sequence specificity.

Ribonucleic acid (RNA) biosynthesis in human cancer

In many respects, the most remarkable chemical substances within the genome of eukaryotic cells are remarkable proteins which are the critical structural and functional units of living cells. The specifications for everything that goes in the cell are natural digital-to-digital decoding process in an archive sequence by deoxyribonucleic acid (DNA) and an articulate construction by ribonucleic acid (RNA). The products of DNA transcription are long polymers of ribonucleotides rather than deoxyribonucleotides and are termed ribonucleic acids. Certain deoxyribonucleotide sequences, or genes, give rise to transfer RNA (tRNA) and other ribosomal RNA (rRNA) when transcribed. The ribonucleotide sequences fold extensively and rRNA is associated with specific proteins to yield the essential cell components, ribosomes. Transcription of other special sequences yields messenger RNAs (mRNAs) that contain ribonucleotide sequences that will be ultimately translated into new types of amino acid sequences of functional cellular protein molecules. This switch to a different variety of cellular molecular sequences is complex, but each sequence of the three ribonucleotides specifies the insertion of one particular amino acid into the polypeptide chain under production. Whilst mRNA is considered the vehicle by which genetic information is transmitted from the genome and allocated in the appropriate cytoplasmic sites for translation into protein via cap-dependent mechanism, the actual translation depends also on the presence of other so-called household and luxury protein molecules. Recent evidence suggests RNA species are required at initiation, because treatment of cells with antibiotics or drugs that inhibit RNA synthesis cause a decrease in protein synthesis. The rRNA is necessary as a structural constituent of the ribosomes upon which translation takes place, whereas tRNA is necessary as an adaptor in amino acid activation and elongation protein chains to ribosomes. In this article, we review malignant tumor, with stem like properties, and recent technical advances into the phenomenon of micro-particles and micro-vesicles containing cell-free nucleic acids that circulate plasma. New areas of research have been opened into screening tumor telomerase progression, prognosis of aptamers targeting cell surface, monitoring the efficacy of anticancer therapies, oncogenic transformation of host cell, and RNA polymerases role in the cell cycle progression and differentiation.

Antibiotic drugs targeting bacterial RNAs

RNAs have diverse structures that include bulges and internal loops able to form tertiary contacts or serve as ligand binding sites. The recent increase in structural and functional information related to RNAs has put them in the limelight as a drug target for small molecule therapy. In addition, the recognition of the marked difference between prokaryotic and eukaryotic rRNA has led to the development of antibiotics that specifically target bacterial rRNA, reduce protein translation and thereby inhibit bacterial growth. To facilitate the development of new antibiotics targeting RNA, we here review the literature concerning such antibiotics, mRNA, riboswitch and tRNA and the key methodologies used for their screening.

Automated identification of RNA 3D modules with discriminative power in RNA structural alignments

Recent progress in predicting RNA structure is moving towards filling the ‘gap’ in 2D RNA structure prediction where, for example, predicted internal loops often form non-canonical base pairs. This is increasingly recognized with the steady increase of known RNA 3D modules. There is a general interest in matching structural modules known from one molecule to other molecules for which the 3D structure is not known yet. We have created a pipeline, metaRNAmodules, which completely automates extracting putative modules from the FR3D database and mapping of such modules to Rfam alignments to obtain comparative evidence. Subsequently, the modules, initially represented by a graph, are turned into models for the RMDetect program, which allows to test their discriminative power using real and randomized Rfam alignments. An initial extraction of 22 495 3D modules in all PDB files results in 977 internal loop and 17 hairpin modules with clear discriminatory power. Many of these modules describe only minor variants of each other. Indeed, mapping of the modules onto Rfam families results in 35 unique locations in 11 different families. The metaRNAmodules pipeline source for the internal loop modules is available at http://rth.dk/resources/mrm.

Oligonucleotide optical switches for intracellular sensing

Fluorescence imaging coupled with nanotechnology is making possible the development of powerful tools in the biological field for applications such as cellular imaging and intracellular messenger RNA monitoring and detection. The delivery of fluorescent probes into cells and tissues is currently receiving growing interest because such molecules, often coupled to nanodimensional materials, can conveniently allow the preparation of small tools to spy on cellular mechanisms with high specificity and sensitivity. The purpose of this review is to provide an exhaustive overview of current research in oligonucleotide optical switches for intracellular sensing with a focus on the engineering methods adopted for these oligonucleotides and the more recent and fascinating techniques for their internalization into living cells. Oligonucleotide optical switches can be defined as specifically designed short nucleic acid molecules capable of turning on or modifying their light emission on molecular interaction with well-defined molecular targets. Molecular beacons, aptamer beacons, hybrid molecular probes, and simpler linear oligonucleotide switches are the most promising optical nanosensors proposed in recent years. The intracellular targets which have been considered for sensing are a plethora of messenger-RNA-expressing cellular proteins and enzymes, or, directly, proteins or small molecules in the case of sensing through aptamer-based switches. Engineering methods, including modification of the oligonucleotide itself with locked nucleic acids, peptide nucleic acids, or L-DNA nucleotides, have been proposed to enhance the stability of nucleases and to prevent false-negative and high background optical signals. Conventional delivery techniques are treated here together with more innovative methods based on the coupling of the switches with nano-objects.

RNAmap2D - calculation, visualization and analysis of contact and distance maps for RNA and protein-RNA complex structures

BACKGROUND

The structures of biological macromolecules provide a framework for studying their biological functions. Three-dimensional structures of proteins, nucleic acids, or their complexes, are difficult to visualize in detail on flat surfaces, and algorithms for their spatial superposition and comparison are computationally costly. Molecular structures, however, can be represented as 2D maps of interactions between the individual residues, which are easier to visualize and compare, and which can be reconverted to 3D structures with reasonable precision. There are many visualization tools for maps of protein structures, but few for nucleic acids.

RESULTS

We developed RNAmap2D, a platform-independent software tool for calculation, visualization and analysis of contact and distance maps for nucleic acid molecules and their complexes with proteins or ligands. The program addresses the problem of paucity of bioinformatics tools dedicated to analyzing RNA 2D maps, given the growing number of experimentally solved RNA structures in the Protein Data Bank (PDB) repository, as well as the growing number of tools for RNA 2D and 3D structure prediction. RNAmap2D allows for calculation and analysis of contacts and distances between various classes of atoms in nucleic acid, protein, and small ligand molecules. It also discriminates between different types of base pairing and stacking.

CONCLUSIONS

RNAmap2D is an easy to use method to visualize, analyze and compare structures of nucleic acid molecules and their complexes with other molecules, such as proteins or ligands and metal ions. Its special features make it a very useful tool for analysis of tertiary structures of RNAs. RNAmap2D for Windows/Linux/MacOSX is freely available for academic users at http://iimcb.genesilico.pl/rnamap2d.html.

Hydrostatic and osmotic pressure study of the RNA hydration

The tertiary structure of nucleic acids results from an equilibrium between electrostatic interactions of phosphates, stacking interactions of bases, hydrogen bonds between polar atoms and water molecules. Water interactions with ribonucleic acid play a key role in its structure formation, stabilization and dynamics. We used high hydrostatic pressure and osmotic pressure to analyze changes in RNA hydration. We analyzed the lead catalyzed hydrolysis of tRNA^Phe from S. cerevisiae as well as hydrolytic activity of leadzyme. Pb(II) induced hydrolysis of the single phosphodiester bond in tRNA^Phe is accompanied by release of 98 water molecules, while other molecule, leadzyme releases 86.

HD-RNAS: An Automated Hierarchical Database of RNA Structures

One of the important goals of most biological investigations is to classify and organize the experimental findings so that they are readily useful for deriving generalized rules. Although there is a huge amount of information on RNA structures in PDB, there are redundant files, ambiguous synthetic sequences etc. Moreover, a systematic hierarchical organization, reflecting RNA classification, is missing in PDB. In this investigation, we have classified all the available RNA structures from PDB through a programmatic approach. Hence, it would be now a simple assignment to regularly update the classification as and when new structures are released. The classification can further determine (i) a non-redundant set of RNA structures and (ii) if available, a set of structures of identical sequence and function, which can highlight structural polymorphism, ligand-induced conformational alterations etc. Presently, we have classified the available structures (2095 PDB entries having RNA chain longer than nine nucleotides solved by X-ray crystallography or NMR spectroscopy) into nine functional classes. The structures of same function and same source are mostly seen to be similar with subtle differences depending on their functional complexation. The web-server is available online at http://www.saha.ac.in/biop/www/HD-RNAS.html and is updated regularly.

Electrochemical nanomaterial-based nucleic acid aptasensors

Recent progress in the development of electrochemical nanomaterial-aptamer-based biosensors is summarized. Aptamers are nucleic acid ligands that can be generated against amino acids, drugs, proteins, and other molecules. They are isolated from a large random library of synthetic nucleic acids by an iterative process of binding, separation, and amplification, called systematic evolution of ligands by exponential enrichment (SELEX). In this review, different methods of integrating aptamers with different nanomaterials and nanoparticles for electrochemical biosensing application are described.

Nucleic acid and peptide aptamers: fundamentals and bioanalytical aspects

In recent years new nucleic acid and protein-based combinatorial molecules have attracted the attention of researchers working in various areas of science, ranging from medicine to analytical chemistry. These molecules, called aptamers, have been proposed as alternatives to antibodies in many different applications. The aim of this Review is to illustrate the peculiarities of these combinatorial molecules which have initially been explored for their importance in molecular medicine, but have enormous potential in other biotechnological fields historically dominated by antibodies, such as bioassays. A description of these molecules is given, and the methods for their selection and production are also summarized. Moreover, critical aspects related to these molecules are discussed.

RNA and protein 3D structure modeling: similarities and differences

In analogy to proteins, the function of RNA depends on its structure and dynamics, which are encoded in the linear sequence. While there are numerous methods for computational prediction of protein 3D structure from sequence, there have been very few such methods for RNA. This review discusses template-based and template-free approaches for macromolecular structure prediction, with special emphasis on comparison between the already tried-and-tested methods for protein structure modeling and the very recently developed “protein-like” modeling methods for RNA. We highlight analogies between many successful methods for modeling of these two types of biological macromolecules and argue that RNA 3D structure can be modeled using “protein-like” methodology. We also highlight the areas where the differences between RNA and proteins require the development of RNA-specific solutions.

Figure

Structural characterization of naturally occurring RNA single mismatches

RNA is known to be involved in several cellular processes; however, it is only active when it is folded into its correct 3D conformation. The folding, bending and twisting of an RNA molecule is dependent upon the multitude of canonical and non-canonical secondary structure motifs. These motifs contribute to the structural complexity of RNA but also serve important integral biological functions, such as serving as recognition and binding sites for other biomolecules or small ligands. One of the most prevalent types of RNA secondary structure motifs are single mismatches, which occur when two canonical pairs are separated by a single non-canonical pair. To determine sequence–structure relationships and to identify structural patterns, we have systematically located, annotated and compared all available occurrences of the 30 most frequently occurring single mismatch-nearest neighbor sequence combinations found in experimentally determined 3D structures of RNA-containing molecules deposited into the Protein Data Bank. Hydrogen bonding, stacking and interaction of nucleotide edges for the mismatched and nearest neighbor base pairs are described and compared, allowing for the identification of several structural patterns. Such a database and comparison will allow researchers to gain insight into the structural features of unstudied sequences and to quickly look-up studied sequences.

SCAN--a high-throughput assay for detecting small molecule binding to RNA targets

A novel optical-based high-throughput screening technology has been developed for increasing the rate of discovering chemical leads against RNA targets. SCAN ( Screen for Compounds with Affinity for Nucleic Acids) is an affinity-based assay that identifies small molecules that bind and recognize structured RNA elements. This technology provides the opportunity to conduct high-throughput screening of a new class of targets-RNA. SCAN offers many attractive features including a simple homogeneous format, low screening costs, and the ability to use common laboratory equipment. A SCAN assay was developed for the HCV IRES Loop IIId RNA domain. A high-throughput screen of our entire compound library resulted in the identification of small molecule ligands that bind to Loop IIId. The Z' values were greater than 0.8, showing this to be a robust high-throughput screening assay. A correlation between SCAN EC50 and KD values is reported suggesting the ability to use the assay for compound optimization.

Methods to crystallize RNA

Preparation of suitably large and well-ordered single crystals is usually the rate-limiting step in the determination of the three-dimensional structure of RNAs and their complexes with proteins by X-ray crystallography. This unit discusses a variety of experimental considerations for obtaining crystals of RNAs and RNA-protein complexes. Topics include design of crystallizable constructs, screening, and optimization of crystallization conditions.

NMR (13)C-relaxation study of base and sugar dynamics in GCAA RNA hairpin tetraloop

Intramolecular dynamics of a 14-mer RNA hairpin including GCAA tetraloop was investigated by (13)C NMR relaxation. R(1) and R(1rho) relaxation rates were measured for all protonated base carbons as well as for C1' carbons of ribose sugars at several magnetic field strengths. The data has been interpreted in the framework of modelfree analysis [G. Lipari and A. Szabo. J Am Chem Soc 104, 4546-4559 (1982); G. Lipari and A. Szabo. J Am Chem Soc 104, 4559-4570 (1982)] characterizing the internal dynamics of the molecule by order parameters and correlation times for fast motions on picosecond to nanosecond time scale and by contributions of the chemical exchange. The fast dynamics reveals a rather rigid stem and a significantly more flexible loop. The cytosine and the last adenine bases in the loop as well as all the loop sugars exhibit a significant contribution of conformational equilibrium on microsecond to millisecond time scale. The high R(1rho) values detected on both base and sugar moieties of the loop indicate coordinated motions in this region. A semiquantitative analysis of the conformational equilibrium suggests the exchange rates on the order of 10(4) s(-1). The results are in general agreement with dynamics studies of GAAA loops by NMR relaxation and fluorescent spectroscopy and support the data on the GCAA loop dynamics obtained by MD simulations.

Riboswitches: small-molecule recognition by gene regulatory RNAs

Riboswitches demonstrate the ability of highly structured RNA molecules to recognize small-molecule metabolites with high specificity and subsequently harness the binding energy for the control of gene expression. Crystal structures have now been determined for the metabolite-binding domains of riboswitches that respond to purines, thiamine pyrophosphate and S-adenosylmethionine, as well as for the glmS ribozyme, a catalytic riboswitch that is activated by the metabolite glucosamine-6-phosphate. In addition to these riboswitch structures, a solution NMR structure has been reported for a ribosensor that regulates heat shock genes in response to changes in temperature. These studies reveal the structural basis of the remarkable selectivity of riboswitches and, in conjunction with biochemical and biophysical measurements, provide a framework for detailed mechanistic understanding of riboswitch-mediated modulation of gene expression.

Crystal structures of the thi-box riboswitch bound to thiamine pyrophosphate analogs reveal adaptive RNA-small molecule recognition

Riboswitches are noncoding mRNA elements that bind small-molecule metabolites with high affinity and specificity, and they regulate the expression of associated genes. The thi-box riboswitch can exhibit a 1000-fold higher affinity for thiamine pyrophosphate over closely related noncognate compounds such as thiamine monophosphate. To understand the chemical basis of thi-box pyrophosphate specificity, we have determined crystal structures of an E. coli thi-box bound to thiamine pyrophosphate, thiamine monophosphate, and the structural analogs benfotiamine and pyrithiamine. When bound to monophosphorylated compounds, the RNA elements that recognize the thiamine and phosphate moieties of the ligand move closer together. This allows the riboswitch to recognize the monophosphate in a manner similar to how it recognizes the beta-phosphate of thiamine pyrophosphate. In the pyrithiamine complex, the pyrophosphate binding site is largely unstructured. These results show how the riboswitch can bind to various metabolites, and why the thi-box preferentially binds thiamine pyrophosphate.

The RNA Ontology Consortium: an open invitation to the RNA community

The aim of the RNA Ontology Consortium (ROC) is to create an integrated conceptual framework-an RNA Ontology (RO)-with a common, dynamic, controlled, and structured vocabulary to describe and characterize RNA sequences, secondary structures, three-dimensional structures, and dynamics pertaining to RNA function. The RO should produce tools for clear communication about RNA structure and function for multiple uses, including the integration of RNA electronic resources into the Semantic Web. These tools should allow the accurate description in computer-interpretable form of the coupling between RNA architecture, function, and evolution. The purposes for creating the RO are, therefore, (1) to integrate sequence and structural databases; (2) to allow different computational tools to interoperate; (3) to create powerful software tools that bring advanced computational methods to the bench scientist; and (4) to facilitate precise searches for all relevant information pertaining to RNA. For example, one initial objective of the ROC is to define, identify, and classify RNA structural motifs described in the literature or appearing in databases and to agree on a computer-interpretable definition for each of these motifs. To achieve these aims, the ROC will foster communication and promote collaboration among RNA scientists by coordinating frequent face-to-face workshops to discuss, debate, and resolve difficult conceptual issues. These meeting opportunities will create new directions at various levels of RNA research. The ROC will work closely with the PDB/NDB structural databases and the Gene, Sequence, and Open Biomedical Ontology Consortia to integrate the RO with existing biological ontologies to extend existing content while maintaining interoperability.

MeRNA: a database of metal ion binding sites in RNA structures

Metal ions are essential for the folding of RNA into stable tertiary structures and for the catalytic activity of some RNA enzymes. To aid in the study of the roles of metal ions in RNA structural biology, we have created MeRNA (Metals in RNA), a comprehensive compilation of all metal binding sites identified in RNA 3D structures available from the PDB and Nucleic Acid Database. Currently, our database contains information relating to binding of 9764 metal ions corresponding to 23 distinct elements, in 256 RNA structures. The metal ion locations were confirmed and ligands characterized using original literature references. MeRNA includes eight manually identified metal-ion binding motifs, which are described in the literature. MeRNA is searchable by PDB identifier, metal ion, method of structure determination, resolution and R-values for X-ray structure and distance from metal to any RNA atom or to water. New structures with their respective binding motifs will be added to the database as they become available. The MeRNA database will further our understanding of the roles of metal ions in RNA folding and catalysis and have applications in structural and functional analysis, RNA design and engineering. The MeRNA database is accessible at .

Conformational transitions in RNA single uridine and adenosine bulge structures: a molecular dynamics free energy simulation study

Extra unmatched nucleotides (single base bulges) are common structural motifs in folded RNA molecules and can participate in RNA-ligand binding and RNA tertiary structure formation. Often these processes are associated with conformational transitions in the bulge region such as flipping out of the bulge base from an intrahelical stacked toward a looped out state. Knowledge of the flexibility of bulge structures and energetics of conformational transitions is an important prerequisite to better understand the function of this RNA motif. Molecular dynamics simulations were performed on single uridine and adenosine bulge nucleotides at the center of eight basepair RNA molecules and indicated larger flexibility of the bulge bases compared to basepaired regions. The umbrella sampling method was applied to study the bulge base looping out process and accompanying conformational and free energy changes. Looping out toward the major groove resulted in partial disruption of adjacent basepairs and was found to be less favorable compared to looping out toward the minor groove. For both uridine and adenosine bulges, a positive free energy change for full looping out was obtained which was approximately 1.5 kcal mol-1 higher in the case of the adenosine compared to the uridine bulge system. The simulations also indicated stable partially looped out states with the bulge bases located in the RNA minor groove and forming base triples with 5'-neighboring basepairs. In the case of the uridine bulge this state was more stable than the intrahelical stacked bulge structure. Induced looping out toward the minor groove involved crossing of an energy barrier of approximately 3.5 kcal mol-1 before reaching the base triple state. A continuum solvent analysis of intermediate bulge states indicated that electrostatic interactions stabilize looped out and base triple states, whereas van der Waals interactions and nonpolar contributions favor the stacked bulge conformation.

Crystal structure of a group I intron splicing intermediate

A recently reported crystal structure of an intact bacterial group I self-splicing intron in complex with both its exons provided the first molecular view into the mechanism of RNA splicing. This intron structure, which was trapped in the state prior to the exon ligation reaction, also reveals the architecture of a complex RNA fold. The majority of the intron is contained within three internally stacked, but sequence discontinuous, helical domains. Here the tertiary hydrogen bonding and stacking interactions between the domains, and the single-stranded joiner segments that bridge between them, are fully described. Features of the structure include: (1) A pseudoknot belt that circumscribes the molecule at its longitudinal midpoint; (2) two tetraloop-tetraloop receptor motifs at the peripheral edges of the structure; (3) an extensive minor groove triplex between the paired and joiner segments, P6-J6/6a and P3-J3/4, which provides the major interaction interface between the intron's two primary domains (P4-P6 and P3-P9.0); (4) a six-nucleotide J8/7 single stranded element that adopts a mu-shaped structure and twists through the active site, making critical contacts to all three helical domains; and (5) an extensive base stacking architecture that realizes 90% of all possible stacking interactions. The intron structure was validated by hydroxyl radical footprinting, where strong correlation was observed between experimental and predicted solvent accessibility. Models of the pre-first and pre-second steps of intron splicing are proposed with full-sized tRNA exons. They suggest that the tRNA undergoes substantial angular motion relative to the intron between the two steps of splicing.

Efficient Synthesis of 2‘,3‘-Dideoxy-2‘-amino-3‘-thiouridine

[reaction: see text] Metal ion rescue experiments provide a powerful approach to establish the presence and role of divalent metal ions in the biological function of RNA. The utility of this approach depends on the availability of suitable nucleoside analogues. To expand the range of this experimental strategy, we describe the first synthesis of 2',3'-dideoxy-2'-amino-3'-thiouridine (12) in 19.5% overall yield starting from 2,2'-anhydrouridine (1).

The N-glycosidase activity of the ribosome-inactivating protein ME1 targets single-stranded regions of nucleic acids independent of sequence or structural motifs

ME(1), a type I ribosome-inactivating protein (RIP), belongs to a family of enzymes long believed to possess rRNA N-glycosidase activity directed solely at the universally conserved residue A4324 in the sarcin/ricin loop of large eukaryotic and prokaryotic rRNAs. We have investigated the effect of modifying the structure of nonribosomal RNA substrates on their interaction with ME(1) and other RIPs. ME(1) was shown to depurinate a variety of partially denatured nucleic acids, randomly removing adenine residues from single-stranded regions and, to a lesser extent, guanine residues from wobble base-pairs in hairpin stems. A defined sequence motif was not required for recognition of non-paired adenosines and cleavage of the N-glycosidic bond. Substrate recognition and ME(1) activity appeared to depend on the physical availability of nucleotides, and denaturation of nucleic acid substrates increased their interaction with ME(1). Pretreatment of mRNA at 75 degrees C rather than 60 degrees C, for example, lowered the apparent K(D) from 87.1 to 73.9 nm, making it more vulnerable to depurination by RIPs. Exposure to ME(1) in vitro completely abolished the infectivity of partially denatured RNA transcripts of the potato spindle tuber viroid, suggesting that RIPs may target invading nucleic acids before they reach host ribosomes in vivo. Our data suggest that the extensive folding of many potential substrates interferes with their ability to interact with RIPs, thereby blocking their inactivation by ME(1) (or other RIPs).

Molecular recognition of aminoglycoside antibiotics by ribosomal RNA and resistance enzymes: an analysis of x-ray crystal structures

The potential of RNA molecules to be used as therapeutic targets by small inhibitors is now well established. In this fascinating wide-open field, aminoglycoside antibiotics constitute the most studied family of RNA binding drugs. Within the last three years, several x-ray crystal structures were solved for aminoglycosides complexed to one of their main natural targets in the bacterial cell, the decoding aminoacyl-tRNA site (A site). Other crystallographic structures have revealed the binding modes of aminoglycosides to the three existing types of resistance-associated enzymes. The present review summarizes the various aspects of the molecular recognition of aminoglycosides by these natural RNA or protein receptors. The analysis and the comparisons of the detailed interactions offer insights that are helpful in designing new generations of antibiotics.

Tools for the automatic identification and classification of RNA base pairs

Three programs have been developed to aid in the classification and visualization of RNA structure. BPViewer provides a web interface for displaying three-dimensional (3D) coordinates of individual base pairs or base pair collections. A web server, RNAview, automatically identifies and classifies the types of base pairs that are formed in nucleic acid structures by various combinations of the three edges, Watson-Crick, Hoogsteen and the Sugar edge. RNAView produces two-dimensional (2D) diagrams of secondary and tertiary structure in either Postscript, VRML or RNAML formats. The application RNAMLview can be used to rearrange various parts of the RNAView 2D diagram to generate a standard representation (like the cloverleaf structure of tRNAs) or any layout desired by the user. A 2D diagram can be rapidly reformatted using RNAMLview since all the parts of RNA (like helices and single strands) are dynamically linked while moving the selected parts. With the base pair annotation and the 2D graphic display, RNA motifs are rapidly identified and classified. A survey has been carried out for 41 unique structures selected from the NDB database. The statistics for the occurrence of each edge and of each of the 12 bp families are given for the combinations of the four bases: A, G, U and C. The program also allows for visualization of the base pair interactions by using a symbolic convention previously proposed for base pairs. The web servers for BPViewer and RNAview are available at http://ndbserver.rutgers.edu/services/. The application RNAMLview can also be downloaded from this site. The 2D diagrams produced by RNAview are available for RNA structures in the Nucleic Acid Database (NDB) at http://ndbserver.rutgers.edu/atlas/.

A thermodynamic framework for the magnesium-dependent folding of RNA

The goal of this review is to present a unified picture of the relationship between ion binding and RNA folding based on recent theoretical and computational advances. In particular, we present a model describing how the association of magnesium ions is coupled to the tertiary structure folding of several well-characterized RNA molecules. This model is developed in terms of the nonlinear Poisson-Boltzmann (NLPB) equation, which provides a rigorous electrostatic description of the interaction between Mg(2+) and specific RNA structures. In our description, most of the ions surrounding an RNA behave as a thermally fluctuating ensemble distributed according to a Boltzmann weighted average of the mean electrostatic potential around the RNA. In some cases, however, individual ions near the RNA may shed some of their surrounding waters to optimize their Coulombic interactions with the negatively charged ligands on the RNA. These chelated ions are energetically distinct from the surrounding ensemble and must be treated explicitly. This model is used to explore several different RNA systems that interact differently with Mg(2+). In each case, the NLPB equation accurately describes the stoichiometric and energetic linkage between Mg(2+) binding and RNA folding without requiring any fitted parameters in the calculation. Based on this model, we present a physical description of how Mg(2+) binds and stabilizes specific RNA structures to promote the folding reaction.

Conformational flexibility of two RNA trimers explored by computational tools and database search

Two RNA sequences, AAA and AUG, were studied by the conformational search program CICADA and by molecular dynamics (MD) in the framework of the AMBER force field, and also via thorough PDB database search. CICADA was used to provide detailed information about conformers and conformational interconversions on the energy surfaces of the above molecules. Several conformational families were found for both sequences. Analysis of the results shows differences, especially between the energy of the single families, and also in flexibility and concerted conformational movement. Therefore, several MD trajectories (altogether 16 ns) were run to obtain more details about both the stability of conformers belonging to different conformational families and about the dynamics of the two systems. Results show that the trajectories strongly depend on the starting structure. When the MD start from the global minimum found by CICADA, they provide a stable run, while MD starting from another conformational family generates a trajectory where several different conformational families are visited. The results obtained by theoretical methods are compared with the thorough database search data. It is concluded that all except for the highest energy conformational families found in theoretical result also appear in experimental data. Registry numbers: adenylyl-(3' --> 5')-adenylyl-(3' --> 5')-adenosine [917-44-2] adenylyl-(3' --> 5')-uridylyl-(3' --> 5')-guanosine [3494-35-7].

RNA as a drug target

In the antiviral and antibacterial area, increasing drug resistance means that there is an ever growing need for novel approaches towards structures and mechanisms which avoid the current problems. The huge increase in high resolution structural data is set to make a dramatic impact on targeting RNA as a drug target. The examples of the RNA binding antibiotics, particularly, the totally synthetic oxazolidinones, should help persuade the skceptics that clinically useful, selective drugs can be obtained from targeting RNA directly.

Catalytic strategies of the hepatitis delta virus ribozymes

The hepatitis delta virus (HDV) ribozymes are self-cleaving RNA sequences critical to the replication of a small RNA genome. A recently determined crystal structure together with biochemical and biophysical studies provides new insight into the possible catalytic mechanism of these ribozymes. The HDV ribozymes are examples of naturally occurring small ribozymes that catalyze cleavage of the RNA backbone with a rate enhancement of 10(6)- to 10(7)-fold over the uncatalyzed rate. To achieve this level of rate enhancement, the HDV ribozymes have been proposed to employ several catalytic strategies that include the use of metal ions, intrinsic binding energy, and a novel example of general acid-base catalysis with a cytosine side chain acting as a proton donor or acceptor.

Activity of HDV ribozymes to trans-cleave HCV RNA

AIM: To explore whether HDV ribozymes have the ability to trans-cleave HCV RNA.

METHODS: Three HDV genomic ribozymes were designed and named RzC1, RzC2 and RzC3. The substrate RNA contained HCV RNA 5’-noncoding region and 5'-fragment of C region (5'-NCR-C). All the ribozymes and HCV RNA 5'-NCR-C were obtained by transcription in vitro from their DNA templates, and HCV RNA 5'-NCR-C was radiolabelled at its 5’-end. Under certain pH, temperature, appropriate concentration of Mg²⁺ and deionized formamide, these ribozymes were respectively or simultaneously mixed with HCV RNA 5'-NCR-C and reacted for a certain time. The trans-cleavage reaction was stopped at different time points, and the products were separated with polyacrylamide gel electrophoresis (PAGE), displayed by autoradiography. Percentage of trans-cleaved products was measured to indicate the activity of HDV ribozymes.

RESULTS: RzC1 and RzC2 could trans-cleave 26% and 21.8% of HCV RNA 5'-NCR-C under our reaction conditions with 2.5 mol•L^-1 deionized formamide respectively. The percentage of HCV RNA 5'-NCR-C trans-cleaved by RzC1, RzC2 or combined usage of the three ribozymes increased with time, up to 24.9%, 20.3% and 37.3% respectively at 90 min point. Almost no product from RzC3 was observed.

CONCLUSION: HDV ribozymes are able to trans-cleave specifically HCV RNA at certain sites under appropriate conditions, and combination of several ribozymes aiming at different target sites can trans-cleave the substrate more efficiently than using only one of them.

RNA simulations: probing hairpin unfolding and the dynamics of a GNRA tetraloop

Simulations of an RNA hairpin containing a GNRA tetraloop were conducted to allow the characterization of its secondary structure formation and dynamics. Ten 10 ns trajectories of the folded hairpin 5'-GGGC[GCAA]GCCU-3' were generated using stochastic dynamics and the GB/SA implicit solvent model at 300 K. Overall, we find the stem to be a very stable subunit of this molecule, whereas multiple loop conformations and transitions between them were observed. These trajectories strongly suggest that extension of the C6 base away from the loop occurs cooperatively with an N-type-->S-type sugar pucker conversion in that residue and that similar pucker transitions are necessary to stabilize other looped-out bases. In addition, a short-lived conformer with an extended fourth loop residue (A8) lacking this stabilizing 2'-endo pucker mode was observed. Results of thermal perturbation at 400 K support this model of loop dynamics. Unfolding trajectories were produced using this same methodology at temperatures of 500 to 700 K. The observed unfolding events display three-state behavior kinetically (including native, globular, and unfolded populations) and, based on these observations, we propose a folding mechanism that consists of three distinct events: (i) collapse of the random unfolded structure and sampling of the globular state; (ii) passage into the folded region of configurational space as stem base-pairs form and gain helicity; and (iii) attainment of proper loop geometry and organization of loop pairing and stacking interactions. These results are considered in the context of current experimental knowledge of this and similar nucleic acid hairpins.

The linkage between magnesium binding and RNA folding

Understanding the linkage between Mg(2+) binding and RNA folding requires a proper theoretical model describing the energetics of Mg(2+) binding to the folded and unfolded states of RNA. Our current understanding of Mg(2+) binding to these different RNA states derives from empirical thermodynamic models that depend on a number of unjustified assumptions. We present a rigorous theoretical model describing the linkage between RNA folding and magnesium ion binding. In this model, based on the non-linear Poisson-Boltzmann (NLPB) equation, the stabilization of RNA by Mg(2+) arises from two distinct binding modes, diffuse binding and site binding. Diffusely bound Mg(2+) are described as an ensemble of hydrated ions that are attracted to the negative charge of the RNA. Site-bound Mg(2+) are partially desolvated ions that are attracted to electronegative pockets on the RNA surface. We explore two systems, yeast tRNA(Phe) and a 58-nucleotide rRNA fragment, with different Mg(2+) binding properties. The NLPB equation accurately describes both the stoichiometric and energetic linkage between Mg(2+) binding and RNA folding for both of these systems without requiring any fitted parameters in the calculation. Moreover, the NLPB model presents a well-defined physical description of how Mg(2+) binding helps fold an RNA. For both of the molecules studied here, the relevant unfolded state is a disordered intermediate state (I) that contains stable helical secondary structure without any tertiary contacts. Diffusely bound Mg(2+) interact with these secondary structure elements to stabilize the I state. The secondary structural elements of the I state fold into a compact, native tertiary structure (the N state). Diffuse binding plays a dominant role in stabilizing the N state for both RNAs studied. However, for the rRNA fragment, site-binding to a location with extraordinarily high electrostatic potential is also coupled to folding. Our results suggest that much experimental data measuring the linkage between Mg(2+) binding and RNA folding must be reinterpreted.

Mg2+-dependent conformational change of RNA studied by fluorescence correlation and FRET on immobilized single molecules

Fluorescence correlation spectroscopy (FCS) of fluorescence resonant energy transfer (FRET) on immobilized individual fluorophores was used to study the Mg2+-facilitated conformational change of an RNA three-helix junction, a structural element that initiates the folding of the 30S ribosomal subunit. Transitions of the RNA junction between open and folded conformations resulted in fluctuations in fluorescence by FRET. Fluorescence fluctuations occurring between two FRET states on the millisecond time scale were found to be dependent on Mg2+ and Na+ concentrations. Correlation functions of the fluctuations were used to determine transition rates between the two conformations as a function of Mg2+ or Na+ concentration. Both the opening and folding rates were found to vary with changing salt conditions. Assuming specific binding of divalent ions to RNA, the Mg2+ dependence of the observed rates cannot be explained by conformational change induced by Mg2+ binding/unbinding, but is consistent with a model in which the intrinsic conformational change of the RNA junction is altered by uptake of Mg2+ ion(s). This version of FCS/FRET on immobilized single molecules is demonstrated to be a powerful technique in the study of conformational dynamics of biomolecules over time scales ranging from microseconds to seconds.

The hydrogen bonding and hydration of 2'-OH in adenosine and adenosine 3'-ethyl phosphate

The 2'-OH group has major structural implications in the recognition, processing, and catalytic properties of RNA. We report here intra- and intermolecular H-bonding of 2'-OH in adenosine 3'-ethyl phosphate (1), 3'-deoxyadenosine (2), and adenosine (3) by both temperature- and concentration-dependent NMR studies, as well as by detailed endo ((3)J(H,H)) and exocyclic ((3)J(H,OH)) coupling constant analyses. We have also examined the nature of hydration and exchange processes of 2'-OH with water by a combination of NOESY and ROESY experiments in DMSO-d(6) containing 2 mol % HOD. The NMR-constrained molecular modeling (by molecular mechanics as well as by ab initio methods both in the gas and solution phase) has been used to characterize the energy minima among the four alternative dihedrals possible from the solution of the Karplus equation for (3)J(H2',OH) and (3)J(H3',OH) to delineate the preferred orientation of 2'-O-H proton in 1 and 2 as well as for 2'/3'-O-H protons in 3. The NMR line shape analysis of 2'-OH gave the DeltaG(H-bond)(298K) of 7.5 kJ mol(-1) for 1 and 8.4 kJ mol(-1) for 3; similar analyses of the methylene protons of 3'-ethyl phosphate moiety in 1 also gave comparable DeltaG(H-bond)(298K) of 7.3 kJ mol(-1). The donor nature of the 2'-OH in the intramolecular H-bonding in 3 is evident from its relatively reduced flexibility [-TDeltaS++](2'-OH) = -17.9(+/-0.5) kJ mol(-1)] because of the loss of conformational freedom owing to the intramolecular 2'O-H...O3' H-bonding, compared to the acceptor 3'-OH in 3 [-TDeltaS++](3'-OH) = -19.8 (+/- 0.6) kJ mol(-1)] at 298 K. The presence of intramolecular 2'-OH...O3' H-bonding in 3 is also corroborated by the existence of weak long-range (4)J(H2',OH3') in 3 (i.e., W conformation of H2'-C2'-C3'-O3'-H) as well as by (3)J(H,OH) dependent orientation of the 2'- and 3'-OH groups. The ROESY spectra for 1 and 3 at 308 K, in DMSO-d(6), show a clear positive ROE contact of both 2'- and 3'-OH with water. The presence of a hydrophilic 3'-phosphate group in 1 causes a much higher water activity in the vicinity of its 2'-OH, which in turn causes the 2'-OH to exchange faster, culminating in a shorter exchange lifetime (tau) for 2'-OH proton with HOD in 1 (tau2'-OH: 489 ms) compared to that in 3 (tau2'-OH: 6897 ms). The activation energy (E(a)) of the exchange with the bound-water for 2'- and 3'-OH in 3 (48.3 and 45.0 kJ mol(-1), respectively) is higher compared to that of 2'-OH in 1 (31.9 kJ mol(-1)), thereby showing that the kinetic availability of hydrated 2'-OH in 1 for any inter- and intramolecular interactions, in general, is owing to the vicinal 3'-phosphate residue. It also suggests that 2'-OH in native RNA can mediate other inter- or intramolecular interactions only in competition with the bound-water, depending upon the specific chemical nature and spatial orientation of other functions with potential for hydrogen bonding in the neighborhood. This availability of the bound water around 2'-OH in RNA would, however, be dictated by whether the vicinal phosphate is exposed to the bulk water or not. This implies that relatively poor hydration around a specific 2'-OH across a polyribonucleotide chain, owing to some hydrophobic microenvironmental pocket around that hydroxyl, may make it more accessible to interact with other donor or acceptor functions for H-bonding interactions, which might then cause the RNA to fold in a specific manner generating a new motif leading to specific recognition and function. Alternatively, a differential hydration of a specific 2'-OH may modulate its nucleophilicity to undergo stereospecific transesterification reaction as encountered in ubiquitous splicing of pre-mRNA to processed RNA or RNA catalysis, in general.

SCOR: a Structural Classification of RNA database

The Structural Classification of RNA (SCOR) database provides a survey of the three-dimensional motifs contained in 259 NMR and X-ray RNA structures. In one classification, the structures are grouped according to function. The RNA motifs, including internal and external loops, are also organized in a hierarchical classification. The 259 database entries contain 223 internal and 203 external loops; 52 entries consist of fully complementary duplexes. A classification of the well-characterized tertiary interactions found in the larger RNA structures is also included along with examples. The SCOR database is accessible at http://scor.lbl.gov.