RNA structure

Dissecting Functional Biological Interactions Using Modular RNA Nanoparticles

Nucleic acid nanoparticles (NANPs) are an exciting and innovative technology in the context of both basic and biomedical research. Made of DNA, RNA, or their chemical analogs, NANPs are programmed for carrying out specific functions within human cells. NANPs are at the forefront of preventing, detecting, and treating disease. Their nucleic acid composition lends them biocompatibility that provides their cargo with enhanced opportunity for coordinated delivery. Of course, the NANP system of targeting specific cells and tissues is not without its disadvantages. Accumulation of NANPs outside of the target tissue and the potential for off-target effects of NANP-mediated cargo delivery present challenges to research and medical professionals and these challenges must be effectively addressed to provide safe treatment to patients. Importantly, development of NANPs with regulated biological activities and immunorecognition becomes a promising route for developing versatile nucleic acid therapeutics. In a basic research context, NANPs can assist investigators in fine-tuning the structure-function relationship of final formulations and in this review, we explore the practical applications of NANPs in laboratory and clinical settings and discuss how we can use established nucleic acid research techniques to design effective NANPs.

Micrometer-size double-helical structures from phospholipid-modified carbon nanotubes

Biomolecular self-assembly plays a key role in the life system. Herein, double-helical phospholipid-modified carbon nanotube structures were constructed via the self-assembly of phospholipids on carbon nanotubes. These micrometer size spring structures may find potential applications in biocompatible microrobots.

Biotechnological and Therapeutic Applications of Natural Nucleic Acid Structural Motifs

Genetic information and the blueprint of life are stored in the form of nucleic acids. The primary sequence of DNA, read from the canonical double helix, provides the code for RNA and protein synthesis. Yet these already-information-rich molecules have higher-order structures which play critical roles in transcription and translation. Uncovering the sequences, parameters, and conditions which govern the formation of these structural motifs has allowed researchers to study them and to utilize them in biotechnological and therapeutic applications in vitro and in vivo. This review covers both DNA and RNA structural motifs found naturally in biological systems including catalytic nucleic acids, non-coding RNA, aptamers, G-quadruplexes, i-motifs, and Holliday junctions. For each category, an overview of the structural characteristics, biological prevalence, and function will be discussed. The biotechnological and therapeutic applications of these structural motifs are highlighted. Future perspectives focus on the addition of proteins and unnatural modifications to enhance structural stability for greater applicability.

Predicting RNA-Metal Ion Binding with Ion Dehydration Effects

Metal ions play essential roles in nucleic acids folding and stability. The interaction between metal ions and nucleic acids can be highly complicated because of the interplay between various effects such as ion correlation, fluctuation, and dehydration. These effects may be particularly important for multivalent ions such as Mg²⁺ ions. Previous efforts to model ion correlation and fluctuation effects led to the development of the Monte Carlo tightly bound ion model. Here, by incorporating ion hydration/dehydration effects into the Monte Carlo tightly bound ion model, we develop a, to our knowledge, new approach to predict ion binding. The new model enables predictions for not only the number of bound ions but also the three-dimensional spatial distribution of the bound ions. Furthermore, the new model reveals several intriguing features for the bound ions such as the mutual enhancement/inhibition in ion binding between the fully hydrated (diffuse) ions, the outer-shell dehydrated ions, and the inner-shell dehydrated ions and novel features for the monovalent-divalent ion interplay due to the hydration effect.

Comparison of Physicochemical Membrane Properties of Vesicles Modified with Guanidinium Derivatives

Bilayer vesicles have garnered considerable research attention as molecular vehicles capable of noncovalent interaction with biomolecules via electrostatic and hydrophobic bonds and van der Waals interactions. Guanidinium strongly interacts with phosphate groups. Thus, guanidinium modification of vesicles helps intensify the interaction between lipid membranes and nucleic acids. Here, two kinds of guanidinium derivatives, stearylguanidinium (SG) and myristoylarginine (MA), were synthesized and incorporated into 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC) vesicles. Differences in their membrane properties were evaluated using Fourier transform infrared spectroscopy, Raman spectroscopy, and the fluorescent probes 1,6-diphenyl-1,3,5-hexatriene (DPH), 6-lauroyl-2-dimethylaminonaphthalene (Laurdan), and 2-p-toluidinylnaphthalene-6-sulfonate (TNS). The increased SG ratio increased overall hydrophobicity and lipid packing density compared to POPC vesicles, and SG-modified vesicles successfully attracted and then denatured negatively charged tRNAs (tRNAs). In contrast, MA-modified vesicles did not affect the stiffness of POPC membranes, wherein no conformational change in tRNAs was observed in the presence of POPC/MA vesicles. Analyses of the pH-dependent fluorescence emission of TNS suggested that SG and MA molecules render the membrane surfaces cationic and anionic, respectively, which was also revealed by zeta potential measurements. Our results enabled the construction of a model of the headgroup orientation of zwitterionic POPC molecules controlled by modification with guanidinium derivatives. The results also indicate the possibility to regulate the interaction and conformation of biological molecules, such as nucleic acid.

Disclosing 3' UTR cis-elements and putative partners involved in gene expression regulation in Leishmania spp

To identify putative cis-elements involved in gene expression regulation in Leishmania, we previously conducted an in silico investigation to find conserved intercoding sequences (CICS) in the genomes of L. major, L. infantum, and L. braziliensis. Here, the CICS databank was explored to search for sequences that were present in the untranslated regions (UTRs) of groups of genes showing similar expression profiles during in vitro differentiation. Using a selectable marker as a reporter gene, flanked by either an intact 3' UTR or a UTR lacking the conserved element, the regulatory role of a CICS was confirmed. We observed that the pattern of modulation of the mRNA levels was altered in the absence of the CICS. We also identified putative CICS RNA-binding proteins. This study suggests that the publicly available CICS database is a useful tool for identifying regulatory cis-elements for Leishmania genes and suggests the existence of post-transcriptional regulons in Leishmania.

