Spectral diffusion and the energy landscape of a protein

We present a novel type of spectral diffusion experiment in the millikelvin range to characterize the energy landscape of a protein as compared with that of a glass. We measure the time evolution of spectral holes for more than 300 hr after well-defined initial nonequilibrium conditions. We show that the model of noninteracting two-level systems can describe spectral diffusion in the glass, but fails for the protein. Our results further demonstrate that randomness in the energy landscape of a protein shows features of organization. There are “deep minimum” states separated by barriers, the heights of which we are able to estimate. The energy landscape of a glass is featureless by comparison.

Independent evolution of the prochlorophyte and green plant chlorophyll a/b light-harvesting proteins

The prochlorophytes are oxygenic prokaryotes differing from other cyanobacteria by the presence of a light-harvesting system containing both chlorophylls (Chls) a and b and by the absence of phycobilins. We demonstrate here that the Chl a/b binding proteins from all three known prochlorophyte genera are closely related to IsiA, a cyanobacterial Chl a-binding protein induced by iron starvation, and to CP43, a constitutively expressed Chl a antenna protein of photosystem II. The prochlorophyte Chl a/b protein (pcb) genes do not belong to the extended gene family encoding eukaryotic Chl a/b and Chl a/c light-harvesting proteins. Although higher plants and prochlorophytes share common pigment complements, their light-harvesting systems have evolved independently.

Directed evolution of antibody fragments with monovalent femtomolar antigen-binding affinity

Single-chain antibody mutants have been evolved in vitro with antigen-binding equilibrium dissociation constant Kd = 48 fM and slower dissociation kinetics (half-time > 5 days) than those for the streptavidin–biotin complex. These mutants possess the highest monovalent ligand-binding affinity yet reported for an engineered protein by over two orders of magnitude. Optimal kinetic screening of randomly mutagenized libraries of 105–107 yeast surface-displayed antibodies enabled a >1,000-fold decrease in the rate of dissociation after four cycles of affinity mutagenesis and screening. The consensus mutations are generally nonconservative by comparison with naturally occurring mouse Fv sequences and with residues that do not contact the fluorescein antigen in the wild-type complex. The existence of these mutants demonstrates that the antibody Fv architecture is not intrinsically responsible for an antigen-binding affinity ceiling during in vivo affinity maturation.

Adaptive evolution of a candidate gene for aging in Drosophila

Examination of the phenotypic effects of specific mutations has been extensively used to identify candidate genes affecting traits of interest. However, such analyses do not reveal anything about the evolutionary forces acting at these loci, or whether standing allelic variation contributes to phenotypic variance in natural populations. The Drosophila gene methuselah (mth) has been proposed as having major effects on organismal stress response and longevity phenotype. Here, we examine patterns of polymorphism and divergence at mth in population level samples of Drosophila melanogaster, D. simulans, and D. yakuba. Mth has experienced an unusually high level of adaptive amino acid divergence concentrated in the intra- and extracellular loop domains of the receptor protein, suggesting the historical action of positive selection on those regions of the molecule that modulate signal transduction. Further analysis of single nucleotide polymorphisms (SNPs) in D. melanogaster provided evidence for contemporary and spatially variable selection at the mth locus. In ten surveyed populations, the most common mth haplotype exhibited a 40% cline in frequency that coincided with population level differences in multiple life-history traits including lifespan. This clinal pattern was not associated with any particular SNP in the coding region, indicating that selection is operating at a closely linked site that may be involved in gene expression. Together, these consistently nonneutral patterns of inter- and intraspecific variation suggest adaptive evolution of a signal transduction pathway that may modulate lifespan in nature.

Purifying selection and birth-and-death evolution in the ubiquitin gene family

Ubiquitin is a highly conserved protein that is encoded by a multigene family. It is generally believed that this gene family is subject to concerted evolution, which homogenizes the member genes of the family. However, protein homogeneity can be attained also by strong purifying selection. We therefore studied the proportion (pS) of synonymous nucleotide differences between members of the ubiquitin gene family from 28 species of fungi, plants, and animals. The results have shown that pS is generally very high and is often close to the saturation level, although the protein sequence is virtually identical for all ubiquitins from fungi, plants, and animals. A small proportion of species showed a low level of pS values, but these values appeared to be caused by recent gene duplication. It was also found that the number of repeat copies of the gene family varies considerably with species, and some species harbor pseudogenes. These observations suggest that the members of this gene family evolve almost independently by silent nucleotide substitution and are subjected to birth-and-death evolution at the DNA level.

