Protein sequence and structure alignments within one framework

An efficient algorithm for protein structure comparison using elastic shape analysis

BACKGROUND

Protein structure comparison play important role in in silico functional prediction of a new protein. It is also used for understanding the evolutionary relationships among proteins. A variety of methods have been proposed in literature for comparing protein structures but they have their own limitations in terms of accuracy and complexity with respect to computational time and space. There is a need to improve the computational complexity in comparison/alignment of proteins through incorporation of important biological and structural properties in the existing techniques.

RESULTS

An efficient algorithm has been developed for comparing protein structures using elastic shape analysis in which the sequence of 3D coordinates atoms of protein structures supplemented by additional auxiliary information from side-chain properties are incorporated. The protein structure is represented by a special function called square-root velocity function. Furthermore, singular value decomposition and dynamic programming have been employed for optimal rotation and optimal matching of the proteins, respectively. Also, geodesic distance has been calculated and used as the dissimilarity score between two protein structures. The performance of the developed algorithm is tested and found to be more efficient, i.e., running time reduced by 80-90 % without compromising accuracy of comparison when compared with the existing methods. Source codes for different functions have been developed in R. Also, user friendly web-based application called ProtSComp has been developed using above algorithm for comparing protein 3D structures and is accessible free.

CONCLUSIONS

The methodology and algorithm developed in this study is taking considerably less computational time without loss of accuracy (Table 2). The proposed algorithm is considering different criteria of representing protein structures using 3D coordinates of atoms and inclusion of residue wise molecular properties as auxiliary information.

The crystal structure of archaeal serine hydroxymethyltransferase reveals idiosyncratic features likely required to withstand high temperatures

Serine hydroxymethyltransferases (SHMTs) play an essential role in one-carbon unit metabolism and are used in biomimetic reactions. We determined the crystal structure of free (apo) and pyridoxal-5'-phosphate-bound (holo) SHMT from Methanocaldococcus jannaschii, the first from a hyperthermophile, from the archaea domain of life and that uses H₄MPT as a cofactor, at 2.83 and 3.0 Å resolution, respectively. Idiosyncratic features were observed that are likely to contribute to structure stabilization. At the dimer interface, the C-terminal region folds in a unique fashion with respect to SHMTs from eubacteria and eukarya. At the active site, the conserved tyrosine does not make a cation-π interaction with an arginine like that observed in all other SHMT structures, but establishes an amide-aromatic interaction with Asn257, at a different sequence position. This asparagine residue is conserved and occurs almost exclusively in (hyper)thermophile SHMTs. This led us to formulate the hypothesis that removal of frustrated interactions (such as the Arg-Tyr cation-π interaction occurring in mesophile SHMTs) is an additional strategy of adaptation to high temperature. Both peculiar features may be tested by designing enzyme variants potentially endowed with improved stability for applications in biomimetic processes.

Fast alignment and comparison of RNA structures

MOTIVATION

To recognize remote relationships between RNA molecules, one must be able to align structures without regard to sequence similarity. We have implemented a method, which is swift [O(n(2))], sensitive and tolerant of large gaps and insertions. Molecules are broken into overlapping fragments, which are characterized by their memberships in a probabilistic classification based on local geometry and H-bonding descriptors. This leads to a probabilistic similarity measure that is used in a conventional dynamic programming method.

RESULTS

Examples are given of database searching, the detection of structural similarities, which would not be found using sequence based methods, and comparisons with a previously published approach.

AVAILABILITY AND IMPLEMENTATION

Source code (C and perl) and binaries for linux are freely available at www.zbh.uni-hamburg.de/fries.

Elucidation of haem-binding sites in the actinobacterial protein HbpS

The extracellular haem-binding protein from Streptomyces reticuli (HbpS) has been shown to be involved in redox sensing and to bind haem. However, the residues involved in haem coordination are unknown. Structural alignments to distantly related haem-binding proteins from Mycobacterium tuberculosis were used to identify a candidate haem-coordinating residue, and site-directed mutagenesis with UV/Vis spectroscopy was used to assess haem binding in vivo and in vitro. We present strong evidence that HbpS belongs to the small set of proteins, which do not use histidine to coordinate the metal in the haem group. Further spectroscopic evidence strongly indicates that threonine 113 is actively involved in coordination of haem. Subsequent protein/haem titration experiments show a 1 : 2, protein/haem stoichiometry. We also present data showing the degradation of haem by HbpS in vivo. Because HbpS is conserved in many Actinobacteria, the presented results are applicable to related species.

Incorporation of local structural preference potential improves fold recognition

Fold recognition, or threading, is a popular protein structure modeling approach that uses known structure templates to build structures for those of unknown. The key to the success of fold recognition methods lies in the proper integration of sequence, physiochemical and structural information. Here we introduce another type of information, local structural preference potentials of 3-residue and 9-residue fragments, for fold recognition. By combining the two local structural preference potentials with the widely used sequence profile, secondary structure information and hydrophobic score, we have developed a new threading method called FR-t5 (fold recognition by use of 5 terms). In benchmark testings, we have found the consideration of local structural preference potentials in FR-t5 not only greatly enhances the alignment accuracy and recognition sensitivity, but also significantly improves the quality of prediction models.

Structural alphabets derived from attractors in conformational space

Background

The hierarchical and partially redundant nature of protein structures justifies the definition of frequently occurring conformations of short fragments as 'states'. Collections of selected representatives for these states define Structural Alphabets, describing the most typical local conformations within protein structures. These alphabets form a bridge between the string-oriented methods of sequence analysis and the coordinate-oriented methods of protein structure analysis.

Results

A Structural Alphabet has been derived by clustering all four-residue fragments of a high-resolution subset of the protein data bank and extracting the high-density states as representative conformational states. Each fragment is uniquely defined by a set of three independent angles corresponding to its degrees of freedom, capturing in simple and intuitive terms the properties of the conformational space. The fragments of the Structural Alphabet are equivalent to the conformational attractors and therefore yield a most informative encoding of proteins. Proteins can be reconstructed within the experimental uncertainty in structure determination and ensembles of structures can be encoded with accuracy and robustness.

Conclusions

The density-based Structural Alphabet provides a novel tool to describe local conformations and it is specifically suitable for application in studies of protein dynamics.

The SALAMI protein structure search server

Protein structures often show similarities to another which would not be seen at the sequence level. Given the coordinates of a protein chain, the SALAMI server atwww.zbh.uni-hamburg.de/salami will search the protein data bank and return a set of similar structures without using sequence information. The results page lists the related proteins, details of the sequence and structure similarity and implied sequence alignments. Via a simple structure viewer, one can view superpositions of query and library structures and finally download superimposed coordinates. The alignment method is very tolerant of large gaps and insertions, and tends to produce slightly longer alignments than other similar programs.