Stabilities of consecutive A.C, C.C, G.G, U.C, and U.U mismatches in RNA internal loops: Evidence for stable hydrogen-bonded U.U and C.C.+ pairs

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

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

Thermodynamic characterization of naturally occurring RNA pentaloops

RNA folding is hierarchical; therefore, predicting RNA secondary structure from sequence is an intermediate step in predicting tertiary structure. Secondary structure prediction is based on a nearest neighbor model using free energy minimization. To improve secondary structure prediction, all types of naturally occurring secondary structure motifs need to be thermodynamically characterized. However, not all secondary structure motifs are well characterized. Pentaloops, the second most abundant hairpin size, is one such uncharacterized motif. In fact, the current thermodynamic model used to predict the stability of pentaloops was derived from a small data set of pentaloops and from data for other hairpins of different sizes. Here, the most commonly occurring pentaloops were identified and optically melted. New experimental data for 22 pentaloop sequences were combined with previously published data for nine pentaloop sequences. Using linear regression, a pentaloop-specific model was derived. This new model is simpler and more accurate than the current model. The new experimental data and improved model can be incorporated into software that is used to predict RNA secondary structure from sequence.

An NMR-based approach reveals the core structure of the functional domain of SINEUP lncRNAs

Long non-coding RNAs (lncRNAs) are attracting widespread attention for their emerging regulatory, transcriptional, epigenetic, structural and various other functions. Comprehensive transcriptome analysis has revealed that retrotransposon elements (REs) are transcribed and enriched in lncRNA sequences. However, the functions of lncRNAs and the molecular roles of the embedded REs are largely unknown. The secondary and tertiary structures of lncRNAs and their embedded REs are likely to have essential functional roles, but experimental determination and reliable computational prediction of large RNA structures have been extremely challenging. We report here the nuclear magnetic resonance (NMR)-based secondary structure determination of the 167-nt inverted short interspersed nuclear element (SINE) B2, which is embedded in antisense Uchl1 lncRNA and upregulates the translation of sense Uchl1 mRNAs. By using NMR 'fingerprints' as a sensitive probe in the domain survey, we successfully divided the full-length inverted SINE B2 into minimal units made of two discrete structured domains and one dynamic domain without altering their original structures after careful boundary adjustments. This approach allowed us to identify a structured domain in nucleotides 31-119 of the inverted SINE B2. This approach will be applicable to determining the structures of other regulatory lncRNAs.

Stability prediction of canonical and non-canonical structures of nucleic acids in various molecular environments and cells

This review provides the biophysicochemical background and recent advances in stability prediction of canonical and non-canonical structures of nucleic acids in various molecular environments and cells.

Divalent cations can control a switch-like behavior in heterotypic and homotypic RNA coacervates

Liquid-liquid phase separation (LLPS) of RNA-protein complexes plays a major role in the cellular function of membraneless organelles (MLOs). MLOs are sensitive to changes in cellular conditions, such as fluctuations in cytoplasmic ion concentrations. To investigate the effect of these changes on MLOs, we studied the influence of divalent cations on the physical and chemical properties of RNA coacervates. Using a model system comprised of an arginine-rich peptide and RNA, we predicted and observed that variations in signaling cations exert interaction-dependent effects on RNA LLPS. Changing the ionic environment has opposing effects on the propensity for heterotypic peptide-RNA and homotypic RNA LLPS, which results in a switch between coacervate types. Furthermore, divalent ion variations continuously tune the microenvironments and fluid properties of heterotypic and homotypic droplets. Our results may provide a general mechanism for modulating the biochemical environment of RNA coacervates in a cellular context.

Modeling Structure, Stability, and Flexibility of Double-Stranded RNAs in Salt Solutions

Double-stranded (ds) RNAs play essential roles in many processes of cell metabolism. The knowledge of three-dimensional (3D) structure, stability, and flexibility of dsRNAs in salt solutions is important for understanding their biological functions. In this work, we further developed our previously proposed coarse-grained model to predict 3D structure, stability, and flexibility for dsRNAs in monovalent and divalent ion solutions through involving an implicit structure-based electrostatic potential. The model can make reliable predictions for 3D structures of extensive dsRNAs with/without bulge/internal loops from their sequences, and the involvement of the structure-based electrostatic potential and corresponding ion condition can improve the predictions for 3D structures of dsRNAs in ion solutions. Furthermore, the model can make good predictions for thermal stability for extensive dsRNAs over the wide range of monovalent/divalent ion concentrations, and our analyses show that the thermally unfolding pathway of dsRNA is generally dependent on its length as well as its sequence. In addition, the model was employed to examine the salt-dependent flexibility of a dsRNA helix, and the calculated salt-dependent persistence lengths are in good accordance with experiments.

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.

Structure and Dynamics of RNA Repeat Expansions That Cause Huntington's Disease and Myotonic Dystrophy Type 1

RNA repeat expansions cause a host of incurable, genetically defined diseases. The most common class of RNA repeats consists of trinucleotide repeats. These long, repeating transcripts fold into hairpins containing 1 × 1 internal loops that can mediate disease via a variety of mechanism(s) in which RNA is the central player. Two of these disorders are Huntington's disease and myotonic dystrophy type 1, which are caused by r(CAG) and r(CUG) repeats, respectively. We report the structures of two RNA constructs containing three copies of a r(CAG) [r(3×CAG)] or r(CUG) [r(3×CUG)] motif that were modeled with nuclear magnetic resonance spectroscopy and simulated annealing with restrained molecular dynamics. The 1 × 1 internal loops of r(3×CAG) are stabilized by one-hydrogen bond (cis Watson-Crick/Watson-Crick) AA pairs, while those of r(3×CUG) prefer one- or two-hydrogen bond (cis Watson-Crick/Watson-Crick) UU pairs. Assigned chemical shifts for the residues depended on the identity of neighbors or next nearest neighbors. Additional insights into the dynamics of these RNA constructs were gained by molecular dynamics simulations and a discrete path sampling method. Results indicate that the global structures of the RNA are A-form and that the loop regions are dynamic. The results will be useful for understanding the dynamic trajectory of these RNA repeats but also may aid in the development of therapeutics.

Advancing viral RNA structure prediction: measuring the thermodynamics of pyrimidine-rich internal loops

