Selection for mutational robustness in finite populations

We investigate the evolutionary dynamics of a finite population of RNA sequences replicating on a neutral network. Despite the lack of differential fitness between viable sequences, we observe typical properties of adaptive evolution, such as increase of mean fitness over time and punctuated-equilibrium transitions, after initial mutation-selection balance has been reached. We find that a product of population size and mutation rate of approximately 30 or larger is sufficient to generate selection pressure for mutational robustness, even if the population size is orders of magnitude smaller than the neutral network on which the population resides. Our results show that quasispecies effects and neutral drift can occur concurrently, and that the relative importance of each is determined by the product of population size and mutation rate.     

Quasi-species evolution in subdivided populations favours maximally deleterious mutations

Most models of quasi-species evolution predict that populations will evolve to occupy areas of sequence space with the greatest concentration of neutral sequences, thus minimizing the deleterious mutation rate and creating mutationally 'robust' genomes. In contrast, empirical studies of the principal model of quasi-species evolution, RNA viruses, suggest that the effects of deleterious mutations are more severe than in similar DNA-based microbes. We demonstrate that populations divided into discrete patches connected by dispersal may favour genotypes where the deleterious effect of non-neutral mutations is maximized. This effect is especially strong in the absence of back mutation and when the amount of time spent in hosts prior to dispersal is intermediate. Our results indicate that RNA viruses that produce acute infections initiated by a small number of virions are expected to evolve fragile genetic architectures when compared with other RNA viruses.     

Statistics of RNA melting kinetics

We present and study the behavior of a simple kinetic model for the melting of RNA secondary structures, given that those structures are known. The model is then used as a map that. assigns structure dependent overall rate constants of melting (or refolding) to a sequence. This induces a landscape of reaction rates, or activation energies, over the space of sequences with fixed length. We study the distribution and the correlation structure of these activation energies. Correspondence to: P. Schuster     

Simulating the Folding Pathway of RNA Secondary Structure Using the Modified Ant Colony Algorithm

<正> A new method for simulating the folding pathway of RNA secondary structure using the modified ant colony algorithmis proposed.For a given RNA sequence,the set of all possible stems is obtained and the energy of each stem iscalculated and stored at the initial stage.Furthermore,a more realistic formula is used to compute the energy ofmulti-branch loop in the following iteration.Then a folding pathway is simulated,including such processes as constructionof the heuristic information,the rule of initializing the pheromone,the mechanism of choosing the initial andnext stem and the strategy of updating the pheromone between two different stems.Finally by testing RNA sequences withknown secondary structures from the public databases,we analyze the experimental data to select appropriate values forparameters.The measure indexes show that our procedure is more consistent with phylogenetically proven structures thansoftware RNAstructure sometimes and more effective than the standard Genetic Algorithm.     

Translational standby sites: how ribosomes may deal with the rapid folding kinetics of mRNA

We have previously shown that stable base-pairing at a translational initiation site in Escherichia coli can inhibit translation by competing with the binding of ribosomes. When the base-pairing is not too strong, this competition is won by the ribosomes, resulting in efficient translation from a structured ribosome binding site (RBS). We now re-examine these results in the light of RNA folding kinetics and find that the window during which a folded RBS is open is generally much too short to recruit a 30S ribosomal subunit from the cytoplasm. We argue that to achieve efficient expression, a 30S subunit must already be in contact with the mRNA while this is still folded, to shift into place as soon as the structure opens. Single-stranded regions flanking the structure may constitute a standby site, to which the 30S subunit can attach non-specifically. We propose a steady-state kinetic model for the early steps of translational initiation and use this to examine various quantitative aspects of standby binding. The kinetic model provides an explanation of why the earlier equilibrium competition model predicted implausibly high 30S-mRNA affinities. Because all RNA is structured to some degree, standby binding is probably a general feature of translational initiation.     

Nanomanipulation of Single RNA Molecules by Optical Tweezers

A large portion of the human genome is transcribed but not translated. In this post genomic era, regulatory functions of RNA have been shown to be increasingly important. As RNA function often depends on its ability to adopt alternative structures, it is difficult to predict RNA three-dimensional structures directly from sequence. Single-molecule approaches show potentials to solve the problem of RNA structural polymorphism by monitoring molecular structures one molecule at a time. This work presents a method to precisely manipulate the folding and structure of single RNA molecules using optical tweezers. First, methods to synthesize molecules suitable for single-molecule mechanical work are described. Next, various calibration procedures to ensure the proper operations of the optical tweezers are discussed. Next, various experiments are explained. To demonstrate the utility of the technique, results of mechanically unfolding RNA hairpins and a single RNA kissing complex are used as evidence. In these examples, the nanomanipulation technique was used to study folding of each structural domain, including secondary and tertiary, independently. Lastly, the limitations and future applications of the method are discussed.     

