Molecular function prediction for a family exhibiting evolutionary tendencies toward substrate specificity swapping: recurrence of tyrosine aminotransferase activity in the Iα subfamily.

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.

Graph analysis of ß2 adrenergic receptor structures: a "social network" of GPCR residues.

PURPOSE: G protein-coupled receptors (GPCRs) are a superfamily of membrane proteins of vast pharmaceutical interest. Here, we describe a graph theory-based analysis of the structure of the ß2 adrenergic receptor (ß2 AR), a prototypical GPCR. In particular, we illustrate the network of direct and indirect interactions that link each amino acid residue to any other residue of the receptor. METHODS: Networks of interconnected amino acid residues in proteins are analogous to social networks of interconnected people. Hence, they can be studied through the same analysis tools typically employed to analyze social networks - or networks in general - to reveal patterns of connectivity, influential members, and dynamicity. We focused on the analysis of closeness-centrality, which is a measure of the overall connectivity distance of the member of a network to all other members. RESULTS: The residues endowed with the highest closeness-centrality are located in the middle of the seven transmembrane domains (TMs). In particular, they are mostly located in the middle of TM2, TM3, TM6 or TM7, while fewer of them are located in the middle of TM1, TM4 or TM5. At the cytosolic end of TM6, the centrality detected for the active structure is markedly lower than that detected for the corresponding residues in the inactive structures. Moreover, several residues acquire centrality when the structures are analyzed in the presence of ligands. Strikingly, there is little overlap between the residues that acquire centrality in the presence of the ligand in the blocker-bound structures and the agonist-bound structures. CONCLUSIONS: Our results reflect the fact that the receptor resembles a bow tie, with a rather tight knot of closely interconnected residues and two ends that fan out in two opposite directions: one toward the extracellular space, which hosts the ligand binding cavity, and one toward the cytosol, which hosts the G protein binding cavity. Moreover, they underscore how interaction network is by the conformational rearrangements concomitant with the activation of the receptor and by the presence of agonists or blockers.

Comparative analysis of mutant tyrosine kinase chemical rescue.

Protein tyrosine kinases are critical cell signaling enzymes. These enzymes have a highly conserved Arg residue in their catalytic loop which is present two residues or four residues downstream from an absolutely conserved Asp catalytic base. Prior studies on protein tyrosine kinases Csk and Src revealed the potential for chemical rescue of catalytically deficient mutant kinases (Arg to Ala mutations) by small diamino compounds, particularly imidazole; however, the potency and efficiency of rescue was greater for Src. This current study further examines the structural and kinetic basis of rescue for mutant Src as compared to mutant Abl tyrosine kinase. An X-ray crystal structure of R388A Src revealed the surprising finding that a histidine residue of the N-terminus of a symmetry-related kinase inserts into the active site of the adjacent Src and mimics the hydrogen-bonding pattern seen in wild-type protein tyrosine kinases. Abl R367A shows potent and efficient rescue more comparable to Src, even though its catalytic loop is more like that of Csk. Various enzyme redesigns of the active sites indicate that the degree and specificity of rescue are somewhat flexible, but the overall properties of the enzymes and rescue agents play an overarching role. The newly discovered rescue agent 2-aminoimidazole is about as efficient as imidazole in rescuing R/A Src and Abl. Rate vs pH studies with these imidazole analogues suggest that the protonated imidazolium is the preferred form for chemical rescue, consistent with structural models. The efficient rescue seen with mutant Abl points to the potential of this approach to be used effectively to analyze Abl phosphorylation pathways in cells.

Recombinant expression of twelve evolutionarily diverse subfamily Ialpha aminotransferases.

Aminotransferases are essential enzymes involved in the central metabolism of all organisms. The Ialpha subfamily of aspartate and tyrosine aminotransferases (AATases and TATases) is the best-characterized 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 Ialpha enzymes from an evolutionarily divergent set of organisms. Nine were purified to homogeneity after heterologous expression in Escherichia coli in native, intein-tagged or His(6)-tagged forms. The two Saccharomyces cerevisiae isoforms were recombinantly produced in yeast. The effects of the C-terminal tags on expression, purification and enzyme activity are discussed.

A lock on phosphotyrosine signaling.

Protein phosphorylation plays a major role in cell signaling and human disease, so understanding the effects of tyrosine phosphorylation on protein structure and function is an area of intense investigation. A new technique allows site-specific incorporation of a non-hydrolyzable phosphotyrosine analogue into recombinant proteins, providing a new strategy for research in this important area.

Protein molecular function prediction by Bayesian phylogenomics.

We present a statistical graphical model to infer specific molecular function for unannotated protein sequences using homology. Based on phylogenomic principles, SIFTER (Statistical Inference of Function Through Evolutionary Relationships) accurately predicts molecular function for members of a protein family given a reconciled phylogeny and available function annotations, even when the data are sparse or noisy. Our method produced specific and consistent molecular function predictions across 100 Pfam families in comparison to the Gene Ontology annotation database, BLAST, GOtcha, and Orthostrapper. We performed a more detailed exploration of functional predictions on the adenosine-5'-monophosphate/adenosine deaminase family and the lactate/malate dehydrogenase family, in the former case comparing the predictions against a gold standard set of published functional characterizations. Given function annotations for 3% of the proteins in the deaminase family, SIFTER achieves 96% accuracy in predicting molecular function for experimentally characterized proteins as reported in the literature. The accuracy of SIFTER on this dataset is a significant improvement over other currently available methods such as BLAST (75%), GeneQuiz (64%), GOtcha (89%), and Orthostrapper (11%). We also experimentally characterized the adenosine deaminase from Plasmodium falciparum, confirming SIFTER's prediction. The results illustrate the predictive power of exploiting a statistical model of function evolution in phylogenomic problems. A software implementation of SIFTER is available from the authors.