Accurate thermodynamic parameters improve RNA structure predictions and thus accelerate understanding of RNA function and the identification of RNA drug binding sites. Many viral RNA structures, such as internal ribosome entry sites, have internal loops and bulges that are potential drug target sites. Current models used to predict internal loops are biased toward small, symmetric purine loops, and thus poorly predict asymmetric, pyrimidine-rich loops with >6 nucleotides (nt) that occur frequently in viral RNA. This article presents new thermodynamic data for 40 pyrimidine loops, many of which can form UU or protonated CC base pairs. Uracil and protonated cytosine base pairs stabilize asymmetric internal loops. Accurate prediction rules are presented that account for all thermodynamic measurements of RNA asymmetric internal loops. New loop initiation terms for loops with >6 nt are presented that do not follow previous assumptions that increasing asymmetry destabilizes loops. Since the last 2004 update, 126 new loops with asymmetry or sizes greater than 2 × 2 have been measured. These new measurements significantly deepen and diversify the thermodynamic database for RNA. These results will help better predict internal loops that are larger, pyrimidine-rich, and occur within viral structures such as internal ribosome entry sites.

Pseudouridine Modification Inhibits Muscleblind-like 1 (MBNL1) Binding to CCUG Repeats and Minimally Structured RNA through Reduced RNA Flexibility

Myotonic dystrophy type 2 is a genetic neuromuscular disease caused by the expression of expanded CCUG repeat RNAs from the non-coding region of the CCHC-type zinc finger nucleic acid-binding protein (CNBP) gene. These CCUG repeats bind and sequester a family of RNA-binding proteins known as Muscleblind-like 1, 2, and 3 (MBNL1, MBNL2, and MBNL3), and sequestration plays a significant role in pathogenicity. MBNL proteins are alternative splicing regulators that bind to the consensus RNA sequence YGCY (Y = pyrimidine). This consensus sequence is found in the toxic RNAs (CCUG repeats) and in cellular RNA substrates that MBNL proteins have been shown to bind. Replacing the uridine in CCUG repeats with pseudouridine (Ψ) resulted in a modest reduction of MBNL1 binding. Interestingly, Ψ modification of a minimally structured RNA containing YGCY motifs resulted in more robust inhibition of MBNL1 binding. The different levels of inhibition between CCUG repeat and minimally structured RNA binding appear to be due to the ability to modify both pyrimidines in the YGCY motif, which is not possible in the CCUG repeats. Molecular dynamic studies of unmodified and pseudouridylated minimally structured RNAs suggest that reducing the flexibility of the minimally structured RNA leads to reduced binding by MBNL1.

Noncanonical registers and base pairs in human 5' splice-site selection

Accurate recognition of splice sites is essential for pre-messenger RNA splicing. Mammalian 5' splice sites are mainly recognized by canonical base-pairing to the 5' end of U1 small nuclear RNA, yet we described multiple noncanonical base-pairing registers by shifting base-pair positions or allowing one-nucleotide bulges. By systematic mutational and suppressor U1 analyses, we prove three registers involving asymmetric loops and show that two-nucleotide bulges but not longer can form in this context. Importantly, we established that a noncanonical uridine-pseudouridine interaction in the 5' splice site/U1 helix contributes to the recognition of certain 5' splice sites. Thermal melting experiments support the formation of noncanonical registers and uridine-pseudouridine interactions. Overall, we experimentally validated or discarded the majority of predicted noncanonical registers, to derive a list of 5' splice sites using such alternative mechanisms that is much different from the original. This study allows not only the mechanistic understanding of the recognition of a wide diversity of mammalian 5' splice sites, but also the future development of better splice-site scoring methods that reliably predict the effects of disease-causing mutations at these sequences.

Watson-Crick-like pairs in CCUG repeats: evidence for tautomeric shifts or protonation

RNA transcripts that include expanded CCUG repeats are associated with myotonic dystrophy type 2. Crystal structures of two CCUG-containing oligomers show that the RNA strands associate into slipped duplexes that contain noncanonical C-U pairs that have apparently undergone tautomeric transition or protonation resulting in an unusual Watson-Crick-like pairing. The overhanging ends of the duplexes interact forming U-U pairs, which also show tautomerism. Duplexes consisting of CCUG repeats are thermodynamically less stable than the trinucleotide repeats involved in the TRED genetic disorders, but introducing LNA residues increases their stability and raises the melting temperature of the studied oligomers by ∼10°C, allowing detailed crystallographic studies. Quantum mechanical calculations were performed to test the possibility of the tautomeric transitions or protonation within the noncanonical pairs. The results indicate that tautomeric or ionic shifts of nucleobases can manifest themselves in biological systems, supplementing the canonical "rules of engagement."

A Tandem Oligonucleotide Approach for SNP-Selective RNA Degradation Using Modified Antisense Oligonucleotides

Antisense oligonucleotides have been studied for many years as a tool for gene silencing. One of the most difficult cases of selective RNA silencing involves the alleles of single nucleotide polymorphisms, in which the allele sequence is differentiated by a single nucleotide. A new approach to improve the performance of allele selectivity for antisense oligonucleotides is proposed. It is based on the simultaneous application of two oligonucleotides. One is complementary to the mutated form of the targeted RNA and is able to activate RNase H to cleave the RNA. The other oligonucleotide, which is complementary to the wild type allele of the targeted RNA, is able to inhibit RNase H cleavage. Five types of SNPs, C/G, G/C, G/A, A/G, and C/U, were analyzed within the sequence context of genes associated with neurodegenerative disorders such as Alzheimer’s disease, Parkinson’s disease, ALS (Amyotrophic Lateral Sclerosis), and Machado-Joseph disease. For most analyzed cases, the application of the tandem approach increased allele-selective RNA degradation 1.5–15 fold relative to the use of a single antisense oligonucleotide. The presented study proves that differentiation between single substitution is highly dependent on the nature of the SNP and surrounding nucleotides. These variables are crucial for determining the proper length of the inhibitor antisense oligonucleotide. In the tandem approach, the comparison of thermodynamic stability of the favorable duplexes WT RNA-inhibitor and Mut RNA-gapmer with the other possible duplexes allows for the evaluation of chances for the allele-selective degradation of RNA. A larger difference in thermodynamic stability between favorable duplexes and those that could possibly form, usually results in the better allele selectivity of RNA degradation.

A U⋅U Pair-to-U⋅C Pair Mutation-Induced RNA Native Structure Destabilisation and Stretching-Force-Induced RNA Misfolding

