SimulFold: simultaneously inferring RNA structures including pseudoknots, alignments, and trees using a Bayesian MCMC framework.
Computational methods for predicting evolutionarily conserved rather than thermodynamic RNA structures have recently attracted increased interest. These methods are indispensable not only for elucidating the regulatory roles of known RNA transcripts, but also for predicting RNA genes. It has been notoriously difficult to devise them to make the best use of the available data and to predict high-quality RNA structures that may also contain pseudoknots. We introduce a novel theoretical framework for co-estimating an RNA secondary structure including pseudoknots, a multiple sequence alignment, and an evolutionary tree, given several RNA input sequences. We also present an implementation of the framework in a new computer program, called SimulFold, which employs a Bayesian Markov chain Monte Carlo method to sample from the joint posterior distribution of RNA structures, alignments, and trees. We use the new framework to predict RNA structures, and comprehensively evaluate the quality of our predictions by comparing our results to those of several other programs. We also present preliminary data that show SimulFold's potential as an alignment and phylogeny prediction method. SimulFold overcomes many conceptual limitations that current RNA structure prediction methods face, introduces several new theoretical techniques, and generates high-quality predictions of conserved RNA structures that may include pseudoknots. It is thus likely to have a strong impact, both on the field of RNA structure prediction and on a wide range of data analyses.
Meyer
Irmtraud M
IM
UBC Bioinformatics Centre, University of British Columbia, Vancouver, British Columbia, Canada. irmtraud.meyer@cantab.net
Miklós
István
I
D000465:Algorithms; D001185:Artificial Intelligence; D001483:Base Sequence; D001499:Bayes Theorem; D003198:Computer Simulation; D017124:Conserved Sequence / Q000235:genetics; D019143:Evolution, Molecular; D008956:Models, Chemical; D008957:Models, Genetic; D008958:Models, Molecular; D008969:Molecular Sequence Data; D009010:Monte Carlo Method; D010363:Pattern Recognition, Automated; D012313:RNA / Q000737:chemistry* / Q000235:genetics* / Q000648:ultrastructure; D016415:Sequence Alignment / Q000379:methods*; D017423:Sequence Analysis, RNA / Q000379:methods*; D012984:Software; D013329:Structure-Activity Relationship
D012313:RNA
Bioinformatics. 2004 Jan 1;20(1):58-66
14693809
Bull Math Biol. 1999 May;61(3):437-67
17883226
Nucleic Acids Res. 2004 Sep 24;32(16):4925-36
15448187
Nucleic Acids Res. 2003 Jul 1;31(13):3497-500
12824352
J Mol Biol. 1999 May 21;288(5):911-40
10329189
J Comput Biol. 1998 Fall;5(3):517-30
9773347
BMC Bioinformatics. 2005 Mar 24;6:73
15790387
BMC Bioinformatics. 2006 Sep 04;7:400
16952317
IEEE/ACM Trans Comput Biol Bioinform. 2004 Apr-Jun;1(2):66-77
17048382
BMC Bioinformatics. 2004 Oct 26;5:166
15507142
Nucleic Acids Res. 1981 Oct 10;9(19):5109-24
6171773
J Mol Biol. 1980 Jun 5;139(4):601-25
6997498
Bioinformatics. 1999 Jun;15(6):446-54
10383470
Nucleic Acids Res. 2005 Apr 28;33(8):2433-9
15860779
J Mol Evol. 1981;17(6):368-76
7288891
BMC Bioinformatics. 2004 Sep 30;5:140
15458580
Bioinformatics. 2005 Mar;21(6):817-20
15479709
Bioinformatics. 2005 May 15;21(10):2246-53
15731207
J Mol Biol. 2002 Mar 22;317(2):191-203
11902836
Nucleic Acids Res. 1994 Jun 11;22(11):2079-88
8029015
Mol Biol Evol. 2006 Feb;23(2):254-67
16221896
BMC Bioinformatics. 2005 Apr 01;6:83
15804354
BMC Mol Biol. 2004 Aug 06;5:10
15298702
BMC Bioinformatics. 2005 Apr 08;6:92
15819989
Proc Natl Acad Sci U S A. 2000 Jun 6;97(12):6515-20
10823910
Bioinformatics. 2004 Jul 10;20(10):1591-602
14962926
J Mol Biol. 1995 Jun 30;250(1):37-51
7541471
Pac Symp Biocomput. 2002;:163-74
11928472
Nucleic Acids Res. 2004 Jul 1;32(Web Server issue):W142-5
15215367
Proc Natl Acad Sci U S A. 2003 Dec 23;100(26):15310-5
14676318
Proc IEEE Comput Soc Bioinform Conf. 2003;2:159-68
16452790
Nucleic Acids Res. 2003 Jul 1;31(13):3429-31
12824340
Nucleic Acids Res. 1981 Jan 10;9(1):133-48
6163133
Bioinformatics. 2005 May 1;21(9):1815-24
15657094
PLoS Biol. 2005 Jun;3(6):e213
15941360
Nucleic Acids Res. 2003 Jul 1;31(13):3406-15
12824337
Nucleic Acids Res. 1991 May 11;19(9):2489-94
1710358
Mol Biol Evol. 2004 Oct;21(10):1913-22
15229291
Bioinformatics. 1998;14(8):691-9
9789095
Bioinformatics. 2000 Apr;16(4):334-40
10869031
J Comput Biol. 2000;7(3-4):409-27
11108471
Nature. 1971 Apr 9;230(5293):362-7
4927725
Bioinformatics. 2003 Aug 12;19(12):1572-4
12912839
Nat New Biol. 1973 Nov 14;246(150):40-1
4519026
Genetics. 2002 Jul;161(3):1307-20
12136032
Comput Appl Biosci. 1994 Jul;10(4):389-99
7528630
Bioinformatics. 2003 Jan;19(1):108-16
12499300
Proc Natl Acad Sci U S A. 2005 Jul 26;102(30):10557-62
16000407
Biopolymers. 1993 Sep;33(9):1389-404
7691201
Proc Natl Acad Sci U S A. 1980 Nov;77(11):6309-13
6161375
BMC Bioinformatics. 2004 Aug 04;5:104
15294028
Nucleic Acids Res. 2005 Jul 1;33(Web Server issue):W650-3
15980555
J Comput Chem. 2003 Oct;24(13):1664-77
12926009
IEEE/ACM Trans Comput Biol Bioinform. 2004 Jan-Mar;1(1):53-62
17048408
J Mol Biol. 1999 Feb 5;285(5):2053-68
9925784
Biopolymers. 1999 Feb;49(2):145-65
10070264
J Mol Biol. 2002 Jun 21;319(5):1059-66
12079347
Bioinformatics. 2003;19 Suppl 1:i66-73
12855439
Nucleic Acids Res. 2003 Jul 1;31(13):3423-8
12824339