Cellular functions are determined by interactions among proteins in the cells. Recognition of these interactions forms an important step in understanding biology at the systems level. Here, we report an interaction network of Escherichia coli, obtained by training a Support Vector Machine on the high quality of interactions in the EcoCyc database, and with the assumption that the periplasmic and cytoplasmic proteins may not interact with each other. The data features included correlation coefficient between bit score phylogenetic profiles, frequency of their co-occurrence in predicted operons, and a new measure-the distance between translational start sites of the genes. The combined genome context methods show a high accuracy of prediction on the test data and predict a total of 78,122 binary interactions. The majority of the interactions identified by high-throughput experimental methods correspond to indirect interaction (interactions through neighbors) in the predicted network. Correlation of the predicted network with the gene essentiality data shows that the essential genes in E. coli exhibit a high linking number, whereas the nonessential genes exhibit a low linking number. Furthermore, our predicted protein-protein interaction network shows that the proteins involved in replication, DNA repair, transcription, translation, and cell wall synthesis are highly connected. We therefore believe that our predicted network will serve as a useful resource in understanding prokaryotic biology.

Abstract. The chaperonins, GroEL and GroES, are present ubiquitously and provide a paradigm in the understanding of assisted protein folding. Due to its essentiality of function, GroEL exhibits high sequence conservation across species. Complete genome sequencing has shown the occurrence of duplicate or multiple copies of groEL genes in bacteria such as Mycobacterium tuberculosis and Corynebacterium glutamicum. Monophyly of each bacterial clade in the phylogenetic tree generated for the GroEL protein suggests a lineage-specific duplication. The duplicated groEL gene in Actinobacteria is not accompanied by the operonic groES despite the presence of upstream regulatory elements. Our analysis suggests that in these bacteria the duplicated groEL genes have undergone rapid evolution and divergence to function in a GroES-independent manner. Evaluation of multiple sequence alignment demonstrates that the duplicated genes have acquired mutations at functionally significant positions including those involved in substrate binding, ATP binding, and GroES binding and those involved in inter-ring and intra-ring interactions. We propose that the duplicate groEL genes in different bacterial clades have evolved independently to meet specific requirements of each clade. We also propose that the groEL gene, although essential and conserved, accumulates nonconservative substitutions to exhibit structural and functional variations.

High throughput structural genomics efforts have been making the structures of proteins available even before their function has been fully characterized. Therefore, methods that exploit the structural knowledge to provide evidence about the functions of proteins would be useful. Such methods would be needed to complement the sequence-based function annotation approaches. The current study describes generation of 3D-structural motifs for metal-binding sites from the known metalloproteins. It then scans all the available protein structures in the PDB database for putative metal-binding sites. Our analysis predicted more than 1000 novel metal-binding sites in proteins using three-residue templates, and more than 150 novel metal-binding sites using four-residue templates. Prediction of metal-binding site in a yeast protein YDR533c led to the hypothesis that it might function as metal-dependent amidopeptidase. The structural motifs identified by our method present novel metal-binding sites that reveal newer mechanisms for a few well-known proteins.

PAR-3D (http://sunserver.cdfd.org.in:8080/protease/PAR_3D/index.html) is a web-based tool that exploits the fact that relative juxtaposition of active site residues is a conserved feature in functionally related protein families. The server uses previously calculated and stored values of geometrical parameters of a set of known proteins (training set) for prediction of active site residues in a query protein structure. PAR-3D stores motifs for different classes of proteases, the ten glycolytic pathway enzymes and metal-binding sites. The server accepts the structures in the pdb format. The first step during the prediction is the extraction of probable active site residues from the query structure. Spatial arrangement of the probable active site residues is then determined in terms of geometrical parameters. These are compared with stored geometries of the different motifs. Its speed and efficiency make it a beneficial tool for structural genomics projects, especially when the biochemical function of the protein has not been characterized.