Volume 13 Supplement 4
Evaluating our ability to predict the structural disruption of RNA by SNPs
© Ritz et al; licensee BioMed Central Ltd. 2012
Published: 18 June 2012
The structure of RiboNucleic Acid (RNA) has the potential to be altered by a Single Nucleotide Polymorphism (SNP). Disease-associated SNPs mapping to non-coding regions of the genome that are transcribed into RiboNucleic Acid (RNA) can potentially affect cellular regulation (and cause disease) by altering the structure of the transcript. We performed a large-scale meta-analysis of Selective 2'-Hydroxyl Acylation analyzed by Primer Extension (SHAPE) data, which probes the structure of RNA. We found that several single point mutations exist that significantly disrupt RNA secondary structure in the five transcripts we analyzed. Thus, every RNA that is transcribed has the potential to be a “RiboSNitch;” where a SNP causes a large conformational change that alters regulatory function. Predicting the SNPs that will have the largest effect on RNA structure remains a contemporary computational challenge. We therefore benchmarked the most popular RNA structure prediction algorithms for their ability to identify mutations that maximally affect structure. We also evaluated metrics for rank ordering the extent of the structural change. Although no single algorithm/metric combination dramatically outperformed the others, small differences in AUC (Area Under the Curve) values reveal that certain approaches do provide better agreement with experiment. The experimental data we analyzed nonetheless show that multiple single point mutations exist in all RNA transcripts that significantly disrupt structure in agreement with the predictions.
RNA (Ribonucleic Acid) is a ubiquitous messenger of genetic information in the cell and plays a central role in the regulation of molecular processes [1–5]. Unlike DNA, RNA is generally single stranded and has a high propensity to fold into functionally important structures [6–10]. These structures can be significantly disrupted by mutations including Single Nucleotide Polymorphisms (SNPs) [11, 12]. Genome-Wide Association Studies (GWAS) regularly identify disease-associated SNPs in non-coding regions of the genome. Disease-associated SNPs do not necessarily directly reveal the molecular cause of the disease and require further analysis [11, 13–15].
A majority of the genome is transcribed into RNA [16, 17]; as a result a majority of genetic mutations will also be transferred to the transcriptome. From a structural perspective, we distinguish two broad classes of RNA; highly structured RNAs (e.g. the Ribosome, tRNAs, self splicing introns, RNAse P) and RNAs that potentially adopt multiple conformations (e.g. mRNAs and non-coding RNAs) [3, 4, 18]. Structured RNAs are under significant evolutionary pressure to adopt a single, functional conformation . However, mRNAs and non-coding RNAs are not necessarily evolved to adopt a single conformation but rather adopt an ensemble of conformations [20–23]. We have recently found specific disease-associated mutations that alter the ensemble partitioning of mRNA affecting gene regulation and thus cause disease . Thus, structure is likely an important functional feature even in RNAs traditionally thought of as “unstructured.”
Algorithms to evaluate the structural and functional consequences of mutations on proteins (e.g. PolyPhen and SIFT) are commonly used to assess the potential deleterious effects of mutations [25–27]. In addition, several groups are actively developing web servers to compute the potential deleterious effects of SNPs on RNA structure and function [28, 29]. The structural basis for deleterious mutations to a structured protein is rationalized through an understanding of protein folding. For example, replacing a hydrophobic residue in the hydrophobic core of a protein with a hydrophobic amino acid will likely cause the protein to misfold [26, 27]. In RNA however, the physico-chemical properties of the four-nucleotides are not as diverse as the amino acids. Furthermore, RNA does not fold through the formation of a hydrophobic core . Instead the structure is a complex network of base-pairing and stacking interactions [3, 8]. To observe a large conformational change in an RNA, the mutation must not only disrupt an existing base-pair, but also favor a completely alternative base-pairing network. The functional consequences of structure disruption depend on whether the affected region is involved in important regulatory interactions. In certain cases, small local changes in the RNA structure may have functional consequences [15, 30]. In this manuscript we are interested in identifying the mutations that globally affect RNA structure and are thus likely to have significant functional consequences.
We initially interrogate high-throughput SHAPE chemical mapping of multiple non-coding RNAs and associated single point mutations [31, 32]. We aim to determine whether single point mutations, like in proteins, can significantly alter the structure of the RNA. We then evaluate the performance of multiple RNA structure prediction algorithms to determine the optimal strategy for identifying the mutations that disrupt RNA structure. As GWAS (Genome Wide Association Studies) continue to focus more on non-coding regions of the genome, it will become increasingly important to have accurate algorithms for assessing the potential deleterious consequences of SNPs on the transcriptome.
Results and discussion
Single mutations disrupt RNA structure
The “on” and “off” conformations of the Riboswitch are present in the Boltzmann ensemble of the WT sequence (Figure 1A, green and magenta clusters, respectively). This is consistent with recent models that suggest that Adenine riboswitching is kinetically controlled at the transcriptional level . Moreover, two other conformations (cyan and red clusters, Figure 1A) are not highly populated in the WT ensemble. If we repeat the Boltzmann sampling procedure for a sequence containing the C77G mutation (Figure 1B), we see a drastic shift in the ensemble favoring the cyan and red conformations. A majority of mutations, however, are like the U39A mutation and have very little effect on the suboptimal ensemble (Figure 1C).
To experimentally validate the prediction made by sub-optimal sampling made in Figures 1A-C, we queried the SNRNASM (Single Nucleotide Resolution Nucleic Acid Structure Mapping) archive as well as the RNA Mapping Database (RMDB, http://rmdb.stanford.edu) for chemical mapping data of the Adenine Riboswitch . We found SHAPE chemical mapping data for the WT, C77G and U39A transcripts under identical solution conditions (10 mM MgCl2 and 100 mM KCl). This data provides single nucleotide resolution measurements of base-pairing in the Riboswitch . A high normalized SHAPE reactivity indicates high flexibility and thus low probability of base-pairing, while low reactivity indicates high likelihood of base-pairing [40, 41]. The data in Figure 1D therefore experimentally validates the predictions made in Figures 1A-C. We see that the C77G (red trace) is significantly different from the black (WT) and blue (U39A) traces, consistent with a large shift in the predominant structures in the ensemble. The significant increase in SHAPE reactivity in residues 32-43 and 62-68 are consistent with the hairpin structure represented by cyan cluster.
We compute the experimental Structure Disruption Coefficient (eSDC) to evaluate the effect of a SNP on the RNA structure as described in the Methods (Equation 1). This value measures the disruptive effect of a SNP on an RNA, the higher it is the greater the structural disruption. In this case it is 2.0 for C77G and 0.1 for U39A. Furthermore, we can use the multiple repeats of the experiments to evaluate the statistical significance (p-value) of these eSDC values, i.e. the probability that we would obtain the value due to noise in the data. For the C77G, the difference is statistically significant (p-value < 0.001) while for U39A it is not (p-value >0.5).
Systematic eSDC analysis of five non-coding RNAs
The results of our analyses are plotted on Figure 2A and reveal that in all cases certain mutations (e.g. U22G.A196G in FTL, U113A in the Glycine Riboswitch) significantly disrupt RNA structure. However, a majority of mutations (e.g. U39A and U32A in the Adenine and Glycine Riboswitches) have very small effects on structure. We plotted representative SHAPE data for structurally disruptive (red) and non-disrupting mutations for the Glycine Riboswitch and P4P6 intron in Figures 2B and C, respectively. To evaluate the significance of the structural disruption, we computed the “within” distribution for multiple repeats (6-fold) of the FTL UTR RNA SHAPE data and plot the resulting histogram to the right of Figure 2A. This allows us to determine the expected eSDC values due to the noise in the experimental data, and evaluate the p-value for any given eSDC. Clearly, single point mutations exist that significantly disrupt RNA structure, however a majority of mutations result in no measurable effect.
The FTL UTR data set is particularly interesting, as this RNA is a “RiboSNitch,” i.e. an RNA in which specific SNPs can alter structure and cause disease [24, 42]. In this case, FTL is associated with Hyperferritinemia Cataract Syndrome, a rare genetic disorder that is characterized by early onset cataracts due to excess ferritin in the retina [43, 44]. We indicate the disease-associated SNPs as magenta text in Figure 2A. All the disease-associated SNPs alter the structure of the RNA significantly (p-value < 0.001).
Three of the RNAs tested in Figure 2A are Riboswitches and undergo a conformational change if ligand is present. We can therefore compute an eSDC value for SHAPE traces in the presence and absence of ligand. We indicate these eSDC values with a green horizontal line in Figure 2A. The reason this result is important is that the structural change caused by ligand binding to a Riboswitch is sufficient to regulate gene expression [37, 45, 46]. Thus the Riboswitch ligand eSDC value (green line Figure 2A) represents a “biological” threshold above which the structure change is likely to affect function. A particularly important result of this analysis is the identification of multiple SNPs with much larger eSDC values compared to ligand binding in the Riboswitches. Thus, it is likely that a majority of these SNPs will have important functional consequences.
Performance of RNA structure prediction algorithms for RiboSNitch detection
RNA is a ubiquitous regulatory molecule in the cell and there is growing evidence that structure is a central component of its function [52, 53]. The Riboswitches studied in this manuscript are one of many examples where RNA structure change regulates bacterial metabolism [46, 54, 55]. In the case of the 5’ UTR, disease-associated SNPs disrupt structure and deregulate Ferritin levels in the eye, resulting in early onset cataracts . The T. thermophila group I intron (P4P6) must fold into its correct three-dimensional structure to self-catalyze its splicing reaction [8, 56]. In these examples, structure change is central to the RNA’s function in the cell.
The data we present in Figures 1 and 2 reveals the extent to which a single point mutation can disrupt RNA structure. Our systematic analysis of 470 mutations on five RNAs reveals that large scale SNP induced structure change is common in RNA and can potentially contribute to disease . Interestingly, all RNA secondary structure prediction algorithms predict that a small subset of mutations will have a large effect on secondary structure. The data we present in Figures 1 and 2 cover a relatively comprehensive set of mutations in each RNA, but are nonetheless limited to five functional molecules. As such, the generalizability of these results will require the analysis of larger experimental data sets as they become available .
The mechanism for this change is best illustrated in Figure 1, where we see how a single mutation (in this case C77G) can completely alter the thermodynamic folding landscape of the RNA, favoring an alternative conformation. The data we present in Figure 2 suggest that the thermodynamic models used to predict RNA structure are sound, as we find mutations experimentally in all RNAs studied that disrupt structure. All RNA structure prediction algorithms predict that certain mutations will significantly disrupt structure. In addition, a recent study of common SNPs in the human genome revealed that these affect local RNA structure .
An important result in our analysis of the Riboswitch SHAPE data is the comparison of the eSDC values for mutations relative to ligand induced conformational change (see green lines, Figure 2A). For all three Riboswitches, multiple mutations exist that result in far larger structural changes (as measured by our SDC metric) than ligand binding. This is highly relevant, as ligand binding induced structure change can completely turn on (or off) gene expression translationally and/or transcriptionally . Thus the mutations above the green lines in Figure 2A have even greater potential to regulate cellular function. This means any functional RNA has the potential to be a “RiboSNitch,” as there exists mutations that can significantly disrupt its structure.
The data we present in Figure 2 are ideal for benchmarking RNA structure prediction algorithms. The analysis we carried out in this manuscript is different from previous secondary structure prediction benchmarks, because we are specifically interested in identifying mutations that globally disrupt a given secondary structure. We developed metrics based on RNA secondary structure prediction algorithms analogous to our eSDC calculations. We can use such an analogy, since SHAPE data is correlated with base-pair probability. The SDC metrics are purposefully global, and we did not evaluate algorithms for their ability to predict the specific local changes in structure, but rather whether they predict that a specific mutation will disrupt structure relative to others. Our reasoning for this approach is that for the analysis of disease-associated SNPs, we are most interested in identifying the most structurally deleterious mutations.
Although RNA structure prediction algorithms correctly predicted that all RNAs are disrupted by certain mutations, it is clear that predicting exactly which mutation will alter structure remains very challenging. Although there is some variation in the relative performance of the different algorithmic and metric combinations we tested, the AUC values reported in Figure 5B remain relatively low. This result is not necessarily surprising, as none of the RNA structure prediction algorithms (other than RNAmutants) have been optimized to predict which mutations disrupt structure. In fact, an algorithm’s sensitivity to point mutations is often viewed as a weakness, favoring methods that are less sensitive to mutation. However, the experimental data clearly show that SNPs can profoundly change an RNA’s folding landscape.
The attempts to constantly refine algorithms so as to have them always converge on a single “correct” RNA structure may not improve their ability to identify RiboSNitches. Although only anecdotal, mFold’s good performance in our benchmark (AUC 0.62, Figure 5B) may indicate that simpler energy functions, which tend to predict more alternative structures, may ultimately perform better for identifying RiboSNitches. Indeed RNAStructure’s relatively low performance in our benchmark is surprising, since it has the most sophisticated and accurate energy function and is most accurate in structure prediction [48, 50]. Improvements in our ability to predict RiboSNitches will likely require a better understanding of the suboptimal ensemble and how mutations affect it in addition to improved energy functions. With the growing number of sequencing efforts revealing ever more single nucleotide variants in the non-coding regions of the genome, accurate algorithms predicting the structural consequences of these mutations are likely to play an important role in genomic interpretation.
Data collection and analysis
where p CC is the WT/mutant pearson correlation coefficient and n is the length of the RNA. The eSDC quantitatively evaluates the effect of a mutation on RNA structure. Prior to the calculation of the eSDC, normalized SHAPE values were capped at one in order to increase the metric's ability to reflect changes in structure identified by differences in the peaking pattern and not minor differences in peak intensity. Significance testing for structure disruption was adjusted using a Bonferroni correction.
PCA analysis of the ensemble of structures and clustering
Principal components were calculated (as described previously) from a total of 10,000 sampled structures generated equally from a WT sequence and mutants of interest . The principal components were generated from the binary representation of these 10,000 structures. These structures were then projected onto the first two principal components and subjected to the k-means clustering algorithm to reveal distinct clusters . The centroid structure of each cluster was identified from the k-means clustering algorithm and then drawn using R2R . Individual mutant structures were then generated (as discussed in Fig. 3) and projected onto the first two principal components. Each structure projection is colored according to their cluster.
Computation of the partition functions from sampled structures and calculation of the Z centroid
Partition functions were generated for each ensemble of structures. Each structure is first transformed to matrix form as described in . This is accomplished by creating an NxN matrix where N is the length of the sequence and placing a 1 at position i,j and j,i if nucleotides i and j are paired and a 0 if they are not paired. When all the matrices representing the structures are summed together and then divided by the total number of structures, the resulting matrix is the partition matrix. This matrix contains the probability of nucleotide i being paired to j. The Z centroid is defined as the structure with all the probability of pairing for each pair greater than 50%.
ROC analysis of prediction performance
Analogously to the eSDC calculation, we compute a pSDC (predicted Structure Disruption Coefficient) by computing the Pearson Correlation Coefficient ( pred CC) between WT and mutant for each RNA structure prediction algorithm. This value is analogous to the eSDC in that it allows us to rank order the disruptive effect of mutations on RNA.
To determine ROC values, the mutations were listed from highest to lowest according to their eSDC value. The top 50% of eSDC values were considered to disrupt the structure while the lowest 50% preserved structure. A second list was generated using the same mutants but using the pSDC values instead. A true positive was defined as having a pSDC value above a cutoff and experimentally disrupting the structure while a true negative was defined as having a pSDC below a cutoff and experimentally preserving a structure. A false positive or false negative is recorded when the predictions contrast with the experimental results. The pSDC cutoff was defined by stepping through the pSDC ranks. The resulting true positive rates and false positive rates were then used to generate an ROC curve. The area under the curve was calculated for each ROC using the trapezoidal method. This process was bootstrapped for each program/metric 5000 times using 20 randomly selected mutants from each set. Due to the fact that each of the RNA data sets has a differing number of mutants, the bootstrapping is done by sampling 20 mutants from each of the other data sets besides FTL, in order to correct for any bias that might come up due to one program/metric favoring one data set over another. This results in the ROC being run on 145 mutants at a time, not the full 470. The average area under the curve was calculated with the standard deviation between runs generating the error. The closer the area under the curve was to one the better the predictive power for a given program/metric.
Precise WT sequences, corresponding mutations (SNPs), eSDC values and normalized SHAPE data are provided as separate excel spreadsheets in the additional files. These data should facilitate further benchmarking efforts for novel algorithms to predict RNA structure change.
List of abbreviations
experimental Structure Disruption Coefficientl
predicted Structure Disruption Coefficient
Selective 2'-Hydroxyl Acylation analyzed by Primer Extension
Single Nucleotide Polymorphism
False Positive Rate
True Positive Rate
Single Nucleotide Resolution Nucleic Acid Structure Mapping.
We would like to thank Rhiju Das and Pablo Cordero (Stanford University) for providing SHAPE data sets in the SNRNASM format through their database (http://rmdb.stanford.edu/). This work is supported by grants R00 GM079953 (NIGMS), R21 MH087336 (NIMH), R01 HL111527 and R01 GM101237 to A.L.
This article has been published as part of BMC Genomics Volume 13 Supplement 4, 2012: SNP-SIG 2011: Identification and annotation of SNPs in the context of structure, function and disease. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/13/S4.
- Preall JB, Sontheimer EJ: RNAi: RISC gets loaded. Cell. 2005, 123 (4): 543-545. 10.1016/j.cell.2005.11.006.View ArticlePubMed
- Pesole G, Liuni S, Grillo G, Ippedico M, Larizza A, Makalowski W, Saccone C: UTRdb: a specialized database of 5' and 3' untranslated regions of eukaryotic mRNAs. Nucleic Acids Res. 1999, 27 (1): 188-191. 10.1093/nar/27.1.188.PubMed CentralView ArticlePubMed
- Laederach A, Chan JM, Schwartzman A, Willgohs E, Altman RB: Coplanar and coaxial orientations of RNA bases and helices. Rna. 2007, 13 (5): 643-650. 10.1261/rna.381407.PubMed CentralView ArticlePubMed
- Laederach A: Informatics challenges in structured RNA. Briefings in bioinformatics. 2007, 8 (5): 294-303. 10.1093/bib/bbm026.PubMed CentralView ArticlePubMed
- Doherty EA, Doudna JA: Ribozyme structures and mechanisms. Annu Rev Biophys Biomol Struct. 2001, 30: 457-475. 10.1146/annurev.biophys.30.1.457.View ArticlePubMed
- Martin JS, Simmons K, Laederach A: Exhaustive Enumeration of Kinetic Model Topologies for the Analysis of Time-Resoved RNA Folding. Algorithms. 2009, 2 (1): 200-214. 10.3390/a2010200.PubMed CentralView ArticlePubMed
- Laederach A, Shcherbakova I, Liang MP, Brenowitz M, Altman RB: Local kinetic measures of macromolecular structure reveal partitioning among multiple parallel pathways from the earliest steps in the folding of a large RNA molecule. J Mol Biol. 2006, 358 (4): 1179-1190. 10.1016/j.jmb.2006.02.075.PubMed CentralView ArticlePubMed
- Laederach A, Shcherbakova I, Jonikas MA, Altman RB, Brenowitz M: Distinct contribution of electrostatics, initial conformational ensemble, and macromolecular stability in RNA folding. Proc Natl Acad Sci U S A. 2007, 104 (17): 7045-7050. 10.1073/pnas.0608765104.PubMed CentralView ArticlePubMed
- Das R, Kwok LW, Millett IS, Bai Y, Mills TT, Jacob J, Maskel GS, Seifert S, Mochrie SG, Thiyagarajan P, et al: The fastest global events in RNA folding: electrostatic relaxation and tertiary collapse of the Tetrahymena ribozyme. J Mol Biol. 2003, 332 (2): 311-319. 10.1016/S0022-2836(03)00854-4.View ArticlePubMed
- Chu VB, Herschlag D: Unwinding RNA's secrets: advances in the biology, physics, and modeling of complex RNAs. Curr Opin Struct Biol. 2008, 18 (3): 305-314. 10.1016/j.sbi.2008.05.002.PubMed CentralView ArticlePubMed
- Zwarts KY, Clee SM, Zwinderman AH, Engert JC, Singaraja R, Loubser O, James E, Roomp K, Hudson TJ, Jukema JW, et al: ABCA1 regulatory variants influence coronary artery disease independent of effects on plasma lipid levels. Clin Genet. 2002, 61 (2): 115-125. 10.1034/j.1399-0004.2002.610206.x.View ArticlePubMed
- Kimchi-Sarfaty C, Oh JM, Kim IW, Sauna ZE, Calcagno AM, Ambudkar SV, Gottesman MM: A "silent" polymorphism in the MDR1 gene changes substrate specificity. Science. 2007, 315 (5811): 525-528. 10.1126/science.1135308.View ArticlePubMed
- Showalter HD, Denny WA: A roadmap for drug discovery and its translation to small molecule agents in clinical development for tuberculosis treatment. Tuberculosis (Edinburgh, Scotland). 2008, 88 (Suppl 1): S3-17.View Article
- Kohnlein T, Welte T: Alpha-1 antitrypsin deficiency: pathogenesis, clinical presentation, diagnosis, and treatment. The American journal of medicine. 2008, 121 (1): 3-9. 10.1016/j.amjmed.2007.07.025.View ArticlePubMed
- Glinsky GV: SNP-guided microRNA maps (MirMaps) of 16 common human disorders identify a clinically accessible therapy reversing transcriptional aberrations of nuclear import and inflammasome pathways. Cell cycle (Georgetown, Tex. 2008, 7 (22): 3564-3576. 10.4161/cc.7.22.7073.View Article
- Weinstock GM: ENCODE: more genomic empowerment. Genome Res. 2007, 17 (6): 667-668. 10.1101/gr.6534207.View ArticlePubMed
- Karolchik D, Kuhn RM, Baertsch R, Barber GP, Clawson H, Diekhans M, Giardine B, Harte RA, Hinrichs AS, Hsu F, et al: The UCSC Genome Browser Database: 2008 update. Nucleic Acids Res. 2008, 36 (Database issue): D773-779.PubMed CentralPubMed
- Reeder J, Hochsmann M, Rehmsmeier M, Voss B, Giegerich R: Beyond Mfold: recent advances in RNA bioinformatics. J Biotechnol. 2006, 124 (1): 41-55. 10.1016/j.jbiotec.2006.01.034.View ArticlePubMed
- Kertesz M, Wan Y, Mazor E, Rinn JL, Nutter RC, Chang HY, Segal E: Genome-wide measurement of RNA secondary structure in yeast. Nature. 2010, 467 (7311): 103-107. 10.1038/nature09322.View ArticlePubMed
- Ying SY, Chang DC, Miller JD, Lin SL: The microRNA: overview of the RNA gene that modulates gene functions. Methods Mol Biol. 2006, 342: 1-18.PubMed
- Tucker BJ, Breaker RR: Riboswitches as versatile gene control elements. Curr Opin Struct Biol. 2005, 15 (3): 342-348. 10.1016/j.sbi.2005.05.003.View ArticlePubMed
- Lemay JF, Lemieux C, St-Andre O, Bachand F: Crossing the borders: poly(A)-binding proteins working on both sides of the fence. RNA Biol. 2010, 7 (3): 291-295. 10.4161/rna.7.3.11649.View ArticlePubMed
- Ding Y, Chan CY, Lawrence CE: Sfold web server for statistical folding and rational design of nucleic acids. Nucleic Acids Res. 2004, 32 (Web Server issue): W135-141.PubMed CentralView ArticlePubMed
- Halvorsen M, Martin JS, Broadaway S, Laederach A: Disease-associated mutations that alter the RNA structural ensemble. PLoS genetics. 2010, 6 (8): e1001074-10.1371/journal.pgen.1001074.PubMed CentralView ArticlePubMed
- Gorlov IP, Gorlova OY, Sunyaev SR, Spitz MR, Amos CI: Shifting paradigm of association studies: value of rare single-nucleotide polymorphisms. American journal of human genetics. 2008, 82 (1): 100-112. 10.1016/j.ajhg.2007.09.006.PubMed CentralView ArticlePubMed
- Ng PC, Henikoff S: SIFT: Predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003, 31 (13): 3812-3814. 10.1093/nar/gkg509.PubMed CentralView ArticlePubMed
- Ng PC, Henikoff S: Predicting deleterious amino acid substitutions. Genome Res. 2001, 11 (5): 863-874. 10.1101/gr.176601.PubMed CentralView ArticlePubMed
- Churkin A, Gabdank I, Barash D: The RNAmute web server for the mutational analysis of RNA secondary structures. Nucleic Acids Res. 2011, 39 (Web Server issue): W92-99.PubMed CentralView ArticlePubMed
- Barash D, Churkin A: Mutational analysis in RNAs: comparing programs for RNA deleterious mutation prediction. Briefings in bioinformatics. 2011, 12 (2): 104-114. 10.1093/bib/bbq059.View ArticlePubMed
- Glinsky GV: Disease phenocode analysis identifies SNP-guided microRNA maps (MirMaps) associated with human "master" disease genes. Cell cycle (Georgetown, Tex. 2008, 7 (23): 3680-3694. 10.4161/cc.7.23.7153.View Article
- Das R, Kudaravalli M, Jonikas M, Laederach A, Fong R, Schwans JP, Baker D, Piccirilli JA, Altman RB, Herschlag D: Structural inference of native and partially folded RNA by high-throughput contact mapping. Proc Natl Acad Sci U S A. 2008, 105 (11): 4144-4149. 10.1073/pnas.0709032105.PubMed CentralView ArticlePubMed
- Das R, Karanicolas J, Baker D: Atomic accuracy in predicting and designing noncanonical RNA structure. Nat Methods. 2010, 7 (4): 291-294. 10.1038/nmeth.1433.PubMed CentralView ArticlePubMed
- Weinberg Z, Barrick JE, Yao Z, Roth A, Kim JN, Gore J, Wang JX, Lee ER, Block KF, Sudarsan N, et al: Identification of 22 candidate structured RNAs in bacteria using the CMfinder comparative genomics pipeline. Nucleic Acids Res. 2007, 35 (14): 4809-4819. 10.1093/nar/gkm487.PubMed CentralView ArticlePubMed
- Lipfert J, Das R, Chu VB, Kudaravalli M, Boyd N, Herschlag D, Doniach S: Structural transitions and thermodynamics of a glycine-dependent riboswitch from Vibrio cholerae. J Mol Biol. 2007, 365 (5): 1393-1406. 10.1016/j.jmb.2006.10.022.PubMed CentralView ArticlePubMed
- Lemay JF, Penedo JC, Mulhbacher J, Lafontaine DA: Molecular basis of RNA-mediated gene regulation on the adenine riboswitch by single-molecule approaches. Methods Mol Biol. 2009, 540: 65-76. 10.1007/978-1-59745-558-9_6.View ArticlePubMed
- Lemay JF, Lafontaine DA: Core requirements of the adenine riboswitch aptamer for ligand binding. Rna. 2007, 13 (3): 339-350. 10.1261/rna.142007.PubMed CentralView ArticlePubMed
- Lemay JF, Penedo JC, Tremblay R, Lilley DM, Lafontaine DA: Folding of the adenine riboswitch. Chem Biol. 2006, 13 (8): 857-868. 10.1016/j.chembiol.2006.06.010.View ArticlePubMed
- Rocca-Serra P, Bellaousov S, Birmingham A, Chen C, Cordero P, Das R, Davis-Neulander L, Duncan CD, Halvorsen M, Knight R, et al: Sharing and archiving nucleic acid structure mapping data. RNA. 2011, 17 (7): 1204-1212. 10.1261/rna.2753211.PubMed CentralView ArticlePubMed
- Wilkinson KA, Merino EJ, Weeks KM: Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE): quantitative RNA structure analysis at single nucleotide resolution. Nature protocols. 2006, 1 (3): 1610-1616. 10.1038/nprot.2006.249.View ArticlePubMed
- Wilkinson KA, Vasa SM, Deigan KE, Mortimer SA, Giddings MC, Weeks KM: Influence of nucleotide identity on ribose 2'-hydroxyl reactivity in RNA. Rna. 2009
- Jones CN, Wilkinson KA, Hung KT, Weeks KM, Spremulli LL: Lack of secondary structure characterizes the 5' ends of mammalian mitochondrial mRNAs. RNA. 2008, 14 (5): 862-871. 10.1261/rna.909208.PubMed CentralView ArticlePubMed
- Yoon S, Kim J, Hum J, Kim H, Park S, Kladwang W, Das R: HiTRACE: high-throughput robust analysis for capillary electrophoresis. Bioinformatics. 2011, 27 (13): 1798-1805. 10.1093/bioinformatics/btr277.View ArticlePubMed
- Ferrari F, Foglieni B, Arosio P, Camaschella C, Daraio F, Levi S, Garcia Erce JA, Beaumont C, Cazzola M, Ferrari M, et al: Microelectronic DNA chip for hereditary hyperferritinemia cataract syndrome, a model for large-scale analysis of disorders of iron metabolism. Human mutation. 2006, 27 (2): 201-208. 10.1002/humu.20294.View ArticlePubMed
- Cremonesi L, Paroni R, Foglieni B, Galbiati S, Fermo I, Soriani N, Belloli S, Ruggeri G, Biasiotto G, Cazzola M, et al: Scanning mutations of the 5'UTR regulatory sequence of L-ferritin by denaturing high-performance liquid chromatography: identification of new mutations. British journal of haematology. 2003, 121 (1): 173-179. 10.1046/j.1365-2141.2003.04253.x.View ArticlePubMed
- Lemay JF, Desnoyers G, Blouin S, Heppell B, Bastet L, St-Pierre P, Masse E, Lafontaine DA: Comparative study between transcriptionally- and translationally-acting adenine riboswitches reveals key differences in riboswitch regulatory mechanisms. PLoS genetics. 2011, 7 (1): e1001278-10.1371/journal.pgen.1001278.PubMed CentralView ArticlePubMed
- Grundy FJ, Henkin TM: From ribosome to riboswitch: control of gene expression in bacteria by RNA structural rearrangements. Crit Rev Biochem Mol Biol. 2006, 41 (6): 329-338. 10.1080/10409230600914294.View ArticlePubMed
- Waldispuhl J, Clote P: Computing the partition function and sampling for saturated secondary structures of RNA, with respect to the Turner energy model. J Comput Biol. 2007, 14 (2): 190-215. 10.1089/cmb.2006.0012.View ArticlePubMed
- Mathews DH: Using an RNA secondary structure partition function to determine confidence in base pairs predicted by free energy minimization. RNA. 2004, 10 (8): 1178-1190. 10.1261/rna.7650904.PubMed CentralView ArticlePubMed
- Hofacker IL, Stadler PF: Memory efficient folding algorithms for circular RNA secondary structures. Bioinformatics. 2006, 22 (10): 1172-1176. 10.1093/bioinformatics/btl023.View ArticlePubMed
- Wilkinson KA, Gorelick RJ, Vasa SM, Guex N, Rein A, Mathews DH, Giddings MC, Weeks KM: High-throughput SHAPE analysis reveals structures in HIV-1 genomic RNA strongly conserved across distinct biological states. PLoS biology. 2008, 6 (4): e96-10.1371/journal.pbio.0060096.PubMed CentralView ArticlePubMed
- Ding Y, Chan CY, Lawrence CE: RNA secondary structure prediction by centroids in a Boltzmann weighted ensemble. Rna. 2005, 11 (8): 1157-1166. 10.1261/rna.2500605.PubMed CentralView ArticlePubMed
- Cloonan N, Grimmond SM: Transcriptome content and dynamics at single-nucleotide resolution. Genome biology. 2008, 9 (9): 234-10.1186/gb-2008-9-9-234.PubMed CentralView ArticlePubMed
- Ashraf SI, Kunes S: A trace of silence: memory and microRNA at the synapse. Curr Opin Neurobiol. 2006, 16 (5): 535-539. 10.1016/j.conb.2006.08.007.View ArticlePubMed
- Lu C, Smith AM, Fuchs RT, Ding F, Rajashankar K, Henkin TM, Ke A: Crystal structures of the SAM-III/S(MK) riboswitch reveal the SAM-dependent translation inhibition mechanism. Nat Struct Mol Biol. 2008, 15 (10): 1076-1083. 10.1038/nsmb.1494.PubMed CentralView ArticlePubMed
- Gilbert SD, Rambo RP, Van Tyne D, Batey RT: Structure of the SAM-II riboswitch bound to S-adenosylmethionine. Nat Struct Mol Biol. 2008, 15 (2): 177-182. 10.1038/nsmb.1371.View ArticlePubMed
- Laederach A, Shcherbakova I, Liang M, Brenowitz M, Altman RB: Local kinetic measures of macromolecular structure reveal partitioning among multiple parallel pathways from the earliest steps in the folding of a large RNA molecule. J Mol Biol. 2006, 358 (358): 1179-1190.PubMed CentralView ArticlePubMed
- Sadee W, Wang D, Papp AC, Pinsonneault JK, Smith RM, Moyer RA, Johnson AD: Pharmacogenomics of the RNA world: structural RNA polymorphisms in drug therapy. Clin Pharmacol Ther. 2011, 89 (3): 355-365. 10.1038/clpt.2010.314.PubMed CentralView ArticlePubMed
- Rocca-Serra P, Bellaousov S, Birmingham A, Chen C, Cordero P, Das R, Davis-Neulander L, Duncan CD, Halvorsen M, Knight R, et al: Sharing and archiving nucleic acid structure mapping data. RNA. 2011,
- Quarrier S, Martin JS, Davis-Neulander L, Beauregard A, Laederach A: Evaluation of the information content of RNA structure mapping data for secondary structure prediction. RNA. 2010, 16 (6): 1108-1117. 10.1261/rna.1988510.PubMed CentralView ArticlePubMed
- Tibshirani RJ, Walther G, Hastie T: Estimating the number of clusters in a data set via the gap statistic. Journal of the Royal Society: Series B (Statistical Methodology). 2001, 63 (2): 411-423. 10.1111/1467-9868.00293.
- Weinberg Z, Breaker RR: R2R--software to speed the depiction of aesthetic consensus RNA secondary structures. BMC bioinformatics. 2011, 12: 3-10.1186/1471-2105-12-3.PubMed CentralView ArticlePubMed
- Sonego P, Kocsor A, Pongor S: ROC analysis: applications to the classification of biological sequences and 3D structures. Briefings in bioinformatics. 2008, 9 (3): 198-209. 10.1093/bib/bbm064.View ArticlePubMed
- Mitra S, Laederach A, Golden BL, Altman RB, Brenowitz M: RNA molecules with conserved catalytic cores but variable peripheries fold along unique energetically optimized pathways. RNA. 2011, 17 (8): 1589-1603. 10.1261/rna.2694811.PubMed CentralView ArticlePubMed
- Tafer H, Ameres SL, Obernosterer G, Gebeshuber CA, Schroeder R, Martinez J, Hofacker IL: The impact of target site accessibility on the design of effective siRNAs. Nature biotechnology. 2008, 26 (5): 578-583. 10.1038/nbt1404.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.