title: Sampled ensemble neutrality as a feature to classify potential structured RNAs.
abstract: BACKGROUND
  Structured RNAs have many biological functions ranging from catalysis of chemical reactions to gene regulation. Yet, many homologous structured RNAs display most of their conservation at the secondary or tertiary structure level. As a result, strategies for structured RNA discovery rely heavily on identification of sequences sharing a common stable secondary structure. However, correctly distinguishing structured RNAs from surrounding genomic sequence remains challenging, especially during de novo discovery. RNA also has a long history as a computational model for evolution due to the direct link between genotype (sequence) and phenotype (structure). From these studies it is clear that evolved RNA structures, like protein structures, can be considered robust to point mutations. In this context, an RNA sequence is considered robust if its neutrality (extent to which single mutant neighbors maintain the same secondary structure) is greater than that expected for an artificial sequence with the same minimum free energy structure.
  
  
  RESULTS
  In this work, we bring concepts from evolutionary biology to bear on the structured RNA de novo discovery process. We hypothesize that alignments corresponding to structured RNAs should consist of neutral sequences. We evaluate several measures of neutrality for their ability to distinguish between alignments of structured RNA sequences drawn from Rfam and various decoy alignments. We also introduce a new measure of RNA structural neutrality, the structure ensemble neutrality (SEN). SEN seeks to increase the biological relevance of existing neutrality measures in two ways. First, it uses information from an alignment of homologous sequences to identify a conserved biologically relevant structure for comparison. Second, it only counts base-pairs of the original structure that are absent in the comparison structure and does not penalize the formation of additional base-pairs.
  
  
  CONCLUSIONS
  We find that several measures of neutrality are effective at separating structured RNAs from decoy sequences, including both shuffled alignments and flanking genomic sequence. Furthermore, as an independent feature classifier to identify structured RNAs, SEN yields comparable performance to current approaches that consider a variety of features including stability and sequence identity. Finally, SEN outperforms other measures of neutrality at detecting mutational robustness in bacterial regulatory RNA structures.
authors: Boston College, 140 Commonwealth Ave., Chestnut Hill, 02467, MA, USA. shermin.pei@bc.edu.
authors: Shermin
authors: S
authors: Pei
authors: Boston College, 140 Commonwealth Ave., Chestnut Hill, 02467, MA, USA. jon.anthony@bc.edu.
authors: Jon S
authors: JS
authors: Anthony
authors: Boston College, 140 Commonwealth Ave., Chestnut Hill, 02467, MA, USA. m.meyer@bc.edu.
authors: Michelle M
authors: MM
authors: Meyer
mesh_terms: D001419:Bacteria / Q000235:genetics; D001483:Base Sequence; D008962:Models, Theoretical; D009154:Mutation; D009690:Nucleic Acid Conformation; D012313:RNA / Q000737:chemistry / Q000145:classification / Q000235:genetics*; D012329:RNA, Bacterial / Q000737:chemistry / Q000235:genetics*; D016415:Sequence Alignment; D017423:Sequence Analysis, RNA; D012689:Sequence Homology, Nucleic Acid; D013329:Structure-Activity Relationship
chemical_list: D012329:RNA, Bacterial; D012313:RNA
references: RNA. 1996 May;2(5):452-62
references: 8665412
references: Cell. 2014 Mar 27;157(1):77-94
references: 24679528
references: Biopolymers. 1999 Feb;49(2):145-65
references: 10070264
references: Nucleic Acids Res. 1999 Apr 1;27(7):1578-84
references: 10075987
references: Proc Natl Acad Sci U S A. 1999 Aug 17;96(17):9716-20
references: 10449760
references: Proc Natl Acad Sci U S A. 2005 Feb 15;102(7):2454-9
references: 15665081
references: RNA. 2005 May;11(5):578-91
references: 15840812
references: Bioinformatics. 2006 Feb 15;22(4):445-52
references: 16357030
references: Proc Natl Acad Sci U S A. 2006 Apr 25;103(17):6593-8
references: 16608911
references: BMC Bioinformatics. 2006;7:221
references: 16638137
references: Nat Rev Microbiol. 2006 Oct;4(10):729-40
references: 16980936
references: PLoS Genet. 2010 Aug;6(8):e1001074
references: 20808897
references: J Bioinform Comput Biol. 2010 Dec;8(6):1013-26
references: 21121024
references: BMC Bioinformatics. 2006;7:404
references: 16956394
references: Mol Biol Evol. 2006 Nov;23(11):2123-30
references: 16901984
references: Nucleic Acids Res. 2007 Jul;35(Web Server issue):W314-9
references: 17567615
references: PLoS Comput Biol. 2007 Aug;3(8):e123
references: 17784778
references: Brief Bioinform. 2007 Nov;8(6):396-414
references: 17483123
references: BMC Bioinformatics. 2008;9:33
references: 18215258
references: BMC Bioinformatics. 2008;9:122
references: 18302738
references: BMC Bioinformatics. 2008;9:222
references: 18445289
references: BMC Bioinformatics. 2008;9:248
references: 18505553
references: PLoS Comput Biol. 2008;4(8):e1000124
references: 18688270
references: BMC Bioinformatics. 2008;9:474
references: 19014431
references: Mol Biol Evol. 2009 Apr;26(4):867-74
references: 19168567
references: PLoS Comput Biol. 2009 Apr;5(4):e1000338
references: 19343219
references: RNA. 2009 Dec;15(12):2075-82
references: 19861422
references: Trends Biotechnol. 2010 Jan;28(1):9-19
references: 19942311
references: Genetics. 2000 Feb;154(2):909-21
references: 10655240
references: Bioinformatics. 2000 Jul;16(7):583-605
references: 11038329
references: Nat Rev Genet. 2001 Dec;2(12):919-29
references: 11733745
references: BMC Bioinformatics. 2001;2:8
references: 11801179
references: Trends Genet. 2004 Jan;20(1):44-50
references: 14698618
references: J Mol Evol. 2004 Jun;58(6):681-91
references: 15461425
references: J Mol Evol. 1975 Oct 3;6(1):61-76
references: 1185796
references: Nucleic Acids Res. 1992 Nov 11;20(21):5785-95
references: 1454539
references: Microbiol Rev. 1994 Mar;58(1):10-26
references: 8177168
references: Proc Biol Sci. 1994 Mar 22;255(1344):279-84
references: 7517565
references: Proc Natl Acad Sci U S A. 1996 Jan 9;93(1):397-401
references: 8552647
references: Biochemistry (Mosc). 2011 Dec;76(13):1450-64
references: 22339598
references: Bioinformatics. 2012 Apr 15;28(8):1093-101
references: 22373787
references: BMC Evol Biol. 2012;12:52
references: 22500888
references: BMC Genomics. 2012;13 Suppl 4:S6
references: 22759654
references: Bioinformatics. 2012 Nov 15;28(22):3006-8
references: 22976082
references: Nucleic Acids Res. 2013 Jan;41(Database issue):D226-32
references: 23125362
references: Nucleic Acids Res. 2013 Jan 7;41(1):44-53
references: 23125360
references: Genome Biol Evol. 2013;5(1):181-6
references: 23275486
references: Hum Mutat. 2013 Apr;34(4):546-56
references: 23315997
references: Nucleic Acids Res. 1997 Sep 15;25(18):3724-32
references: 9278497