Prediction of sequentially optimal RNA secondary structures

A rigorous mathematical modeling of the RNA sequential folding process during transcription is proposed. It is based, at each transcription step, on a homogeneous markovian jump process, the state space of which is the set of structures constructible on the part of the RNA already transcribed. A theoretical formula permitting the computation of the structures probabilities at the end of the RNA transcription is derived. Successive approximations, aimed at reducing the size of the state space, permit the design of a prediction algorithm. The algorithm is tested on some structural RNAs (tRNA, 5S, 16S, hammerhead, ...), results are discussed and possible improvements are proposed.