The Robustness of Naturally and Artificially Selected Nucleic Acid Secondary Structures

Thermodynamic stability and mutational robustness of secondary structure are critical to the function and evolutionary longevity of RNA molecules. We hypothesize that natural and artificial selection for functional molecules favors the formation of structures that are stable to both thermal and mutational perturbation. There is little direct evidence, however, that functional RNA molecules have been selected for their stability. Here we use thermodynamic secondary structure prediction algorithms to compare the thermal and mutational robustness of over 1000 naturally and artificially evolved molecules. Although we find evidence for the evolution of both types of stability in both sets of molecules, the naturally evolved functional RNA molecules were significantly more stable than those selected in vitro, and artificially evolved catalysts (ribozymes) were more stable than artificially evolved binding species (aptamers). The thermostability of RNA molecules bred in the laboratory is probably not constrained by a lack of suitable variation in the sequence pool but, rather, by intrinsic biases in the selection process.     

RNA as a catalyst: Natural and designed ribozymes

RNA can catalyse chemical reactions through its ability to fold into complex three-dimensional structures and to specifically bind small molecules and divalent metal ions. The 2′-hydroxyl groups of the ribose moieties contribute to this exceptional reactivity of RNA, compared to DNA. RNA is not only able to catalyse phosphate ester transfer reactions in ribonucleic acids, but can also show aminoacyl esterase activity, and is probably able to promote peptide bond formation. Bearing its potential for functioning both as a genome and as a gene product, RNA is suitable for in vitro evolution experiments enabling the selection of molecules with new properties. The growing repertoire of RNA catalysed reactions will establish RNA as a primordial molecule in the evolution of life.     

Competitive speciation

A new mode of speciation, competitive speciation, is suggested. It assumes that fitness is depressed by the density of a phenotype's competitors, and that the adaptive landscape of phenotypes is complex. From this it follows that some intermediate forms may be fit if and only if some extreme forms are rare or absent. Subsequent to the evolution and population growth of both extreme forms, the intermediate may disappear and homogamy evolve among each of the extremes because of disruptive selection If so, sympatric speciation has occurred and niche space has been rendered into discrete segments.
The limitations of the forces leading to competitive speciation are explored. Competitive speciation is discussed in relation to stasipatric speciation and host race formation. It may be responsible for both. Finally the rates of geographical speciation and polyploidy are compared to those of competitive speciation. The latter should be almost as fast as polyploidy and may be at the root of adaptive radiation. Unlike either polyploidy or geographical speciation, competitive speciation accelerates when species diversity declines.     

Low Selection Pressure Aids the Evolution of Cooperative Ribozyme Mutations in Cells

Understanding the evolution of functional RNA molecules is important for our molecular understanding of biology. Here we tested experimentally how two evolutionary parameters, selection pressure and recombination, influenced the evolution of an evolving RNA population. This was done using four parallel evolution experiments that employed low or gradually increasing selection pressure, and recombination events either at the end or dispersed throughout the evolution. As model system, a trans-splicing group I intron ribozyme was evolved in Escherichia coli cells over 12 rounds of selection and amplification, including mutagenesis and recombination. The low selection pressure resulted in higher efficiency of the evolved ribozyme populations, whereas differences in recombination did not have a strong effect. Five mutations were responsible for the highest efficiency. The first mutation swept quickly through all four evolving populations, whereas the remaining four mutations accumulated later and more efficiently under low selection pressure. To determine why low selection pressure aided this evolution, all evolutionary intermediates between the wild type and the 5-mutation variant were constructed, and their activities at three different selection pressures were determined. The resulting fitness profiles showed a high cooperativity among the four late mutations, which can explain why high selection pressure led to inefficient evolution. These results show experimentally how low selection pressure can benefit the evolution of cooperative mutations in functional RNAs.     

Rate of polymer formation and entropy production during competitive replication

