Manzourolajdad and Arnold BMC Bioinformatics (2015) 16:133 DOI 10.1186/s12859-015-0523-2 METHODOLOGY ARTICLE Open Access Secondary structural entropy in RNA switch (Riboswitch) identification Amirhossein Manzourolajdad1,2* and Jonathan Arnold1,3 *Correspondence:
[email protected] Abstract 1 Institute of Bioinformatics, Background: RNA regulatory elements play a significant role in gene regulation. University of Georgia, Davison Life Sciences Bldg, Room B118B, 120 Riboswitches, a widespread group of regulatory RNAs, are vital components of many Green St, 30602 Athens, USA bacterial genomes. These regulatory elements generally function by forming a 2 National Center for Biotechnology ligand-induced alternative fold that controls access to ribosome binding sites or other Information (NCBI), NIH, Building 38A, RM 6S614K, 8600 Rockville regulatory sites in RNA. Riboswitch-mediated mechanisms are ubiquitous across Pike, 20894 Bethesda, USA bacterial genomes. A typical class of riboswitch has its own unique structural and Full list of author information is available at the end of the article biological complexity, making de novo riboswitch identification a formidable task. Traditionally, riboswitches have been identified through comparative genomics based on sequence and structural homology. The limitations of structural-homology-based approaches, coupled with the assumption that there is a great diversity of undiscovered riboswitches, suggests the need for alternative methods for riboswitch identification, possibly based on features intrinsic to their structure. As of yet, no such reliable method has been proposed. Results: We used structural entropy of riboswitch sequences as a measure of their secondary structural dynamics. Entropy values of a diverse set of riboswitches were compared to that of their mutants, their dinucleotide shuffles, and their reverse complement sequences under different stochastic context-free grammar folding models.