PubMed Health⌕ Search

Biomedical subjects

Patricia C Babbitt

Publications and source records attributed to Patricia C Babbitt.

At least 19 recordsLinked to original sources

Mechanisms of protein evolution and their application to protein engineering.

Protein engineering holds great promise for the development of new biosensors, diagnostics, therapeutics, and agents for bioremediation. Despite some remarkable successes in experimental and computational protein design, engineered proteins rarely achieve the efficiency or specificity of natural enzymes. Current protein design methods utilize evolutionary concepts, including mutation, recombination, and selection, but the inability to fully recapitulate the success of natural evolution suggests that some evolutionary principles have not been fully exploited. One aspect of protein engineering that has received little attention is how to select the most promising proteins to serve as templates, or scaffolds, for engineering. Two evolutionary concepts that could provide a rational basis for template selection are the conservation of catalytic mechanisms and functional promiscuity. Knowledge of the catalytic motifs responsible for conserved aspects of catalysis in mechanistically diverse superfamilies could be used to identify promising templates for protein engineering. Second, protein evolution often proceeds through promiscuous intermediates, suggesting that templates which are naturally promiscuous for a target reaction could enhance protein engineering strategies. This review explores these ideas and alternative hypotheses concerning protein evolution and engineering. Future research will determine if application of these principles will lead to a protein engineering methodology governed by predictable rules for designing efficient, novel catalysts.

Binding Sites↗

Comparison of methods for genomic localization of gene trap sequences.

BACKGROUND: Gene knockouts in a model organism such as mouse provide a valuable resource for the study of basic biology and human disease. Determining which gene has been inactivated by an untargeted gene trapping event poses a challenging annotation problem because gene trap sequence tags, which represent sequence near the vector insertion site of a trapped gene, are typically short and often contain unresolved residues. To understand better the localization of these sequences on the mouse genome, we compared stand-alone versions of the alignment programs BLAT, SSAHA, and MegaBLAST. A set of 3,369 sequence tags was aligned to build 34 of the mouse genome using default parameters for each algorithm. Known genome coordinates for the cognate set of full-length genes (1,659 sequences) were used to evaluate localization results. RESULTS: In general, all three programs performed well in terms of localizing sequences to a general region of the genome, with only relatively subtle errors identified for a small proportion of the sequence tags. However, large differences in performance were noted with regard to correctly identifying exon boundaries. BLAT correctly identified the vast majority of exon boundaries, while SSAHA and MegaBLAST missed the majority of exon boundaries. SSAHA consistently reported the fewest false positives and is the fastest algorithm. MegaBLAST was comparable to BLAT in speed, but was the most susceptible to localizing sequence tags incorrectly to pseudogenes. CONCLUSION: The differences in performance for sequence tags and full-length reference sequences were surprisingly small. Characteristic variations in localization results for each program were noted that affect the localization of sequence at exon boundaries, in particular.

Algorithms↗

Evolution of enzyme superfamilies.

Enzyme evolution is often constrained by aspects of catalysis. Sets of homologous proteins that catalyze different overall reactions but share an aspect of catalysis, such as a common partial reaction, are called mechanistically diverse superfamilies. The common mechanistic steps and structural characteristics of several of these superfamilies, including the enolase, Nudix, amidohydrolase, and haloacid dehalogenase superfamilies have been characterized. In addition, studies of mechanistically diverse superfamilies are helping to elucidate mechanisms of functional diversification, such as catalytic promiscuity. Understanding how enzyme superfamilies evolve is vital for accurate genome annotation, predicting protein functions, and protein engineering.

Binding Sites↗

Evolution of structure and function in the o-succinylbenzoate synthase/N-acylamino acid racemase family of the enolase superfamily.