Little is known about how a non-Watson-Crick pair affects the RNA folding dynamics. We studied the effects of a U⋅U-to-U⋅C pair mutation on the folding of a hairpin in human telomerase RNA. The ensemble thermal melting of the hairpins shows an on-pathway intermediate with the disruption of the internal loop structure containing the U⋅U/U⋅C pairs. By using optical tweezers, we applied a stretching force on the terminal ends of the hairpins to probe directly the non-nearest-neighbour effects upon the mutations. The single U⋅U to U⋅C mutations are observed to 1) lower the mechanical unfolding force by approximately 1 picoNewton (pN) per mutation without affecting the unfolding reaction transition-state position (thus suggesting that removing a single hydrogen bond affects the structural dynamics at least two base pairs away), 2) result in more frequent misfolding into a small hairpin at approximately 10 pN and 3) shift the folding reaction transition-state position towards the native hairpin structure and slightly increase the mechanical folding kinetics (thus suggesting that untrapping from the misfolded state is not the rate-limiting step).

Structural insights of non-canonical U*U pair and Hoogsteen interaction probed with Se atom

Unlike DNA, in addition to the 2′-OH group, uracil nucleobase and its modifications play essential roles in structure and function diversities of non-coding RNAs. Non-canonical U•U base pair is ubiquitous in non-coding RNAs, which are highly diversified. However, it is not completely clear how uracil plays the diversifing roles. To investigate and compare the uracil in U-A and U•U base pairs, we have decided to probe them with a selenium atom by synthesizing the novel 4-Se-uridine (^SeU) phosphoramidite and Se-nucleobase-modified RNAs (^SeU-RNAs), where the exo-4-oxygen of uracil is replaced by selenium. Our crystal structure studies of U-A and U•U pairs reveal that the native and Se-derivatized structures are virtually identical, and both U-A and U•U pairs can accommodate large Se atoms. Our thermostability and crystal structure studies indicate that the weakened H-bonding in U-A pair may be compensated by the base stacking, and that the stacking of the trans-Hoogsteen U•U pairs may stabilize RNA duplex and its junction. Our result confirms that the hydrogen bond (O4^…H-C5) of the Hoogsteen pair is weak. Using the Se atom probe, our Se-functionalization studies reveal more insights into the U•U interaction and U-participation in structure and function diversification of nucleic acids.

Genome-wide analyses of Epstein-Barr virus reveal conserved RNA structures and a novel stable intronic sequence RNA

BACKGROUND

Epstein-Barr virus (EBV) is a human herpesvirus implicated in cancer and autoimmune disorders. Little is known concerning the roles of RNA structure in this important human pathogen. This study provides the first comprehensive genome-wide survey of RNA and RNA structure in EBV.

RESULTS

Novel EBV RNAs and RNA structures were identified by computational modeling and RNA-Seq analyses of EBV. Scans of the genomic sequences of four EBV strains (EBV-1, EBV-2, GD1, and GD2) and of the closely related Macacine herpesvirus 4 using the RNAz program discovered 265 regions with high probability of forming conserved RNA structures. Secondary structure models are proposed for these regions based on a combination of free energy minimization and comparative sequence analysis. The analysis of RNA-Seq data uncovered the first observation of a stable intronic sequence RNA (sisRNA) in EBV. The abundance of this sisRNA rivals that of the well-known and highly expressed EBV-encoded non-coding RNAs (EBERs).

CONCLUSION

This work identifies regions of the EBV genome likely to generate functional RNAs and RNA structures, provides structural models for these regions, and discusses potential functions suggested by the modeled structures. Enhanced understanding of the EBV transcriptome will guide future experimental analyses of the discovered RNAs and RNA structures.

Thermodynamic characterization of RNA 2 × 3 nucleotide internal loops

To better elucidate RNA structure-function relationships and to improve the design of pharmaceutical agents that target specific RNA motifs, an understanding of RNA primary, secondary, and tertiary structure is necessary. The prediction of RNA secondary structure from sequence is an intermediate step in predicting RNA three-dimensional structure. RNA secondary structure is typically predicted using a nearest neighbor model based on free energy parameters. The current free energy parameters for 2 × 3 nucleotide loops are based on a 23-member data set of 2 × 3 loops and internal loops of other sizes. A database of representative RNA secondary structures was searched to identify 2 × 3 nucleotide loops that occur in nature. Seventeen of the most frequent 2 × 3 nucleotide loops in this database were studied by optical melting experiments. Fifteen of these loops melted in a two-state manner, and the associated experimental ΔG°(37,2×3) values are, on average, 0.6 and 0.7 kcal/mol different from the values predicted for these internal loops using the predictive models proposed by Lu, Turner, and Mathews [Lu, Z. J., Turner, D. H., and Mathews, D. H. (2006) Nucleic Acids Res. 34, 4912-4924] and Chen and Turner [Chen, G., and Turner, D. H. (2006) Biochemistry 45, 4025-4043], respectively. These new ΔG°(37,2×3) values can be used to update the current algorithms that predict secondary structure from sequence. To improve free energy calculations for duplexes containing 2 × 3 nucleotide loops that still do not have experimentally determined free energy contributions, an updated predictive model was derived. This new model resulted from a linear regression analysis of the data reported here combined with 31 previously studied 2 × 3 nucleotide internal loops. Most of the values for the parameters in this new predictive model are within experimental error of those of the previous models, suggesting that approximations and assumptions associated with the derivation of the previous nearest neighbor parameters were valid. The updated predictive model predicts free energies of 2 × 3 nucleotide internal loops within 0.4 kcal/mol, on average, of the experimental free energy values. Both the experimental values and the updated predictive model can be used to improve secondary structure prediction from sequence.

Folding of the hammerhead ribozyme: pyrrolo-cytosine fluorescence separates core folding from global folding and reveals a pH-dependent conformational change

The catalytic activity of the hammerhead ribozyme is limited by its ability to fold into the native tertiary structure. Analysis of folding has been hampered by a lack of assays that can independently monitor the environment of nucleobases throughout the ribozyme-substrate complex in real time. Here, we report the development and application of a new folding assay in which we use pyrrolo-cytosine (pyC) fluorescence to (1) probe active-site formation, (2) examine the ability of peripheral ribozyme domains to support native folding, (3) identify a pH-dependent conformational change within the ribozyme, and (4) explore its influence on the equilibrium between the folded and unfolded core of the hammerhead ribozyme. We conclude that the natural ribozyme folds in two distinct noncooperative steps and the pH-dependent correlation between core folding and activity is linked to formation of the G8-C3 base pair.

Crystal structures of CGG RNA repeats with implications for fragile X-associated tremor ataxia syndrome

