Sequence-structure relations of biopolymers.
DNA data is transcribed into single-stranded RNA, which folds into specific molecular structures. In this paper we pose the question to what extent sequence- and structure-information correlate. We view this correlation as structural semantics of sequence data that allows for a different interpretation than conventional sequence alignment. Structural semantics could enable us to identify more general embedded ‘patterns’ in DNA and RNA sequences.



We compute the partition function of sequences with respect to a fixed structure and connect this computation to the mutual information of a sequence–structure pair for RNA secondary structures. We present a Boltzmann sampler and obtain the a priori probability of specific sequence patterns. We present a detailed analysis for the three PDB-structures, 2JXV (hairpin), 2N3R (3-branch multi-loop) and 1EHZ (tRNA). We localize specific sequence patterns, contrast the energy spectrum of the Boltzmann sampled sequences versus those sequences that refold into the same structure and derive a criterion to identify native structures. We illustrate that there are multiple sequences in the partition function of a fixed structure, each having nearly the same mutual information, that are nevertheless poorly aligned. This indicates the possibility of the existence of relevant patterns embedded in the sequences that are not discoverable using alignments.



The source code is freely available at http://staff.vbi.vt.edu/fenixh/Sampler.zip



duckcr@vbi.vt.edu



Supplementary data are available at Bioinformatics online.
Biocomplexity Institute of Virginia Tech, Virginia Tech University, Blacksburg, VA, USA.
Christopher
C
Barrett
Biocomplexity Institute of Virginia Tech, Virginia Tech University, Blacksburg, VA, USA.
Fenix W
FW
Huang
Biocomplexity Institute of Virginia Tech, Virginia Tech University, Blacksburg, VA, USA.
Christian M
CM
Reidys
D000465:Algorithms; D019295:Computational Biology / Q000379:methods*; D009690:Nucleic Acid Conformation; D011336:Probability; D012313:RNA / Q000737:chemistry* / Q000378:metabolism; D017423:Sequence Analysis, RNA / Q000379:methods*; D012984:Software
D012313:RNA