Understanding how proteins evolve to provide both exquisite specificity and proficient activity is a fundamental problem in biology that has implications for protein function prediction and protein engineering. To study this problem, we analyzed the evolution of structure and function in the o-succinylbenzoate synthase/N-acylamino acid racemase (OSBS/NAAAR) family, part of the mechanistically diverse enolase superfamily. Although all characterized members of the family catalyze the OSBS reaction, this family is extraordinarily divergent, with some members sharing <15% identity. In addition, a member of this family, Amycolatopsis OSBS/NAAAR, is promiscuous, catalyzing both dehydration and racemization. Although the OSBS/NAAAR family appears to have a single evolutionary origin, no sequence or structural motifs unique to this family could be identified; all residues conserved in the family are also found in enolase superfamily members that have different functions. Based on their species distribution, several uncharacterized proteins similar to Amycolatopsis OSBS/NAAAR appear to have been transmitted by lateral gene transfer. Like Amycolatopsis OSBS/NAAAR, these might have additional or alternative functions to OSBS because many are from organisms lacking the pathway in which OSBS is an intermediate. In addition to functional differences, the OSBS/NAAAR family exhibits surprising structural variations, including large differences in orientation between the two domains. These results offer several insights into protein evolution. First, orthologous proteins can exhibit significant structural variation, and specificity can be maintained with little conservation of ligand-contacting residues. Second, the discovery of a set of proteins similar to Amycolatopsis OSBS/NAAAR supports the hypothesis that new protein functions evolve through promiscuous intermediates. Finally, a combination of evolutionary, structural, and sequence analyses identified characteristics that might prime proteins, such as Amycolatopsis OSBS/NAAAR, for the evolution of new activities.

Actinobacteria↗

Leveraging enzyme structure-function relationships for functional inference and experimental design: the structure-function linkage database.

The study of mechanistically diverse enzyme superfamilies-collections of enzymes that perform different overall reactions but share both a common fold and a distinct mechanistic step performed by key conserved residues-helps elucidate the structure-function relationships of enzymes. We have developed a resource, the structure-function linkage database (SFLD), to analyze these structure-function relationships. Unique to the SFLD is its hierarchical classification scheme based on linking the specific partial reactions (or other chemical capabilities) that are conserved at the superfamily, subgroup, and family levels with the conserved structural elements that mediate them. We present the results of analyses using the SFLD in correcting misannotations, guiding protein engineering experiments, and elucidating the function of recently solved enzyme structures from the structural genomics initiative. The SFLD is freely accessible at http://sfld.rbvi.ucsf.edu.

Animals↗

D-Ribulose 5-phosphate 3-epimerase: functional and structural relationships to members of the ribulose-phosphate binding (beta/alpha)8-barrel superfamily.

The "ribulose phosphate binding" superfamily defined by the Structural Classification of Proteins (SCOP) database is considered the result of divergent evolution from a common (beta/alpha)(8)-barrel ancestor. The superfamily includes d-ribulose 5-phosphate 3-epimerase (RPE), orotidine 5'-monophosphate decarboxylase (OMPDC), and 3-keto-l-gulonate 6-phosphate decarboxylase (KGPDC), members of the OMPDC suprafamily, as well as enzymes involved in histidine and tryptophan biosynthesis that utilize phosphorylated metabolites as substrates. We now report studies of the functional and structural relationships of RPE to the members of the superfamily. As suggested by the results of crystallographic studies of the RPEs from rice [Jelakovic, S., Kopriva, S., Suss, K. H., and Schulz, G. E. (2003) J. Mol. Biol. 326, 127-35] and Plasmodium falciparum [Caruthers, J., Bosch, J., Bucker, F., Van Voorhis, W., Myler, P., Worthey, E., Mehlin, C., Boni, E., De Titta, G., Luft, J., Kalyuzhniy, O., Anderson, L., Zucker, F., Soltis, M., and Hol, W. G. J. (2006) Proteins 62, 338-42], the RPE from Streptococcus pyogenes is activated by Zn(2+) which binds with a stoichiometry of one ion per polypeptide. Although wild type RPE has a high affinity for Zn(2+) and inactive apoenzyme cannot be prepared, the affinity for Zn(2+) is decreased by alanine substitutions for the two histidine residues that coordinate the Zn(2+) ion (H34A and H67A); these mutant proteins can be prepared in an inactive, metal-free form and activated by exogenous Zn(2+). The crystal structure of the RPE was solved at 1.8 A resolution in the presence of d-xylitol 5-phosphate, an inert analogue of the d-xylulose 5-phosphate substrate. This structure suggests that the 2,3-enediolate intermediate in the 1,1-proton transfer reaction is stabilized by bidentate coordination to the Zn(2+) that also is liganded to His 34, Asp 36, His 67, and Asp 176; the carboxylate groups of the Asp residues are positioned also to function as the acid/base catalysts. Although the conformation of the bound analogue resembles those of ligands bound in the active sites of OMPDC and KGPDC, the identities of the active site residues that coordinate the essential Zn(2+) and participate as acid/base catalysts are not conserved. We conclude that only the phosphate binding motif located at the ends of the seventh and eighth beta-strands of the (beta/alpha)(8)-barrel is functionally conserved among RPE, OMPDC, and KGPDC, consistent with the hypothesis that the members of the "ribulose phosphate binding" (beta/alpha)(8)-barrel "superfamily" as defined by SCOP have not evolved by evolutionary processes involving the intact (beta/alpha)(8)-barrel. Instead, this "superfamily" may result from assembly from smaller modules, including the conserved phosphate binding motif associated with the C-terminal (beta/alpha)(2)-quarter barrel.

