An unbiased adaptive sampling algorithm for the exploration of RNA mutational landscapes under evolutionary pressure.
The analysis of the relationship between sequences and structures (i.e., how mutations affect structures and reciprocally how structures influence mutations) is essential to decipher the principles driving molecular evolution, to infer the origins of genetic diseases, and to develop bioengineering applications such as the design of artificial molecules. Because their structures can be predicted from the sequence data only, RNA molecules provide a good framework to study this sequence-structure relationship. We recently introduced a suite of algorithms called RNAmutants which allows a complete exploration of RNA sequence-structure maps in polynomial time and space. Formally, RNAmutants takes an input sequence (or seed) to compute the Boltzmann-weighted ensembles of mutants with exactly k mutations, and sample mutations from these ensembles. However, this approach suffers from major limitations. Indeed, since the Boltzmann probabilities of the mutations depend of the free energy of the structures, RNAmutants has difficulties to sample mutant sequences with low G+C-contents. In this article, we introduce an unbiased adaptive sampling algorithm that enables RNAmutants to sample regions of the mutational landscape poorly covered by classical algorithms. We applied these methods to sample mutations with low G+C-contents. These adaptive sampling techniques can be easily adapted to explore other regions of the sequence and structural landscapes which are difficult to sample. Importantly, these algorithms come at a minimal computational cost. We demonstrate the insights offered by these techniques on studies of complete RNA sequence structures maps of sizes up to 40 nucleotides. Our results indicate that the G+C-content has a strong influence on the size and shape of the evolutionary accessible sequence and structural spaces. In particular, we show that low G+C-contents favor the apparition of internal loops and thus possibly the synthesis of tertiary structure motifs. On the other hand, high G+C-contents significantly reduce the size of the evolutionary accessible mutational landscapes.
School of Computer Science & McGill Center for Bioinformatics, McGill University, Montreal, Canada. jeromew@cs.mcgill.ca
Jérôme
J
Waldispühl
Yann
Y
Ponty
D000465:Algorithms; D001482:Base Composition; D020029:Base Pairing; D001483:Base Sequence; D003198:Computer Simulation; D003627:Data Interpretation, Statistical; D019143:Evolution, Molecular; D008957:Models, Genetic; D008958:Models, Molecular; D008969:Molecular Sequence Data; D009154:Mutation; D009690:Nucleic Acid Conformation; D012313:RNA / Q000737:chemistry / Q000235:genetics*; D020871:RNA Stability; D012641:Selection, Genetic; D013816:Thermodynamics
D012313:RNA