Summary The rate of increase in the mean polymer formation rate constant during competitive replication byQ RNA variants (Kramer et al., 1974) has been shown to agree statistically with the variance in their formation rate constants. This result demonstrates that Fisher's fundamental theorem of natural selection (Fisher, 1930) can define time variations in the mean rate of synthesis for a heterogeneous population of replicating polymers. It was also revealed that RNA replication, far from equilibrium, accompanied a progressive decrease in the order of the entropy production derivative, with respect to time, that reached a maximum (with the next higher order being zero). Maximization of entropy at equilibrium, in compliance with the second law of thermodynamics, therefore appears as a natural extension of the earlier non-equilibrium pattern of entropy production within the system. The order of the zero-valued entropy production derivative was shown to be determined by the chemical affinity, and its rate of decrease was specified by the mean polymer formation rate constant.     

A test of alternative hypotheses for the evolution of reproductive isolation between spadefoot toads: support for the reinforcement hypothesis

Abstract How do species that interbreed become reproductively isolated? If hybrids are less fit than parental types, natural selection should promote reproductive isolation by favoring the evolution of premating mechanisms that prevent hybridization (a process termed reinforcement). Although reinforcement should generate a decline in hybridization over time, countervailing forces of gene flow and recombination are thought to preclude natural selection from enhancing and finalizing reproductive isolation. Here, I present recent estimates of hybridization frequency between two species of spadefoot toad, Spea multiplicata and S. bombifrons. I compare these recent measures of hybrid frequency with previously published estimates and show that hybridization between these species has declined precipitously over the past 27 years. Although previous studies suggest that reinforcement possibly accounts for this decline in hybrids over time, three alternative hypotheses also can explain the observed decrease in hybridization. First, if one of the two interacting species becomes rare, opportunities for and incidence of hybridization may decrease. Second, if one of the two interacting species is initially rare, hybridization may be initially common if the rare species has difficulty locating conspecific mates. Third, if hybrids are produced only in particular environments, hybrid frequency may decline if habitat changes result in loss of those environments that promote hybrid formation. I found no support for these three alternative explanations of the decline in hybrids. Instead, reinforcement appears to best account for the evolution of enhanced reproductive isolation between these species. Moreover, the finding that hybridization declined precipitously in only 27 years suggests that many systems that have undergone reinforcement may be overlooked because reproductive isolation between the interacting populations or species may already be complete.     

Darwinian aspects of molecular evolution at sublinear propagation rates

Mean fitness is non-decreasing in the symmetry sector of the frequency trajectory followed in competitive replication at sublinear propagation rates (parabolic time course). This sector contains the pairwise symmetric distribution of species frequencies and its neighboring states, and represents at least half the possible states of an evolving sublinear system. States in the non-symmetry sector produce a negative rate of change in mean fitness. The heterogeneous steady state attained in a finite sublinear system is destabilized by formation of a variant with above-threshold fitness. Evolution in the post-steady-state interval elevates the fitness threshold for coexistence. Contrary to the proposition that ‘parabolic growth invariably results in the survival of all competing species’, only species with sufficient fitness to avoid subthreshold fitness survive. An erratum to this article is available at .     

A combined in vitro / in vivo selection for polymerases with novel promoter specificities

Background  

The DNA-dependent RNA polymerase from T7 bacteriophage (T7 RNAP) has been extensively characterized, and like other phage RNA polymerases it is highly specific for its promoter. A combined in vitro / in vivo selection method has been developed for the evolution of T7 RNA polymerases with altered promoter specificities. Large (10³ – 10⁶) polymerase libraries were made and cloned downstream of variant promoters. Those polymerase variants that can recognize variant promoters self-amplify both themselves and their attendent mRNAs in vivo. Following RT / PCR amplification in vitro, the most numerous polymerase genes are preferentially cloned and carried into subsequent rounds of selection.     

Ancient RNA world

Based on the modern concept of multiplicity of RNA functions, certain new mechanisms, which could have played a key role in the origin and evolution of the ancient RNA world, are discussed. In particular, the reaction of spontaneous transesterification of polyribonucleotides, which was discovered by A.B. Chetverin and colleagues, could result in elongation of short initial oligoribonucleotides and produce sequence variants for further natural selection of accidentally arising functionally active molecules. Further, the formation of mixed molecular colonies of RNA on moist solid media, such as clays, could have provided compartmentation of functional RNA ensembles in the absence of envelopes or membranes, which were necessary for further evolution of the RNA world. The systematic exponential enrichment of the RNA population with functionally the best molecules because of alternating dissolution of colonies due to flooding and formation of new colonies due to drying in ancient pools (“primordial natural SELEX”) could be critical for the evolutionary process in the RNA world.     

