Statistical mechanics of RNA folding: importance of alphabet size

We construct a minimalist model of RNA secondary-structure formation and use it to study the mapping from sequence to structure. There are strong, qualitative differences between two-letter and four or six-letter alphabets. With only two kinds of bases, there are many alternate folding configurations, yielding thermodynamically stable ground-states only for a small set of structures of high designability, i.e., total number of associated sequences. In contrast, sequences made from four bases, as found in nature, or six bases have far fewer competing folding configurations, resulting in a much greater average stability of the ground state.Comment: 7 figures; uses revtex

Evolving new genetic codes

Although the genetic code is almost universal, natural variations exist that have caused evolutionary biologists to speculate about codon evolution. There are two predominant hypotheses that specify either a gradual (ambiguous intermediate) or stochastic (codon capture) change in the code. These hypotheses are similar to two biotechnology techniques that have been used to engineer the genetic code: a 'top down' approach, in which the whole organism is evolved for the ability to incorporate unnatural amino acids, and a 'bottom up' approach, in which aminoacyl-tRNA synthetases and their cognate tRNAs are engineered. The biotechnology experiments provide insights into natural codon evolution, and a combination of these approaches should enable the evolution of organisms that can incorporate unnatural amino acids throughout their proteomes

AptaCluster – A Method to Cluster HT-SELEX Aptamer Pools and Lessons from Its Application

Systematic Evolution of Ligands by EXponential Enrichment (SELEX) is a well established experimental procedure to identify aptamers - synthetic single-stranded (ribo)nucleic molecules that bind to a given molecular target. Recently, new sequencing technologies have revolutionized the SELEX protocol by allowing for deep sequencing of the selection pools after each cycle. The emergence of High Throughput SELEX (HT-SELEX) has opened the field to new computational opportunities and challenges that are yet to be addressed. To aid the analysis of the results of HT-SELEX and to advance the understanding of the selection process itself, we developed AptaCluster. This algorithm allows for an efficient clustering of whole HT-SELEX aptamer pools; a task that could not be accomplished with traditional clustering algorithms due to the enormous size of such datasets. We performed HT-SELEX with Interleukin 10 receptor alpha chain (IL-10RA) as the target molecule and used AptaCluster to analyze the resulting sequences. AptaCluster allowed for the first survey of the relationships between sequences in different selection rounds and revealed previously not appreciated properties of the SELEX protocol. As the first tool of this kind, AptaCluster enables novel ways to analyze and to optimize the HT-SELEX procedure. Our AptaCluster algorithm is available as a very fast multiprocessor implementation upon request