Bistable secondary structures of small RNAs and their structural probing by comparative imino proton NMR spectroscopy

We investigate 25-34 nucleotide RNA sequences, that have been rationally designed to adopt two different secondary structures that are in thermodynamic equilibrium. Experimental evidence for the co-existence of the two conformers results from the NH...N 1H NMR spectra. When compared to the NH...N 1H NMR spectra of appropriate reference sequences the equilibrium position is easily quantifiable even without the assignment of the individual NH resonances. The reference sequences represent several Watson-Crick base-paired double helical segments, each encountered in either of the two conformers of the bistable target sequence. In addition, we rationalize the influence of nucleotide mutations on the equilibrium position of one of the bistable RNA sequences. The approach further allows a detailed thermodynamic analysis and the evaluation of secondary structure predictions for multistable RNAs obtained by computational methods.     

ProbKnot: Fast prediction of RNA secondary structure including pseudoknots

It is a significant challenge to predict RNA secondary structures including pseudoknots. Here, a new algorithm capable of predicting pseudoknots of any topology, ProbKnot, is reported. ProbKnot assembles maximum expected accuracy structures from computed base-pairing probabilities in O(N²) time, where N is the length of the sequence. The performance of ProbKnot was measured by comparing predicted structures with known structures for a large database of RNA sequences with fewer than 700 nucleotides. The percentage of known pairs correctly predicted was 69.3%. Additionally, the percentage of predicted pairs in the known structure was 61.3%. This performance is the highest of four tested algorithms that are capable of pseudoknot prediction. The program is available for download at: http://rna.urmc.rochester.edu/RNAstructure.html.     

Epistasis can lead to fragmented neutral spaces and contingency in evolution

In evolution, the effects of a single deleterious mutation can sometimes be compensated for by a second mutation which recovers the original phenotype. Such epistatic interactions have implications for the structure of genome space--namely, that networks of genomes encoding the same phenotype may not be connected by single mutational moves. We use the folding of RNA sequences into secondary structures as a model genotype-phenotype map and explore the neutral spaces corresponding to networks of genotypes with the same phenotype. In most of these networks, we find that it is not possible to connect all genotypes to one another by single point mutations. Instead, a network for a phenotypic structure with n bonds typically fragments into at least 2(n) neutral components, often of similar size. While components of the same network generate the same phenotype, they show important variations in their properties, most strikingly in their evolvability and mutational robustness. This heterogeneity implies contingency in the evolutionary process.     

正在崛起的RNA纳米技术

RNA纳米技术得益于纽约大学西曼(Nadrian C.Seeman)教授开创的DNA纳米技术,RNA是由腺嘌呤(A)、尿嘧啶(U)、鸟嘌呤(G)和胞嘧啶(C)构成的一种核糖核酸高分子,与DNA的Watson-Crick碱基配对(A-T,G-C)的双螺旋链的结构不完全一样,RNA的二级结构里经常出现一些非传统的碱基配对如环环相互作用,这些非传统配对促使RNA分子折叠成刚性结构。本文综述了正在崛起的RNA纳米技术,列举了一些著名的实验,如郭培宣(Peixuan Guo)等从自然界的phi29噬菌体中发现的pRNA纳米马达是由六个小RNA分子构成的六环结构,Jaeger等发展了RNA构造术(RNA-tectonics),根据已知的RNA分子的碱基和非传统配对,他们设计利用小RNA分子构造二聚体、一维线性多聚体、和二维网状的七巧板迷宫(jigsaw puzzle)等图案,用tRNA分子或设计用几条RNA分子来构建多面体如立方体和八面体等立体结构等。RNA纳米技术正在崛起,它将在医学、生物技术、合成生物学和纳米技术领域扮演重要的角色。     

Structure of the pyrimidine-rich internal loop in the poliovirus 3'-UTR: the importance of maintaining pseudo-2-fold symmetry in RNA helices containing two adjacent non-canonical base-pairs

Formation of non-canonical base-pairs in RNA often plays a very important functional role. In addition they frequently serve as factors in stabilizing the secondary structure elements that provide the frame of large compact RNA structures. Here we describe the structure of an internal loop containing a 5'CU3'/5'UU3' non-canonical tandem base-pair motif, which is conserved within the 3'-UTR of poliovirus-like enteroviruses. Structural details reveal striking regularities of the local helix geometry, resulting from alternating geometrical adjustments, which are important for understanding and predicting stabilities and configurations of tandem non-canonical base-pairs. The C-U and U-U base-pairs severely contract the minor groove of the sugar-phosphate backbone, which might be important for protein recognition or binding to other RNA elements.