Novel folded protein domains generated by combinatorial shuffling of polypeptide segments

It has been proposed that the architecture of protein domains has evolved by the combinatorial assembly and/or exchange of smaller polypeptide segments. To investigate this proposal, we fused DNA encoding the N-terminal half of a β-barrel domain (from cold shock protein CspA) with fragmented genomic Escherichia coli DNA and cloned the repertoire of chimeric polypeptides for display on filamentous bacteriophage. Phage displaying folded polypeptides were selected by proteolysis; in most cases the protease-resistant chimeric polypeptides comprised genomic segments in their natural reading frames. Although the genomic segments appeared to have no sequence homologies with CspA, one of the originating proteins had the same fold as CspA, but another had a different fold. Four of the chimeric proteins were expressed as soluble polypeptides; they formed monomers and exhibited cooperative unfolding. Indeed, one of the chimeric proteins contained a set of very slowly exchanging amides and proved more stable than CspA itself. These results indicate that native-like proteins can be generated directly by combinatorial segment assembly from nonhomologous proteins, with implications for theories of the evolution of new protein folds, as well as providing a means of creating novel domains and architectures in vitro.

Evolution of an enzyme active site: The structure of a new crystal form of muconate lactonizing enzyme compared with mandelate racemase and enolase

Muconate lactonizing enzyme (MLE), a component of the β-ketoadipate pathway of Pseudomonas putida, is a member of a family of related enzymes (the “enolase superfamily”) that catalyze the abstraction of the α-proton of a carboxylic acid in the context of different overall reactions. New untwinned crystal forms of MLE were obtained, one of which diffracts to better than 2.0-Å resolution. The packing of the octameric enzyme in this crystal form is unusual, because the asymmetric unit contains three subunits. The structure of MLE presented here contains no bound metal ion, but is very similar to a recently determined Mn2+-bound structure. Thus, absence of the metal ion does not perturb the structure of the active site. The structures of enolase, mandelate racemase, and MLE were superimposed. A comparison of metal ligands suggests that enolase may retain some characteristics of the ancestor of this enzyme family. Comparison of other residues involved in catalysis indicates two unusual patterns of conservation: (i) that the position of catalytic atoms remains constant, although the residues that contain them are located at different points in the protein fold; and (ii) that the positions of catalytic residues in the protein scaffold are conserved, whereas their identities and roles in catalysis vary.

A method for directed evolution and functional cloning of enzymes

A general scheme is described for the in vitro evolution of protein catalysts in a biologically amplifiable system. Substrate is covalently and site specifically attached by a flexible tether to the pIII coat protein of a filamentous phage that also displays the catalyst. Intramolecular conversion of substrate to product provides a basis for selecting active catalysts from a library of mutants, either by release from or attachment to a solid support. This methodology has been developed with the enzyme staphylococcal nuclease as a model. An analysis of factors influencing the selection efficiency is presented, and it is shown that phage displaying staphylococcal nuclease can be enriched 100-fold in a single step from a library-like ensemble of phage displaying noncatalytic proteins. Additionally, this approach should allow one to functionally clone natural enzymes, based on their ability to catalyze specific reactions (e.g., glycosyl transfer, sequence-specific proteolysis or phosphorylation, polymerization, etc.) rather than their sequence- or structural homology to known enzymes.

Protein fluctuations are sensed by stimulated infrared echoes of the vibrations of carbon monoxide and azide probes

The correlation functions of the fluctuations of vibrational frequencies of azide ions and carbon monoxide in proteins are determined directly from stimulated photon echoes generated with femtosecond infrared pulses. The asymmetric stretching vibration of azide bound to carbonic anhydrase II exhibits a pronounced evolution of its vibrational frequency distribution on the time scale of a few picoseconds, which is attributed to modifications of the ligand structure through interactions with the nearby Thr-199. When azide is bound in hemoglobin, a more complex evolution of the protein structure is required to interchange the different ligand configurations, as evidenced by the much slower relaxation of the frequency distribution in this case. The time evolution of the distribution of frequencies of carbon monoxide bound in hemoglobin occurs on the ≈10-ps time scale and is very nonexponential. The correlation functions of the frequency fluctuations determine the evolution of the protein structure local to the probe and the extent to which the probe can navigate those parts of the energy landscape where the structural configurations are able to modify the local potential energy function of the probe.

