Statistical significance of protein structure prediction by threading

In this study, we estimate the statistical significance of structure prediction by threading. We introduce a single parameter ɛ that serves as a universal measure determining the probability that the best alignment is indeed a native-like analog. Parameter ɛ takes into account both length and composition of the query sequence and the number of decoys in threading simulation. It can be computed directly from the query sequence and potential of interactions, eliminating the need for sequence reshuffling and realignment. Although our theoretical analysis is general, here we compare its predictions with the results of gapless threading. Finally we estimate the number of decoys from which the native structure can be found by existing potentials of interactions. We discuss how this analysis can be extended to determine the optimal gap penalties for any sequence-structure alignment (threading) method, thus optimizing it to maximum possible performance.

Energy-based de novo protein folding by conformational space annealing and an off-lattice united-residue force field: Application to the 10-55 fragment of staphylococcal protein A and to apo calbindin D9K

The conformational space annealing (CSA) method for global optimization has been applied to the 10-55 fragment of the B-domain of staphylococcal protein A (protein A) and to a 75-residue protein, apo calbindin D9K (PDB ID code 1CLB), by using the UNRES off-lattice united-residue force field. Although the potential was not calibrated with these two proteins, the native-like structures were found among the low-energy conformations, without the use of threading or secondary-structure predictions. This is because the CSA method can find many distinct families of low-energy conformations. Starting from random conformations, the CSA method found that there are two families of low-energy conformations for each of the two proteins, the native-like fold and its mirror image. The CSA method converged to the same low-energy folds in all cases studied, as opposed to other optimization methods. It appears that the CSA method with the UNRES force field, which is based on the thermodynamic hypothesis, can be used in prediction of protein structures in real time.

An iterative method for extracting energy-like quantities from protein structures.

We present a method (ENERGI) for extracting energy-like quantities from a data base of protein structures. In this paper, we use the method to generate pairwise additive amino acid "energy" scores. These scores are obtained by iteration until they correctly discriminate a set of known protein folds from decoy conformations. The method succeeds in lattice model tests and in the gapless threading problem as defined by Maiorov and Crippen [Maiorov, V. N. & Crippen, G. M. (1992) J. Mol. Biol. 227, 876-888]. A more challenging test of threading a larger set of test proteins derived from the representative set of Hobohm and Sander [Hobohm, U. & Sander, C. (1994) Protein Sci. 3, 522-524] is used as a "workbench" for exploring how the ENERGI scores depend on their parameter sets.

Replication infidelity during a single cycle of Ty1 retrotransposition.

Retroviruses undergo a high frequency of genetic alterations during the process of copying their RNA genomes. However, little is known about the replication fidelity of other elements that transpose via reverse transcription of an RNA intermediate. The complete sequence of 29 independently integrated copies of the yeast retrotransposon Ty1 (173,043 nt) was determined, and the mutation rate during a single cycle of replication was calculated. The observed base substitution rate of 2.5 x 10(-5) bp per replication cycle suggests that this intracellular element can mutate as rapidly as retroviruses. The pattern and distribution of errors in the Ty1 genome is nonrandom and provides clues to potential in vivo molecular mechanisms of reverse transcriptase-mediated error generation, including heterogeneous RNase H cleavage of Ty1 RNA, addition of terminal nontemplated bases, and transient dislocation and realignment of primer-templates. Overall, analysis of errors generated during Ty1 replication underscores the utility of a genetically tractable model system for the study of reverse transcriptase fidelity.