Amino Acid Motifs↗

A gold standard set of mechanistically diverse enzyme superfamilies.

Superfamily and family analyses provide an effective tool for the functional classification of proteins, but must be automated for use on large datasets. We describe a 'gold standard' set of enzyme superfamilies, clustered according to specific sequence, structure, and functional criteria, for use in the validation of family and superfamily clustering methods. The gold standard set represents four fold classes and differing clustering difficulties, and includes five superfamilies, 91 families, 4,887 sequences and 282 structures.

Conserved Sequence↗

Automated discovery of 3D motifs for protein function annotation.

MOTIVATION: Function inference from structure is facilitated by the use of patterns of residues (3D motifs), normally identified by expert knowledge, that correlate with function. As an alternative to often limited expert knowledge, we use machine-learning techniques to identify patterns of 3-10 residues that maximize function prediction. This approach allows us to test the assumption that residues that provide function are the most informative for predicting function. RESULTS: We apply our method, GASPS, to the haloacid dehalogenase, enolase, amidohydrolase and crotonase superfamilies and to the serine proteases. The motifs found by GASPS are as good at function prediction as 3D motifs based on expert knowledge. The GASPS motifs with the greatest ability to predict protein function consist mainly of known functional residues. However, several residues with no known functional role are equally predictive. For four groups, we show that the predictive power of our 3D motifs is comparable with or better than approaches that use the entire fold (Combinatorial-Extension) or sequence profiles (PSI-BLAST). AVAILABILITY: Source code is freely available for academic use by contacting the authors. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Algorithms↗

The International Gene Trap Consortium Website: a portal to all publicly available gene trap cell lines in mouse.

Gene trapping is a method of generating murine embryonic stem (ES) cell lines containing insertional mutations in known and novel genes. A number of international groups have used this approach to create sizeable public cell line repositories available to the scientific community for the generation of mutant mouse strains. The major gene trapping groups worldwide have recently joined together to centralize access to all publicly available gene trap lines by developing a user-oriented Website for the International Gene Trap Consortium (IGTC). This collaboration provides an impressive public informatics resource comprising approximately 45 000 well-characterized ES cell lines which currently represent approximately 40% of known mouse genes, all freely available for the creation of knockout mice on a non-collaborative basis. To standardize annotation and provide high confidence data for gene trap lines, a rigorous identification and annotation pipeline has been developed combining genomic localization and transcript alignment of gene trap sequence tags to identify trapped loci. This information is stored in a new bioinformatics database accessible through the IGTC Website interface. The IGTC Website (www.genetrap.org) allows users to browse and search the database for trapped genes, BLAST sequences against gene trap sequence tags, and view trapped genes within biological pathways. In addition, IGTC data have been integrated into major genome browsers and bioinformatics sites to provide users with outside portals for viewing this data. The development of the IGTC Website marks a major advance by providing the research community with the data and tools necessary to effectively use public gene trap resources for the large-scale characterization of mammalian gene function.

Animals↗

Divergent evolution in the enolase superfamily: the interplay of mechanism and specificity.

