PSiFR: an integrated resource for prediction of protein structure and function

In the post-genomic era, the annotation of protein function facilitates the understanding of various biological processes. To extend the range of function annotation methods to the twilight zone of sequence identity, we have developed approaches that exploit both protein tertiary structure and/or protein sequence evolutionary relationships. To serve the scientific community, we have integrated the structure prediction tools, TASSER, TASSER-Lite and METATASSER, and the functional inference tools, FINDSITE, a structure-based algorithm for binding site prediction, Gene Ontology molecular function inference and ligand screening, EFICAz(2), a sequence-based approach to enzyme function inference and DBD-hunter, an algorithm for predicting DNA-binding proteins and associated DNA-binding residues, into a unified web resource, Protein Structure and Function prediction Resource (PSiFR).

AVAILABILITY AND IMPLEMENTATION

PSiFR is freely available for use on the web at http://psifr.cssb.biology.gatech.edu/

EFICAz2: enzyme function inference by a combined approach enhanced by machine learning

BACKGROUND

We previously developed EFICAz, an enzyme function inference approach that combines predictions from non-completely overlapping component methods. Two of the four components in the original EFICAz are based on the detection of functionally discriminating residues (FDRs). FDRs distinguish between member of an enzyme family that are homofunctional (classified under the EC number of interest) or heterofunctional (annotated with another EC number or lacking enzymatic activity). Each of the two FDR-based components is associated to one of two specific kinds of enzyme families. EFICAz exhibits high precision performance, except when the maximal test to training sequence identity (MTTSI) is lower than 30%. To improve EFICAz's performance in this regime, we: i) increased the number of predictive components and ii) took advantage of consensual information from the different components to make the final EC number assignment.

RESULTS

We have developed two new EFICAz components, analogs to the two FDR-based components, where the discrimination between homo and heterofunctional members is based on the evaluation, via Support Vector Machine models, of all the aligned positions between the query sequence and the multiple sequence alignments associated to the enzyme families. Benchmark results indicate that: i) the new SVM-based components outperform their FDR-based counterparts, and ii) both SVM-based and FDR-based components generate unique predictions. We developed classification tree models to optimally combine the results from the six EFICAz components into a final EC number prediction. The new implementation of our approach, EFICAz2, exhibits a highly improved prediction precision at MTTSI < 30% compared to the original EFICAz, with only a slight decrease in prediction recall. A comparative analysis of enzyme function annotation of the human proteome by EFICAz2 and KEGG shows that: i) when both sources make EC number assignments for the same protein sequence, the assignments tend to be consistent and ii) EFICAz2 generates considerably more unique assignments than KEGG.

CONCLUSION

Performance benchmarks and the comparison with KEGG demonstrate that EFICAz2 is a powerful and precise tool for enzyme function annotation, with multiple applications in genome analysis and metabolic pathway reconstruction. The EFICAz2 web service is available at: http://cssb.biology.gatech.edu/skolnick/webservice/EFICAz2/index.html.

Identification of metabolites with anticancer properties by computational metabolomics

Background

Certain endogenous metabolites can influence the rate of cancer cell growth. For example, diacylglycerol, ceramides and sphingosine, NAD⁺ and arginine exert this effect by acting as signaling molecules, while carrying out other important cellular functions. Metabolites can also be involved in the control of cell proliferation by directly regulating gene expression in ways that are signaling pathway-independent, e.g. by direct activation of transcription factors or by inducing epigenetic processes. The fact that metabolites can affect the cancer process on so many levels suggests that the change in concentration of some metabolites that occurs in cancer cells could have an active role in the progress of the disease.

Results

CoMet, a fully automated Computational Metabolomics method to predict changes in metabolite levels in cancer cells compared to normal references has been developed and applied to Jurkat T leukemia cells with the goal of testing the following hypothesis: Up or down regulation in cancer cells of the expression of genes encoding for metabolic enzymes leads to changes in intracellular metabolite concentrations that contribute to disease progression. All nine metabolites predicted to be lowered in Jurkat cells with respect to lymphoblasts that were examined (riboflavin, tryptamine, 3-sulfino-L-alanine, menaquinone, dehydroepiandrosterone, α-hydroxystearic acid, hydroxyacetone, seleno-L-methionine and 5,6-dimethylbenzimidazole), exhibited antiproliferative activity that has not been reported before, while only two (bilirubin and androsterone) of the eleven tested metabolites predicted to be increased or unchanged in Jurkat cells displayed significant antiproliferative activity.

