Testing ancient RNA–protein interactions

The past decade in molecular biology has seen remarkable advances in the study of the origin and early evolution of life. The mathematical tools for analyzing DNA and protein sequences, coupled with the availability of complete microbial genome sequences, provide insight almost as far back as the age of the nucleic acids themselves. Experimental evolution in the laboratory and especially in vitro evolution of RNA provide insight into a hypothetical world where RNA, or a close relative, may have debuted as a primary functional and informational molecule. The ability to isolate new functional RNAs from random sequences now ultimately makes the world of possible primitive chemical interactions accessible even when the molecules or reactions are no longer present in modern species. Thus we can at last form direct experimental tests of specific models for the origin of RNA–protein associations, such as those that influenced the genetic code. This marks a turning point for probing the origin and early history of life at the molecular level.

Active barnase variants with completely random hydrophobic cores.

The central structural feature of natural proteins is a tightly packed and highly ordered hydrophobic core. If some measure of exquisite, native-like core packing is necessary for enzymatic function, this would constitute a significant obstacle to the development of novel enzymes, either by design or by natural or experimental evolution. To test the minimum requirements for a core to provide sufficient structural integrity for enzymatic activity, we have produced mutants of the ribonuclease barnase in which 12 of the 13 core residues have together been randomly replaced by hydrophobic alternatives. Using a sensitive biological screen, we find that a strikingly high proportion of these mutants (23%) retain enzymatic activity in vivo. Further substitution at the 13th core position shows that a similar proportion of completely random hydrophobic cores supports enzyme function. Of the active mutants produced, several have no wild-type core residues. These results imply that hydrophobicity is nearly a sufficient criterion for the construction of a functional core and, in conjunction with previous studies, that refinement of a crudely functional core entails more stringent sequence constraints than does the initial attainment of crude core function. Since attainment of crude function is the critical initial step in evolutionary innovation, the relatively scant requirements contributed by the hydrophobic core would greatly reduce the initial hurdle on the evolutionary pathway to novel enzymes. Similarly, experimental development of novel functional proteins might be simplified by limiting core design to mere specification of hydrophobicity and using iterative mutation-selection to optimize core structure.

Flower color variation: A model for the experimental study of evolution

We review the study of flower color polymorphisms in the morning glory as a model for the analysis of adaptation. The pathway involved in the determination of flower color phenotype is traced from the molecular and genetic levels to the phenotypic level. Many of the genes that determine the enzymatic components of flavonoid biosynthesis are redundant, but, despite this complexity, it is possible to associate discrete floral phenotypes with individual genes. An important finding is that almost all of the mutations that determine phenotypic differences are the result of transposon insertions. Thus, the flower color diversity seized on by early human domesticators of this plant is a consequence of the rich variety of mobile elements that reside in the morning glory genome. We then consider a long history of research aimed at uncovering the ecological fate of these various flower phenotypes in the southeastern U.S. A large body of work has shown that insect pollinators discriminate against white phenotypes when white flowers are rare in populations. Because the plant is self-compatible, pollinator bias causes an increase in self-fertilization in white maternal plants, which should lead to an increase in the frequency of white genes, according to modifier gene theory. Studies of geographical distributions indicate other, as yet undiscovered, disadvantages associated with the white phenotype. The ultimate goal of connecting ecology to molecular genetics through the medium of phenotype is yet to be attained, but this approach may represent a model for analyzing the translation between these two levels of biological organization.

Calibrating bacterial evolution

Attempts to calibrate bacterial evolution have relied on the assumption that rates of molecular sequence divergence in bacteria are similar to those of higher eukaryotes, or to those of the few bacterial taxa for which ancestors can be reliably dated from ecological or geological evidence. Despite similarities in the substitution rates estimated for some lineages, comparisons of the relative rates of evolution at different classes of nucleotide sites indicate no basis for their universal application to all bacteria. However, there is evidence that bacteria have a constant genome-wide mutation rate on an evolutionary time scale but that this rate differs dramatically from the rate estimated by experimental methods.

Stochastic processes strongly influence HIV-1 evolution during suboptimal protease-inhibitor therapy

It has long been assumed that HIV-1 evolution is best described by deterministic evolutionary models because of the large population size. Recently, however, it was suggested that the effective population size (Ne) may be rather small, thereby allowing chance to influence evolution, a situation best described by a stochastic evolutionary model. To gain experimental evidence supporting one of the evolutionary models, we investigated whether the development of resistance to the protease inhibitor ritonavir affected the evolution of the env gene. Sequential serum samples from five patients treated with ritonavir were used for analysis of the protease gene and the V3 domain of the env gene. Multiple reverse transcription–PCR products were cloned, sequenced, and used to construct phylogenetic trees and to calculate the genetic variation and Ne. Genotypic resistance to ritonavir developed in all five patients, but each patient displayed a unique combination of mutations, indicating a stochastic element in the development of ritonavir resistance. Furthermore, development of resistance induced clear bottleneck effects in the env gene. The mean intrasample genetic variation, which ranged from 1.2% to 5.7% before treatment, decreased significantly (P < 0.025) during treatment. In agreement with these findings, Ne was estimated to be very small (500–15,000) compared with the total HIV-1 RNA copy number. This study combines three independent observations, strong population bottlenecking, small Ne, and selection of different combinations of protease-resistance mutations, all of which indicate that HIV-1 evolution is best described by a stochastic evolutionary model.

