  Conformational entropy for atomic-level, three dimensional biomolecules is
known experimentally to play an important role in protein-ligand
discrimination, yet reliable computation of entropy remains a difficult
problem. Here we describe the first two accurate and efficient algorithms to
compute the conformational entropy for RNA secondary structures, with respect
to the Turner energy model, where free energy parameters are determined from UV
aborption experiments. An algorithm to compute the derivational entropy for RNA
secondary structures had previously been introduced, using stochastic context
free grammars (SCFGs). However, the numerical value of derivational entropy
depends heavily on the chosen context free grammar and on the training set used
to estimate rule probabilities. Using data from the Rfam database, we determine
that both of our thermodynamic methods, which agree in numerical value, are
substantially faster than the SCFG method. Thermodynamic structural entropy is
much smaller than derivational entropy, and the correlation between
length-normalized thermodynamic entropy and derivational entropy is moderately
weak to poor. In applications, we plot the structural entropy as a function of
temperature for known thermoswitches, determine that the correlation between
hammerhead ribozyme cleavage activity and total free energy is improved by
including an additional free energy term arising from conformational entropy,
and plot the structural entropy of windows of the HIV-1 genome. Our software
RNAentropy can compute structural entropy for any user-specified temperature,
and supports both the Turner'99 and Turner'04 energy parameters. It follows
that RNAentropy is state-of-the-art software to compute RNA secondary structure
conformational entropy. The software is available at
http://bioinformatics.bc.edu/clotelab/RNAentropy.