Conclusion

These results: a) demonstrate that CoMet is a valuable method to identify potential compounds for experimental validation, b) indicate that cancer cell metabolism may be regulated to reduce the intracellular concentration of certain antiproliferative metabolites, leading to uninhibited cellular growth and c) suggest that many other endogenous metabolites with important roles in carcinogenesis are awaiting discovery.

Aerobic uranium (VI) bioprecipitation by metal-resistant bacteria isolated from radionuclide- and metal-contaminated subsurface soils

In this study, the immobilization of toxic uranium [U(VI)] mediated by the intrinsic phosphatase activities of naturally occurring bacteria isolated from contaminated subsurface soils was examined. The phosphatase phenotypes of strains belonging to the genera, Arthrobacter, Bacillus and Rahnella, previously isolated from subsurface soils at the US Department of Energy's (DOE) Oak Ridge Field Research Center (ORFRC), were determined. The ORFRC represents a unique, extreme environment consisting of highly acidic soils with co-occurring heavy metals, radionuclides and high nitrate concentrations. Isolates exhibiting phosphatase-positive phenotypes indicative of constitutive phosphatase activity were subsequently tested in U(VI) bioprecipitation assays. When aerobically grown in synthetic groundwater (pH 5.5) amended with 10 mM glycerol-3-phosphate (G3P), phosphatase-positive Bacillus and Rahnella spp. strains Y9-2 and Y9602 liberated sufficient phosphate to precipitate 73% and 95% of total soluble U added as 200 microM uranyl acetate respectively. In contrast, an Arthrobacter sp. X34 exhibiting a phosphatase-negative phenotype did not liberate phosphate from G3P or promote U(VI) precipitation. This study provides the first evidence of U(VI) precipitation via the phosphatase activity of naturally occurring Bacillus and Rahnella spp. isolated from the acidic subsurface at the DOE ORFRC.

TASSER: an automated method for the prediction of protein tertiary structures in CASP6

The recently developed TASSER (Threading/ASSembly/Refinement) method is applied to predict the tertiary structures of all CASP6 targets. TASSER is a hierarchical approach that consists of template identification by the threading program PROSPECTOR_3, followed by tertiary structure assembly via rearranging continuous template fragments. Assembly occurs using parallel hyperbolic Monte Carlo sampling under the guide of an optimized, reduced force field that includes knowledge-based statistical potentials and spatial restraints extracted from threading alignments. Models are automatically selected from the Monte Carlo trajectories in the low-temperature replicas using the clustering program SPICKER. For all 90 CASP targets/domains, PROSPECTOR_3 generates initial alignments with an average root-mean-square deviation (RMSD) to native of 8.4 A with 79% coverage. After TASSER reassembly, the average RMSD decreases to 5.4 A over the same aligned residues; the overall cumulative TM-score increases from 39.44 to 52.53. Despite significant improvements over the PROSPECTOR_3 template alignment observed in all target categories, the overall quality of the final models is essentially dictated by the quality of threading templates: The average TM-scores of TASSER models in the three categories are, respectively, 0.79 [comparative modeling (CM), 43 targets/domains], 0.47 [fold recognition (FR), 37 targets/domains], and 0.30 [new fold (NF), 10 targets/domains]. This highlights the need to develop novel (or improved) approaches to identify very distant targets as well as better NF algorithms.

On the origin and highly likely completeness of single-domain protein structures

