The block spectrum of RNA pseudoknot structures

• Published in: Journal of mathematical biology 2019-08, Vol.79 (3), p.791-822
• Main Authors: Li, Thomas J. X., Burris, Christie S., Reidys, Christian M.
• Format: Article
• Language: English
• Subjects: Applications of Mathematics; Binomial distribution; Building components; Event-related potentials; Mathematical and Computational Biology; Mathematics; Mathematics and Statistics; Nesting; Polynomials; Rainbows; Ribonucleic acid; RNA; Substructures
• ISSN: 0303-6812; 1432-1416
• DOI: 10.1007/s00285-019-01379-8

In this paper we analyze the length-spectrum of blocks in γ -structures. γ -structures are a class of RNA pseudoknot structures that play a key role in the context of polynomial time RNA folding. A γ -structure is constructed by nesting and concatenating specific building components having topological genus at most γ . A block is a substructure enclosed by crossing maximal arcs with respect to the partial order induced by nesting. We show that, in uniformly generated γ -structures, there is a significant gap in this length-spectrum, i.e., there asymptotically almost surely exists a unique longest block of length at least n - O ( n 1 / 2 ) and that with high probability any other block has finite length. For fixed γ , we prove that the length of the complement of the longest block converges to a discrete limit law, and that the distribution of short blocks of given length tends to a negative binomial distribution in the limit of long sequences. We refine this analysis to the length spectrum of blocks of specific pseudoknot types, such as H-type and kissing hairpins. Our results generalize the rainbow spectrum on secondary structures by the first and third authors and are being put into context with the structural prediction of long non-coding RNAs.