How universal are fold recognition parameters. A comprehensive study of alignment and scoring function parameters influence on recognition of distant folds.

A choice of sequence-structure similarity scoring function parameters can significantly alter results of the performance of the recognition of distantly related folds. It therefore constitutes a critical part of fold recognition process. In order to increase an understanding of the influence of parameter choice, a comprehensive benchmark of very hard (SFOLD) and medium hard (SFAM) fold recognition examples has been derived from the SCOP database of protein structure families. These benchmarks have subsequently been used to optimize, validate and analyze dependence of recognition sensitivity on alignment and fold similarity score parameters for different scoring functions. Significant variation of the common parameters has been observed for different functions, leading to the conclusion that optimal parameter sets are not universal. The scope of solutions common to any pair of scoring function is relatively small, hence, using jury method for fold prediction seems not appropriate. Also, using a redundant version of fold libraries significantly increases odds of identification of distantly related fold.

From fold recognition to homology modeling: an analysis of protein modeling challenges at different levels of prediction complexity.

An analysis of different approaches to protein structure prediction is presented based solely on the range of models submitted to the third Critical Assessment of Protein Structure Prediction (CASP3) conference. CASP conferences evaluate the current state of the art of protein structure prediction by comparing blind prediction efforts of many groups for the same set of target sequences. Target sequences may be highly similar to those with known structure or can be totally (at least superficially) sequentially dissimilar. Techniques applied to those blind predictions (over 40 targets) ranges from a detailed homology prediction to the detection of remote homologues well below a twilight zone of protein sequence similarity. For the CASP3 conference, we have submitted predictions, totaling 35, with various levels of difficulty and complexity. For ten submitted homology targets, eight of them were determined by experiment so far. The RMSD of C-alpha atoms are 1.2-1.7, 2.3, and 4.6-17.9 A for the three easy targets, two hard targets, and three very hard homology targets, respectively. Out of 18-fold recognition predictions available for analysis, we got six correct predictions, five near misses, three tough near misses and four far misses. Here we analyze successes and failures of those predictions in an attempt to identify common problems and common achievements.

Folding simulations and computer redesign of protein A three-helix bundle motifs.

In solution, the B domain of protein A from Staphylococcus aureus (B domain) possesses a three-helix bundle structure. This simple motif has been previously reproduced by Kolinski and Skolnick (Proteins 18: 353-366, 1994) using a reduced representation lattice model of proteins with a statistical interaction scheme. In this paper, an improved version of the potential has been used, and the robustness of this result has been tested by folding from the random state a set of three-helix bundle proteins that are highly homologous to the B domain of protein A. Furthermore, an attempt to redesign the B domain native structure to its topological mirror image fold has been made by multiple mutations of the hydrophobic core and the turn region between helices I and II. A sieve method for scanning a large set of mutations to search for this desired property has been proposed. It has been shown that mutations of native B domain hydrophobic core do not introduce significant changes in the protein motif. Mutations in the turn region were also very conservative; nevertheless, a few mutants acquired the desired topological mirror image motif. A set of all atom models of the most probable mutant was reconstructed from the reduced models and refined using a molecular dynamics algorithm in the presence of water. The packing of all atom structures obtained corroborates the lattice model results. We conclude that the change in the handedness of the turn induced by the mutations, augmented by the repacking of hydrophobic core and the additional burial of the second helix N-cap side chain, are responsible for the predicted preferential adoption of the mirror image structure.

Does a backwardly read protein sequence have a unique native state?

Amino acid sequences of native proteins are generally not palindromic. Nevertheless, the protein molecule obtained as a result of reading the sequence backwards, i.e. a retro-protein, obviously has the same amino acid composition and the same hydrophobicity profile as the native sequence. The important questions which arise in the context of retro-proteins are: does a retro-protein fold to a well defined native-like structure as natural proteins do and, if the answer is positive, does a retro-protein fold to a structure similar to the native conformation of the original protein? In this work, the fold of retro-protein A, originated from the retro-sequence of the B domain of Staphylococcal protein A, was studied. As a result of lattice model simulations, it is conjectured that the retro-protein A also forms a three-helix bundle structure in solution. It is also predicted that the topology of the retro-protein A three-helix bundle is that of the native protein A, rather than that corresponding to the mirror image of native protein A. Secondary structure elements in the retro-protein do not exactly match their counterparts in the original protein structure; however, the amino acid side chain contract pattern of the hydrophobic core is partly conserved.