H-Bond Self-Assembly: Folding versus Duplex Formation

Linear oligomers equipped with complementary H-bond donor (D) and acceptor (A) sites can interact via intermolecular H-bonds to form duplexes or fold via intramolecular H-bonds. These competing equilibria have been quantified using NMR titration and dilution experiments for seven systems featuring different recognition sites and backbones. For all seven architectures, duplex formation is observed for homo-sequence 2-mers (AA·DD) where there are no competing folding equilibria. The corresponding hetero-sequence AD 2-mers also form duplexes, but the observed self-association constants are strongly affected by folding equilibria in the monomeric states. When the backbone is flexible (five or more rotatable bonds separating the recognition sites), intramolecular H-bonding is favored, and the folded state is highly populated. For these systems, the stability of the AD·AD duplex is 1-2 orders of magnitude lower than that of the corresponding AA·DD duplex. However, for three architectures which have more rigid backbones (fewer than five rotatable bonds), intramolecular interactions are not observed, and folding does not compete with duplex formation. These systems are promising candidates for the development of longer, mixed-sequence synthetic information molecules that show sequence-selective duplex formation.

Binding interaction of phenothiazinium dyes with double stranded RNAs: Spectroscopic and calorimetric investigation

RNA targeting through small molecules is an emerging and promising therapeutic route that necessitates identification of small molecules that can selectively target specific RNA structures. In this context a comparative study of the interaction of two phenothiazinium dyes thionine (THN) and toluidine blue O (TBO) with three double stranded RNA polynucleotides (ds RNAs) viz. poly(I).poly(C), poly(A).poly(U) and poly(C).poly(G) was conducted by various biophysical techniques. A higher binding of THN with poly(I).poly(C) over poly(A).poly(U) and poly(C).poly(G) was observed. The intercalative binding and RNA induced fluorescence quenching of the dyes through a static mechanism was confirmed by viscosity studies and steady state polarization anisotropy experiments. Binding induced structural perturbation in the RNA polynucleotides was confirmed from circular dichroism spectroscopy. DSC and thermal melting experiments confirmed that the binding resulted in strong thermal stabilization. The binding affinity of THN with poly(I).poly(C) was the highest followed by that to poly(A).poly(U) and poly(C).poly(G). The trend was the same for TBO also, but THN bound stronger than TBO. The binding of the dyes was characterized by strong negative enthalpy changes with minimum positive entropy changes indicating typical intercalative interaction. The results presented here may be useful to design new types of RNA binding antitumor, antibacterial and anticancer agents.

Evidence for widespread association of mammalian splicing and conserved long-range RNA structures

Pre-mRNA structure impacts many cellular processes, including splicing in genes associated with disease. The contemporary paradigm of RNA structure prediction is biased toward secondary structures that occur within short ranges of pre-mRNA, although long-range base-pairings are known to be at least as important. Recently, we developed an efficient method for detecting conserved RNA structures on the genome-wide scale, one that does not require multiple sequence alignments and works equally well for the detection of local and long-range base-pairings. Using an enhanced method that detects base-pairings at all possible combinations of splice sites within each gene, we now report RNA structures that could be involved in the regulation of splicing in mammals. Statistically, we demonstrate strong association between the occurrence of conserved RNA structures and alternative splicing, where local RNA structures are generally more frequent at alternative donor splice sites, while long-range structures are more associated with weak alternative acceptor splice sites. As an example, we validated the RNA structure in the human SF1 gene using minigenes in the HEK293 cell line. Point mutations that disrupted the base-pairing of two complementary boxes between exons 9 and 10 of this gene altered the splicing pattern, while the compensatory mutations that reestablished the base-pairing reverted splicing to that of the wild-type. There is statistical evidence for a Dscam-like class of mammalian genes, in which mutually exclusive RNA structures control mutually exclusive alternative splicing. In sum, we propose that long-range base-pairings carry an important, yet unconsidered part of the splicing code, and that, even by modest estimates, there must be thousands of such potentially regulatory structures conserved throughout the evolutionary history of mammals.

Evidence for a thermodynamically distinct Mg2+ ion associated with formation of an RNA tertiary structure

A folding strategy adopted by some RNAs is to chelate cations in pockets or cavities, where the ions neutralize charge from solvent-inaccessible phosphate. Although such buried Mg(2+)-RNA chelates could be responsible for a significant fraction of the Mg(2+)-dependent stabilization free energy of some RNA tertiary structures, direct measurements have not been feasible because of the difficulty of finding conditions under which the free energy of Mg(2+) chelation is uncoupled from RNA folding and from unfavorable interactions with Mg(2+) ions in other environments. In a 58mer rRNA fragment, we have used a high-affinity thermophilic ribosomal protein to trap the RNA in a structure nearly identical to native; Mg(2+)- and protein-stabilized structures differ in the solvent exposure of a single nucleotide located at the chelation site. Under these conditions, titration of a high affinity chelation site takes place in a micromolar range of Mg(2+) concentration, and is partially resolved from the accumulation of Mg(2+) in the ion atmosphere. From these experiments, we estimate the total and site-specific Mg(2+)-RNA interaction free energies over the range of accessed Mg(2+) concentrations. At 0.1 mM Mg(2+) and 60 mM K(+), specific site binding contributes ∼-3 kcal/mol of the total Mg(2+) interaction free energy of ∼-13 kcal/mol from all sources; at higher Mg(2+) concentrations the site-binding contribution becomes a smaller proportion of the total (-4.5 vs -33 kcal/mol). Under approximately physiological ionic conditions, the specific binding site will be saturated but will provide only a fraction of the total free energy of Mg(2+)-RNA interactions.

