Structure-based functional classification of hypothetical protein MTH538 from Methanobacterium thermoautotrophicum

Characterization of a thaumarchaeal symbiont that drives incomplete nitrification in the tropical sponge Ianthella basta

Marine sponges represent one of the few eukaryotic groups that frequently harbour symbiotic members of the Thaumarchaeota, which are important chemoautotrophic ammonia-oxidizers in many environments. However, in most studies, direct demonstration of ammonia-oxidation by these archaea within sponges is lacking, and little is known about sponge-specific adaptations of ammonia-oxidizing archaea (AOA). Here, we characterized the thaumarchaeal symbiont of the marine sponge Ianthella basta using metaproteogenomics, fluorescence in situ hybridization, qPCR and isotope-based functional assays. 'Candidatus Nitrosospongia ianthellae' is only distantly related to cultured AOA. It is an abundant symbiont that is solely responsible for nitrite formation from ammonia in I. basta that surprisingly does not harbour nitrite-oxidizing microbes. Furthermore, this AOA is equipped with an expanded set of extracellular subtilisin-like proteases, a metalloprotease unique among archaea, as well as a putative branched-chain amino acid ABC transporter. This repertoire is strongly indicative of a mixotrophic lifestyle and is (with slight variations) also found in other sponge-associated, but not in free-living AOA. We predict that this feature as well as an expanded and unique set of secreted serpins (protease inhibitors), a unique array of eukaryotic-like proteins, and a DNA-phosporothioation system, represent important adaptations of AOA to life within these ancient filter-feeding animals.

Structure of the archaeal chemotaxis protein CheY in a domain-swapped dimeric conformation

Archaea are motile by the rotation of the archaellum. The archaellum switches between clockwise and counterclockwise rotation, and movement along a chemical gradient is possible by modulation of the switching frequency. This modulation involves the response regulator CheY and the archaellum adaptor protein CheF. In this study, two new crystal forms and protein structures of CheY are reported. In both crystal forms, CheY is arranged in a domain-swapped conformation. CheF, the protein bridging the chemotaxis signal transduction system and the motility apparatus, was recombinantly expressed, purified and subjected to X-ray data collection.

Exploring dynamics of protein structure determination and homology-based prediction to estimate the number of superfamilies and folds

Background

As tertiary structure is currently available only for a fraction of known protein families, it is important to assess what parts of sequence space have been structurally characterized. We consider protein domains whose structure can be predicted by sequence similarity to proteins with solved structure and address the following questions. Do these domains represent an unbiased random sample of all sequence families? Do targets solved by structural genomic initiatives (SGI) provide such a sample? What are approximate total numbers of structure-based superfamilies and folds among soluble globular domains?

Results

To make these assessments, we combine two approaches: (i) sequence analysis and homology-based structure prediction for proteins from complete genomes; and (ii) monitoring dynamics of the assigned structure set in time, with the accumulation of experimentally solved structures. In the Clusters of Orthologous Groups (COG) database, we map the growing population of structurally characterized domain families onto the network of sequence-based connections between domains. This mapping reveals a systematic bias suggesting that target families for structure determination tend to be located in highly populated areas of sequence space. In contrast, the subset of domains whose structure is initially inferred by SGI is similar to a random sample from the whole population. To accommodate for the observed bias, we propose a new non-parametric approach to the estimation of the total numbers of structural superfamilies and folds, which does not rely on a specific model of the sampling process. Based on dynamics of robust distribution-based parameters in the growing set of structure predictions, we estimate the total numbers of superfamilies and folds among soluble globular proteins in the COG database.

Conclusion

The set of currently solved protein structures allows for structure prediction in approximately a third of sequence-based domain families. The choice of targets for structure determination is biased towards domains with many sequence-based homologs. The growing SGI output in the future should further contribute to the reduction of this bias. The total number of structural superfamilies and folds in the COG database are estimated as ~4000 and ~1700. These numbers are respectively four and three times higher than the numbers of superfamilies and folds that can currently be assigned to COG proteins.

Automated prediction of protein function and detection of functional sites from structure

