This Article Full Text Full Text (PDF) Supplement All Versions of this Article: biophysj.107.113415v1 94/11/4202    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Google Scholar Articles by Fulle, S. Articles by Gohlke, H. PubMed PubMed Citation Articles by Fulle, S. Articles by Gohlke, H.

Analyzing the Flexibility of RNA Structures by Constraint Counting

Department of Biological Sciences, Molecular Bioinformatics Group, J. W. Goethe-University, Frankfurt, Germany

Correspondence: Address reprint requests to Holger Gohlke, Max-von-Laue-Str. 9, 60438 Frankfurt, Germany. Tel.: 49-69-798-29411; Fax: 49-69-798-29527; E-mail: gohlke{at}bioinformatik.uni-frankfurt.de.

RNA requires conformational dynamics to undergo its diverse functional roles. Here, a new topological network representation of RNA structures is presented that allows analyzing RNA flexibility/rigidity based on constraint counting. The method extends the FIRST approach, which identifies flexible and rigid regions in atomic detail in a single, static, three-dimensional molecular framework. Initially, the network rigidity of a canonical A-form RNA is analyzed by counting on constraints of network elements of increasing size. These considerations demonstrate that it is the inclusion of hydrophobic contacts into the RNA topological network that is crucial for an accurate flexibility prediction. The counting also explains why a protein-based parameterization results in overly rigid RNA structures. The new network representation is then validated on a tRNA^ASP structure and all NMR-derived ensembles of RNA structures currently available in the Protein Data Bank (with chain length 40). The flexibility predictions demonstrate good agreement with experimental mobility data, and the results are superior compared to predictions based on two previously used network representations. Encouragingly, this holds for flexibility predictions as well as mobility predictions obtained by constrained geometric simulations on these networks. Potential applications of the approach to analyzing the flexibility of DNA and RNA/protein complexes are discussed.