Conformational change of single-stranded RNAs induced by liposome binding

The interaction between single-stranded RNAs and liposomes was studied using UV, Fourier Transform Infrared spectroscopy (FTIR) and Circular Dichroism spectroscopy (CD). The effect of the surface characteristics of liposomes, which were composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and modified with cholesterol (Ch) or 1,2-dioleoyl-3-trimethylammonium propane (DOTAP), on the liposome–RNA interaction was investigated. The fluorescence of 6-(p-toluidino)naphthalene-2-sulfonate (TNS) embedded in the liposome surface (ε = 30–40) was decreased in the presence of tRNA, suggesting that single-stranded tRNA could bind onto the liposome. The dehydration of –PO₂⁻ –, guanine (G) and cytosine (C) of tRNA molecules in the presence of liposomes suggested both an electrostatic interaction (phosphate backbone of tRNA and trimethylammonium group of POPC, DOTAP) and a hydrophobic interaction (guanine or cytosine of tRNA and aliphatic tail of lipid). The tRNA conformation on the liposome was determined by CD spectroscopy. POPC/Ch (70/30) maintained tRNA conformation without any denaturation, while POPC/DOTAP(70/30) drastically denatured it. The mRNA translation was evaluated in an Escherichia coli cell-free translation system. POPC/Ch(70/30) enhanced expression of green fluorescent protein (GFP) (116%) while POPC/DOTAP(70/30) inhibited (37%), suggesting that the conformation of RNAs was closely related to the translation efficiency. Therefore, single-stranded RNAs could bind to liposomal membranes through electrostatic and hydrophobic attraction, after which conformational changes were induced depending on the liposome characteristics.

Liposomes destabilize tRNA during heat stress

Biomembranes play an important role in cellular response to heat stress. In this study, we focus on the interaction between liposomes and tRNA. Upon heat treatment we determined circular dichroism spectra of tRNA in presence of liposomes prepared from POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and cholesterol (Ch). To compare thermal stability, midpoint temperature (T(m)) of tRNA was calculated from normalized theta(208). Addition of POPC/Ch liposomes decreased the T(m) value of tRNA from 48 degrees C to 38 degrees C. We conclude that POPC/Ch liposomes interact with tRNA and destabilize its conformation under heat stress.

The PKR-binding domain of adenovirus VA RNAI exists as a mixture of two functionally non-equivalent structures

VA RNA_I is a non-coding adenoviral transcript that counteracts the host cell anti-viral defenses such as immune responses mediated via PKR. We investigated potential alternate secondary structure conformations within the PKR-binding domain of VA RNA_I using site-directed mutagenesis, RNA UV-melting analysis and enzymatic RNA secondary structure probing. The latter data clearly indicated that the wild-type VA RNA_I apical stem can adopt two different conformations and that it exists as a mixed population of these two structures. In contrast, in two sequence variants we designed to eliminate one of the possible structures, while leaving the other intact, each formed a unique secondary structure. This clarification of the apical stem pairing also suggests a small alteration to the apical stem–loop secondary structure. The relative ability of the two apical stem conformations to bind PKR and inhibit kinase activity was measured by isothermal titration calorimetry and PKR autophosphorylation inhibition assay. We found that the two sequence variants displayed markedly different activities, with one being a significantly poorer binder and inhibitor of PKR. Whether the presence of the VA RNA_I conformation with reduced PKR inhibitory activity is directly beneficial to the virus in the cell for some other function requires further investigation.

Self-assembly of single-stranded RNA on carbon nanotube: polyadenylic acid to form a duplex structure

All messenger-RNA (mRNA) molecules in eukaryotic cells have a polyadenylic acid [poly(rA)] tail at the 3'-end and human poly(rA) polymerase (PAP) has been considered as a tumor-specific target. A ligand that is capable of recognizing and binding to the poly(rA) tail of mRNA might interfere with the full processing of mRNA by PAP and can be a potential therapeutic agent. We report here for the first time that single-walled carbon nanotubes (SWNTs) can cause single-stranded poly(rA) to self-structure and form a duplex structure, which is studied by UV melting, atomic force microscopy, circular dichroism spectroscopy, and NMR spectrometry. SWNTs have shown potential applications that range from nanodevices, gene therapy, and drug delivery to membrane separations. Our studies may provide new insights into the application of SWNTs under physiological conditions, possibly being used as probes that target specific gene sequences.

Applications of ESI-MS in drug discovery: interrogation of noncovalent complexes

For many years, analytical mass spectrometry has had numerous supporting roles in the drug development process, including the assessment of compound purity; quantitation of absorption, distribution, metabolism and excretion; and compound-specific pharmacokinetic analyses. More recently, mass spectrometry has emerged as an effective technique for identifying lead compounds on the basis of the characterization of noncovalent ligand-macromolecular target interactions. This approach offers several attractive properties for screening applications in drug discovery compared with other strategies, including the small quantities of target and ligands required, and the capacity to study ligands or targets without having to label them. Here, we review the application of electrospray ionization mass spectrometry to the interrogation of noncovalent complexes, highlighting examples from drug discovery efforts aimed at a range of target classes.

RNA aptamers as potential pharmaceuticals against infections with African trypanosomes