The size and origin of the protein fold universe is of fundamental and practical importance. Analyzing randomly generated, compact sticky homopolypeptide conformations constructed in generic simplified and all-atom protein models, all have similar folds in the library of solved structures, the Protein Data Bank, and conversely, all compact, single-domain protein structures in the Protein Data Bank have structural analogues in the compact model set. Thus, both sets are highly likely complete, with the protein fold universe arising from compact conformations of hydrogen-bonded, secondary structures. Because side chains are represented by their Cbeta atoms, these results also suggest that the observed protein folds are insensitive to the details of side-chain packing. Sequence specificity enters both in fine-tuning the structure and thermodynamically stabilizing a given fold with respect to the set of alternatives. Scanning the models against a three-dimensional active-site library, close geometric matches are frequently found. Thus, the presence of active-site-like geometries also seems to be a consequence of the packing of compact, secondary structural elements. These results have significant implications for the evolution of protein structure and function.

EFICAz: a comprehensive approach for accurate genome-scale enzyme function inference

EFICAz (Enzyme Function Inference by Combined Approach) is an automatic engine for large-scale enzyme function inference that combines predictions from four different methods developed and optimized to achieve high prediction accuracy: (i) recognition of functionally discriminating residues (FDRs) in enzyme families obtained by a Conservation-controlled HMM Iterative procedure for Enzyme Family classification (CHIEFc), (ii) pairwise sequence comparison using a family specific Sequence Identity Threshold, (iii) recognition of FDRs in Multiple Pfam enzyme families, and (iv) recognition of multiple Prosite patterns of high specificity. For FDR (i.e. conserved positions in an enzyme family that discriminate between true and false members of the family) identification, we have developed an Evolutionary Footprinting method that uses evolutionary information from homofunctional and heterofunctional multiple sequence alignments associated with an enzyme family. The FDRs show a significant correlation with annotated active site residues. In a jackknife test, EFICAz shows high accuracy (92%) and sensitivity (82%) for predicting four EC digits in testing sequences that are <40% identical to any member of the corresponding training set. Applied to Escherichia coli genome, EFICAz assigns more detailed enzymatic function than KEGG, and generates numerous novel predictions.

Large-scale assessment of the utility of low-resolution protein structures for biochemical function assignment

MOTIVATION

Several protein function prediction methods employ structural features captured in three-dimensional (3D) descriptors of biologically relevant sites. These methods are successful when applied to high-resolution structures, but their detection ability in lower resolution predicted structures has only been tested for a few cases.

RESULTS

A method that automatically generates a library of 3D functional descriptors for the structure-based prediction of enzyme active sites (automated functional templates, 593 in total for 162 different enzymes), based on functional and structural information automatically extracted from public databases, has been developed and evaluated using decoy structures. The applicability to predicted structures was investigated by analyzing decoys of varying quality, derived from enzyme native structures. For 35% of decoy structures, our method identifies the active site in models having 3-4 A coordinate root mean square deviation from the native structure, a quality that is reachable using state of the art protein structure prediction algorithms.

AVAILABILITY

See http://www.bioinformatics.buffalo.edu/resources/aft/

TOUCHSTONE: A unified approach to protein structure prediction

We have applied the TOUCHSTONE structure prediction algorithm that spans the range from homology modeling to ab initio folding to all protein targets in CASP5. Using our threading algorithm PROSPECTOR that does not utilize input from metaservers, one threads against a representative set of PDB templates. If a template is significantly hit, Generalized Comparative Modeling designed to span the range from closely to distantly related proteins from the template is done. This involves freezing the aligned regions and relaxing the remaining structure to accommodate insertions or deletions with respect to the template. For all targets, consensus predicted side chain contacts from at least weakly threading templates are pooled and incorporated into ab initio folding. Often, TOUCHSTONE performs well in the CM to FR categories, with PROSPECTOR showing significant ability to identify analogous templates. When ab initio folding is done, frequently the best models are closer to the native state than the initial template. Among the particularly good predictions are T0130 in the CM/FR category, T0138 in the FR(H) category, T0135 in the FR(A) category, T0170 in the FR/NF category and T0181 in the NF category. Improvements in the approach are needed in the FR/NF and NF categories. Nevertheless, TOUCHSTONE was one of the best performing algorithms over all categories in CASP5.

Multimeric threading-based prediction of protein-protein interactions on a genomic scale: application to the Saccharomyces cerevisiae proteome

