Comparison between CUUG and UUCG tetraloops: thermodynamic stability and structural features analyzed by UV absorption and vibrational spectroscopy

Prediction of 3D RNA Structures from Sequence Using Energy Landscapes of RNA Dimers: Application to RNA Tetraloops

Access to the three-dimensional structure of RNA enables an ability to gain a more profound understanding of its biological mechanisms, as well as the ability to design RNA-targeting drugs, which can take advantage of the unique chemical environment imposed by a folded RNA structure. Due to the dynamic and structurally complex properties of RNA, both experimental and traditional computational methods have difficulty in determining RNA's 3D structure. Herein, we introduce TAPERSS (Theoretical Analyses, Prediction, and Evaluation of RNA Structures from Sequence), a physics-based fragment assembly method for predicting 3D RNA structures from sequence. Using a fragment library created using discrete path sampling calculations of RNA dinucleoside monophosphates, TAPERSS can sample the physics-based energy landscapes of any RNA sequence with relatively low computational complexity. We have benchmarked TAPERSS on 21 RNA tetraloops, using a combinatorial algorithm as a proof-of-concept. We show that TAPERSS was successfully able to predict the apo-state structures of all 21 RNA hairpins, with 16 of those structures also having low predicted energies as well. We demonstrate that TAPERSS performs most accurately on GNRA-like tetraloops with mostly stacked loop-nucleotides, while having limited success with more dynamic UNCG and CUYG tetraloops, most likely due to the influence of the RNA force field used to create the fragment library. Moreover, we show that TAPERSS can successfully predict the majority of the experimental non-apo states, highlighting its potential in anticipating biologically significant yet unobserved states. This holds great promise for future applications in drug design and related studies. With discussed improvements and implementation of more efficient sampling algorithms, we believe TAPERSS may serve as a useful tool for a physics-based conformational sampling of large RNA structures.

Single Functionalized pRNA/Gold Nanoparticle for Ultrasensitive MicroRNA Detection Using Electrochemical Surface-Enhanced Raman Spectroscopy

Controlling the selective one-to-one conjugation of RNA with nanoparticles is vital for future applications of RNA nanotechnology. Here, the monofunctionalization of a gold nanoparticle (AuNP) with a single copy of RNA is developed for ultrasensitive microRNA-155 quantification using electrochemical surface-enhanced Raman spectroscopy (EC-SERS). A single AuNP is conjugated with one copy of the packaging RNA (pRNA) three-way junction (RNA 3WJ). pRNA 3WJ containing one strand of the 3WJ is connected to a Sephadex G100 aptamer and a biotin group at each arm (SEPapt/3WJ/Bio) which is then immobilized to the Sephadex G100 resin. The resulting complex is connected to streptavidin-coated AuNP (STV/AuNP). Next, the STV/AuNP-Bio/3WJa is purified and reassembled with another 3WJ to form a single-labeled 3WJ/AuNP. Later, the monoconjugate is immobilized onto the AuNP-electrodeposited indium tin oxide coated substrate for detecting microRNA-155 based on EC-SERS. Application of an optimum potential of +0.2 V results in extraordinary amplification (≈7 times) of methylene blue (reporter) SERS signal compared to the normal SERS signal. As a result, a highly sensitive detection of 60 × 10-18 m microRNA-155 in 1 h in serum based on monoconjugated AuNP/RNA is achieved. Thus, the monofunctionalization of RNA onto nanoparticle can provide a new methodology for biosensor construction and diverse RNA nanotechnology development.

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

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

Biophysical properties, thermal stability and functional impact of 8-oxo-7,8-dihydroguanine on oligonucleotides of RNA-a study of duplex, hairpins and the aptamer for preQ1 as models

A better understanding of the effects that oxidative lesions have on RNA is of importance to understand their role in the development/progression of disease. 8-oxo-7,8-dihydroguanine was incorporated into RNA to understand its structural and functional impact on RNA:RNA and RNA:DNA duplexes, hairpins and pseudoknots. One to three modifications were incorporated into dodecamers of RNA [AAGAGGGAUGAC] resulting in thermal destabilization (ΔT_m – 10°C per lesion). Hairpins with tetraloops c-UUCG*-g* (8-10), a-ACCG-g* (11-12), c-UUG*G*-g* (13-16) and c-ACG*G*-g* (17-20) were modified and used to determine thermal stabilities, concluding that: (i) modifying the stem leads to destabilization unless adenosine is the opposing basepair of 8-oxoGua; (ii) modification at the loop is position- and sequence-dependent and varies from slight stabilization to large destabilization, in some cases leading to formation of other secondary structures (hairpin→duplex). Functional effects were established using the aptamer for preQ₁ as model. Modification at G5 disrupted the stem P1 and inhibited recognition of the target molecule 7-methylamino-7-deazaguanine (preQ₁). Modifying G11 results in increased thermal stability, albeit with a K_d 4-fold larger than its canonical analog. These studies show the capability of 8-oxoG to affect structure and function of RNA, resulting in distinct outcomes as a function of number and position of the lesion.