The members of the mechanistically diverse enolase superfamily catalyze different overall reactions. Each shares a partial reaction in which an active site base abstracts the alpha-proton of the carboxylate substrate to generate an enolate anion intermediate that is stabilized by coordination to the essential Mg(2+) ion; the intermediates are then directed to different products in the different active sites. In this minireview, our current understanding of structure/function relationships in the divergent members of the superfamily is reviewed, and the use of this knowledge for our future studies is proposed.

Amino Acid Sequence↗

Introduction to informatics approaches in structural genomics: modeling and representation of function from macromolecular structure.

Despite the advantages provided by the enormous recent increases in the availability of structural information, functional assignment for the large number of proteins represented in the sequence and structural genomics projects remains a pressing problem for genomic era biology. This section describes work relevant to this problem from several perspectives, including new approaches that take advantage of combined structure and sequence-based classification. Leveraging of genomic context and evolutionary information to improve classification and predictive power is a second prominent theme in the papers represented here. Finally, issues in building a database for linking sequence, structural, and functional information are explored.

Genomics↗

Representing structure-function relationships in mechanistically diverse enzyme superfamilies.

The prediction of protein function from structure or sequence data remains a problem best addressed by leveraging information available from previously determined structure-function relationships. In the case of enzymes, the study of mechanistically diverse superfamilies can provide a rich source of structure-function information useful in functional determination and enzyme engineering. To access these relationships using a computational resource, several issues must be addressed regarding the representation of enzyme function, the organization of structure-function relationships in the superfamily context, the handling of misannotations, and reliability of classifications and evidence. We discuss here our approaches to solving these problems in the development of a Structure-Function Linkage Database (SFLD) (online at http://sfld.rbvi.ucsf.edu).

Amino Acid Sequence↗

Divergence of function in the thioredoxin fold suprafamily: evidence for evolution of peroxiredoxins from a thioredoxin-like ancestor.

The thioredoxin fold is found in proteins that serve a wide variety of functions. Among these are peroxiredoxins, which catalyze the reduction of hydrogen peroxide and alkyl peroxides. Although the common structural fold shared by thioredoxins and peroxiredoxins suggests the possibility that they have evolved from a common progenitor, it has been difficult to examine this hypothesis in depth because pairwise sequence identities between proteins in these two superfamilies are statistically insignificant. Using the Shotgun program, we have found that sequences of reductases involved in maturation of cytochromes in certain bacteria bridge the sequences of thioredoxins and peroxiredoxins. Analysis of motifs found in a divergent set of thioredoxins, cytochrome maturation proteins, and peroxiredoxins provides further support for an evolutionary relationship between these proteins. Within the conserved motifs are specific residues that are characteristic of individual protein classes, and therefore are likely to be involved in the specific functions of those classes. We have used this information, in combination with existing structural and functional information, to gain new insight into the structure-function relationships in these proteins and to construct a model for the emergence of peroxiredoxins from a thioredoxin-like ancestor.

Amino Acid Motifs↗

Isoleucine 69 and valine 325 form a specificity pocket in human muscle creatine kinase.

Creatine kinase (CK) catalyzes the reversible phosphorylation of creatine by ATP. From a structural perspective, the enzyme utilizes two flexible loop regions to sequester and position the substrates for catalysis. There has been debate over the specific roles of the flexible loops in substrate specificity and catalysis in CK and other related phosphagen kinases. In CK, two hydrophobic loop residues, I69 and V325, make contacts with the N-methyl group of creatine. In this study, we report the alteration of the substrate specificity of CK through the mutagenesis of V325. The V325 to glutamate mutation results in a more than 100-fold preference for glycocyamine, while mutation of V325 to alanine results in a slight preference of the enzyme for cyclocreatine (1-carboxymethyl-2-iminoimidazolidine). This study enhances our understanding of how the active sites of phosphagen kinases have evolved to recognize their respective substrates and catalyze their reactions.

Amino Acid Sequence↗

Superfamily active site templates.

We show that three-dimensional signatures consisting of only a few functionally important residues can be diagnostic of membership in superfamilies of enzymes. Using the enolase superfamily as a model system, we demonstrate that such a signature, or template, can identify superfamily members in structural databases with high sensitivity and specificity. This is remarkable because superfamilies can be highly diverse, with members catalyzing many different overall reactions; the unifying principle can be a conserved partial reaction or chemical capability. Our definition of a superfamily thus hinges on the disposition of residues involved in a conserved function, rather than on fold similarity alone. A clear advantage of basing structure searches on such active site templates rather than on fold similarity is the specificity with which superfamilies with distinct functional characteristics can be identified within a large set of proteins with the same fold, such as the (beta/alpha)8 barrels. Preliminary results are presented for an additional group of enzymes with a different fold, the haloacid dehalogenase superfamily, suggesting that this approach may be generally useful for assigning reading frames of unknown function to specific superfamilies and thereby allowing inference of some of their functional properties.

Binding Sites↗

Intersect: identification and visualization of overlaps in database search results.

UNLABELLED: The determination of distant evolutionary relationships remains an important biological problem, and distant homologs often appear in statistically insignificant regions of sequence similarity searches. Intersect is a computer program designed to identify and visualize the overlaps between sets of sequences reported by multiple database searches. This capability extends the usefulness of database search results and aids researchers in identifying the individual sequences that best bridge sequence families and superfamilies. AVAILABILITY: The Intersect program is available from the Babbitt laboratory website at http://www.babbittlab.ucsf.edu/software/intersect

Algorithms↗

Evolutionary potential of (beta/alpha)8-barrels: functional promiscuity produced by single substitutions in the enolase superfamily.

The members of the mechanistically diverse, (beta/alpha)(8)-barrel fold-containing enolase superfamily evolved from a common progenitor but catalyze different reactions using a conserved partial reaction. The molecular pathway for natural divergent evolution of function in the superfamily is unknown. We have identified single-site mutants of the (beta/alpha)(8)-barrel domains in both the l-Ala-d/l-Glu epimerase from Escherichia coli (AEE) and the muconate lactonizing enzyme II from Pseudomonas sp. P51 (MLE II) that catalyze the o-succinylbenzoate synthase (OSBS) reaction as well as the wild-type reaction. These enzymes are members of the MLE subgroup of the superfamily, share conserved lysines on opposite sides of their active sites, but catalyze acid- and base-mediated reactions with different mechanisms. A comparison of the structures of AEE and the OSBS from E. coli was used to design the D297G mutant of AEE; the E323G mutant of MLE II was isolated from directed evolution experiments. Although neither wild-type enzyme catalyzes the OSBS reaction, both mutants complement an E. coli OSBS auxotroph and have measurable levels of OSBS activity. The analogous mutations in the D297G mutant of AEE and the E323G mutant of MLE II are each located at the end of the eighth beta-strand of the (beta/alpha)(8)-barrel and alter the ability of AEE and MLE II to bind the substrate of the OSBS reaction. The substitutions relax the substrate specificity, thereby allowing catalysis of the mechanistically diverse OSBS reaction with the assistance of the active site lysines. The generation of functionally promiscuous and mechanistically diverse enzymes via single-amino acid substitutions likely mimics the natural divergent evolution of enzymatic activities and also highlights the utility of the (beta/alpha)(8)-barrel as a scaffold for new function.

Amino Acid Sequence↗

BayGenomics: a resource of insertional mutations in mouse embryonic stem cells.

The BayGenomics gene-trap resource (http://baygenomics.ucsf.edu) provides researchers with access to thousands of mouse embryonic stem (ES) cell lines harboring characterized insertional mutations in both known and novel genes. Each cell line contains an insertional mutation in a specific gene. The identity of the gene that has been interrupted can be determined from a DNA sequence tag. Approximately 75% of our cell lines contain insertional mutations in known mouse genes or genes that share strong sequence similarities with genes that have been identified in other organisms. These cell lines readily transmit the mutation to the germline of mice and many mutant lines of mice have already been generated from this resource. BayGenomics provides facile access to our entire database, including sequence tags for each mutant ES cell line, through the World Wide Web. Investigators can browse our resource, search for specific entries, download any portion of our database and BLAST sequences of interest against our entire set of cell line sequence tags. They can then obtain the mutant ES cell line for the purpose of generating knockout mice.

Animals↗