Experimental surgery to create subgenomes of Bacillus subtilis 168

The 4,188-kb circular genome of Bacillus subtilis 168 was artificially dissected into two stable circular chromosomes in vivo, one being the 3,878-kb main genome and the other the 310-kb subgenome that was recovered as covalently closed circular DNA in CsCl-ethidium bromide ultracentrifugation. The minimal requirements to physically separate the 310-kb DNA segment out of the genome were two interrepeat homologous sequences and an origin of DNA replication between them. The subgenome originated from the 1,255–1,551-kb region of the B. subtilis genome was essential for the cell to survive because the subgenome was not lost from the cell. The finding that the B. subtilis genome has a potential to be divided and the resulting two replicons stably maintained may shed light on origins and formation mechanisms of giant plasmids or second chromosomes present in many bacteria. Similar excision or its reversal process, i.e., integration of large sized covalently closed circular DNA pieces into the main genome, implies significant roles of subgenomes in the exchange of genetic information and size variation of bacterial genomes in bacterial evolution.

Dynamic regulation of the tryptophan operon: A modeling study and comparison with experimental data

A mathematical model for regulation of the tryptophan operon is presented. This model takes into account repression, feedback enzyme inhibition, and transcriptional attenuation. Special attention is given to model parameter estimation based on experimental data. The model's system of delay differential equations is numerically solved, and the results are compared with experimental data on the temporal evolution of enzyme activity in cultures of Escherichia coli after a nutritional shift (minimal + tryptophan medium to minimal medium). Good agreement is obtained between the numeric simulations and the experimental results for wild-type E. coli, as well as for two different mutant strains.

Computational method to reduce the search space for directed protein evolution

We introduce a computational method to optimize the in vitro evolution of proteins. Simulating evolution with a simple model that statistically describes the fitness landscape, we find that beneficial mutations tend to occur at amino acid positions that are tolerant to substitutions, in the limit of small libraries and low mutation rates. We transform this observation into a design strategy by applying mean-field theory to a structure-based computational model to calculate each residue's structural tolerance. Thermostabilizing and activity-increasing mutations accumulated during the experimental directed evolution of subtilisin E and T4 lysozyme are strongly directed to sites identified by using this computational approach. This method can be used to predict positions where mutations are likely to lead to improvement of specific protein properties.

The evolution of RNA viruses: A population genetics view

RNA viruses are excellent experimental models for studying evolution under the theoretical framework of population genetics. For a proper justification of this thesis we have introduced some properties of RNA viruses that are relevant for studying evolution. On the other hand, population genetics is a reductionistic theory of evolution. It does not consider or make simplistic assumptions on the transformation laws within and between genotypic and phenotypic spaces. However, such laws are minimized in the case of RNA viruses because the phenotypic space maps onto the genotypic space in a much more linear way than on higher DNA-based organisms. Under experimental conditions, we have tested the role of deleterious and beneficial mutations in the degree of adaptation of vesicular stomatitis virus (VSV), a nonsegmented virus of negative strand. We also have studied how effective population size, initial genetic variability in populations, and environmental heterogeneity shapes the impact of mutations in the evolution of vesicular stomatitis virus. Finally, in an integrative attempt, we discuss pros and cons of the quasispecies theory compared with classic population genetics models for haploid organisms to explain the evolution of RNA viruses.

Computer simulations of enzyme catalysis: Finding out what has been optimized by evolution

The origin of the catalytic power of enzymes is discussed, paying attention to evolutionary constraints. It is pointed out that enzyme catalysis reflects energy contributions that cannot be determined uniquely by current experimental approaches without augmenting the analysis by computer simulation studies. The use of energy considerations and computer simulations allows one to exclude many of the popular proposals for the way enzymes work. It appears that the standard approaches used by organic chemists to catalyze reactions in solutions are not used by enzymes. This point is illustrated by considering the desolvation hypothesis and showing that it cannot account for a large increase in kcat relative to the corresponding kcage for the reference reaction in a solvent cage. The problems associated with other frequently invoked mechanisms also are outlined. Furthermore, it is pointed out that mutation studies are inconsistent with ground state destabilization mechanisms. After considering factors that were not optimized by evolution, we review computer simulation studies that reproduced the overall catalytic effect of different enzymes. These studies pointed toward electrostatic effects as the most important catalytic contributions. The nature of this electrostatic stabilization mechanism is far from being obvious because the electrostatic interaction between the reacting system and the surrounding area is similar in enzymes and in solution. However, the difference is that enzymes have a preorganized dipolar environment that does not have to pay the reorganization energy for stabilizing the relevant transition states. Apparently, the catalytic power of enzymes is stored in their folding energy in the form of the preorganized polar environment.