Protozoal pathogens cause symptomatic as well as asymptomatic infections. They have a worldwide impact, which in part is reflected in the long-standing search for antiprotozoal chemotherapy. Unfortunately, effective treatments for the different diseases are by and large not available. This is especially true for African trypanosomiasis, also known as sleeping sickness. The disease is an increasing problem in many parts of sub-Saharan Africa, which is due to the lack of new therapeutics and the increasing resistance against traditional drugs such as melarsoprol, berenil and isometamidium. Considerable progress has been made over the past 10 years in the development of nucleic acid-based drug molecules using a variety of different technologies. One approach is a combinatorial technology that involves an iterative Darwinian-type in vitro evolution process, which has been termed SELEX for "systematic evolution of ligands by exponential enrichment". The procedure is a highly efficient method of identifying rare ligands from combinatorial nucleic acid libraries of very high complexity. It allows the selection of nucleic acid molecules with desired functions, and it has been instrumental in the identification of a number of synthetic DNA and RNA molecules, so-called aptamers that recognize ligands of different chemical origin. Aptamers typically bind their target with high affinity and high specificity and have successfully been converted into pharmaceutically active compounds. Here we summarize the recent examples of the SELEX technique within the context of identifying high-affinity RNA ligands against the surface of the protozoan parasite Trypanosoma brucei, which is the causative agent of sleeping sickness.

Molecular aspects on the specific interaction of cytotoxic plant alkaloid palmatine to poly(A)

The interaction of the protoberberine alkaloid palmatine with single and double stranded structures of poly(A) was studied by various biophysical techniques. Comparative binding studies were also performed with double stranded DNA, t-RNA, poly(C).poly(G), poly(U) and poly(C). The results of competition dialysis, fluorescence, and absorption spectral studies converge to reveal the molecular aspects of the strong and specific binding of palmatine to single stranded poly(A). The binding affinity of palmatine to natural DNA, t-RNA and double stranded poly(A) was weaker while no binding was apparent with single stranded poly(U), poly(C) and double stranded poly(C).poly(G). The strong affinity of the alkaloid to single stranded poly(A) in comparison to the double stranded structure was also revealed from circular dichroic and viscometric studies. The effect of [Na+] ion concentration on the binding process revealed the significant role of electrostatic forces in the complexation. The presence of bound alkaloid also remarkably affected denaturation-renaturation of stacked helical poly(A). The energetics of the strong binding to poly(A) was studied from thermodynamic estimation from van Hoff' analysis of the temperature dependent binding constants and ultra sensitive isothermal titration calorimertry, both suggesting the binding to be exothermic and enthalpy driven. This study provides detailed insight into the binding specificity of the natural alkaloid to single stranded poly(A) over several other single and double stranded nucleic acid structures suggesting its potential as a lead compound for RNA based drug targeting.

RNA specific molecules: cytotoxic plant alkaloid palmatine binds strongly to poly(A)

The cytotoxic plant alkaloid palmatine was found to bind strongly by partial intercalation to single stranded poly(A) structure with binding affinity (Ka) of (8.36+/-0.26) x 10(5) M(-1). The binding of palmatine was characterized to be exothermic and enthalpy driven with one palmatine for every two adenine residues. On the other hand, the binding to the double stranded poly(A) has been found to be significantly weak. This study identifies poly(A) as a potential bio-target for the alkaloid palmatine and its use as a lead compound in antitumor drug screening.

The binding of DNA intercalating and non-intercalating compounds to A-form and protonated form of poly(rC).poly(rG): spectroscopic and viscometric study

Polymorphic RNA conformations may serve as potential targets for structure specific antiviral agents. As an initial step in the development of such drugs, the interaction of a wide variety of compounds which are characterized to bind to DNA through classical or partial intercalation or by mechanism of groove binding, with the A-form and the protonated form of poly(rC).poly(rG), been evaluated by multifaceted spectroscopic and viscometric techniques. Results of this study suggest that (i) ethidium intercalates to the A-form of RNA, but does not intercalate to the protonated form, (ii) methylene blue intercalates to the protonated form of the RNA but does not intercalate to the A-form, (iii) actinomycin D does not bind to either conformations of the RNA, and (iv) berberine binds to the protonated form by partial intercalation process, while its binding to the A-form is very weak. The DNA groove binder distamycin A has much higher affinity to the protonated form of the RNA compared to the A-form and binds to both structures by non-intercalative mechanism. We conclude that the binding affinity characteristics of these DNA binding molecules to the RNA conformations are vastly different and may serve as data for the development of RNA based antiviral drugs.

Specific and non-specific purine trap in the T-loop of normal and suppressor tRNAs

To elucidate the general constraints imposed on the structure of the D and T-loops in functional tRNAs, active suppressor tRNAs were selected in vivo from a combinatorial tRNA gene library in which several nucleotide positions in these loops were randomized. Analysis of the nucleotide sequences of the selected clones demonstrates that most of them contain combination U54-A58 allowing the formation of the standard reverse-Hoogsteen base-pair 54-58 in the T-loop. With only one exception, all these clones fall into two groups, each characterized by a distinct sequence pattern. Analysis of these two groups has allowed us to suggest two different types of nucleotide arrangement in the DT region. The first type, the so-called specific purine trap, is found in 12 individual tRNA clones and represents a generalized version of the standard D-T loop interaction. It consists of purine 18 sandwiched between the reverse-Hoogsteen base-pair U54-A58 and purine 57. The identity of purine 18 is restricted by the specific base-pairing with nucleotide 55. Depending on whether nucleotide 55 is U or G, purine 18 should be, respectively, G or A. The second structural type, the so-called non-specific purine trap, corresponds to the nucleotide sequence pattern found in 16 individual tRNA clones and is described here for the first time. It consists of purine 18 sandwiched between two reverse-Hoogsteen base-pairs U54-A58 and A55-C57 and, unlike the specific purine trap, requires the T-loop to contain an extra eighth nucleotide. Since purine 18 does not form a base-pair in the non-specific purine trap, both purines, G18 and A18, fit to the structure equally well. The important role of both the specific and non-specific purine traps in the formation of the tRNA L-shape is discussed.

Determination of thermodynamic parameters for HIV DIS type loop-loop kissing complexes