Hemispherical platinum : silver core : shell nanoparticles for miRNA detection

Novel electrochemical detection based on regioselective functionalised electrocatalytic nanoparticles suitable for the detection of low-abundance molecular biomarkers, miR-132 of neuroblastoma.

Nuclear magnetic resonance structure of the III-IV-V three-way junction from the Varkud satellite ribozyme and identification of magnesium-binding sites using paramagnetic relaxation enhancement

The VS ribozyme is a catalytic RNA found within some natural isolates of Neurospora that is being used as a model system to improve our understanding of RNA structure, catalysis, and engineering. The catalytic domain contains five helical domains (SLII-SLVI) that are organized by two three-way junctions. The III-IV-V junction is required for high-affinity binding of the substrate domain (SLI) through formation of a kissing loop interaction with SLV. Here, we determine the high-resolution nuclear magnetic resonance (NMR) structure of a 47-nucleotide RNA containing the III-IV-V junction (J345). The J345 RNA adopts a Y-shaped fold typical of the family C three-way junctions, with coaxial stacking between stems III and IV and an acute angle between stems III and V. The NMR structure reveals that the core of the III-IV-V junction contains four stacked base triples, a U-turn motif, a cross-strand stacking interaction, an A-minor interaction, and a ribose zipper. In addition, the NMR structure shows that the cCUUGg tetraloop used to stabilize stem IV adopts a novel RNA tetraloop fold, different from the known gCUUGc tetraloop structure. Using Mn(2+)-induced paramagnetic relaxation enhancement, we identify six Mg(2+)-binding sites within J345, including one associated with the cCUUGg tetraloop and two with the junction core. The NMR structure of J345 likely represents the conformation of the III-IV-V junction in the context of the active VS ribozyme and suggests that this junction functions as a dynamic hinge that contributes to substrate recognition and catalysis. Moreover, this study highlights a new role for family C three-way junctions in long-range tertiary interactions.

Investigating the thermodynamics of UNCG tetraloops using infrared spectroscopy

Using infrared (IR) absorption spectroscopy, we have explored the folding thermodynamics of the UNCG class of RNA hairpin tetraloops (N = U, A, C, or G). Without the need to introduce non-native probes, IR spectroscopy makes it possible to distinguish specific structural elements such as base pairing versus base stacking or loop versus stem motions. Our results show that different structural components exhibit different thermodynamics. Specifically, we have found that tetraloop melting proceeds in a thermally sequential fashion, where base pairing in the stem is disrupted before (i.e., at a lower temperature) base stacking along the entire chain. In addition, for N = A, our data argue that the structure immediately surrounding adenine is particularly stable and melts at a higher temperature than either base-pairing or base-stacking interactions. Taken together, these results suggest that hairpin loop formation is not a simple two-state process, even in the equilibrium limit.

Time-resolved events on the reaction pathway of transcript initiation by a single-subunit RNA polymerase: Raman crystallographic evidence

The nucleotidyl transfer reaction leading to formation of the first phosphodiester bond has been followed in real time by Raman microscopy, as it proceeds in single crystals of the N4 phage virion RNA polymerase (RNAP). The reaction is initiated by soaking nucleoside triphosphate (NTP) substrates and divalent cations into the RNAP and promoter DNA complex crystal, where the phosphodiester bond formation is completed in about 40 min. This slow reaction allowed us to monitor the changes of the RNAP and DNA conformations as well as bindings of substrate and metal through Raman spectra taken every 5 min. Recently published snapshot X-ray crystal structures along the same reaction pathway assisted the spectroscopic assignments of changes in the enzyme and DNA, while isotopically labeled NTP substrates allowed differentiation of the Raman spectra of bases in substrates and DNA. We observed that substrates are bound at 2-7 min after soaking is commenced, the O-helix completes its conformational change, and binding of both divalent metals required for catalysis in the active site changes the conformation of the ribose triphosphate at position +1. These are followed by a slower decrease of NTP triphosphate groups due to phosphodiester bond formation that reaches completion at about 15 min and even slower complete release of the divalent metals at about 40 min. We have also shown that the O-helix movement can be driven by substrate binding only. The kinetics of the in crystallo nucleotidyl transfer reaction revealed in this study suggest that soaking the substrate and metal into the RNAP-DNA complex crystal for a few minutes generates novel and uncharacterized intermediates for future X-ray and spectroscopic analysis.

