Dynamics in Sequence Space for RNA Secondary Structure Design

J Chem Theory Comput. 2012 Oct 9;8(10):3663-70. doi: 10.1021/ct300267j. Epub 2012 Aug 23.

Abstract

We have implemented a method for the design of RNA sequences that should fold to arbitrary secondary structures. A popular energy model allows one to take the derivative with respect to composition, which can then be interpreted as a force and used for Newtonian dynamics in sequence space. Combined with a negative design term, one can rapidly sample sequences which are compatible with a desired secondary structure via simulated annealing. Results for 360 structures were compared with those from another nucleic acid design program using measures such as the probability of the target structure and an ensemble-weighted distance to the target structure.