The HIV-1 type dimerization initiation signal (DIS) loop was used as a starting point for the analysis of the stability of Watson-Crick (WC) base pairs in a tertiary structure context. We used ultraviolet melting to determine thermodynamic parameters for loop-loop tertiary interactions and compared them with regular secondary structure RNA helices of the same sequences. In 1 M Na+ the loop-loop interaction of a HIV-1 DIS type pairing is 4 kcal/mol more stable than its sequence in an equivalent regular and isolated RNA helix. This difference is constant and sequence independent, suggesting that the rules governing the stability of WC base pairs in the secondary structure context are also valid for WC base pairs in the tertiary structure context. Moreover, the effect of ion concentration on the stability of loop-loop tertiary interactions differs considerably from that of regular RNA helices. The stabilization by Na+ and Mg2+ is significantly greater if the base pairing occurs within the context of a loop-loop interaction. The dependence of the structural stability on salt concentration was defined via the slope of a T(m)/log [ion] plot. The short base-paired helices are stabilized by 8 degrees C/log [Mg2+] or 11 degrees C/log [Na+], whereas base-paired helices forming tertiary loop-loop interactions are stabilized by 16 degrees C/log [Mg2+] and 26 degrees C/log [Na+]. The different dependence on ionic strength that is observed might reflect the contribution of specific divalent ion binding to the preformation of the hairpin loops poised for the tertiary kissing loop-loop contacts.

Study of the Functional Interaction of the 900 Tetraloop of 16S Ribosomal RNA with Helix 24 within the Bacterial Ribosome

The 900 tetraloop that caps helix 27 of 16S ribosomal RNA (rRNA) is amongst the most conserved regions of rRNA. This tetraloop forms a GNRA motif that docks into the minor groove of three base-pairs at the bottom of helix 24 of 16S rRNA in the 30S subunit. Both the tetraloop and its receptor in helix 24 contact the 23S rRNA, forming the intersubunit bridge B2c. Here, we investigated the interaction between the 900 tetraloop and its receptor by genetic complementation. We used a specialized ribosome system in combination with an in vivo instant evolution approach to select mutations in helix 24 compensating for a mutation in the 900 tetraloop (A900G) that severely decreases ribosomal activity, impairing subunit association and translational fidelity. We selected two mutants where the G769-C810 base-pair of helix 24 was substituted with either U-A or C x A. When these mutations in helix 24 were investigated in the context of a wild-type 900 tetraloop, the C x A but not the U-A mutation severely impaired ribosome activity, interfering with subunit association and decreasing translational fidelity. In the presence of the A900G mutation, both mutations in helix 24 increased the ribosome activity to the same extent. Subunit association and translational fidelity were increased to the same level. Computer modeling was used to analyze the effect of the mutations in helix 24 on the interaction between the tetraloop and its receptor. This study demonstrates the functional importance of the interaction between the 900 tetraloop and helix 24.

Roles of protein subunits in RNA-protein complexes: lessons from ribonuclease P

Ribonucleoproteins (RNP) are involved in many essential processes in life. However, the roles of RNA and protein subunits in an RNP complex are often hard to dissect. In many RNP complexes, including the ribosome and the Group II introns, one main function of the protein subunits is to facilitate RNA folding. However, in other systems, the protein subunits may perform additional functions, and can affect the biological activities of the RNP complexes. In this review, we use ribonuclease P (RNase P) as an example to illustrate how the protein subunit of this RNP affects different aspects of catalysis. RNase P plays an essential role in the processing of the precursor to transfer RNA (pre-tRNA) and is found in all three domains of life. While every cell has an RNase P (ribonuclease P) enzyme, only the bacterial and some of the archaeal RNase P RNAs (RNA component of RNase P) are active in vitro in the absence of the RNase P protein. RNase P is a remarkable enzyme in the fact that it has a conserved catalytic core composed of RNA around which a diverse array of protein(s) interact to create the RNase P holoenzyme. This combination of highly conserved RNA and altered protein components is a puzzle that allows the dissection of the functional roles of protein subunits in these RNP complexes.

In vitro selection of high-affinity nucleic acid ligands to parasite target molecules

The logic of using nucleic acids as pharmaceutical reagents is in part based on their capacity to interact with high affinity and specificity with other biological components. Considerable progress has been made over the past 10 years in the development of nucleic acid-based drug molecules using a variety of different technologies. One approach is a combinatorial technology that involves an iterative Darwinian-type in vitro evolution process, which has been termed SELEX for 'systematic evolution of ligands by exponential enrichment'. The procedure is a highly efficient method of identifying rare ligands from combinatorial nucleic acid libraries of very high complexity. It allows the selection of nucleic acid molecules with desired functions and it has been instrumental in the identification of a number of synthetic DNA and RNA molecules, so-called aptamers that recognise ligands of different chemical origin. The method is fast, it does not require special equipment and the selected aptamers typically bind their target with high affinity and high specificity. Here we summarise the recent examples of the SELEX technique within the context of identifying high-affinity ligands against parasite target molecules.

RNA labeling and immobilization for nano-displacement measurement: probing three-dimensional RNA structure

RNA molecules play essential roles in many biological processes, including the storage and transfer of information in the cell. These events are mediated via RNA-protein interactions or by catalytic RNA molecules. It is now recognized that unique RNA folds are associated with biological functions. Therefore, to study the intrinsic structural changes and dynamics which regulate the various functions of RNA, it is necessary to probe its three-dimensional structure in solution. In this respect, using single-molecule methodologies may allow study of native RNA molecules independent of their size and in real time. However, this may require the immobilization of RNA on a surface. Here, we report a novel approach to immobilize RNA on a glass. The procedures involve both chemical and enzymatic modifications of long RNA molecules. In addition, we demonstrate the application of an optical tweezers apparatus to measure the length and, hence, the dynamics of immobilized intact ribosomal RNA molecules as a function of different solution conditions.

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].