Current structural genomics projects are yielding structures for proteins whose functions are unknown. Accordingly, there is a pressing requirement for computational methods for function prediction. Here we present PHUNCTIONER, an automatic method for structure-based function prediction using automatically extracted functional sites (residues associated to functions). The method relates proteins with the same function through structural alignments and extracts 3D profiles of conserved residues. Functional features to train the method are extracted from the Gene Ontology (GO) database. The method extracts these features from the entire GO hierarchy and hence is applicable across the whole range of function specificity. 3D profiles associated with 121 GO annotations were extracted. We tested the power of the method both for the prediction of function and for the extraction of functional sites. The success of function prediction by our method was compared with the standard homology-based method. In the zone of low sequence similarity (approximately 15%), our method assigns the correct GO annotation in 90% of the protein structures considered, approximately 20% higher than inheritance of function from the closest homologue.

Thermophile-specific proteins: the gene product of aq_1292 from Aquifex aeolicus is an NTPase

Background

To identify thermophile-specific proteins, we performed phylogenetic patterns searches of 66 completely sequenced microbial genomes. This analysis revealed a cluster of orthologous groups (COG1618) which contains a protein from every thermophile and no sequence from 52 out of 53 mesophilic genomes. Thus, COG1618 proteins belong to the group of thermophile-specific proteins (THEPs) and therefore we here designate COG1618 proteins as THEP1s. Since no THEP1 had been analyzed biochemically thus far, we characterized the gene product of aq_1292 which is THEP1 from the hyperthermophilic bacterium Aquifex aeolicus (aaTHEP1).

Results

aaTHEP1 was cloned in E. coli, expressed and purified to homogeneity. At a temperature optimum between 70 and 80°C, aaTHEP1 shows enzymatic activity in hydrolyzing ATP to ADP + P_i with k_cat = 5 × 10^-3 s^-1 and K_m = 5.5 × 10^-6 M. In addition, the enzyme exhibits GTPase activity (k_cat = 9 × 10^-3 s^-1 and K_m= 45 × 10^-6 M). aaTHEP1 is inhibited competitively by CTP, UTP, dATP, dGTP, dCTP, and dTTP. As shown by gel filtration, aaTHEP1 in its purified state appears as a monomer. The enzyme is resistant to limited proteolysis suggesting that it consists of a single domain. Although THEP1s are annotated as "predicted nucleotide kinases" we could not confirm such an activity experimentally.

Conclusion

Since aaTHEP1 is the first member of COG1618 that is characterized biochemically and functional information about one member of a COG may be transferred to the entire COG, we conclude that COG1618 proteins are a family of thermophilic NTPases.

Annotation in three dimensions. PINTS: Patterns in Non-homologous Tertiary Structures