Origin of a chloroplast protein importer

During evolution, chloroplasts have relinquished the majority of their genes to the nucleus. The products of transferred genes are imported into the organelle with the help of an import machinery that is distributed across the inner and outer plastid membranes. The evolutionary origin of this machinery is puzzling because, in the putative predecessors, the cyanobacteria, the outer two membranes, the plasma membrane, and the lipopolysaccharide layer lack a functionally similar protein import system. A 75-kDa protein-conducting channel in the outer envelope of pea chloroplasts, Toc75, shares ≈22% amino acid identity to a similarly sized protein, designated SynToc75, encoded in the Synechocystis PCC6803 genome. Here we show that SynToc75 is located in the outer membrane (lipopolysaccharide layer) of Synechocystis PCC6803 and that SynToc75 forms a voltage-gated, high conductance channel with a high affinity for polyamines and peptides in reconstituted liposomes. These findings suggest that a component of the chloroplast protein import system, Toc75, was recruited from a preexisting channel-forming protein of the cyanobacterial outer membrane. Furthermore, the presence of a protein in the chloroplastic outer envelope homologous to a cyanobacterial protein provides support for the prokaryotic nature of this chloroplastic membrane.

Molecular evolution of two vertebrate aryl hydrocarbon (dioxin) receptors (AHR1 and AHR2) and the PAS family

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor through which halogenated aromatic hydrocarbons such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) cause altered gene expression and toxicity. The AHR belongs to the basic helix–loop–helix/Per-ARNT-Sim (bHLH-PAS) family of transcriptional regulatory proteins, whose members play key roles in development, circadian rhythmicity, and environmental homeostasis; however, the normal cellular function of the AHR is not yet known. As part of a phylogenetic approach to understanding the function and evolutionary origin of the AHR, we sequenced the PAS homology domain of AHRs from several species of early vertebrates and performed phylogenetic analyses of these AHR amino acid sequences in relation to mammalian AHRs and 24 other members of the PAS family. AHR sequences were identified in a teleost (the killifish Fundulus heteroclitus), two elasmobranch species (the skate Raja erinacea and the dogfish Mustelus canis), and a jawless fish (the lamprey Petromyzon marinus). Two putative AHR genes, designated AHR1 and AHR2, were found both in Fundulus and Mustelus. Phylogenetic analyses indicate that the AHR2 genes in these two species are orthologous, suggesting that an AHR gene duplication occurred early in vertebrate evolution and that multiple AHR genes may be present in other vertebrates. Database searches and phylogenetic analyses identified four putative PAS proteins in the nematode Caenorhabditis elegans, including possible AHR and ARNT homologs. Phylogenetic analysis of the PAS gene family reveals distinct clades containing both invertebrate and vertebrate PAS family members; the latter include paralogous sequences that we propose have arisen by gene duplication early in vertebrate evolution. Overall, our analyses indicate that the AHR is a phylogenetically ancient protein present in all living vertebrate groups (with a possible invertebrate homolog), thus providing an evolutionary perspective to the study of dioxin toxicity and AHR function.

Rapid evolution of the DNA-binding site in LAGLIDADG homing endonucleases

Sequence analysis of chloroplast and mitochondrial large subunit rRNA genes from over 75 green algae disclosed 28 new group I intron-encoded proteins carrying a single LAGLIDADG motif. These putative homing endonucleases form four subfamilies of homologous enzymes, with the members of each subfamily being encoded by introns sharing the same insertion site. We showed that four divergent endonucleases from the I-CreI subfamily cleave the same DNA substrates. Mapping of the 66 amino acids that are conserved among the members of this subfamily on the 3-dimensional structure of I-CreI bound to its recognition sequence revealed that these residues participate in protein folding, homodimerization, DNA recognition and catalysis. Surprisingly, only seven of the 21 I-CreI amino acids interacting with DNA are conserved, suggesting that I-CreI and its homologs use different subsets of residues to recognize the same DNA sequence. Our sequence comparison of all 45 single-LAGLIDADG proteins identified so far suggests that these proteins share related structures and that there is a weak pressure in each subfamily to maintain identical protein–DNA contacts. The high sequence variability we observed in the DNA-binding site of homologous LAGLIDADG endonucleases provides insight into how these proteins evolve new DNA specificity.

Identification of a sequence element directing a protein to nuclear speckles