Regulation of splicing enhancer activities by RNA secondary structures

In this report, we studied the effect of RNA structures on the activity of exonic splicing enhancers on the SMN1 minigene model by engineering known ESEs into different positions of stable hairpins. We found that as short as 7-bp stem is sufficient to abolish the enhancer activity. When placing ESEs in the loop region, AG-rich ESEs are fully active, but a UCG-rich ESE is not because of additional structural constraints. ESEs placed adjacent to the 3' end of the hairpin structure display high enhancer activity, regardless of their sequence identities. These rules explain the suppression of multiple ESEs by point mutations that result in a stable RNA structure, and provide an additional mechanism for the C6T mutation in SMN2.

Heterochrony in the germ ring closure and tail bud formation in embryonic development of rainbow trout (Oncorhynchus mykiss)

Due to their large yolk size, salmonid embryos take a longer time for epiboly movements and germ ring closure compared with most other teleost species. Here we analyzed the germ ring closure, tail bud formation and development of the notochord and somites in rainbow trout using live embryo imaging and in situ hybridization with the rt-ntl probe. Rt-ntl is expressed in the germ ring (blastula, gastrula and somitogenesis stage), notochord, tail bud and somites (somitogenesis stage). When epiboly covers half the yolk, a tail bud-like structure is formed and somitogenesis starts. By the time epiboly is completed, the yolk covered and the germ ring closed, the embryo has already reached the 20 somite stage. Therefore, the timing of germ ring closure and tail bud formation is reversed in trout embryos compared with zebrafish and other small model fish embryos (heterochrony). Based on this result, we re-examined the definition of tail bud formation.

Arrangement of 3D structural motifs in ribosomal RNA

Structural 3D motifs in RNA play an important role in the RNA stability and function. Previous studies have focused on the characterization and discovery of 3D motifs in RNA secondary and tertiary structures. However, statistical analyses of the distribution of 3D motifs along the RNA appear to be lacking. Herein, we present a novel strategy for evaluating the distribution of 3D motifs along the RNA chain and those motifs whose distributions are significantly non-random are identified. By applying it to the X-ray structure of the large ribosomal subunit from Haloarcula marismortui, helical motifs were found to cluster together along the chain and in the 3D structure, whereas the known tetraloops tend to be sequentially and spatially dispersed. That the distribution of key structural motifs such as tetraloops differ significantly from a random one suggests that our method could also be used to detect novel 3D motifs of any size in sufficiently long/large RNA structures. The motif distribution type can help in the prediction and design of 3D structures of large RNA molecules.

Finding 3D motifs in ribosomal RNA structures

The identification of small structural motifs and their organization into larger subassemblies is of fundamental interest in the analysis, prediction and design of 3D structures of large RNAs. This problem has been studied only sparsely, as most of the existing work is limited to the characterization and discovery of motifs in RNA secondary structures. We present a novel geometric method for the characterization and identification of structural motifs in 3D rRNA molecules. This method enables the efficient recognition of known 3D motifs, such as tetraloops, E-loops, kink-turns and others. Furthermore, it provides a new way of characterizing complex 3D motifs, notably junctions, that have been defined and identified in the secondary structure but have not been analyzed and classified in three dimensions. We demonstrate the relevance and utility of our approach by applying it to the Haloarcula marismortui large ribosomal unit. Pending the implementation of a dedicated web server, the code accompanying this article, written in JAVA, is available upon request from the contact author.

Raman and Raman optical activity (ROA) analysis of RNA structural motifs in Domain I of the EMCV IRES

Raman and Raman optical activity (ROA) spectra were collected for four RNA oligonucleotides based on the EMCV IRES Domain I to assess the contributions of helix, GNRA tetraloop, U.C mismatch base pair and pyrimidine-rich bulge structures to each. Both Raman and ROA spectra show overall similarities for all oligonucleotides, reflecting the presence of the same base paired helical regions and GNRA tetraloop in each. Specific bands are sensitive to the effect of the mismatch and asymmetric bulge on the structure of the RNA. Raman band changes are observed that reflect the structural contexts of adenine residues, disruption of A-form helical structure, and incorporation of pyrimidine bases in non-helical regions. The ROA spectra are also sensitive to conformational mobility of ribose sugars, and verify a decrease in A-type helix content upon introduction of the pyrimidine-rich bulge. Several Raman and ROA bands also clearly show cooperative effects between the mismatch and pyrimidine-rich bulge motifs on the structure of the RNA. The complementary nature of Raman and ROA spectra provides detailed and highly sensitive information about the local environments of bases, and secondary and tertiary structures, and has the potential to yield spectral signatures for a wide range of RNA structural motifs.