Evolution in vitro: analysis of a lineage of ribozymes

Background: Catalytic RNAs, or ribozymes, possessing both a genotype and a phenotype, are ideal molecules for evolution experiments in vitro. A large, heterogeneous pool of RNAs can be subjected to multiple rounds of selection, amplification and mutation, leading to the development of variants that have some desired phenotype. Such experiments allow the investigator to correlate specific genetic changes with quantifiable alterations of the catalytic properties of the RNA. In addition, patterns of evolutionary change can be discerned through a detailed examination of the genotypic composition of the evolving RNA population. Results: Beginning with a pool of 10(13) variants of the Tetrahymena ribozyme, we carried out in vitro evolution experiments that led to the generation of ribozymes with the ability to cleave an RNA substrate in the presence of Ca2+ ions, an activity that does not exist for the wild-type molecule. Over the course of 12 generations, a seven-error variant emerged that has substantial Ca(2+)-dependent RNA-cleavage activity. Advantageous mutations increased in frequency in the population according to three distinct dynamics--logarithmic, linear and transient. Through a comparative analysis of 31 individual variants, we infer how certain mutations influence the catalytic properties of the ribozyme. Conclusions: In vitro evolution experiments make it possible to elucidate important aspects of both evolutionary biology and structural biochemistry on a reasonable short time scale.     

Analysis of energetically biased transcripts of viruses and transposable elements

RNA interference (RNAi) is a natural endogenous process by which double-stranded RNA molecules trigger potent and specific gene silencing in eukaryotic cells and is characterized by target RNA cleavage. In mammals, small interfering RNAs (siRNAs) are the trigger molecules of choice and constitute a new class of RNA-based antiviral agents. In an efficient RNAi response, the antisense strand of siRNAs must enter the RNA-induced silencing complex (RISC) in a process mediated by thermodynamic features. In this report, we hypothesize that silent mutations capable of inverting thermodynamic properties can promote resistance to siRNAs. Extensive computational analyses were used to assess whether continuous selective pressure that promotes such mutations could lead to the emergence of viral strains completely resistant to RNAi (i.e., prone to transfer only the sense strands to RISC). Based on our findings, we propose that, although synonymous mutations may produce functional resistance, this strategy cannot be systematically adopted by viruses since the longest RNAi-refractory sequence is only 10 nt long. This finding also suggests that all mRNAs display fluctuating thermodynamic landscapes and that, in terms of thermodynamic features, RNAi is a very efficient antiviral system since there will always be sites susceptible to siRNAs.     

Kinetics of rapid RNA evolution in vitro

Summary A rapidly acquired partial resistance to the replicase antagonist, ethidium bromide (EB), seen by Spiegelman and coresearchers in Q RNA variants competitively replicating under defined conditions in vitro, reflected existence of a pool of mutant RNA molecules, preadapted to EB, and their cross-propagation from the pre-EB optimum species, MDV-1, and from other kindred variants, some of which remained undetected, according to this quantitative analysis of midivariant RNA replication kinetics. DNAlike features of their evolution, such as the cloning of variants from an MDV-1 subtype and a complicance with the fundamental theorem of natural selection, resulted from the suppression, both real and apparent, of intrinsic RNA heterogeneity through sampling and detection methods, and also by the ascendency of self-propagation over cross-propagation with advancement of a superior variant. The deficit in mean polymer fitness, compared with optimum levels, determines the lower limit of this heterogeneity. Stability conditions for frequency equilibrium and strategies for counteracting viral drug resistance have been considered.     

A Thermodynamic Perspective on Natural Selection

A novel thermodynamic perspective on natural selection is presented. In the case that life continuity is optimized in an ideal system, where relatively constant and homogeneous selective pressures favour a given competing species, natural selection leads that system to a stationary state of maximum genotypic uniformity of life and maximum sustainable consumption of available energy by life (competitive equilibrium). Structurally and functionally, this optimizing tendency towards competitive equilibrium looks similar to the optimizing tendency towards thermodynamic equilibrium of classical thermodynamics (maximum energetic uniformity and maximum degradation of available energy). The principle of competitive exclusion may thus be conceptually viewed as an ecological manifestation of the second law of (classical) thermodynamics. On the other hand, the novel thermodynamic perspective on natural selection is discussed with regard to the open and nonequilibrium system of nature, where selective pressures vary in space and time. In this case, natural selection can induce diversity instead of uniformity, though an optimizing tendendcy towards maximum sustainable consumption of resources (optimization of life continuity) always remains. Overall, it is concluded that the action of natural selection favours the maximization of the sustainable consumption of energy at the level of individual organism.     

