Genetic definition of a protein-splicing domain: Functional mini-inteins support structure predictions and a model for intein evolution

Inteins are protein-splicing elements, most of which contain conserved sequence blocks that define a family of homing endonucleases. Like group I introns that encode such endonucleases, inteins are mobile genetic elements. Recent crystallography and computer modeling studies suggest that inteins consist of two structural domains that correspond to the endonuclease and the protein-splicing elements. To determine whether the bipartite structure of inteins is mirrored by the functional independence of the protein-splicing domain, the entire endonuclease component was deleted from the Mycobacterium tuberculosis recA intein. Guided by computer modeling studies, and taking advantage of genetic systems designed to monitor intein function, the 440-aa Mtu recA intein was reduced to a functional mini-intein of 137 aa. The accuracy of splicing of several mini-inteins was verified. This work not only substantiates structure predictions for intein function but also supports the hypothesis that, like group I introns, mobile inteins arose by an endonuclease gene invading a sequence encoding a small, functional splicing element.

On roads not taken in the evolution of protein catalysts: Antibody steroid isomerases that use an enamine mechanism

Reactive immunization has emerged as a new tool for the study of biological catalysis. A powerful application resulted in catalytic antibodies that use an enamine mechanism akin to that used by the class I aldolases. With regard to the evolution of enzyme mechanisms, we investigated the utility of an enamine pathway for the allylic rearrangement exemplified by Δ5-3-ketosteroid isomerase (KSI; EC 5.3.3.1). Our aldolase antibodies were found to catalyze the isomerization of both steroid model compounds and steroids. The kinetic and chemical studies showed that the antibodies afforded rate accelerations up to a factor of 104 by means of an enamine mechanism in which imine formation was the rate-determining step. In light of our observations and the enzyme studies by other workers, we suggest that an enamine pathway could have been an early, viable KSI mechanism. Although this pathway is amenable to optimization for increased catalytic power, it appears that certain factors precluded its evolution in known KSI enzymes.

Directed evolution of a thermostable esterase

We have used in vitro evolution to probe the relationship between stability and activity in a mesophilic esterase. Previous studies of these properties in homologous enzymes evolved for function at different temperatures have suggested that stability at high temperatures is incompatible with high catalytic activity at low temperatures through mutually exclusive demands on enzyme flexibility. Six generations of random mutagenesis, recombination, and screening stabilized Bacillus subtilis p-nitrobenzyl esterase significantly (>14°C increase in Tm) without compromising its catalytic activity at lower temperatures. Furthermore, analysis of the stabilities and activities of large numbers of random mutants indicates that these properties are not inversely correlated. Although enhanced thermostability does not necessarily come at the cost of activity, the process by which the molecule adapts is important. Mutations that increase thermostability while maintaining low-temperature activity are very rare. Unless both properties are constrained (by natural selection or screening) the evolution of one by the accumulation of single amino acid substitutions typically comes at the cost of the other, regardless of whether the two properties are inversely correlated or not correlated at all.

In vitro evolution of horse heart myoglobin to increase peroxidase activity