The RNA helicase DbpA exhibits a markedly different conformation in the ADP-bound state when compared with the ATP- or RNA-bound states

The motor enzymes that belong to the family of RNA helicases catalyze the strand separation of duplex RNA via ATP hydrolysis. Among these enzymes, Escherichia coli DbpA is a unique RNA helicase because it possesses ATPase-specific activity toward the peptidyl transferase center in 23 S ribosomal RNA. For this reason, it has been the subject of numerous biochemical and structure-function studies. The ATP-stimulated unwinding activity of DbpA toward specific and nonspecific RNA duplexes has been demonstrated. However, the underlying molecular and structural basis, which facilitates its helicase activities, is presently not known. We combined time-dependent limited proteolysis digestion, fluorescence spectroscopy, and three-dimensional structural homology modeling techniques to study the structural conformations of DbpA with respect to its binding to stoichiometric ratios of RNA and cofactors. We show that the conformational state of DbpA is markedly different in the ADP-bound state than in any other state (ATP- or RNA-bound). These results, together with structural homology studies, suggest that a hinge region located in the core domain of DbpA mediates such conformational changes.

Functional studies of the 900 tetraloop capping helix 27 of 16S ribosomal RNA

The 900 tetraloop (positions 898-901) of Escherichia coli 16S rRNA caps helix 27, which is involved in a conformational switch crucial for the decoding function of the ribosome. This tetraloop forms a GNRA motif involved in intramolecular RNA-RNA interactions with its receptor in helix 24 of 16S rRNA. It is involved also in an intersubunit bridge, via an interaction with helix 67 in domain IV of 23S rRNA. Using a specialized ribosome system and an instant-evolution procedure, the four nucleotides of this loop were randomized and 15 functional mutants were selected in vivo. Positions 899 and 900, responsible for most of the tetraloop/receptor interactions, were found to be the most critical for ribosome activity. Functional studies showed that mutations in the 900 tetraloop impair subunit association and decrease translational fidelity. Computer modeling of the mutations allows correlation of the effect of mutations with perturbations of the tetraloop/receptor interactions.

Molecular dynamics simulations of the denaturation and refolding of an RNA tetraloop

Tetraloops are very abundant structural elements of RNA that are formed by four nucleotides in a hairpin loop which is closed by a double stranded helical stem with some Watson-Crick base pairs. A tetraloop r(GCGAAGGC) was identified from the crystal structure of the central domain of 16S rRNA (727-730) in the Thermus thermophilus 30S ribosomal complex. The crystal structure of the 30S complex includes a total of 104 nucleotides from the central domain of the 16S rRNA and three ribosomal proteins S6, S15 and S18. Independent biochemical experiments have demonstrated that protein S15 plays the role in initiating the formation of the central domain of this complex. In the crystal, the tetraloop interacts with the protein S15 at two sites: one of them is associated with hydrogen bond interactions between residue His50 and nucleotide G730, and the other is associated with the occurrence of residue Arg53 beside A728. This paper uses molecular dynamics (MD) simulations to investigate the protein-dependent structural stability of the tetraloop and demonstrates the folding pathway of this tetraloop via melting denaturation and its subsequent refolding. Three important results are derived from these simulations: (i) The stability of nucleotide A728 appears to be protein dependent. Without the interaction with S15, A728 flips away from stacking with A729. (ii) The melting temperature demonstrated by the simulations is analogous to the results of thermodynamic experiments. In addition, the simulated folding of the tetraloop is stepwise: the native shape of the backbone is formed first; this is then followed by the formation of the Watson- Crick base pairs in the stem; and finally the hydrogen bonds and base stacking in the loop are formed. (iii) The tetraloop structure is similar to the crystal structure at salt concentrations of 0.1 M and 1.0 M used for the simulations, but the refolded structure at 0.1 M salt is more comparable to the crystal structure than at 1.0 M. The results from the simulations using both the Generalized Born continuum model and explicit solvent model (Particle Mesh Ewald) generate a similar pathway for unfolding/refolding of the tetraloop.

RNA conformation and folding studied with fluorescence resonance energy transfer

Fluorescence resonance energy transfer (FRET) results from nonradiative coupling of two fluorophores and reports on distances in the range 10-100 A. It is therefore a suitable probe to determine distances in RNA molecules and define their global structure, to follow kinetics of RNA conformational changes during folding in real time, to monitor ion binding, or to analyze conformational equilibria and assess the thermodynamic stability of tertiary structure conformers. Along with the basic principles of steady-state and time-resolved fluorescence resonance energy transfer measurements, approaches to investigate RNA conformational transitions and folding are described and illustrated with selected examples. The versatility of FRET-based techniques has recently been demonstrated by implementations of FRET in high-throughput screening of potential drugs as well as studies of energy transfer that monitor RNA conformational changes on the single-molecule level.

HyPaLib: a database of RNAs and RNA structural elements defined by hybrid patterns

The database, called HyPaLib (for Hybrid Pattern Library), contains annotated structural elements characteristic for certain classes of structural and/or functional RNAs. These elements are described in a language specifically designed for this purpose. The language allows convenient specification of hybrid patterns, i.e. motifs consisting of sequence features and structural elements together with sequence similarity and thermodynamic constraints. We are currently developing software tools that allow a user to search sequence databases for any pattern in HyPaLib, thus providing functionality which is similar to PROSITE, but dedicated to the more complex patterns in RNA sequences. HyPaLib is available at http://bibiserv. techfak.uni-bielefeld.de/HyPa/.

Quantitative studies of Mn(2+)-promoted specific and non-specific cleavages of a large RNA: Mn(2+)-GAAA ribozymes and the evolution of small ribozymes