Evolutionary biotechnology – theory,facts and perspectives

Molecular evolution has recently been applied in biotechnology which consist of the development of evolutionary strategies in the design of biopolymers with predefined properties and functions. At the heart of this new technology are the in vitro replication and random synthesis of RNA or DNA molecules, producing large libraries of genotypes that are subjected to selection techniques following DARWIN's principle. By means of these evolutionary methods, RNA molecules were derived which specifically bind to predefined target molecules. Ribozymes with new catalytic functions were obtained as well as RNA molecules that are resistant to cleavage by specific RNases. In addition, the catalytic specificities of group I introns, a special class of ribozymes, were modified by variation and selection. Efficient applications of molecular evolution to problems in biotechnology require a fundamental and detailed understanding of the evolutionary process. Two basic questions are of primary importance: (i) How can evolutionary methods be successful as the numbers of possible genotypes are so large that the chance of obtaining a particular sequence by random processes is practically zero, and (ii) how can populations avoid being caught in evolutionary traps corresponding to local fitness optima? This review is therefore concerned with an abridged account of the theory of molecular evolution, as well as its application to biotechnology. We add a brief discussion of new techniques for the massively parallel handling and screening of very small probes as is required for the spatial separation and selection of genotypes. Finally, some imminent prospects concerning the evolutionary design of biopolymers are presented.     

Interaction of human serum albumin with 10-hydroxycamptothecin: spectroscopic and molecular modeling studies

In this work, fluorescence spectroscopy in combination with circular dichroism spectroscopy and molecular modeling was employed to investigate the binding of 10-hydroxycamptothecin (HCPT) to human serum albumin (HSA) under simulative physiological conditions. The experiment results showed that the fluorescence quenching of HSA by HCPT was a result of the formation of HCPT–HSA complex. The corresponding association constants (K _a) between HCPT and HSA at four different temperatures were determined according to the modified Stern–Volmer equation. The results of thermodynamic parameters ΔG, ΔH, and ΔS indicated that hydrogen bonds and van der Waals forces played major roles for HCPT–HSA association. Site marker competitive displacement experiment indicated that the binding of HCPT to HSA primarily took place in sub-domain IIA (site I). Molecular docking study further confirmed the binding mode and the binding site obtained by fluorescence and site marker competitive experiments. The conformational investigation showed that the presence of HCPT decreased the α-helical content of HSA and induced the slight unfolding of the polypeptides of protein, which confirmed some micro-environmental and conformational changes of HSA molecules.     

Product analysis of RNA generated de novo by Qβ replicase

Q_β replicase synthesizes self-replicating RNA in the absence of exogenous template after a certain lag time (“template-free synthesis”). The products of the template-free RNA synthesis have been investigated by gel electrophoresis and fingerprinting techniques. It has been found that a multitude of self-replicating RNA species appears in the early phases of reaction with variable lengths and sequences. Template-free synthesis in different samples under completely identical conditions yields RNA products with very different and unrelated fingerprints. The early products rapidly undergo an evolution process that alters the phenotypic properties of the self-replicating RNA, and leads to a concomitant increase of replication efficiency. Fingerprints and electrophoretic mobilities of the self-replicating RNA species are altered discontinuously during the evolution process. The evolution process ends with the selection of optimized self-replicating RNA species, whose phenotypes are conserved even after many serial transfers. Some optimized RNA species and midivariant RNA apparently have related sequences, since they contain many identical spots in their fingerprints. The properties of the RNA species produced by template-free synthesis match those of 6 S RNA found in Q_β-infected Escherichia coli cells. The results are in full agreement with the finding of Sumper & Luce (1975), who have presented evidence that Q_β replicase synthesized RNA de novo in the absence of exogenous template.