MULTIPROSPECTOR, a multimeric threading algorithm for the prediction of protein-protein interactions, is applied to the genome of Saccharomyces cerevisiae. Each possible pairwise interaction among more than 6000 encoded proteins is evaluated against a dimer database of 768 complex structures by using a confidence estimate of the fold assignment and the magnitude of the statistical interfacial potentials. In total, 7321 interactions between pairs of different proteins are predicted, based on 304 complex structures. Quality estimation based on the coincidence of subcellular localizations and biological functions of the predicted interactors shows that our approach ranks third when compared with all other large-scale methods. Unlike other in silico methods, MULTIPROSPECTOR is able to identify the residues that participate directly in the interaction. Three hundred seventy-four of our predictions can be found by at least one of the other studies, which is compatible with the overlap between two different other methods. From the analysis of the mRNA abundance data, our method does not bias towards proteins with high abundance. Finally, several relevant predictions involved in various functions are presented. In summary, we provide a novel approach to predict protein-protein interactions on a genomic scale that is a useful complement to experimental methods.

Involvement of the flavin si-face tyrosine on the structure and function of ferredoxin-NADP+ reductases

In ferredoxin-NADP(+) reductase (FNR), FAD is bound outside of an anti-parallel beta-barrel with the isoalloxazine lying in a two-tyrosine pocket. To elucidate the function of the flavin si-face tyrosine (Tyr-89 in pea FNR) on the enzyme structure and catalysis, we performed ab initio molecular orbital calculations and site-directed mutagenesis. Our results indicate that the position of Tyr-89 in pea FNR is mainly governed by the energetic minimum of the pairwise interaction between the phenol ring and the flavin. Moreover, most of FNR-like proteins displayed geometries for the si-face tyrosine phenol and the flavin, which correspond to the more negative free energy theoretical value. FNR mutants were obtained replacing Tyr-89 by Phe, Trp, Ser, or Gly. Structural and functional features of purified FNR mutants indicate that aromaticity on residue 89 is essential for FAD binding and proper folding of the protein. Moreover, hydrogen bonding through the Tyr-89 hydroxyl group may be responsible of the correct positioning of FAD and the substrate NADP(+)

Interaction of the targeting sequence of chloroplast precursors with Hsp70 molecular chaperones

We have analyzed the interaction of DnaK and plant Hsp70 proteins with the wild-type ferredoxin-NADP+ reductase precursor (preFNR) and mutants containing amino-acid replacements in the targeting sequence. Using an algorithm already developed [Rüdiger, S., Germeroth, L., Schneider-Mergener, J. & Bukau, B. (1997) EMBO J. 16, 1501-1507] we observed that 75% of the 727 plastid precursor proteins analyzed contained at least one site with high likelihood of DnaK binding in their transit peptides. Statistical analysis showed a decrease of DnaK binding site frequency within the first 15 amino-acid residues of the transit peptides. Using fusion proteins we detected the interaction of DnaK with the transit peptide of the folded preFNR but not with the mature region of the protein. Discharge of DnaK from the presequence was favored by addition of MgATP. When a putative DnaK binding site was artificially added at the N-terminus of the mature protein, we observed formation of complexes with bacterial and plant Hsp70 molecular chaperones. Reducing the likelihood of DnaK binding by directed mutagenesis of the presequence increased the release of bound DnaK. The Hsp70 proteins from plastids and plant cell cytosol also interacted with the preFNR transit peptide. Overall results are discussed in the context of the proposed models to explain the organelle protein import.

Competition between C-terminal tyrosine and nicotinamide modulates pyridine nucleotide affinity and specificity in plant ferredoxin-NADP(+) reductase

