Deciphering globular protein sequence–structure relationships: from observation to prediction

Careful comparison of proteins sharing a same fold but only low or no sequence identity should allow a better understanding of the coding of three-dimensional structures by amino acid sequences. It has already been shown that positions of a given fold occupied mainly by hydrophobic residues in the different proteins of a structural family share very specific physical properties and participate in stabilization of the protein domain. They probably also play a crucial role in the very first steps of folding [ Poupon A, Mornon J.-P (1999) FEBS Lett. 452: 283–289; Mirny LA, Shaknovich EI (1999) J. Mol. Biol. 291: 177–196]. To further understand the sequence–structure relationship, we studied the correlation between allowed mutations at a given three-dimensional position and some of its physical properties. The different amino acids were divided in three groups (hydrophobic, nonpolar or weakly polar and polar or charged), and a correlation was established between the occupation rate of each group at a given position in the fold and the burying, the side-chain dispersion, the interposition distances and the ability to form a network of directly interacting residues. The results are then applied to predict some solvent accessibility. We show that this property can be accurately predicted for about 70% of the residues, providing precious information concerning the corresponding three-dimensional structures. The results are used to predict other structural features, as secondary structures, compactness or long-range interactions between residues remote in sequence. This information will allow the number of possible structures for a given sequence to be reduced considerably, simplifying the ab initio modelling problem to a level where it might be solved by computing methods.