The CGG repeats are present in the 5'-untranslated region (5'-UTR) of the fragile X mental retardation gene FMR1 and are associated with two diseases: fragile X-associated tremor ataxia syndrome (FXTAS) and fragile X syndrome (FXS). FXTAS occurs when the number of repeats is 55-200 and FXS develops when the number exceeds 200. FXTAS is an RNA-mediated disease in which the expanded CGG tracts form stable structures and sequester important RNA binding proteins. We obtained and analysed three crystal structures of double-helical CGG repeats involving unmodified and 8-Br modified guanosine residues. Despite the presence of the non-canonical base pairs, the helices retain an A-form. In the G-G pairs one guanosine is always in the syn conformation, the other is anti. There are two hydrogen bonds between the Watson-Crick edge of G(anti) and the Hoogsteen edge of G(syn): O6·N1H and N7·N2H. The G(syn)-G(anti) pair shows affinity for binding ions in the major groove. G(syn) causes local unwinding of the helix, compensated elsewhere along the duplex. CGG helical structures appear relatively stable compared with CAG and CUG tracts. This could be an important factor in the RNA's ligand binding affinity and specificity.

Thermodynamic examination of the pyrophosphate sensor helix in the thiamine pyrophosphate riboswitch

Riboswitches are functional mRNA that control gene expression. Thiamine pyrophosphate (TPP) binds to thi-box riboswitch RNA and allosterically inhibits genes that code for proteins involved in the biosynthesis and transport of thiamine. Thiamine binding to the pyrimidine sensor helix and pyrophosphate binding to the pyrophosphate sensor helix cause changes in RNA conformation that regulate gene expression. Here we examine the thermodynamic properties of the internal loop of the pyrophosphate binding domain by comparing the wild-type construct (RNA WT) with six modified 2 x 2 bulged RNA and one 2 x 2 bulged DNA. The wild-type construct retains five conserved bases of the pyrophosphate sensor domain, two of which are in the 2 x 2 bulge (C65 and G66). The RNA WT construct was among the most stable (ΔG°₃₇ = -7.7 kcal/mol) in 1 M KCl at pH 7.5. Breaking the A•G mismatch of the bulge decreases the stability of the construct ~0.5-1 kcal/mol, but does not affect magnesium binding to the RNA WT. Guanine at position 48 is important for RNA-Mg²+ interactions of the TPP-binding riboswitch at pH 7.5. In the presence of 9.5 mM magnesium at pH 5.5, the bulged RNA constructs gained an average of 1.1 kcal/mol relative to 1 M salt. Formation of a single A+•C mismatch base pair contributes about 0.5 kcal/mol at pH 5.5, whereas two tandem A+•C mismatch base pairs together contribute about 2 kcal/mol.

In silico tools for qPCR assay design and data analysis

qPCR instruments are supplied with basic software packages that enable the measurement of fluorescent changes, calculations of quantification cycle (Cq) values, the generation of standard curves and subsequent relative target nucleic acid quantity determination. However, detailed assessments of the technical parameters underlying Cq values and their translation into biological meaningful results require validation of these basic calculations through further analyses such as qPCR efficiency correction, normalization to multiple reference genes, averaging and statistical tests. Some instruments incorporate some of these features, while others offer additional tools to complement the basic running software, in many cases providing those that are described below. In this chapter, there is a detailed description of some of these programs and recommended strategies for the design of robust qPCR assays. Some of the packages available for validation of the resulting Cq data and detailed statistical analysis are described.

Computational approaches for RNA energy parameter estimation

Methods for efficient and accurate prediction of RNA structure are increasingly valuable, given the current rapid advances in understanding the diverse functions of RNA molecules in the cell. To enhance the accuracy of secondary structure predictions, we developed and refined optimization techniques for the estimation of energy parameters. We build on two previous approaches to RNA free-energy parameter estimation: (1) the Constraint Generation (CG) method, which iteratively generates constraints that enforce known structures to have energies lower than other structures for the same molecule; and (2) the Boltzmann Likelihood (BL) method, which infers a set of RNA free-energy parameters that maximize the conditional likelihood of a set of reference RNA structures. Here, we extend these approaches in two main ways: We propose (1) a max-margin extension of CG, and (2) a novel linear Gaussian Bayesian network that models feature relationships, which effectively makes use of sparse data by sharing statistical strength between parameters. We obtain significant improvements in the accuracy of RNA minimum free-energy pseudoknot-free secondary structure prediction when measured on a comprehensive set of 2518 RNA molecules with reference structures. Our parameters can be used in conjunction with software that predicts RNA secondary structures, RNA hybridization, or ensembles of structures. Our data, software, results, and parameter sets in various formats are freely available at http://www.cs.ubc.ca/labs/beta/Projects/RNA-Params.

Positional and neighboring base pair effects on the thermodynamic stability of RNA single mismatches

Many naturally occurring RNA structures contain single mismatches, many of which occur near the ends of helices. However, previous thermodynamic studies have focused their efforts on thermodynamically characterizing centrally placed single mismatches. Additionally, algorithms currently used to predict secondary structure from sequence are based on two assumptions for predicting the stability of RNA duplexes containing this motif. It has been assumed that the thermodynamic contribution of small RNA motifs is independent of both its position in the duplex and the identity of the non-nearest neighbors. Thermodynamically characterizing single mismatches three nucleotides from both the 3' and 5' ends (i.e., off-center) of an RNA duplex and comparing these results to those of the same single mismatch-nearest neighbor combination centrally located have allowed for the investigation of these effects. The thermodynamic contributions of 13 single mismatch-nearest neighbor combinations are reported, but only nine combinations are studied at all three duplex positions and are used to determine trends and patterns. In general, the 5'- and 3'-shifted single mismatches are relatively similar, on average, and more favorable in free energy than centrally placed single mismatches. However, close examination and comparison shows there are several associated idiosyncrasies with these identified general trends. These peculiarities may be due, in part, to the identities of the single mismatch, the nearest neighbors, and the non-nearest neighbors, along with the effects of the single mismatch position in the duplex. The prediction algorithm recently proposed by Davis and Znosko [Davis, A. R., and Znosko, B. M. (2008) Biochemistry 47, 10178-10187] is used to predict the thermodynamic parameters of single mismatch contribution, and those values are compared to the measured values presented here. This comparison suggests the proposed model is a good approximation but could be improved by the addition of parameters that account for positional and/or non-nearest neighbor effects. However, more data are required to improve our understanding of these effects and to accurately account for them.

Nuclear magnetic resonance structure of the prohead RNA E-loop hairpin