Manganese (Mn(2+)) promotes specific cleavage at two major (I and III) and four minor (II, IV, V and VI) sites, in addition to slow non-specific cleavage, in a 659-nucleotide RNA containing the Cr.LSU group I intron. The specific cleavages occurred between G and AAA sequences and thus can be considered Mn(2+)-GAAA ribozymes. We have estimated rates of specific and non-specific cleavages under different conditions. Comparisons of the rates of major-specific and background cleavages gave a maximal specificity of approximately 900 for GAAA cleavage. Both specific and non-specific cleavages showed hyperbolic kinetics and there was no evidence of cooperativity with Mn(2+) concentration. Interestingly, at site III, Mg(2+) alone promoted weak, but the same specific cleavage as Mn(2+). When added with Mn(2+), Mg(2+) had a synergistic effect on cleavage at site III, but inhibited cleavage at the other sites. Mn(2+) cleavage at site III also exhibited lower values of K (Mn(2+) requirement), pH-dependency and activation energy than did cleavage at the other sites. In contrast, the pH-dependency and activation energy for cleavage at site I was similar to non-specific cleavage. These results increase our understanding of the Mn(2+)-GAAA ribozyme. The implications for evolution of small ribozymes are also discussed.

Tertiary core rearrangements in a tight binding transfer RNA aptamer

Guided by an in vitro selection experiment designed to obtain tight binding aptamers of Escherichia coli glutamine specific tRNA (tRNAGln) for glutaminyl-tRNA synthetase (GlnRS), we have engineered a tRNA mutant in which the five-nucleotide variable loop sequence 5'-44CAUUC48-3' is replaced by 5'-44AGGU48-3'. This mutant tRNA binds to GlnRS with 30-fold improved affinity compared to the wild type. The 2.7 A cocrystal structure of the RNA aptamer-GlnRS complex reveals major rearrangements in the central tertiary core of the tRNA, while maintaining an RNA-protein interface identical to the wild type. The repacked RNA core features a novel hydrogen bonding arrangement of the trans Levitt pair G15-U48, a new sulfate binding pocket in the major groove, and increased hydrophobic stacking interactions among the bases. These data suggest that enhanced protein binding to a mutant globular RNA can arise from stabilization of RNA tertiary interactions rather than optimization of RNA-protein contacts.

Co-conservation of rRNA tetraloop sequences and helix length suggests involvement of the tetraloops in higher-order interactions

Terminal loops containing four nucleotides (tetraloops) are common in structural RNAs, and they frequently conform to one of three sequence motifs, GNRA, UNCG, or CUUG. Here we compare available sequences and secondary structures for rRNAs from bacteria, and we show that helices capped by phylogenetically conserved GNRA loops display a strong tendency to be of conserved length. The simplest interpretation of this correlation is that the conserved GNRA loops are involved in higher-order interactions, intramolecular or intermolecular, resulting in a selective pressure for maintaining the lengths of these helices. A small number of conserved UNCG loops were also found to be associated with conserved length helices, consistent with the possibility that this type of tetraloop also takes part in higher-order interactions.

ConStruct: a tool for thermodynamic controlled prediction of conserved secondary structure

A tool for prediction of conserved secondary structure of a set of homologous single-stranded RNAs is presented. For each RNA of the set the structure distribution is calculated and stored in a base pair probability matrix. Gaps, resulting from a multiple sequence alignment of the RNA set, are introduced into the individual probability matrices. These 'aligned' probability matrices are summed up to give a consensus probability matrix emphasizing the conserved structural elements of the RNA set. Because the multiple sequence alignment is independent of any structural constraints, such an alignment may result in introduction of gaps into the homologous probability matrices that disrupt a common consensus structure. By use of its graphical user interface the presented tool allows the removal of such misalignments, which are easily recognized, from the individual probability matrices by optimizing the sequence alignment with respect to a structural alignment. From the consensus probability matrix a consensus structure is extracted, which is viewable in three different graphical representations. The functionality of the tool is demonstrated using a small set of U7 RNAs, which are involved in 3'-end processing of histone mRNA precursors. Supplementary Material lists further results obtained. Advantages and drawbacks of the tool are discussed in comparison to several other algorithms.

Identification of a cis-acting dendritic targeting element in MAP2 mRNAs

In neurons, a limited number of mRNAs have been identified in dendritic processes, whereas other transcripts are restricted to the cell soma. Here we have investigated the molecular mechanisms underlying extrasomatic localization of mRNAs encoding microtubule-associated protein 2 (MAP2) in primary neuronal cultures. Vectors expressing recombinant mRNAs were introduced into hippocampal and sympathetic neurons using DNA transfection and microinjection protocols, respectively. Chimeric mRNAs containing the entire 3' untranslated region of MAP2 transcripts fused to a nondendritic reporter mRNA are detected in dendrites. In contrast, RNAs containing MAP2 coding and 5' untranslated regions or tubulin sequences are restricted to the cell soma. Moreover, 640 nucleotides from the MAP2 3' untranslated region (UTR) are both sufficient and essential for extrasomatic localization of chimeric mRNAs in hippocampal and sympathetic neurons. Thus, a cis-acting dendritic targeting element that is effective in two distinct neuronal cell types is contained in the 3' UTR of MAP2 transcripts. The observation of RNA granules in dendrites implies that extrasomatic transcripts seem to assemble into multimolecular complexes that may function as transport units.

How RNA folds

