oskar is the only gene in the animal kingdom necessary and sufficient for specifying functional germ cells. However, oskar has only been identified in holometabolous ("higher") insects that specify their germline using specialized cytoplasm called germ plasm. Here we show that oskar evolved before the divergence of higher insects and provide evidence that its germline role is a recent evolutionary innovation. We identify an oskar ortholog in a basally branching insect, the cricket Gryllus bimaculatus. In contrast to Drosophila oskar, Gb-oskar is not required for germ cell formation or axial patterning. Instead, Gb-oskar is expressed in neuroblasts of the brain and CNS and is required for neural development. Taken together with reports of a neural role for Drosophila oskar, our data demonstrate that oskar arose nearly 50 million years earlier in insect evolution than previously thought, where it may have played an ancestral neural role, and was co-opted to its well-known essential germline role in holometabolous insects.

The Statistical Inference of Function Through Evolutionary Relationships (SIFTER) framework uses a statistical graphical model that applies phylogenetic principles to automate precise protein function prediction. Here we present a revised approach (SIFTER version 2.0) that enables annotations on a genomic scale. SIFTER 2.0 produces equivalently precise predictions compared to the earlier version on a carefully studied family and on a collection of 100 protein families. We have added an approximation method to SIFTER 2.0 and show a 500-fold improvement in speed with minimal impact on prediction results in the functionally diverse sulfotransferase protein family. On the Nudix protein family, previously inaccessible to the SIFTER framework because of the 66 possible molecular functions, SIFTER achieved 47.4% accuracy on experimental data (where BLAST achieved 34.0%). Finally, we used SIFTER to annotate all of the Schizosaccharomyces pombe proteins with experimental functional characterizations, based on annotations from proteins in 46 fungal genomes. SIFTER precisely predicted molecular function for 45.5% of the characterized proteins in this genome, as compared with four current function prediction methods that precisely predicted function for 62.6%, 30.6%, 6.0%, and 5.7% of these proteins. We use both precision-recall curves and ROC analyses to compare these genome-scale predictions across the different methods and to assess performance on different types of applications. SIFTER 2.0 is capable of predicting protein molecular function for large and functionally diverse protein families using an approximate statistical model, enabling phylogenetics-based protein function prediction for genome-wide analyses. The code for SIFTER and protein family data are available at http://sifter.berkeley.edu. [Supplemental material is available for this article.]Automated protein function prediction is an important challenge for computational biology because protein function is difficult to describe and represent, protein databases are littered with annotation errors, and our understanding of how molecular functions arise and mutate over evolutionary time is far from complete. Because biologists depend on protein function annotations for insight and analysis, automated methods have been used extensively to compensate for the relative dearth of experimental characterizations. Although there are 10 7 protein sequences in the comprehensive UniProt database (The UniProt Consortium 2010), <5% have annotations from the Gene Ontology Annotation (GOA) database (Barrell et al. 2009). Far fewer (0.2%) have been manually annotated, and only 0.25% of those manual annotations are from the molecular function ontology in Gene Ontology (GO) (The Gene Ontology Consortium 2010) and are based on experimental evidence. Because of the need for so many annotations, function prediction methods are often assessed based on annotation quantity rather than quality, increasing the number of false positive function annotations and polluting databases (Galper...

The Nudix homology clan encompasses over 80,000 protein domains from all three domains of life, defined by homology to each other. Proteins with a domain from this clan fall into four general functional classes: pyrophosphohydrolases, isopentenyl diphosphate isomerases (IDIs), adenine/guanine mismatch-specific adenine glycosylases (A/G-specific adenine glycosylases), and non-enzymatic activities such as protein/protein interaction and transcriptional regulation. The largest group, pyrophosphohydrolases, encompasses more than 100 distinct hydrolase specificities. To understand the evolution of this vast number of activities, we assembled and analyzed experimental and structural data for 205 Nudix proteins collected from the literature. We corrected erroneous functions or provided more appropriate descriptions for 53 annotations described in the Gene Ontology Annotation database in this family, and propose 275 new experimentally-based annotations. We manually constructed a structure-guided sequence alignment of 78 Nudix proteins. Using the structural alignment as a seed, we then made an alignment of 347 “select” Nudix homology domains, curated from structurally determined, functionally characterized, or phylogenetically important Nudix domains. Based on our review of Nudix pyrophosphohydrolase structures and specificities, we further analyzed a loop region downstream of the Nudix hydrolase motif previously shown to contact the substrate molecule and possess known functional motifs. This loop region provides a potential structural basis for the functional radiation and evolution of substrate specificity within the hydrolase family. Finally, phylogenetic analyses of the 347 select protein domains and of the complete Nudix homology clan revealed general monophyly with regard to function and a few instances of probable homoplasy.