Chloroplast ferredoxin-NADP(+) reductase has a 32,000-fold preference for NADPH over NADH, consistent with its main physiological role of NADP(+) photoreduction for de novo carbohydrate biosynthesis. Although it is distant from the 2'-phosphoryl group of NADP(+), replacement of the C-terminal tyrosine (Tyr(308) in the pea enzyme) by Trp, Phe, Gly, and Ser produced enzyme forms in which the preference for NADPH over NADH was decreased about 2-, 10-, 300-, and 400-fold, respectively. Remarkably, in the case of the Y308S mutant, the k(cat) value for the NADH-dependent activity approached that of the NADPH-dependent activity of the wild-type enzyme. Furthermore, difference spectra of the NAD(+) complexes revealed that the nicotinamide ring of NAD(+) binds at nearly full occupancy in the active site of both the Y308G and Y308S mutants. These results correlate well with the k(cat) values obtained with these mutants in the NADH-ferricyanide reaction. The data presented support the hypothesis that specific recognition of the 2'-phosphate group of NADP(H) is required but not sufficient to ensure a high degree of discrimination against NAD(H) in ferredoxin-NADP(+) reductase. Thus, the C-terminal tyrosine enhances the specificity of the reductase for NADP(H) by destabilizing the interaction of a moiety common to both coenzymes, i.e. the nicotinamide.

A productive NADP+ binding mode of ferredoxin-NADP + reductase revealed by protein engineering and crystallographic studies

The flavoenzyme ferredoxin-NADP+ reductase (FNR) catalyzes the production of NADPH during photosynthesis. Whereas the structures of FNRs from spinach leaf and a cyanobacterium as well as many of their homologs have been solved, none of these studies has yielded a productive geometry of the flavin-nicotinamide interaction. Here, we show that this failure occurs because nicotinamide binding to wild type FNR involves the energetically unfavorable displacement of the C-terminal Tyr side chain. We used mutants of this residue (Tyr 308) of pea FNR to obtain the structures of productive NADP+ and NADPH complexes. These structures reveal a unique NADP+ binding mode in which the nicotinamide ring is not parallel to the flavin isoalloxazine ring, but lies against it at an angle of approximately 30 degrees, with the C4 atom 3 A from the flavin N5 atom.

Cooperation of the DnaK and GroE chaperone systems in the folding pathway of plant ferredoxin-NADP+ reductase expressed in Escherichia coli

The DnaK system is required for the productive folding of pea chloroplast ferredoxin-NADP+ reductase (FNR) expressed in Escherichia coli. The formation of a mature active enzyme was severely impaired in E. coli dnaK, dnaJ or grpE mutants expressing either the cytosolic precursor of the reductase (preFNR) or the mature apoenzyme, and these forms aggregated extensively in these cells. Coexpression of dnaK from a multicopy plasmid in the dnaK-null mutants restored preFNR processing and folding of FNR, rendering a mature-sized active enzyme. Overexpression of GroESL chaperonins failed to prevent preFNR aggregation, but it restored productive folding of FNR in dnaK-null mutants expressing the mature enzyme. Expression of preFNR in OmpT-protease-deficient E. coli cells resulted in the accumulation of the unprocessed precursor in the soluble fraction of the cells. The interaction of this soluble preFNR, but not the mature reductase, with DnaK and GroEL was evidenced by immunoprecipitation studies. We conclude that, in addition to the GroE chaperonins [Carrillo, N., Ceccarelli, E. A., Krapp, A. R., Boggio, S., Ferreyra, R. G. & Viale, A. M. (1992) J. Biol. Chem. 267, 15537-15541], the DnaK chaperone system plays a crucial role in the folding pathway of FNR.

Plant-type ferredoxin-NADP+ reductases: a basal structural framework and a multiplicity of functions

Ferredoxin-NADP+ (oxido)reductase (EC 1.18.1.2, FNR) is an FAD-containing enzyme that catalyzes the reversible electron transfer between NADP(H) and electron carrier proteins such as ferredoxin and flavodoxin. Isoforms of this flavoprotein are present in chloroplasts, mitochondria, and bacteria in which they participate in a wide variety of redox metabolic pathways. Although ferredoxin-NADP+ reductases have been thoroughly investigated and their properties reviewed on several occasions, considerable advances in the understanding of these flavoenzymes have occurred in the last few years, including the characterization of cDNA and genomic clones encoding FNR proteins from plants, algae, vertebrates, and bacteria, determination of the atomic structure of a plant FNR at high resolution, and the expression of functional reductases in microorganisms like Escherichia coli and Saccharomyces cerevisiae. The aim of this article is to summarize information gained through these recent developments, including the phylogenetic relationships among ferredoxin reductases and the key structural features of the plant FNR family. Other aspects such as the catalytic mechanism of FNR and the molecular events underlying biogenesis, intracellular sorting, folding, and holoenzyme assembly of this important flavoenzyme are also discussed in some detail. Ferredoxin-NADP+ reductases display several outstanding properties that make them excellent model proteins to address broad biological questions.