Random mutagenesis and screening for enzymatic activity has been used to engineer horse heart myoglobin to enhance its intrinsic peroxidase activity. A chemically synthesized gene encoding horse heart myoglobin was subjected to successive cycles of PCR random mutagenesis. The mutated myoglobin gene was expressed in Escherichia coli LE392, and the variants were screened for peroxidase activity with a plate assay. Four cycles of mutagenesis and screening produced a series of single, double, triple, and quadruple variants with enhanced peroxidase activity. Steady-state kinetics analysis demonstrated that the quadruple variant T39I/K45D/F46L/I107F exhibits peroxidase activity significantly greater than that of the wild-type protein with k1 (for H2O2 oxidation of metmyoglobin) of 1.34 × 104 M−1 s−1 (≈25-fold that of wild-type myoglobin) and k3 [for reducing the substrate (2, 2′-azino-di-(3-ethyl)benzthiazoline-6-sulfonic acid] of 1.4 × 106 M−1 s−1 (1.6-fold that of wild-type myoglobin). Thermal stability of these variants as measured with circular dichroism spectroscopy demonstrated that the Tm of the quadruple variant is decreased only slightly compared with wild-type (74.1°C vs. 76.5°C). The rate constants for binding of dioxygen exhibited by the quadruple variant are identical to the those observed for wild-type myoglobin (kon, 22.2 × 10−6 M−1 s−1 vs. 22.3 × 10−6 M−1 s−1; koff, 24.3 s−1 vs. 24.2 s−1; KO2, 0.91 × 10−6 M−1 vs. 0.92 × 10−6 M−1). The affinity of the quadruple variant for CO is increased slightly (kon, 0.90 × 10−6 M−1s−1 vs. 0.51 × 10−6 M−1s−1; koff, 5.08 s−1 vs. 3.51 s−1; KCO, 1.77 × 10−7 M−1 vs. 1.45 × 10−7 M−1). All four substitutions are in the heme pocket and within 5 Å of the heme group.

Thermus thermophilus: A link in evolution of the tRNA-dependent amino acid amidation pathways

Thermus thermophilus possesses an aspartyl-tRNA synthetase (AspRS2) able to aspartylate efficiently tRNAAsp and tRNAAsn. Aspartate mischarged on tRNAAsn then is converted into asparagine by an ω amidase that differs structurally from all known asparagine synthetases. However, aspartate is not misincorporated into proteins because the binding capacity of aminoacylated tRNAAsn to elongation factor Tu is only conferred by conversion of aspartate into asparagine. T. thermophilus additionally contains a second aspartyl-tRNA synthetase (AspRS1) able to aspartylate tRNAAsp and an asparaginyl-tRNA synthetase able to charge tRNAAsn with free asparagine, although the organism does not contain a tRNA-independent asparagine synthetase. In contrast to the duplicated pathway of tRNA asparaginylation, tRNA glutaminylation occurs in the thermophile via the usual pathway by using glutaminyl-tRNA synthetase and free glutamine synthesized by glutamine synthetase that is unique. T. thermophilus is able to ensure tRNA aminoacylation by alternative routes involving either the direct pathway or by conversion of amino acid mischarged on tRNA. These findings shed light on the interrelation between the tRNA-dependent and tRNA-independent pathways of amino acid amidation and on the processes involved in fidelity of the aminoacylation systems.

Molding a peptide into an RNA site by in vivo peptide evolution

Short peptides corresponding to the arginine-rich domains of several RNA-binding proteins are able to bind to their specific RNA sites with high affinities and specificities. In the case of the HIV-1 Rev-Rev response element (RRE) complex, the peptide forms a single α-helix that binds deeply in a widened, distorted RNA major groove and makes a substantial set of base-specific and backbone contacts. Using a reporter system based on antitermination by the bacteriophage λ N protein, it has been possible to identify novel arginine-rich peptides from combinatorial libraries that recognize the RRE with affinities and specificities similar to Rev but that appear to bind in nonhelical conformations. Here we have used codon-based mutagenesis to evolve one of these peptides, RSG-1, into an even tighter binder. After two rounds of evolution, RSG-1.2 bound the RRE with 7-fold higher affinity and 15-fold higher specificity than the wild-type Rev peptide, and in vitro competition experiments show that RSG-1.2 completely displaces the intact Rev protein from the RRE at low peptide concentrations. By fusing RRE-binding peptides to the activation domain of HIV-1 Tat, we show that the peptides can deliver Tat to the RRE site and activate transcription in mammalian cells, and more importantly, that the fusion proteins can inhibit the activity of Rev in chloramphenicol acetyltransferase reporter assays. The evolved peptides contain proline and glutamic acid mutations near the middle of their sequences and, despite the presence of a proline, show partial α-helix formation in the absence of RNA. These directed evolution experiments illustrate how readily complex peptide structures can be evolved within the context of an RNA framework, perhaps reflecting how early protein structures evolved in an “RNA world.”

Convergent evolution of apolipoprotein(a) in primates and hedgehog

Apolipoprotein(a) [apo(a)] is the distinguishing protein component of lipoprotein(a), a major inherited risk factor for atherosclerosis. Human apo(a) is homologous to plasminogen. It contains from 15 to 50 repeated domains closely related to plasminogen kringle four, plus single kringle five-like and inactive protease-like domains. This expressed gene is confined to a subset of primates. Although most mammals lack apo(a), hedgehogs produce an apo(a)-like protein composed of highly repeated copies of a plasminogen kringle three-like domain, with complete absence of protease domain sequences. Both human and hedgehog apo(a)-like proteins form covalently linked lipoprotein particles that can bind to fibrin and other substrates shared with plasminogen. DNA sequence comparisons and phylogenetic analysis indicate that the human type of apo(a) evolved from a duplicated plasminogen gene during recent primate evolution. In contrast, the kringle three-based type of apo(a) evolved from an independent duplication of the plasminogen gene approximately 80 million years ago. In a type of convergent evolution, the plasminogen gene has been independently remodeled twice during mammalian evolution to produce similar forms of apo(a) in two widely divergent groups of species.

Evolution of gilled mushrooms and puffballs inferred from ribosomal DNA sequences

Homobasidiomycete fungi display many complex fruiting body morphologies, including mushrooms and puffballs, but their anatomical simplicity has confounded efforts to understand the evolution of these forms. We performed a comprehensive phylogenetic analysis of homobasidiomycetes, using sequences from nuclear and mitochondrial ribosomal DNA, with an emphasis on understanding evolutionary relationships of gilled mushrooms and puffballs. Parsimony-based optimization of character states on our phylogenetic trees suggested that strikingly similar gilled mushrooms evolved at least six times, from morphologically diverse precursors. Approximately 87% of gilled mushrooms are in a single lineage, which we call the “euagarics.” Recently discovered 90 million-year-old fossil mushrooms are probably euagarics, suggesting that (i) the origin of this clade must have occurred no later than the mid-Cretaceous and (ii) the gilled mushroom morphology has been maintained in certain lineages for tens of millions of years. Puffballs and other forms with enclosed spore-bearing structures (Gasteromycetes) evolved at least four times. Derivation of Gasteromycetes from forms with exposed spore-bearing structures (Hymenomycetes) is correlated with repeated loss of forcible spore discharge (ballistospory). Diverse fruiting body forms and spore dispersal mechanisms have evolved among Gasteromycetes. Nevertheless, it appears that Hymenomycetes have never been secondarily derived from Gasteromycetes, which suggests that the loss of ballistospory has constrained evolution in these lineages.

Origin and evolution of the slime molds (Mycetozoa)

The Mycetozoa include the cellular (dictyostelid), acellular (myxogastrid), and protostelid slime molds. However, available molecular data are in disagreement on both the monophyly and phylogenetic position of the group. Ribosomal RNA trees show the myxogastrid and dictyostelid slime molds as unrelated early branching lineages, but actin and β-tubulin trees place them together as a single coherent (monophyletic) group, closely related to the animal–fungal clade. We have sequenced the elongation factor-1α genes from one member of each division of the Mycetozoa, including Dictyostelium discoideum, for which cDNA sequences were previously available. Phylogenetic analyses of these sequences strongly support a monophyletic Mycetozoa, with the myxogastrid and dictyostelid slime molds most closely related to each other. All phylogenetic methods used also place this coherent Mycetozoan assemblage as emerging among the multicellular eukaryotes, tentatively supported as more closely related to animals + fungi than are green plants. With our data there are now three proteins that consistently support a monophyletic Mycetozoa and at least four that place these taxa within the “crown” of the eukaryote tree. We suggest that ribosomal RNA data should be more closely examined with regard to these questions, and we emphasize the importance of developing multiple sequence data sets.

Evolution of plastid gene rps2 in a lineage of hemiparasitic and holoparasitic plants: Many losses of photosynthesis and complex patterns of rate variation

The plastid genomes of some nonphotosynthetic parasitic plants have experienced an extreme reduction in gene content and an increase in evolutionary rate of remaining genes. Nothing is known of the dynamics of these events or whether either is a direct outcome of the loss of photosynthesis. The parasitic Scrophulariaceae and Orobanchaceae, representing a continuum of heterotrophic ability ranging from photosynthetic hemiparasites to nonphotosynthetic holoparasites, are used to investigate these issues. We present a phylogenetic hypothesis for parasitic Scrophulariaceae and Orobanchaceae based on sequences of the plastid gene rps2, encoding the S2 subunit of the plastid ribosome. Parasitic Scrophulariaceae and Orobanchaceae form a monophyletic group in which parasitism can be inferred to have evolved once. Holoparasitism has evolved independently at least five times, with certain holoparasitic lineages representing single species, genera, and collections of nonphotosynthetic genera. Evolutionary loss of the photosynthetic gene rbcL is limited to a subset of holoparasitic lineages, with several holoparasites retaining a full length rbcL sequence. In contrast, the translational gene rps2 is retained in all plants investigated but has experienced rate accelerations in several hemi- as well as holoparasitic lineages, suggesting that there may be substantial molecular evolutionary changes to the plastid genome of parasites before the loss of photosynthesis. Independent patterns of synonymous and nonsynonymous rate acceleration in rps2 point to distinct mechanisms underlying rate variation in different lineages. Parasitic Scrophulariaceae (including the traditional Orobanchaceae) provide a rich platform for the investigation of molecular evolutionary process, gene function, and the evolution of parasitism.

Evolution of foraging behavior in Drosophila by density-dependent selection

One of the rare examples of a single major gene underlying a naturally occurring behavioral polymorphism is the foraging locus of Drosophila melanogaster. Larvae with the rover allele, forR, have significantly longer foraging path lengths on a yeast paste than do those homozygous for the sitter allele, fors. These variants do not differ in general activity in the absence of food. The evolutionary significance of this polymorphism is not as yet understood. Here we examine the effect of high and low animal rearing densities on the larval foraging path-length phenotype and show that density-dependent natural selection produces changes in this trait. In three unrelated base populations the long path (rover) phenotype was selected for under high-density rearing conditions, whereas the short path (sitter) phenotype was selected for under low-density conditions. Genetic crosses suggested that these changes resulted from alterations in the frequency of the fors allele in the low-density-selected lines. Further experiments showed that density-dependent selection during the larval stage rather than the adult stage of development was sufficient to explain these results. Density-dependent mechanisms may be sufficient to maintain variation in rover and sitter behavior in laboratory populations.

Evolution of the Friedreich’s ataxia trinucleotide repeat expansion: Founder effect and premutations

Friedreich’s ataxia, the most frequent inherited ataxia, is caused, in the vast majority of cases, by large GAA repeat expansions in the first intron of the frataxin gene. The normal sequence corresponds to a moderately polymorphic trinucleotide repeat with bimodal size distribution. Small normal alleles have approximately eight to nine repeats whereas a more heterogeneous mode of large normal alleles ranges from 16 to 34 GAA. The latter class accounts for ≈17% of normal alleles. To identify the origin of the expansion mutation, we analyzed linkage disequilibrium between expansion mutations or normal alleles and a haplotype of five polymorphic markers within or close to the frataxin gene; 51% of the expansions were associated with a single haplotype, and the other expansions were associated with haplotypes that could be related to the major one by mutation at a polymorphic marker or by ancient recombination. Of interest, the major haplotype associated with expansion is also the major haplotype associated with the larger alleles in the normal size range and was almost never found associated with the smaller normal alleles. The results indicate that most if not all large normal alleles derive from a single founder chromosome and that they represent a reservoir for larger expansion events, possibly through “premutation” intermediates. Indeed, we found two such alleles (42 and 60 GAA) that underwent cataclysmic expansion to pathological range in a single generation. This stepwise evolution to large trinucleotide expansions already was suggested for myotonic dystrophy and fragile X syndrome and may relate to a common mutational mechanism, despite sequence motif differences.

Striking sequence similarity in inter- and intra-specific comparisons of class I SLG alleles from Brassica oleracea and Brassica campestris: Implications for the evolution and recognition mechanism

Self-incompatibility in Brassica is controlled by a single multi-allelic locus (S locus), which contains at least two highly polymorphic genes expressed in the stigma: an S glycoprotein gene (SLG) and an S receptor kinase gene (SRK). The putative ligand-binding domain of SRK exhibits high homology to the secretory protein SLG, and it is believed that SLG and SRK form an active receptor kinase complex with a self-pollen ligand, which leads to the rejection of self-pollen. Here, we report 31 novel SLG sequences of Brassica oleracea and Brassica campestris. Sequence comparisons of a large number of SLG alleles and SLG-related genes revealed the following points. (i) The striking sequence similarity observed in an inter-specific comparison (95.6% identity between SLG14 of B. oleracea and SLG25 of B. campestris in deduced amino acid sequence) suggests that SLG diversification predates speciation. (ii) A perfect match of the sequences in hypervariable regions, which are thought to determine S specificity in an intra-specific comparison (SLG8 and SLG46 of B. campestris) and the observation that the hypervariable regions of SLG and SRK of the same S haplotype were not necessarily highly similar suggests that SLG and SRK bind different sites of the pollen ligand and that they together determine S specificity. (iii) Comparison of the hypervariable regions of SLG alleles suggests that intragenic recombination, together with point mutations, has contributed to the generation of the high level of sequence variation in SLG alleles. Models for the evolution of SLG/SRK are presented.

Ecological factors rather than temporal factors dominate the evolution of vesicular stomatitis virus

Vesicular stomatitis New Jersey virus (VSV-NJ) is a rhabdovirus that causes economically important disease in cattle and other domestic animals in endemic areas from southeastern United States to northern South America. Its negatively stranded RNA genome is capable of undergoing rapid evolution, which allows phylogenetic analysis and molecular epidemiology studies to be performed. Previous epidemiological studies in Costa Rica showed the existence of at least two distinct ecological zones of high VSV-NJ activity, one located in the highlands (premontane tropical moist forest) and the other in the lowlands (tropical dry forest). We wanted to test the hypothesis that the viruses circulating in these ecological zones were genetically distinct. For this purpose, we sequenced the hypervariable region of the phosphoprotein gene for 50 VSV-NJ isolates from these areas. Phylogenetic analysis showed that viruses from each ecological zone had distinct genotypes. These genotypes were maintained in each area for periods of up to 8 years. This evolutionary pattern of VSV-NJ suggests an adaptation to ecological factors that could exert selective pressure on the virus. As previous data indicated an absence of virus adaptation to factors related to the bovine host (including immunological pressure), it appears that VSV genetic divergence represents positive selection to adapt to specific vectors and/or reservoirs at each ecological zone.

The carbamate kinase-like carbamoyl phosphate synthetase of the hyperthermophilic archaeon Pyrococcus furiosus, a missing link in the evolution of carbamoyl phosphate biosynthesis

Microbial carbamoyl phosphate synthetases (CPS) use glutamine as nitrogen donor and are composed of two subunits (or domains), one exhibiting glutaminase activity, the other able to synthesize carbamoyl phosphate (CP) from bicarbonate, ATP, and ammonia. The pseudodimeric organization of this synthetase suggested that it has evolved by duplication of a smaller kinase, possibly a carbamate kinase (CK). In contrast to other prokaryotes the hyperthermophilic archaeon Pyrococcus furiosus was found to synthesize CP by using ammonia and not glutamine. We have purified the cognate enzyme and found it to be a dimer of two identical subunits of Mr 32,000. Its thermostability is considerable, 50% activity being retained after 1 h at 100°C or 3 h at 95°C. The corresponding gene was cloned by PCR and found to present about 50% amino acid identity with known CKs. The stoichiometry of the reaction (two ATP consumed per CP synthesized) and the ability of the enzyme to catalyze at high rate a bicarbonate-dependent ATPase reaction however clearly distinguish P. furiosus CPS from ordinary CKs. Thus the CPS of P. furiosus could represent a primeval step in the evolution of CPS from CK. Our results suggest that the first event in this evolution was the emergence of a primeval synthetase composed of subunits able to synthesize both carboxyphosphate and CP; this step would have preceded the duplication assumed to have generated the two subdomains of modern CPSs. The gene coding for this CK-like CPS was called cpkA.