We describe the RNA folding problem and contrast it with the much more difficult protein folding problem. RNA has four similar monomer units, whereas proteins have 20 very different residues. The folding of RNA is hierarchical in that secondary structure is much more stable than tertiary folding. In RNA the two levels of folding (secondary and tertiary) can be experimentally separated by the presence or absence of Mg2+. Secondary structure can be predicted successfully from experimental thermodynamic data on secondary structure elements: helices, loops, and bulges. Tertiary interactions can then be added without much distortion of the secondary structure. These observations suggest a folding algorithm to predict the structure of an RNA from its sequence. However, to solve the RNA folding problem one needs thermodynamic data on tertiary structure interactions, and identification and characterization of metal-ion binding sites. These data, together with force versus extension measurements on single RNA molecules, should provide the information necessary to test and refine the proposed algorithm.

Thermodynamics of stabilization of RNA pseudoknots by cobalt(III) hexaammine

Equilibrium unfolding (folding) studies reveal that the autoregulatory RNA pseudoknots derived from the bacteriophage T2 and T4 gene 32 mRNAs exhibit significant stabilization by increasing concentrations of divalent metal ions in solution. In this report, the apparent affinities of exchange inert trivalent Co(NH(3))(3+)(6) have been determined, relative to divalent Mg(2+), for the folded, partially folded (K(f)), and fully unfolded (K(u)) conformations of these molecules. A general nonspecific, delocalized ion binding model was developed and applied to the analysis of the metal ion concentration dependence of individual two-state unfolding transitions. Trivalent Co(NH(3))(3+)(6) was found to associate with the fully folded and partially unfolded pseudoknotted forms of these RNAs with a K(f) of 5-8 x 10(4) M(-1) in a background of 0.10 M K(+), or 3- to 5-fold larger than the K(f) obtained for two model RNA hairpins and hairpin unfolding intermediates, and approximately 40-50-fold larger than K(f) for Mg(2+). The magnitude of K(f) was found to be strongly dependent on the monovalent salt concentration in a manner qualitatively consistent with polyelectrolyte theory, with K(f) reaching 1.2 x 10(5) M(-1) in 50 mM K(+). Two RNA hairpins were found to have affinities for Co(NH(3))(3+)(6) and Ru(NH(3))(3+)(6) of 1-2 x10(4) M(-1), or approximately 15-fold larger than the K(f) of approximately 1000 M(-1) observed for Mg(2+). Additionally, the K(u) of 4,800 M(-1) for the trivalent ligands is approximately 8-fold larger than the K(u) of 600 M(-1) observed for Mg(2+). These findings suggest that the T2 and T4 gene 32 mRNA pseudoknots possess a site(s) for Mg(2+) and Co(NH(3))(3+)(6) binding of significantly higher affinity than a "duplexlike" delocalized ion binding site that is strongly linked to the thermodynamic stability of these molecules. Imino proton perturbation nmr spectroscopy suggests that this site(s) lies near the base of the pseudoknot stem S2, near a patch of high negative electrostatic potential associated with the region where the single loop L1 adenosine crosses the major groove of stem S2.

Direct identification of NH...N hydrogen bonds in non-canonical base pairs of RNA by NMR spectroscopy

It is shown that the recently developed quantitative J(NN)HNN-COSY experiment can be used for the direct identification of hydrogen bonds in non-canonical base pairs in RNA. Scalar(2h)J(NN)couplings across NH.N hydrogen bonds are observed in imino hydrogen bonded GA base pairs of the hpGA RNA molecule, which contains a tandem GA mismatch, and in the reverse Hoogsteen AU base pairs of the E-loop of Escherichia coli 5S rRNA. These scalar couplings correlate the imino donor(15)N nucleus of guanine or uridine with the acceptor N1 or N7 nucleus of adenine. The values of the corresponding(2h)J(NN)coupling constants are similar in size to those observed in Watson-Crick base pairs. The reverse Hoogsteen base pairs could be directly detected for the E-loop of E.coli 5S rRNA both in the free form and in a complex with the ribosomal protein L25. This supports the notion that the E-loop is a pre-folded RNA recognition site that is not subject to significant induced conformational changes. Since Watson-Crick GC and AU base pairs are also readily detected the HNN-COSY experiment provides a useful and sensitive tool for the rapid identification of RNA secondary structure elements.

Structure and thermodynamics of metal binding in the P5 helix of a group I intron ribozyme

The solution structure of an RNA hairpin modelling the P5 helix of a group I intron, complexed with Co(NH3)63+, has been determined by nuclear magnetic resonance. Co(NH3)63+, which possesses a geometry very close to Mg(H2O)62+, was used to identify and characterize a Mg2+binding site in the RNA. Strong and positive intermolecular nuclear Overhauser effect (NOE) cross-peaks define a specific complex in which the Co(NH3)63+molecule is in the major groove of tandem G.U base-pairs. The structure of the RNA is characterized by a very low twist angle between the two G.U base-pairs, providing a flat and narrowed major groove. The Co(NH3)63+, although highly localized, is free to rotate to hydrogen bond in several ways to the O4 atoms of the uracil bases and to N7 and O6 of the guanine bases. Negative and small NOE cross-peaks to other protons in the sequence reveal a non-specific or delocalized interaction, characterized by a high mobility of the cobalt ion. Mn2+titrations of P5 show specific broadening of protons of the G.U base-pairs that form the metal ion binding site, in agreement with the NOE data from Co(NH3)63+. Binding constants for the interaction of Co(NH3)63+and of Mg2+to P5 were determined by monitoring imino proton chemical shifts during titration of the RNA with the metal ions. Dissociation constants are on the order of 0.1 mM for Co(NH3)63+and 1 mM for Mg2+. Binding studies were done on mutants with sequences corresponding to the three orientations of tandem G.U base-pairs. The affinities of Co(NH3)63+and Mg2+for the tandem G.U base-pairs depend strongly on their sequences; the differences can be understood in terms of the different structures of the corresponding metal ion-RNA complexes. Substitution of G.C or A.U for G.U pairs also affected the binding, as expected. These structural and thermodynamic results provide systematic new information about major groove metal ion binding in RNA.