Expression, assembly and secretion of a fully active plant ferredoxin-NADP+ reductase by Saccharomyces cerevisiae

The flavoprotein ferredoxin-NADP+ reductase catalyzes the final step of the photosynthetic electron transport i.e., the reduction of NADP+ by ferredoxin. Expression and secretion of this enzyme was examined in Saccharomyces cerevisiae using a cDNA cloned from a pea library [Newman, B. J. & Gray, J. C. (1988) Plant Mol. Biol. 10, 511-520]. Two pea library cDNA sequences were employed, one corresponding to the mature enzyme and the other containing, in addition, the sequence of the transit peptide that directs ferredoxin-NADP+ reductase to the chloroplast. These sequences were introduced into a yeast shuttle vector in frame with the mating factor alpha 1 secretion-signal coding region under the control of its natural mating factor alpha 1 promoter. Saccharomyces cerevisiae cells transformed with the recombinant plasmids were able to synthesize and secrete fully active pea ferredoxin-NADP+ reductase. In both cases, a 35-kDa polypeptide was the major product. N-terminal sequencing of the secreted proteins indicates processing at position -1 with respect to the N-terminus of the pea mature enzyme. Yeast cells transformed with plasmid encoding the ferredoxin-NADP+ reductase precursor secrete four-times more ferredoxin-NADP+ reductase to the medium than cells transformed with the plasmid encoding the mature form of the enzyme. Ferredoxin-NADP+ reductases purified from culture medium showed structural and enzymatic properties that were identical, within the experimental error, to those of native plant ferredoxin-NADP+ reductase. The overall results indicate that pea ferredoxin-NADP+ reductase can be properly folded and its prosthetic group assembled in the yeast endoplasmic reticulum, and that its natural transit peptide favors its secretion.

Evolutionary relationships among eubacterial groups as inferred from GroEL (chaperonin) sequence comparisons

The essential GroEL proteins represent a subset of molecular chaperones ubiquitously distributed among species of the eubacterial lineage, as well as in eukaryote organelles. We employed these highly conserved proteins to infer eubacterial phylogenies. GroEL from the species analyzed clustered in distinct groups in evolutionary trees drawn by either the distance or the parsimony method, which were in general agreement with those found by 16S rRNA comparisons (i.e., proteobacteria, chlamydiae, bacteroids, spirochetes, firmicutes [gram-positive bacteria], and cyanobacteria-chloroplasts). Moreover, the analysis indicated specific relationships between some of the aforementioned groups which appeared not to be clearly defined or controversial in rRNA-based phylogenetic studies. For instance, a monophyletic origin for the low-G+C and high-G+C subgroups among the firmicutes, as well as their specific relationship to the cyanobacteria-chloroplasts, was inferred. The general observations suggest that GroEL proteins provide valuable evolutionary tools for defining evolutionary relationships among the eubacterial lineage of life.

The chaperone connection to the origins of the eukaryotic organelles

The heat-shock 60 proteins (Hsp60) constitute a subset of molecular chaperones essential for the survival of the cell, present in eubacteria as well as in eukaryotic organelles. Here, we have employed these highly conserved proteins for the inferences of the origins of the organelles. Hsp60s present in mitochondria from different eukaryotic lineages formed a clade, which showed the closest relationship to that of the Ehrlichia/Rickettsia cluster among the alpha-Proteobacteria. This, in addition to phenotypic characteristics, suggests that these obligate intracellular parasites and the lineage that generated the mitochondrion shared last common ancestry. In turn, Hsp60s present in chloroplasts from plants and a red alga, respectively, clustered specifically with those of the cyanobacteria, suggesting that all plastids derive exclusively from this eubacterial lineage.