How do proteins evolve novel functions? To address this question, we are studying the evolution of a mammalian toxin, the serine protease BLTX [1], from the salivary glands of the North American shrew Blarina brevicauda. Here, we examine the molecular changes responsible for promoting BLTX toxicity. First, we show that regulatory loops surrounding the BLTX active site have evolved adaptively via acquisition of small insertions and subsequent accelerated sequence evolution. Second, these mutations introduce a novel chemical environment into the catalytic cleft of BLTX. Third, molecular-dynamic simulations show that the observed changes create a novel chemical and physical topology consistent with increased enzyme catalysis. Finally, we show that a toxic serine protease from the Mexican beaded lizard (GTX) [2] has evolved convergently through almost identical functional changes. Together, these results suggest that the evolution of toxicity might be predictable-arising via adaptive structural modification of analogous labile regulatory loops of an ancestral serine protease-and thus might aid in the identification of other toxic proteins.

No abstract

The nearly 50,000 known Nudix proteins have a diverse array of functions, of which the most extensively studied is the catalyzed hydrolysis of aberrant nucleotide triphosphates. The functions of 171 Nudix proteins have been characterized to some degree, although physiological relevance of the assayed activities has not always been conclusively demonstrated. We investigated substrate specificity for eight structurally characterized Nudix proteins, whose functions were unknown. These proteins were screened for hydrolase activity against a 74‐compound library of known Nudix enzyme substrates. We found substrates for four enzymes with k cat/K m values >10,000 M−1 s−1: Q92EH0_LISIN of Listeria innocua serovar 6a against ADP‐ribose, Q5LBB1_BACFN of Bacillus fragilis against 5‐Me‐CTP, and Q0TTC5_CLOP1 and Q0TS82_CLOP1 of Clostridium perfringens against 8‐oxo‐dATP and 3'‐dGTP, respectively. To ascertain whether these identified substrates were physiologically relevant, we surveyed all reported Nudix hydrolytic activities against NTPs. Twenty‐two Nudix enzymes are reported to have activity against canonical NTPs. With a single exception, we find that the reported k cat/K m values exhibited against these canonical substrates are well under 105 M−1 s−1. By contrast, several Nudix enzymes show much larger k cat/K m values (in the range of 105 to >107 M−1 s−1) against noncanonical NTPs. We therefore conclude that hydrolytic activities exhibited by these enzymes against canonical NTPs are not likely their physiological function, but rather the result of unavoidable collateral damage occasioned by the enzymes' inability to distinguish completely between similar substrate structures. Proteins 2016; 84:1810–1822. © 2016 The Authors Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.

Aminotransferases are essential enzymes involved in the central metabolism of all organisms. The Iα subfamily of aspartate and tyrosine aminotransferases (AATases and TATases) is the bestcharacterized grouping, but only eight enzymes from this subfamily, representing relatively little sequence diversity, have been experimentally characterized for substrate specificity (i.e., AATase vs. TATase). Genome annotation, based on this limited dataset, provides tentative assignments for all sequenced members of this subfamily. This procedure is, however, subject to error, particularly when the experimental basis set is limited. To address this problem we cloned twelve additional subfamily Iα enzymes from an evolutionarily divergent set of organisms. Nine were purified to homogeneity after heterologous expression in E. coli in native, intein-tagged or His 6 -tagged forms and the two S. cerevisiae isoforms were recombinantly produced in yeast. The effects of the Cterminal tags on expression, purification and enzyme activity are discussed.

The subfamily Iα aminotransferases are typically categorized as having narrow specificity toward carboxylic amino acids (AATases), or broad specificity that includes aromatic amino acid substrates (TATases). Because of their general role in central metabolism and, more specifically, their association with liver-related diseases in humans, this subfamily is biologically interesting. The substrate specificities for only a few members of this subfamily have been reported, and the reliable prediction of substrate specificity from protein sequence has remained elusive. In this study, a diverse set of aminotransferases was chosen for characterization based on a scoring system that measures the sequence divergence of the active site. The enzymes that were experimentally characterized include both narrow-specificity AATases and broad-specificity TATases, as well as AATases with broader-specificity and TATases with narrower-specificity than the previously known family members. Molecular function and phylogenetic analyses underscored the complexity of this family's evolution as the TATase function does not follow a single evolutionary thread, but rather appears independently multiple times during the evolution of the subfamily. The additional functional characterizations described in this article, alongside a detailed sequence and phylogenetic analysis, provide some novel clues to understanding the evolutionary mechanisms at work in this family.