FR3D: finding local and composite recurrent structural motifs in RNA 3D structures.
New methods are described for finding recurrent three-dimensional (3D) motifs in RNA atomic-resolution structures. Recurrent RNA 3D motifs are sets of RNA nucleotides with similar spatial arrangements. They can be local or composite. Local motifs comprise nucleotides that occur in the same hairpin or internal loop. Composite motifs comprise nucleotides belonging to three or more different RNA strand segments or molecules. We use a base-centered approach to construct efficient, yet exhaustive search procedures using geometric, symbolic, or mixed representations of RNA structure that we implement in a suite of MATLAB programs, "Find RNA 3D" (FR3D). The first modules of FR3D preprocess structure files to classify base-pair and -stacking interactions. Each base is represented geometrically by the position of its glycosidic nitrogen in 3D space and by the rotation matrix that describes its orientation with respect to a common frame. Base-pairing and base-stacking interactions are calculated from the base geometries and are represented symbolically according to the Leontis/Westhof basepairing classification, extended to include base-stacking. These data are stored and used to organize motif searches. For geometric searches, the user supplies the 3D structure of a query motif which FR3D uses to find and score geometrically similar candidate motifs, without regard to the sequential position of their nucleotides in the RNA chain or the identity of their bases. To score and rank candidate motifs, FR3D calculates a geometric discrepancy by rigidly rotating candidates to align optimally with the query motif and then comparing the relative orientations of the corresponding bases in the query and candidate motifs. Given the growing size of the RNA structure database, it is impossible to explicitly compute the discrepancy for all conceivable candidate motifs, even for motifs with less than ten nucleotides. The screening algorithm that we describe finds all candidate motifs whose geometric discrepancy with respect to the query motif falls below a user-specified cutoff discrepancy. This technique can be applied to RMSD searches. Candidate motifs identified geometrically may be further screened symbolically to identify those that contain particular basepair types or base-stacking arrangements or that conform to sequence continuity or nucleotide identity constraints. Purely symbolic searches for motifs containing user-defined sequence, continuity and interaction constraints have also been implemented. We demonstrate that FR3D finds all occurrences, both local and composite and with nucleotide substitutions, of sarcin/ricin and kink-turn motifs in the 23S and 5S ribosomal RNA 3D structures of the H. marismortui 50S ribosomal subunit and assigns the lowest discrepancy scores to bona fide examples of these motifs. The search algorithms have been optimized for speed to allow users to search the non-redundant RNA 3D structure database on a personal computer in a matter of minutes.
Sarver
Michael
M
Department of Mathematics and Statistics, Bowling Green State University, Bowling Green, OH 43403, USA.
Zirbel
Craig L
CL
Stombaugh
Jesse
J
Mokdad
Ali
A
Leontis
Neocles B
NB
D000465:Algorithms; D019295:Computational Biology / Q000379:methods*; D009690:Nucleic Acid Conformation; D012313:RNA / Q000737:chemistry*; D012338:RNA, Ribosomal, 23S / Q000737:chemistry; D012341:RNA, Ribosomal, 5S / Q000737:chemistry; D012276:Ricin / Q000737:chemistry; D012984:Software
D012338:RNA, Ribosomal, 23S; D012341:RNA, Ribosomal, 5S; D012313:RNA; D012276:Ricin
Nucleic Acids Res. 2003 Nov 1;31(21):6249-57
14576313
Nucleic Acids Res. 2000 Jan 1;28(1):235-42
10592235
Structure. 2005 Mar;13(3):381-8
15766539
Nucleic Acids Res. 2003 Jul 1;31(13):3450-60
12824344
Cell. 2001 Nov 30;107(5):679-88
11733066
Curr Opin Struct Biol. 2003 Jun;13(3):300-8
12831880
J Mol Biol. 1998 Dec 18;284(5):1465-78
9878364
Science. 2000 Aug 11;289(5481):905-20
10937989
Nucleic Acids Res. 2004;32(8):2342-52
15121895
Bioinformatics. 2005 Aug 1;21(15):3320-1
15905274
RNA. 2006 Apr;12(4):533-41
16484377
Nucleic Acids Res. 2003 Aug 15;31(16):4755-61
12907716
RNA. 2001 Apr;7(4):499-512
11345429
Curr Opin Struct Biol. 2006 Jun;16(3):279-87
16713707
Chem Rev. 1999 Nov 10;99(11):3247-76
11749516
J Mol Biol. 2001 Oct 12;313(1):229-37
11601858
Curr Opin Struct Biol. 2005 Jun;15(3):302-8
15963891
Science. 1991 Sep 13;253(5025):1255-60
1716375
J Comput Aided Mol Des. 2003 Aug;17(8):537-49
14703124
J Mol Biol. 1994 Mar 18;237(1):125-56
8133513
Nature. 2000 Sep 21;407(6802):327-39
11014182
EMBO J. 2001 Aug 1;20(15):4214-21
11483524
Mol Cell Biol. 2005 Jun;25(11):4752-66
15899876
Bioinformatics. 2005 Sep 1;21 Suppl 2:ii47-53
16204124
Science. 2005 Nov 4;310(5749):827-34
16272117
J Chem Inf Model. 2005 May-Jun;45(3):696-707
15921459
Nucleic Acids Res. 2005;33(8):2395-409
15860776
Nucleic Acids Res. 2004;32(5):1666-77
15016910
Nucleic Acids Res. 2002 Aug 15;30(16):3497-531
12177293
J Mol Biol. 2001 May 18;308(5):919-36
11352582
Proc Natl Acad Sci U S A. 2005 Sep 20;102(38):13392-7
16157868
Proc Natl Acad Sci U S A. 2003 Nov 25;100(24):13904-9
14612579
Nat Struct Mol Biol. 2005 Jan;12(1):82-9
15580277
Nucleic Acids Res. 2004;32(22):6650-9
15608296
Nucleic Acids Res. 2005 Jan 1;33(Database issue):D233-7
15608185
RNA. 2005 Apr;11(4):412-23
15769871
Biochem Soc Trans. 2005 Jun;33(Pt 3):485-7
15916548
Nucleic Acids Res. 2004 Jan 1;32(Database issue):D223-5
14681399
Proc Natl Acad Sci U S A. 2003 Jun 10;100(12):7003-8
12782785
Biochimie. 2002 Sep;84(9):961-73
12458088
Nucleic Acids Res. 2005;33(17):5677-90
16214802
Nature. 2004 Jul 1;430(6995):45-50
15175762