The Bacillus subtilis phage phi29 packaging motor requires prohead RNA for genome encapsidation. The nuclear magnetic resonance structure of the prohead RNA E-loop hairpin, r(5'AUUGAGUU), is presented and compared to predictions from MC-SYM. The prohead RNA E-loop hairpins contain sequences similar to rRNA hairpins. Comparison of predicted and experimentally determined prohead and ribosomal hairpin structures reveals that sequence similarity is a stronger determinant of hairpin structural similarity than grouping similar types of RNA. All the hairpins contain a U-turn motif but differ in the first noncanonical pair and backbone orientation. These structures provide benchmarks for further improvements in RNA structure predictions from sequence.

Driving forces for nucleic acid pK(a) shifting in an A(+).C wobble: effects of helix position, temperature, and ionic strength

Secondary structure plays critical roles in nucleic acid function. Mismatches in DNA can lead to mutation and disease, and some mismatches involve a protonated base. Among protonated mismatches, A(+).C wobble pairs form near physiological pH and have relatively minor effects on helix geometry, making them especially important in biology. Herein, we investigate effects of helix position, temperature, and ionic strength on pK(a) shifting in A(+).C wobble pairs in DNA. We observe that pK(a) shifting is favored by internal A(+).C wobbles, which have low cooperativities of folding and make large contributions to stability, and disfavored by external A(+).C wobbles, which have high folding cooperativities but make small contributions to stability. An inverse relationship between pK(a) shifting and temperature is also found, which supports a model in which protonation is enthalpically favored overall and entropically correlated with cooperativity of folding. We also observe greater pK(a) shifts as the ionic strength decreases, consistent with anticooperativity between proton binding and counterion-condensed monovalent cation. Under the most favorable temperature and ionic strength conditions tested, a pK(a) of 8.0 is observed for the A(+).C wobble pair, which represents an especially large shift ( approximately 4.5 pK(a) units) from the unperturbed pK(a) value of adenosine. This study suggests that protonated A(+).C wobble pairs exist in DNA under biologically relevant conditions, where they can drive conformational changes and affect replication and transcription fidelity.

Thermodynamic characterization of naturally occurring RNA tetraloops

Although tetraloops are one of the most frequently occurring secondary structure motifs in RNA, less than one-third of the 30 most frequently occurring RNA tetraloops have been thermodynamically characterized. Therefore, 24 stem-loop sequences containing common tetraloops were optically melted, and the thermodynamic parameters DeltaH degrees , DeltaS degrees , DeltaG degrees (37,) and T(M) for each stem-loop were determined. These new experimental values, on average, are 0.7 kcal/mol different from the values predicted for these tetraloops using the model proposed by Vecenie CJ, Morrow CV, Zyra A, Serra MJ. 2006. Biochemistry 45: 1400-1407. The data for the 24 tetraloops reported here were then combined with the data for 28 tetraloops that were published previously. A new model, independent of terminal mismatch data, was derived to predict the free energy contribution of previously unmeasured tetraloops. The average absolute difference between the measured values and the values predicted using this proposed model is 0.4 kcal/mol. This new experimental data and updated predictive model allow for more accurate calculations of the free energy of RNA stem-loops containing tetraloops and, furthermore, should allow for improved prediction of secondary structure from sequence. It was also shown that tetraloops within the sequence 5'-GCCNNNNGGC-3' are, on average, 0.6 kcal/mol more stable than the same tetraloop within the sequence 5'-GGCNNNNGCC-3'. More systemic studies are required to determine the full extent of non-nearest-neighbor effects on tetraloop stability.

Optical melting measurements of nucleic acid thermodynamics

Optical melting experiments provide measurements of thermodynamic parameters for nucleic acids. These thermodynamic parameters are widely used in RNA structure prediction programs and DNA primer design software. This review briefly summarizes the theory and underlying assumptions of the method and provides practical details for instrument calibration, experimental design, and data interpretation.

A CA(+) pair adjacent to a sheared GA or AA pair stabilizes size-symmetric RNA internal loops

RNA internal loops are often important sites for folding and function. Residues in internal loops can have pK_a values shifted close to neutral pH because of the local structural environment. A series of RNA internal loops were studied at different pH by UV absorbance versus temperature melting experiments and imino proton nuclear magnetic resonance (NMR). A stabilizing CA pair forms at pH 7 in the and nearest neighbors when the CA pair is the first noncanonical pair (loop-terminal pair) in 3 × 3 nucleotide and larger size-symmetric internal loops. These and nearest neighbors, with CA adjacent to a closing Watson−Crick pair, are further stabilized when the pH is lowered from 7 to 5.5. The results are consistent with a significantly larger fraction (from ∼20% at pH 7 to ∼90% at pH 5.5) of adenines being protonated at the N1 position to form stabilizing wobble CA⁺ pairs adjacent to a sheared GA or AA pair. The noncanonical pair adjacent to the GA pair in can either stabilize or destabilize the loop, consistent with the sequence-dependent thermodynamics of GA pairs. No significant pH-dependent stabilization is found for most of the other nearest neighbor combinations involving CA pairs (e.g., and ), which is consistent with the formation of various nonwobble pairs observed in different local sequence contexts in crystal and NMR structures. A revised free-energy model, including stabilization by wobble CA⁺ pairs, is derived for predicting stabilities of medium-size RNA internal loops.

Correlation of RNA secondary structure statistics with thermodynamic stability and applications to folding

The accurate prediction of the secondary and tertiary structure of an RNA with different folding algorithms is dependent on several factors, including the energy functions. However, an RNA higher-order structure cannot be predicted accurately from its sequence based on a limited set of energy parameters. The inter- and intramolecular forces between this RNA and other small molecules and macromolecules, in addition to other factors in the cell such as pH, ionic strength, and temperature, influence the complex dynamics associated with transition of a single stranded RNA to its secondary and tertiary structure. Since all of the factors that affect the formation of an RNAs 3D structure cannot be determined experimentally, statistically derived potential energy has been used in the prediction of protein structure. In the current work, we evaluate the statistical free energy of various secondary structure motifs, including base-pair stacks, hairpin loops, and internal loops, using their statistical frequency obtained from the comparative analysis of more than 50,000 RNA sequences stored in the RNA Comparative Analysis Database (rCAD) at the Comparative RNA Web (CRW) Site. Statistical energy was computed from the structural statistics for several datasets. While the statistical energy for a base-pair stack correlates with experimentally derived free energy values, suggesting a Boltzmann-like distribution, variation is observed between different molecules and their location on the phylogenetic tree of life. Our statistical energy values calculated for several structural elements were utilized in the Mfold RNA-folding algorithm. The combined statistical energy values for base-pair stacks, hairpins and internal loop flanks result in a significant improvement in the accuracy of secondary structure prediction; the hairpin flanks contribute the most.

Thermodynamic characterization of tandem mismatches found in naturally occurring RNA

Although all sequence symmetric tandem mismatches and some sequence asymmetric tandem mismatches have been thermodynamically characterized and a model has been proposed to predict the stability of previously unmeasured sequence asymmetric tandem mismatches [Christiansen,M.E. and Znosko,B.M. (2008) Biochemistry, 47, 4329–4336], experimental thermodynamic data for frequently occurring tandem mismatches is lacking. Since experimental data is preferred over a predictive model, the thermodynamic parameters for 25 frequently occurring tandem mismatches were determined. These new experimental values, on average, are 1.0 kcal/mol different from the values predicted for these mismatches using the previous model. The data for the sequence asymmetric tandem mismatches reported here were then combined with the data for 72 sequence asymmetric tandem mismatches that were published previously, and the parameters used to predict the thermodynamics of previously unmeasured sequence asymmetric tandem mismatches were updated. The average absolute difference between the measured values and the values predicted using these updated parameters is 0.5 kcal/mol. This updated model improves the prediction for tandem mismatches that were predicted rather poorly by the previous model. This new experimental data and updated predictive model allow for more accurate calculations of the free energy of RNA duplexes containing tandem mismatches, and, furthermore, should allow for improved prediction of secondary structure from sequence.

The structure of the human tRNALys3 anticodon bound to the HIV genome is stabilized by modified nucleosides and adjacent mismatch base pairs

Replication of human immunodeficiency virus (HIV) requires base pairing of the reverse transcriptase primer, human tRNA^Lys3, to the viral RNA. Although the major complementary base pairing occurs between the HIV primer binding sequence (PBS) and the tRNA's 3′-terminus, an important discriminatory, secondary contact occurs between the viral A-rich Loop I, 5′-adjacent to the PBS, and the modified, U-rich anticodon domain of tRNA^Lys3. The importance of individual and combined anticodon modifications to the tRNA/HIV-1 Loop I RNA's interaction was determined. The thermal stabilities of variously modified tRNA anticodon region sequences bound to the Loop I of viral sub(sero)types G and B were analyzed and the structure of one duplex containing two modified nucleosides was determined using NMR spectroscopy and restrained molecular dynamics. The modifications 2-thiouridine, s²U₃₄, and pseudouridine, Ψ₃₉, appreciably stabilized the interaction of the anticodon region with the viral subtype G and B RNAs. The structure of the duplex results in two coaxially stacked A-form RNA stems separated by two mismatched base pairs, U₁₆₂•Ψ₃₉ and G₁₆₃•A₃₈, that maintained a reasonable A-form helix diameter. The tRNA's s²U₃₄ stabilized the interaction between the A-rich HIV Loop I sequence and the U-rich anticodon, whereas the tRNA's Ψ₃₉ stabilized the adjacent mismatched pairs.

3' terminal nucleotides determine thermodynamic stabilities of mismatches at the ends of RNA helices

The thermodynamic stabilities of consecutive mismatches at the ends of RNA helices are determined by the 3' terminal nucleotides. More than 40 RNA duplexes containing terminal motifs of 3 or more nucleotides were studied by optical melting experiments. Up to three noncanonical pairs of nucleotides at the end of RNA helices provide additional thermodynamic stability. 3' nucleotides contribute more stability than 5' nucleotides, and purines contribute more stability than pyrimidines. The additional stability of a second or third 3' nucleotide stacking on a purine is the same for both dangling ends and consecutive terminal mismatches. Current predictions underestimate RNA duplex stabilities with terminal motifs by 1.4 kcal/mol on average, which is an order of magnitude in a binding constant at 37 degrees C. Accurate thermodynamic parameters for these terminal motifs will contribute to improvements in RNA secondary structure predictions, identification of microRNA targets, and design of siRNA therapeutics with fewer off-target effects.

Thermodynamic characterization of naturally occurring RNA single mismatches with G-U nearest neighbors

Due to their prevalence and roles in biological systems, single mismatches adjacent to G-U pairs are important RNA structural elements. Since there are only limited experimental values for the stability of single mismatches adjacent to G-U pairs, current algorithms using free energy minimization to predict RNA secondary structure from sequence assign predicted thermodynamic values to these types of single mismatches. Here, thermodynamic data are reported for frequently occurring single mismatches adjacent to at least one G-U pair. This experimental data can be used in place of predicted thermodynamic values in algorithms that predict secondary structure from sequence using free energy minimization. When predicting the thermodynamic contributions of previously unmeasured single mismatches, most algorithms apply the same thermodynamic penalty for an A-U pair adjacent to a single mismatch and a G-U pair adjacent to a single mismatch. A recent study, however, suggests that the penalty for a G-U pair adjacent to a tandem mismatch should be 1.2 +/- 0.1 kcal/mol, and the penalty for an A-U pair adjacent to a tandem mismatch should be 0.5 +/- 0.2 kcal/mol [Christiansen, M. E. and Znosko, B. M. (2008) Biochemistry 47, 4329-4336]. Therefore, the data reported here are combined with the existing thermodynamic dataset of single mismatches, and nearest neighbor parameters are derived for an A-U pair adjacent to a single mismatch (1.1 +/- 0.1 kcal/mol) and a G-U pair adjacent to a single mismatch (1.4 +/- 0.1 kcal/mol).

Synthesis of polynucleotide analogs containing a polyvinyl alcohol backbone

Water soluble homo-base polynucleotide analogues were synthesized in which polyvinyl alcohol and partially phosphonated polyvinyl alcohol constituted the backbones,onto which were grafted uracil or adenine via 1,3-dioxane spacers formed by acetal formation with the 1,3-diol moieties in PVA. The resulting adenine-PVA polynucleotide analogs exhibited hyperchromic effects, which was not the case for the corresponding uracil compounds. Mixtures of the adenine- and aracil PVA-phosphate polynucleotide analogs in solutions exhibited characteristic S-shaped UV-absorbance vs temperature and melting curves with melting points at approximately 40 degrees C.

Thermodynamic characterization of the complete set of sequence symmetric tandem mismatches in RNA and an improved model for predicting the free energy contribution of sequence asymmetric tandem mismatches

Because of the availability of an abundance of RNA sequence information, the ability to rapidly and accurately predict the secondary structure of RNA from sequence is becoming increasingly important. A common method for predicting RNA secondary structure from sequence is free energy minimization. Therefore, accurate free energy contributions for every RNA secondary structure motif are necessary for accurate secondary structure predictions. Tandem mismatches are prevalent in naturally occurring sequences and are biologically important. A common method for predicting the stability of a sequence asymmetric tandem mismatch relies on the stabilities of the two corresponding sequence symmetric tandem mismatches [Mathews, D. H., Sabina, J., Zuker, M., and Turner, D. H. (1999) J. Mol. Biol. 288, 911-940]. To improve the prediction of sequence asymmetric tandem mismatches, the experimental thermodynamic parameters for the 22 previously unmeasured sequence symmetric tandem mismatches are reported. These new data, however, do not improve prediction of the free energy contributions of sequence asymmetric tandem mismatches. Therefore, a new model, independent of sequence symmetric tandem mismatch free energies, is proposed. This model consists of two penalties to account for destabilizing tandem mismatches, two bonuses to account for stabilizing tandem mismatches, and two penalties to account for A-U and G-U adjacent base pairs. This model improves the prediction of asymmetric tandem mismatch free energy contributions and is likely to improve the prediction of RNA secondary structure from sequence.

Structure and thermodynamics of a conserved U2 snRNA domain from yeast and human

The spliceosome is a dynamic ribonucleoprotein complex responsible for the removal of intron sequences from pre-messenger RNA. The highly conserved 5' end of the U2 small nuclear RNA (snRNA) makes key base-pairing interactions with the intron branch point sequence and U6 snRNA. U2 stem I, a stem-loop located in the 5' region of U2, has been implicated in spliceosome assembly and may modulate the folding of the U2 and U6 snRNAs in the spliceosome active site. Here we present the NMR structures of U2 stem I from human and Saccharomyces cerevisiae. These sequences represent the two major classes of U2 stem I, distinguished by the identity of tandem wobble pairs (UU/UU in yeast and CA/GU in human) and the presence of post-transcriptional modifications (four 2'-O-methyl groups and two pseudouracils in human). The structures reveal that the UU/UU and CA/GU tandem wobble pairs are nearly isosteric. The tandem wobble pairs separate two thermodynamically distinct regions of Watson-Crick base pairs, with the modified nucleotides in human stem I conferring a significant increase in stability. We hypothesize that the separate thermodynamic stabilities of U2 stem I exist to allow the structure to transition through different folded conformations during spliceosome assembly and catalysis.

NMR structures of (rGCUGAGGCU)2 and (rGCGGAUGCU)2: probing the structural features that shape the thermodynamic stability of GA pairs

The NMR structures of [see text] and [see text] are reported. The internal loop, [see text], is about 2 kcal/mol more stable than [see text] at 37 degrees C. The duplexes assemble into similar global folds characterized by the formation of tandem sheared GA pairs. The different stabilities of the loops are accompanied by differences in the local structure of the closing GU pairs. In the [see text] internal loop, the GU pairs form canonical wobble configurations with two hydrogen bonds, whereas in [see text], the GU pairs form a single hydrogen bond involving the amino group, GH22, and the carbonyl group, UO4. This pairing is similar to the GU closing pair of the 690 hairpin loop found in E. coli 16S rRNA. The [see text] and [see text] structures reveal how the subtle interplay between stacking and hydrogen bonding determines sequence dependent conformation and thermodynamic stability. Thus, this work provides structural and thermodynamic benchmarks for theoreticians in the ongoing effort to understand the sequence dependence of RNA physicochemical properties.

tRNA's wobble decoding of the genome: 40 years of modification

The genetic code is degenerate, in that 20 amino acids are encoded by 61 triplet codes. In 1966, Francis Crick hypothesized that the cell's limited number of tRNAs decoded the genome by recognizing more than one codon. The ambiguity of that recognition resided in the third base-pair, giving rise to the Wobble Hypothesis. Post-transcriptional modifications at tRNA's wobble position 34, especially modifications of uridine 34, enable wobble to occur. The Modified Wobble Hypothesis proposed in 1991 that specific modifications of a tRNA wobble nucleoside shape the anticodon architecture in such a manner that interactions were restricted to the complementary base plus a single wobble pairing for amino acids with twofold degenerate codons. However, chemically different modifications at position 34 would expand the ability of a tRNA to read three or even four of the fourfold degenerate codons. One foundation of Crick's Wobble Hypothesis was that a near-constant geometry of canonical base-pairing be maintained in forming all three base-pairs between the tRNA anticodon and mRNA codon on the ribosome. In accepting an aminoacyl-tRNA, the ribosome requires maintenance of a specific geometry for the anticodon-codon base-pairing. However, it is the post-transcriptional modifications at tRNA wobble position 34 and purine 37, 3'-adjacent to the anticodon, that pre-structure the anticodon domain to ensure the correct codon binding. The modifications create both the architecture and the stability needed for decoding through restraints on anticodon stereochemistry and conformational space, and through selective hydrogen bonding. A physicochemical understanding of modified nucleoside contributions to the tRNA anticodon domain architecture and its decoding of the genome has advanced RNA world evolutionary theory, the principles of RNA chemistry, and the application of this knowledge to the introduction of new amino acids to proteins.

A set of nearest neighbor parameters for predicting the enthalpy change of RNA secondary structure formation

A complete set of nearest neighbor parameters to predict the enthalpy change of RNA secondary structure formation was derived. These parameters can be used with available free energy nearest neighbor parameters to extend the secondary structure prediction of RNA sequences to temperatures other than 37°C. The parameters were tested by predicting the secondary structures of sequences with known secondary structure that are from organisms with known optimal growth temperatures. Compared with the previous set of enthalpy nearest neighbor parameters, the sensitivity of base pair prediction improved from 65.2 to 68.9% at optimal growth temperatures ranging from 10 to 60°C. Base pair probabilities were predicted with a partition function and the positive predictive value of structure prediction is 90.4% when considering the base pairs in the lowest free energy structure with pairing probability of 0.99 or above. Moreover, a strong correlation is found between the predicted melting temperatures of RNA sequences and the optimal growth temperatures of the host organism. This indicates that organisms that live at higher temperatures have evolved RNA sequences with higher melting temperatures.

The AUUCU repeats responsible for spinocerebellar ataxia type 10 form unusual RNA hairpins

Spinocerebellar ataxia type 10 is an autosomal dominant disorder caused by expansion of a pentameric repeat tract (ATTCT.AGAAT)(n) in intron 9 of the gene that encodes ataxin-10. We have analyzed duplex DNA containing the repeat, the individual DNA strands, and the RNA that would be generated by transcription of the repeat. Circular dichroism and UV absorbance measurements suggest that the previously reported tendency of the repeat to unpair when supercoiled is probably related simply to GC content rather than reflecting any unusual property of the duplex. DNA containing d(ATTCT)9 forms a folded structure at relatively low temperatures, whereas the antisense strand, d(AGAAT)9, does not form a structure even at 0 degrees C. In contrast r(AUUCU)9 forms a folded structure under physiologically reasonable conditions. S1 nuclease analysis reveals a single region of hypersensitivity in the middle of the repeat tract, whereas V1 digestion is consistent with a hydrogen bonded or well stacked structure. CD spectroscopy shows that the structure is unimolecular and hydrogen bonded and has a significant amount of A-form helix. NMR spectroscopy demonstrates that these hydrogen bonds comprise an equal number of A.U and U.U base pairs. Our data thus suggest that the repeat forms an unusual RNA hairpin. Thus the ability to form an RNA hairpin seems to be a common property of those Repeat Expansion Diseases that are not recessively inherited and are caused by repeats that are transcribed but not translated.

A hairpin-like structure within an AU-rich mRNA-destabilizing element regulates trans-factor binding selectivity and mRNA decay kinetics

In mammals, rapid mRNA turnover directed by AU-rich elements (AREs) is mediated by selective association of cellular ARE-binding proteins. These trans-acting factors display overlapping RNA substrate specificities and may act to either stabilize or destabilize targeted transcripts; however, the mechanistic features of AREs that promote preferential binding of one trans-factor over another are not well understood. Here, we describe a hairpin-like structure adopted by the ARE from tumor necrosis factor alpha (TNFalpha) mRNA that modulates its affinity for selected ARE-binding proteins. In particular, association of the mRNA-destabilizing factor p37(AUF1) was strongly inhibited by adoption of the higher order ARE structure, whereas binding of the inducible heat shock protein Hsp70 was less severely compromised. By contrast, association of the mRNA-stabilizing protein HuR was only minimally affected by changes in ARE folding. Consistent with the inverse relationship between p37(AUF1) binding affinity and the stability of ARE folding, mutations that stabilized the ARE hairpin also inhibited its ability to direct rapid mRNA turnover in transfected cells. Finally, phylogenetic analyses and structural modeling indicate that TNFalpha mRNA sequences flanking the ARE are highly conserved and may stabilize the hairpin fold in vivo. Taken together, these data suggest that local higher order structures involving AREs may function as potent regulators of mRNA turnover in mammalian cells by modulating trans-factor binding selectivity.

Diversity of base-pair conformations and their occurrence in rRNA structure and RNA structural motifs

In addition to the canonical base-pairs comprising the standard Watson-Crick (C:G and U:A) and wobble U:G conformations, an analysis of the base-pair types and conformations in the rRNAs in the high-resolution crystal structures of the Thermus thermophilus 30S and Haloarcula marismortui 50S ribosomal subunits has identified a wide variety of non-canonical base-pair types and conformations. However, the existing nomenclatures do not describe all of the observed non-canonical conformations or describe them with some ambiguity. Thus, a standardized system is required to classify all of these non-canonical conformations appropriately. Here, we propose a new, simple and systematic nomenclature that unambiguously classifies base-pair conformations occurring in base-pairs, base-triples and base-quadruples that are associated with secondary and tertiary interactions. This system is based on the topological arrangement of the two bases and glycosidic bonds in a given base-pair. Base-pairs in the internal positions of regular secondary structure helices usually form with canonical base-pair groups (C:G, U:A, and U:G) and canonical conformations (C:G WC, U:A WC, and U:G Wb). In contrast, non-helical base-pairs outside of regular structure helices usually have non-canonical base-pair groups and conformations. In addition, many non-helical base-pairs are involved in RNA motifs that form a defined set of non-canonical conformations. Thus, each rare non-canonical conformation may be functionally and structurally important. Finally, the topology-based isostericity of base-pair conformations can rationalize base-pair exchanges in the evolution of RNA molecules.

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.

Closing loop base pairs in RNA loop-loop complexes: structural behavior, interaction energy and solvation analysis through molecular dynamics simulations

Nanosecond molecular dynamics using the Ewald summation method have been performed to elucidate the structural and energetic role of the closing base pair in loop-loop RNA duplexes neutralized by Mg2+ counterions in aqueous phases. Mismatches GA, CU and Watson-Crick GC base pairs have been considered for closing the loop of an RNA in complementary interaction with HIV-1 TAR. The simulations reveal that the mismatch GA base, mediated by a water molecule, leads to a complex that presents the best compromise between flexibility and energetic contributions. The mismatch CU base pair, in spite of the presence of an inserted water molecule, is too short to achieve a tight interaction at the closing-loop junction and seems to force TAR to reorganize upon binding. An energetic analysis has allowed us to quantify the strength of the interactions of the closing and the loop-loop pairs throughout the simulations. Although the water-mediated GA closing base pair presents an interaction energy similar to that found on fully geometry-optimized structure, the water-mediated CU closing base pair energy interaction reaches less than half the optimal value.

ChipCheck--a program predicting total hybridization equilibria for DNA binding to small oligonucleotide microarrays

Presented here is the program ChipCheck that allows the computation of total hybridization equilibria for hybridization experiments involving small oligonucleotide arrays. The calculation requires the free energies of binding for all pairs of probes and targets as well as total strand concentrations and probe molecule numbers. ChipCheck has been tested computationally on microarrays with up to 100 spots and 42 target strands (4200 binding equilibria). It arrives at solutions through iterations employing the multidimensional Newton method. While currently running in simulation mode only, an extension of the approach to the exhaustive analysis of chip results is being outlined and may be implemented in the future. The output displays the extent of correct and cross hybridization both graphically and numerically. In principle, calculating total hybridization equilibria allows for eliminating noise from DNA chip results and thus an improvement in sensitivity and accuracy.

Rational design of DNA sequences for nanotechnology, microarrays and molecular computers using Eulerian graphs

Nucleic acids are molecules of choice for both established and emerging nanoscale technologies. These technologies benefit from large functional densities of 'DNA processing elements' that can be readily manufactured. To achieve the desired functionality, polynucleotide sequences are currently designed by a process that involves tedious and laborious filtering of potential candidates against a series of requirements and parameters. Here, we present a complete novel methodology for the rapid rational design of large sets of DNA sequences. This method allows for the direct implementation of very complex and detailed requirements for the generated sequences, thus avoiding 'brute force' filtering. At the same time, these sequences have narrow distributions of melting temperatures. The molecular part of the design process can be done without computer assistance, using an efficient 'human engineering' approach by drawing a single blueprint graph that represents all generated sequences. Moreover, the method eliminates the necessity for extensive thermodynamic calculations. Melting temperature can be calculated only once (or not at all). In addition, the isostability of the sequences is independent of the selection of a particular set of thermodynamic parameters. Applications are presented for DNA sequence designs for microarrays, universal microarray zip sequences and electron transfer experiments.