Conformational analysis of the telomerase RNA pseudoknot hairpin by Raman spectroscopy

We have measured the temperature-dependent Raman spectra of two 30-mer ribonucleotides that represent the wild-type (WT) and dyskeratosis congenita (DKC) mutant (MT) GC (107-108) --> AG structures of the pseudoknot hairpin region of human telomerase RNA. We have used these structures, previously characterized by UV-melting and NMR, as a model system for our Raman investigation. We observe that Raman hypochromism of vibrational bands, previously assigned to specific bases or conformational RNA markers, reflect temperature-dependent alterations in the pentaloop and stem structures of these two oligonucleotides. We also observe that the intense nu(s)(O-P-O) band at 812 cm(-1) indicates the presence of A-form backbone structure at relatively low temperatures in both the WT and MT RNA sequences. The mutation induces a decrease in the intensity of the uridine (rU) band at 1244 cm(-1) associated with C2'-endo/anti ribose conformation in the pentaloop. Two transition temperatures (T(m) ) were determined from the analysis of Raman difference intensity-temperature profiles of the 1256 cm(-1) band, which is associated with vibrations of cytidine (rC) residues, in particular, the C2'-endo/anti ribose conformation (T(m) 1 = 23.6 +/- 1.6 degrees C for WT and 19.7 +/- 2.8 degrees C for MT; T(m) 2 = 68.9 +/- 1.8 degrees C for WT and 70.9 +/- 1.1 degrees C for MT). From these results we can conclude that the DKC mutant 30-mer exhibits a lower stability in the pentaloop region and a slightly higher stability in the stem region than the WT 30-mer. This demonstrates that Raman bands, previously assigned to specific bases or conformational RNA markers, can be used to probe local structural features of the telomerase pseudoknot hairpin sequence.

Modulating role of RNA structure in alternative splicing of a critical exon in the spinal muscular atrophy genes

Humans have two nearly identical copies of the survival motor neuron (SMN ) gene, SMN1 and SMN2. Homozygous loss of SMN1 causes spinal muscular atrophy (SMA). SMN2 is unable to prevent the disease due to skipping of exon 7. Using a systematic approach of in vivo selection, we have previously demonstrated that a weak 5' splice site (ss) serves as the major cause of skipping of SMN2 exon 7. Here we show the inhibitory impact of RNA structure on the weak 5' ss of exon 7. We call this structure terminal stem-loop 2 (TSL2). Confirming the inhibitory nature of TSL2, point mutations that destabilize TSL2 promote exon 7 inclusion in SMN2, whereas strengthening of TSL2 promotes exon 7 skipping even in SMN1. We also demonstrate that TSL2 negatively affects the recruitment of U1snRNP at the 5' ss of exon 7. Using enzymatic structure probing, we confirm that the sequence at the junction of exon 7/intron 7 folds into TSL2 and show that mutations in TSL2 cause predicted structural changes in this region. Our findings reveal for the first time the critical role of RNA structure in regulation of alternative splicing of human SMN.

Single nucleotide RNA choreography

New structural analysis methods, and a tree formalism re-define and expand the RNA motif concept, unifying what previously appeared to be disparate groups of structures. We find RNA tetraloops at high frequencies, in new contexts, with unexpected lengths, and in novel topologies. The results, with broad implications for RNA structure in general, show that even at this most elementary level of organization, RNA tolerates astounding variation in conformation, length, sequence and context. However the variation is not random; it is well-described by four distinct modes, which are 3-2 switches (backbone topology variations), insertions, deletions and strand clips.

Secondary conformation of short lysine- and leucine-rich peptides assessed by optical spectroscopies: effect of chain length, concentration, solvent, and time

