Network Properties of the Ensemble of RNA Structures

PLoS One. 2015 Oct 21;10(10):e0139476. doi: 10.1371/journal.pone.0139476. eCollection 2015.

Abstract

We describe the first dynamic programming algorithm that computes the expected degree for the network, or graph G = (V, E) of all secondary structures of a given RNA sequence a = a1, …, an. Here, the nodes V correspond to all secondary structures of a, while an edge exists between nodes s, t if the secondary structure t can be obtained from s by adding, removing or shifting a base pair. Since secondary structure kinetics programs implement the Gillespie algorithm, which simulates a random walk on the network of secondary structures, the expected network degree may provide a better understanding of kinetics of RNA folding when allowing defect diffusion, helix zippering, and related conformation transformations. We determine the correlation between expected network degree, contact order, conformational entropy, and expected number of native contacts for a benchmarking dataset of RNAs. Source code is available at http://bioinformatics.bc.edu/clotelab/RNAexpNumNbors.

Publication types

  • Research Support, U.S. Gov't, Non-P.H.S.

MeSH terms

  • Algorithms*
  • Base Pairing
  • Base Sequence
  • Computational Biology / methods*
  • Kinetics
  • Models, Molecular
  • Molecular Sequence Data
  • Nucleic Acid Conformation*
  • RNA / chemistry*
  • RNA / genetics
  • RNA Folding

Substances

  • RNA

Grants and funding

PC received funding from the National Science Foundation under grant DBI-1262439 (www.nsf.gov). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.