SF3b155 is an essential spliceosomal protein, highly conserved during evolution. It has been identified as a subunit of splicing factor SF3b, which, together with a second multimeric complex termed SF3a, interacts specifically with the 12S U2 snRNP and converts it into the active 17S form. The protein displays a characteristic intranuclear localization. It is diffusely distributed in the nucleoplasm but highly concentrated in defined intranuclear structures termed “speckles,” a subnuclear compartment enriched in small ribonucleoprotein particles and various splicing factors. The primary sequence of SF3b155 suggests a multidomain structure, different from those of other nuclear speckles components. To identify which part of SF3b155 determines its specific intranuclear localization, we have constructed expression vectors encoding a series of epitope-tagged SF3b155 deletion mutants as well as chimeric combinations of SF3b155 sequences with the soluble cytoplasmic protein pyruvate kinase. Following transfection of cultured mammalian cells, we have identified (i) a functional nuclear localization signal of the monopartite type (KRKRR, amino acids 196–200) and (ii) a molecular segment with multiple threonine-proline repeats (amino acids 208–513), which is essential and sufficient to confer a specific accumulation in nuclear speckles. This latter sequence element, in particular amino acids 208–440, is required for correct subcellular localization of SF3b155 and is also sufficient to target a reporter protein to nuclear speckles. Moreover, this “speckle-targeting sequence” transfers the capacity for interaction with other U2 snRNP components.

Systemic spread of an RNA insect virus in plants expressing plant viral movement protein genes

Flock house virus (FHV), a single-stranded RNA insect virus, has previously been reported to cross the kingdom barrier and replicate in barley protoplasts and in inoculated leaves of several plant species [Selling, B. H., Allison, R. F. & Kaesberg, P. (1990) Proc. Natl. Acad. Sci. USA 87, 434–438]. There was no systemic movement of FHV in plants. We tested the ability of movement proteins (MPs) of plant viruses to provide movement functions and cause systemic spread of FHV in plants. We compared the growth of FHV in leaves of nontransgenic and transgenic plants expressing the MP of tobacco mosaic virus or red clover necrotic mosaic virus (RCNMV). Both MPs mobilized cell-to-cell and systemic movement of FHV in Nicotiana benthamiana plants. The yield of FHV was more than 100-fold higher in the inoculated leaves of transgenic plants than in the inoculated leaves of nontransgenic plants. In addition, FHV accumulated in the noninoculated upper leaves of both MP-transgenic plants. RCNMV MP was more efficient in mobilizing FHV to noninoculated upper leaves. We also report here that FHV replicates in inoculated leaves of six additional plant species: alfalfa, Arabidopsis, Brassica, cucumber, maize, and rice. Our results demonstrate that plant viral MPs cause cell-to-cell and long-distance movement of an animal virus in plants and offer approaches to the study of the evolution of viruses and mechanisms governing mRNA trafficking in plants as well as to the development of promising vectors for transient expression of foreign genes in plants.

Dynamic evolution of plant mitochondrial genomes: Mobile genes and introns and highly variable mutation rates

We summarize our recent studies showing that angiosperm mitochondrial (mt) genomes have experienced remarkably high rates of gene loss and concomitant transfer to the nucleus and of intron acquisition by horizontal transfer. Moreover, we find substantial lineage-specific variation in rates of these structural mutations and also point mutations. These findings mostly arise from a Southern blot survey of gene and intron distribution in 281 diverse angiosperms. These blots reveal numerous losses of mt ribosomal protein genes but, with one exception, only rare loss of respiratory genes. Some lineages of angiosperms have kept all of their mt ribosomal protein genes whereas others have lost most of them. These many losses appear to reflect remarkably high (and variable) rates of functional transfer of mt ribosomal protein genes to the nucleus in angiosperms. The recent transfer of cox2 to the nucleus in legumes provides both an example of interorganellar gene transfer in action and a starting point for discussion of the roles of mechanistic and selective forces in determining the distribution of genetic labor between organellar and nuclear genomes. Plant mt genomes also acquire sequences by horizontal transfer. A striking example of this is a homing group I intron in the mt cox1 gene. This extraordinarily invasive mobile element has probably been acquired over 1,000 times separately during angiosperm evolution via a recent wave of cross-species horizontal transfers. Finally, whereas all previously examined angiosperm mtDNAs have low rates of synonymous substitutions, mtDNAs of two distantly related angiosperms have highly accelerated substitution rates.