Solution secondary structures of three synthetic cationic peptides, currently used in antisense oligonucleotide delivery into living cells, have been analyzed by means of circular dichroism (CD) and Raman scattering in different buffers as a function of concentration and time. All three peptides are of minimalist conception, i.e., formed by only two types of amino acids (leucine: L and lysine: K). Two of these peptides contain 15 aminoacids: N(ter)- KLLKLLLKLLLKLLK (L(10)K(5)), N(ter)-KLKLKLKLKLKLKLK (L(7)K(8)), and the third one has only 9 residues: N(ter)-KLKLKLKLK (L(4)K(5)). The conformational behavior of the 15-mers in pure water differs considerably one from another. Although both of them are initially disordered in the 50-350 microM range, L(10)K(5) gradually undergoes a disordered to alpha-helix transition for molecular concentrations above 100 microM. In all other solvents used, L(10)K(5) adopts a stable alpha-helical conformation. In methanol and methanol/Tris mixture, nonnative alpha-helices can be induced in both KL-alternating peptides, i.e., L(7)K(8) and L(4)K(5). However, in major cases and with a time delay depending on peptide concentration, beta-like structures can be gradually formed in both solutions. In PBS and methanol/PBS mixture, the tendency for L(7)K(8) and L(4)K(5) is to form structures belonging to beta-family. A discussion has been undertaken on the effect of counterions as well as their nature in the stabilization of ordered structures in both KL-alternating peptides.

Thermal stability, structural features, and B-to-Z transition in DNA tetraloop hairpins as determined by optical spectroscopy in d(CG)(3)T(4)(CG)(3) and d(CG)(3)A(4)(CG)(3) oligodeoxynucleotides

NMR and CD data have previously shown the formation of the T(4) tetraloop hairpin in aqueous solutions, as well as the possibility of the B-to-Z transition in its stem in high salt concentration conditions. It has been shown that the stem B-to-Z transition in T(4) hairpins leads to S (south)- to N (north)-type conformational changes in the loop sugars, as well as anti to syn orientations in the loop bases. In this article, we have compared by means of UV absorption, CD, Raman, and Fourier transform infrared (FTIR), the thermodynamic and structural properties of the T(4) and A(4) tetraloop hairpins formed in 5'-d(CGCGCG-TTTT-CGCGCG)-3' and 5'-d(CGCGCG-AAAA-CGCGCG)-3', respectively. In presence of 5M NaClO(4), a complete B-to-Z transition of the stems is first proved by CD spectra. UV melting profiles are consistent with a higher thermal stability of the T(4) hairpin compared to the A(4) hairpin. Order-to-disorder transition of both hairpins has also been analyzed by means of Raman spectra recorded as a function of temperature. A clear Z-to-B transition of the stem has been confirmed in the T(4) hairpin, and not in the A(4) hairpin. With a right-handed stem, Raman and FTIR spectra have confirmed the C2'-endo/anti conformation for all the T(4) loop nucleosides. With a left-handed stem, a part of the T(4) loop sugars adopt a N-type (C3'-endo) conformation, and the C3'-endo/syn conformation seems to be the preferred one for the dA residues involved in the A(4) tetraloop.

NMR structure of the 3' stem-loop from human U4 snRNA

The NMR structure of the 3' stem-loop (3'SL) from human U4 snRNA was determined to gain insight into the structural basis for conservation of this stem-loop sequence from vertebrates. 3'SL sequences from human, rat, mouse and chicken U4 snRNA each consist of a 7 bp stem capped by a UACG tetraloop. No high resolution structure has previously been reported for a UACG tetraloop. The UACG tetraloop portion of the 3'SL was especially well defined by the NMR data, with a total of 92 NOE-derived restraints (about 15 per residue), including 48 inter-residue restraints (about 8 per residue) for the tetraloop and closing C-G base pair. Distance restraints were derived from NOESY spectra using MARDIGRAS with random error analysis. Refinement of the 20mer RNA hairpin structure was carried out using the programs DYANA and miniCarlo. In the UACG tetraloop, U and G formed a base pair stabilized by two hydrogen bonds, one between the 2'-hydroxyl proton of U and carbonyl oxygen of G, another between the imino proton of G and carbonyl oxygen O2 of U. In addition, the amino group of C formed a hydrogen bond with the phosphate oxygen of A. G adopted a syn orientation about the glycosidic bond, while the sugar puckers of A and C were either C2'-endo or flexible. The conformation of the UACG tetraloop was, overall, similar to that previously reported for UUCG tetraloops, another member of the UNCG class of tetraloops. The presence of an A, rather than a U, at the variable position, however, presents a distinct surface for interaction of the 3'SL tetraloop with either RNA or protein residues that may stabilize interactions important for active spliceosome formation. Such tertiary interactions may explain the conservation of the UACG tetraloop motif in 3'SL sequences from U4 snRNA in vertebrates.