The detection of local structural patterns in proteins (e.g. active sites) can provide insights into protein function in the absence of sequence or fold similarity. Methods to detect such similarities are key during structural annotation, for example with results from Structural Genomics initiatives. PINTS (Patterns in Non-homologous Tertiary Structures, http://pints.embl.de) performs database searches for such patterns and most importantly provides a measure of statistical significance for any similarity uncovered. To aid functional annotation of proteins, we allow comparisons of pre-defined patterns against databases of complete structures and of entire structures to databases of particular residues likely to be functionally important.

Structural proteomics: lessons learnt from the early case studies

The genomics efforts have identified a large number of novel genes and thus provided a pool of interesting but not functionally characterized target proteins. It has been suggested that structural proteomics will significantly impact the success rate of functional characterization of such identified genes and proteins by providing structure-function hypotheses by fold and feature recognition and analysis. Structural proteomics initiatives, both in academic and industrial settings, are today generating protein structures at an unprecedented rate although relatively few large-scale efforts have been displayed in the public domain. However, a number of individual studies have provided a 'road-map' for selected approaches that hold the promise to significantly impact the process of deriving function from structure.

Breaking the singleton of germination protease

Germination protease (GPR) plays an important role in the germination of spores of Bacillus and Clostridium species. A few very similar GPRs form a singleton group without significant sequence similarities to any other proteins. Their active site locations and catalytic mechanisms are unclear, despite the recent 3-D structure determination of Bacillus megaterium GPR. Using structural comparison and sequence analysis, we show that GPR is homologous to bacterial hydrogenase maturation protease (HybD). HybD's activity relies on the recognition and binding of metal ions in Ni-Fe hydrogenase, its substrate. Two highly conserved motifs are shared among GPRs, hydrogenase maturation proteases, and another group of hypothetical proteins. Conservation of two acidic residues in all these homologs indicates that metal binding is important for their function. Our analysis helps localize the active site of GPRs and provides insight into the catalytic mechanisms of a superfamily of putative metal-regulated proteases.

Structural proteomics: developments in structure-to-function predictions

The major challenge for post-genomic research is to functionally assign and validate a large number of novel target genes and their corresponding proteins. Functional genomics approaches have, therefore, gained considerable attention in the quest to convert this massive data set into useful information. One of the crucial components for the functional understanding of unassigned proteins is the analysis of their experimental or modeled 3D structures. Structural proteomics initiatives are generating protein structures at an unprecedented rate but our current knowledge of 3D-structural space is still limited. Estimates on the completeness of the 3D-structural coverage of proteins vary but it is generally accepted that only a minority of the structural proteome has a template structure from which reliable conclusions can be drawn. Thus, structural proteomics has set out to build a map of protein structures that will represent all protein folds included in the 'global proteome'.

A tour of structural genomics

Structural genomics projects aim to provide an experimental or computational three-dimensional model structure for all of the tractable macromolecules that are encoded by complete genomes. To this end, pilot centres worldwide are now exploring the feasibility of large-scale structure determination. Their experimental structures and computational models are expected to yield insight into the molecular function and mechanism of thousands of proteins. The pervasiveness of this information is likely to change the use of structure in molecular biology and biochemistry.

Prediction of functionally important residues based solely on the computed energetics of protein structure

Catalytic and other functionally important residues in proteins can often be mutated to yield more stable proteins. Many of these residues are charged residues that are located in electrostatically unfavorable environments. Here it is demonstrated that because continuum electrostatics methods can identify these destabilizing residues, the same methods can also be used to identify functionally important residues in otherwise uncharacterized proteins. To establish this point, detailed calculations are performed on six proteins for which good structural and mutational data are available from experiments. In all cases it is shown that functionally important residues known to be destabilizing experimentally are among the most destabilizing residues found in the calculations. A larger scale analysis performed on 216 different proteins demonstrates the existence of a general relationship between the calculated electrostatic energy of a charged residue and its degree of evolutionary conservation. This relationship becomes obscured when electrostatic energies are calculated using Coulomb's law instead of the more complete continuum electrostatics method. Finally, in a first predictive application of the method, calculations are performed on three proteins whose structures have recently been reported by a structural genomics consortium.

Solution structure and function of a conserved protein SP14.3 encoded by an essential Streptococcus pneumoniae gene

Streptococcus pneumoniae is a major human pathogen that causes high mortality and morbidity rates and has developed resistance to many antibiotics. The genome of S. pneumoniae has recently been completely sequenced revealing many genes encoding hypothetical proteins of unknown function. We have found that the gene encoding one such conserved protein, SP14.3, is essential for growth of S. pneumonia. Since it is essential, SP14.3 represents a potential target for drug discovery. Here, we describe the three-dimensional solution structure of SP14.3 as determined by NMR spectroscopy. The structure consists of two domains each with an alpha/beta-fold. The N-terminal domain contains two alpha-helices and a three-stranded beta-sheet, while the C-terminal domain is composed of one alpha-helix and a five-stranded beta-sheet. The N-terminal domain of the protein contains a highly negatively charged surface and resembles the fold of the N-terminal domain of Thermus thermophilus ribosomal protein S3. The C-terminal domain has a protein fold similar to human small nuclear ribonucleoprotein Sm D3 and Haloarcula marismortui ribosomal protein L21E. The two domains of the protein tumble in solution overall as a whole with an overall molecular rotational correlation time (tau(m)) of 12.9 ns at 25 degrees C. The relative orientation of the two domains is not defined by the nuclear Overhauser effect data. Indeed, residual dipolar couplings and the structure calculations indicate that the relative orientation of the two domains is not rigidly oriented with respect to one another in solution.

Structural genomics: opportunities and challenges

Following the complete genome sequencing of an increasing number of organisms, structural biology is engaging in a systematic approach of high-throughput structure determination called structural genomics to create a complete inventory of protein folds/structures that will help predict functions for all proteins. First results show that structural genomics will be highly effective in finding functional annotations for proteins of unknown function.

Determination of protein function, evolution and interactions by structural genomics

The genome sequencing projects and knowledge of the entire protein repertoires of many organisms have prompted new procedures and techniques for the large-scale determination of protein structure, function and interactions. Recently, new work has been carried out on the determination of the function and evolutionary relationships of proteins by experimental structural genomics, and the discovery of protein-protein interactions by computational structural genomics.