Detection of distant structural similarities in a set of proteins using a fast graph-based method.
We introduce a method for finding weak structural similarities in a set of protein structures. Proteins are considered at their secondary structure level. The method uses a rigorous graph-theoretical algorithm which finds all structural similarities. Protein structures are modelled as undirected labelled graphs, the so-called protein graphs. We suggest that for detecting the similarities between two protein structures it is sufficient to find similarities in the protein core which consists of tightly packed secondary structure elements. Therefore, we can restrict ourselves to solving the maximal common connected subgraph problem instead of the maximal common subgraph problem. We have modified the algorithm by Bron and Kerbosch for solving that problem. The speed of the algorithm increases drastically. After calculating all maximal common connected substructures for all pairwise comparisons in a set of protein graphs the common substructure in all proteins can be calculated by intersecting them. In this paper we characterize the method briefly and explain the modelling of the protein structure in detail. For the pairwise alignment the similarity of porin (1OMF) with bacteriochlorophyll a (3BCL) and BirA protein (1BIB) with DNA polymerase III (2POL) will be discussed. In the case of the multiple structure alignment the similarity in variants of four phosphatases and in subtilisin Carlsberg, carboxypeptidase, elongation factor Tu, and flavodoxin will be represented. Our first experiments show that the method works correctly and fast. The method can be used for arbitrary graphs. Thus, different graph-theoretical models of protein structures can be examined.
Koch
I
I
Institute for Algorithms and Scientific Computing, GMD-German National Research Center for Information Technology, St. Augustin, Germany. koch@pluto.chem.uni-potsdam.de
Lengauer
T
T
D000465:Algorithms; D001426:Bacterial Proteins / Q000737:chemistry / Q000235:genetics; D001429:Bacteriochlorophylls / Q000737:chemistry / Q000235:genetics; D019731:Carbon-Nitrogen Ligases / Q000737:chemistry / Q000235:genetics; D003198:Computer Simulation; D004258:DNA Polymerase III / Q000737:chemistry / Q000235:genetics; D016208:Databases, Factual; D029968:Escherichia coli Proteins; D008956:Models, Chemical; D008958:Models, Molecular; D010744:Phosphoric Monoester Hydrolases / Q000737:chemistry / Q000235:genetics; D018272:Porins / Q000737:chemistry / Q000235:genetics; D011487:Protein Conformation; D017433:Protein Structure, Secondary; D012097:Repressor Proteins; D016415:Sequence Alignment / Q000379:methods / Q000706:statistics & numerical data; D014157:Transcription Factors
D001426:Bacterial Proteins; D001429:Bacteriochlorophylls; D029968:Escherichia coli Proteins; D018272:Porins; D012097:Repressor Proteins; D014157:Transcription Factors; D004258:DNA Polymerase III; D010744:Phosphoric Monoester Hydrolases; D019731:Carbon-Nitrogen Ligases; C058125:birA protein, E coli