Blockade of T-cell costimulation prevents development of experimental chronic renal allograft rejection.

Blocking CD28-B7 T-cell costimulation by systemic administration of CTLA4Ig, a fusion protein which binds B7 molecules on the surface of antigen-presenting cells, prevents rejection and induces tolerance in experimental acute allograft rejection models. We tested the effect of CTLA4Ig therapy on the process of chronic renal allograft rejection using an established experimental transplantation model. F344 kidneys were transplanted orthotopically into bilaterally nephrectomized LEW recipients. Control animals received low dose cyclosporine for 10 days posttransplantation. Administration of a single injection of CTLA4Ig on day 2 posttransplant alone or in addition to the low dose cyclosporine protocol resulted in improvement of long-term graft survival as compared with controls. More importantly, control recipients which received cyclosporine only developed progressive proteinuria by 8-12 weeks, and morphological evidence of chronic rejection by 16-24 weeks, including widespread transplant arteriosclerosis and focal and segmental glomerulosclerosis, while animals treated with CTLA4Ig alone or in addition to cyclosporine did not. Competitive reverse transcriptase-PCR and immunohistological analysis of allografts at 8, 16, and 24 weeks showed attenuation of lymphocyte and macrophage infiltration and activation in the CTLA4Ig-treated animals, as compared with cyclosporine-alone treated controls. These data confirm that early blockade of the CD28-B7 T-cell costimulatory pathway prevents later development and evolution of chronic renal allograft rejection. Our results indicate that T-cell recognition of alloantigen is a central event in initiating the process of chronic rejection, and that strategies targeted at blocking T-cell costimulation may prove to be a valuable clinical approach to preventing development of the process.

Origin and evolution of circular waves and spirals in Dictyostelium discoideum territories.

Randomly distributed Dictyostelium discoideum cells form cooperative territories by signaling to each other with cAMP. Cells initiate the process by sending out pulsatile signals, which propagate as waves. With time, circular and spiral patterns form. We show that by adding spatial and temporal noise to the levels of an important regulator of external cAMP levels, the cAMP phosphodiesterase inhibitor, we can explain the natural progression of the system from randomly firing cells to circular waves whose symmetries break to form double- and single- or multi-armed spirals. When phosphodiesterase inhibitor is increased with time, mimicking experimental data, the wavelength of the spirals shortens, and a proportion of them evolve into pairs of connected spirals. We compare these results to recent experiments, finding that the temporal and spatial correspondence between experiment and model is very close.

Microevolutionary responses in experimental populations of plants to CO2-enriched environments: parallel results from two model systems.

Despite the critical role that terrestrial vegetation plays in the Earth's carbon cycle, very little is known about the potential evolutionary responses of plants to anthropogenically induced increases in concentrations of atmospheric CO2. We present experimental evidence that rising CO2 concentration may have a direct impact on the genetic composition and diversity of plant populations but is unlikely to result in selection favoring genotypes that exhibit increased productivity in a CO2-enriched atmosphere. Experimental populations of an annual plant (Abutilon theophrasti, velvetleaf) and a temperate forest tree (Betula alleghaniensis, yellow birch) displayed responses to increased CO2 that were both strongly density-dependent and genotype-specific. In competitive stands, a higher concentration of CO2 resulted in pronounced shifts in genetic composition, even though overall CO2-induced productivity enhancements were small. For the annual species, quantitative estimates of response to selection under competition were 3 times higher at the elevated CO2 level. However, genotypes that displayed the highest growth responses to CO2 when grown in the absence of competition did not have the highest fitness in competitive stands. We suggest that increased CO2 intensified interplant competition and that selection favored genotypes with a greater ability to compete for resources other than CO2. Thus, while increased CO2 may enhance rates of selection in populations of competing plants, it is unlikely to result in the evolution of increased CO2 responsiveness or to operate as an important feedback in the global carbon cycle. However, the increased intensity of selection and drift driven by rising CO2 levels may have an impact on the genetic diversity in plant populations.

Founder-effect speciation theory: failure of experimental corroboration.

The theory of founder-effect speciation proposes that colonization by very few individuals of an empty habitat favors rapid genetic changes and the evolution of a new species. We report here the results obtained in a 10-year-long and large-scale experiment with Drosophila pseudoobscura designed to test the theory. In our experimental protocol, populations are established with variable numbers of very few individuals and allowed to expand greatly for several generations until conditions of severe competition for resources are reached and the population crashes. A few random survivors are then taken to start a new population expansion and thus initiate a new cycle of founding events, population flushes, and crashes. Our results provide no support for the theories proposing that new species are very likely to appear as by-products of founder events.