Skip to main content

Residue contact-count potentials are as effective as residue-residue contact-type potentials for ranking protein decoys



For over 30 years potentials of mean force have been used to evaluate the relative energy of protein structures. The most commonly used potentials define the energy of residue-residue interactions and are derived from the empirical analysis of the known protein structures. However, single-body residue 'environment' potentials, although widely used in protein structure analysis, have not been rigorously compared to these classical two-body residue-residue interaction potentials. Here we do not try to combine the two different types of residue interaction potential, but rather to assess their independent contribution to scoring protein structures.


A data set of nearly three thousand monomers was used to compare pairwise residue-residue 'contact-type' propensities to single-body residue 'contact-count' propensities. Using a large and standard set of protein decoys we performed an in-depth comparison of these two types of residue interaction propensities. The scores derived from the contact-type and contact-count propensities were assessed using two different performance metrics and were compared using 90 different definitions of residue-residue contact. Our findings show that both types of score perform equally well on the task of discriminating between near-native protein decoys. However, in a statistical sense, the contact-count based scores were found to carry more information than the contact-type based scores.


Our analysis has shown that the performance of either type of score is very similar on a range of different decoys. This similarity suggests a common underlying biophysical principle for both types of residue interaction propensity. However, several features of the contact-count based propensity suggests that it should be used in preference to the contact-type based propensity. Specifically, it has been shown that contact-counts can be predicted from sequence information alone. In addition, the use of a single-body term allows for efficient alignment strategies using dynamic programming, which is useful for fold recognition, for example. These facts, combined with the relative simplicity of the contact-count propensity, suggests that contact-counts should be studied in more detail in the future.


Accurate descriptions of the different non-covalent interactions involved in protein folding and stability are essential for a number of related problems. Potential energy functions based on such terms have been widely used to facilitate: fold recognition [13], homology modelling [4, 5], docking [6], ab-initio structure prediction [79], sequence design [10] and the analysis of protein folding kinetics [11, 12]. In each case, the purpose of the potential function is to discriminate between a variety of alternative conformations, selecting the most energetically favourable (assumed to be the most native) for further analysis [13]. Different potential energy functions have been defined at different levels of structural resolution [14]. At the atomic level, various pairwise inter-atom potentials (force-fields) have been developed from the detailed analysis of small, protein-like compounds. These include: ECEPP [15, 16], MM [17, 18], AMBER [19, 20], CHARMM [2123] and GROMOS [24]. Potential functions between distinct groups of atoms have also been defined, typically between pairs of residues [8, 2528] or idealised elements of secondary structure [9, 2934]. These 'potentials of mean force' (mean-fields) have the nature of free energies [27, 35], and may be derived by conformational averaging [7] or, more commonly, by empirical methods as described below.

There are two commonly used methods for deriving empirical potential energy functions [36]. The first method employs a statistical analysis of the observed 'interactions' [8, 25, 26, 37, 38]. In this method, the observed occurrence of a particular interaction is weighted by its expected occurrence in a given reference state [27, 39]. The resulting statistical interaction propensities can be either converted into energies using the Boltzmann distribution [8, 25, 26, 38] or log-odds scores [40, 41]. However, it has been shown that these two types of propensity are essentially the same [36]. In the second method, a potential function can be directly optimised in order to discriminate between native and near-native (decoy) structures [42]. This technique resembles machine learning, and has been applied in a variety of different ways, usually by maximising the discrimination between an average decoy and the native structure [4346]. Either of the above two methods may be applied to any feature of the protein structure that can be parameterised [9]. In the current work, we focus on the statistical analysis of residue interaction propensities. Previously, a variety of different methods have been applied to derive empirical residue-residue interaction potentials, often yielding remarkably consistent results [27]. However, the physical basis of the empirically derived potentials remains ambiguous [47]. Specifically, it has been shown that protein structures are inconsistent with the assumptions that underlie the use of the Boltzmann distribution [28, 48].

The major criticism of empirical residue-residue interaction potentials is that they ignore the protein/solvent boundary [27, 28, 48]. Consequently, there is an apparent attractive force between residues that co-segregate into the protein surface or core regions [28]. To address this, several groups have developed residue-specific environment potentials. These residue-specific environment potentials are usually correlated with hydrophobicity, measuring the extent to which each residue is buried in the protein core. In this way these single-body environment potentials capture information about the protein/solvent boundary. Such potentials have been combined with residue-residue interaction potentials: as a 'solvent correction factor' [49, 50], as an ad-hoc repulsive term [38], and using a Bayesian framework to avoid over-counting [40].

The above combination of two-body, residue-residue interaction potentials with single-body, residue-specific environment potentials raises the question as to which type of potential is the most specific for the native protein structure. To address this question, we separated statistical residue interaction propensities into two different types of score: a two-body, residue-residue 'contact-type' score, and a single-body, residue 'contact-count' score.

These two types of score can be expected to capture qualitatively different kinds of residue interaction propensities. The resulting propensities can be understood in terms of biophysical properties of protein structure. For example, the contact-type score can encode the fact that hydrophobic residues tend to interact with other hydrophobic residues in preference to hydrophilic residues. In contrast, the contact-count score can encode the fact that bulky hydrophobic residues tend to have more residue-residue interactions than small hydrophilic residues.

Here we report a comparison of two-body, residue-residue 'contact-type' scores and single-body, residue 'contact-count' score, as described below.


Two different types of residue interaction propensity are studied here, contact-type and contact-count. The 'two-body' residue contact-type propensities are based on the distinct amino acid types of a pair of contacting residues. The 'single-body' residue contact-count propensities are based on the discrete number of residue-residue contacts made by each distinct residue type. These two different interaction propensities are captured by the contact-type and contact-count scoring matrices, respectively. An example contact-type scoring matrix is given in Table 1, and the scores for some residues in an example contact-count scoring matrix are shown in Figure 1. The scores in these matrices reflect the observed residue interaction propensities in a set of native structures, and are defined in comparison to simple, random models of residue interaction. (For details see the Methods Section.)

Table 1 An example of data from a contact-type scoring matrix
Figure 1

An example of data from a contact-count scoring matrix. An example of some data from one of the contact-count scoring matrices used in this work. The scores are defined as described in the Methods Section. In this matrix residue-residue contacts were defined using a 12 Å distance threshold between C β atoms without filtering for short-range interactions. Scores are shown for selected residues as trends across the range of observed 'number of contacts'. Missing values for a residue indicate cells of the matrix that were removed due to lack of data.

Both types of scoring matrix were constructed using several different definitions of residue-residue interaction. Three different structural criteria were used to define residue-residue interaction. Firstly, we tested the effect of the choice of atomic interaction site, representing each residue by either the C α atom or the C β atom, or both. Secondly, the contact distance threshold was varied between 6 and 20 Å in increments of 1 Å, giving a total of 15 different distance cutoffs. Thirdly, we applied a sequence separation filter, either considering all interactions or only the long-range interactions. Long range interactions were defined as interactions between residues that are more than 10 residues apart in the protein sequence [51, 52]. The combination of these criteria gave a total of 3 (C α , C β or both) ×15 (distance cutoffs) ×2 (all or long range contacts) = 90 different residue-residue contact definitions.

The following three sections present the different contact-type and contact-count scoring matrices. First, the scoring matrices themselves are described, as they provide information on the nature of the captured residue interaction propensities. Second, the results of scoring native and 'fully-random' protein structures are presented. Third, the matrices are used to evaluate several sets of protein decoy structures.

In summary the results show that; i) the contact-count scores are much more specific than the contact type scores compared to random models of residue-residue interaction, ii) the C β -C β interaction captures the most specific residue interaction information compared to other atomic interaction sites, iii) both scores can identify 'unusual' proteins in the training set, iv) in contrast to point i, both scores perform equally well on the task of discriminating between decoy structures. The apparent contradiction between point i) and iv) will be returned to in the Discussion.

The magnitude of the scoring matrices

The 'mean absolute score' of a scoring matrix (MAS) was defined as the mean of the absolute value of the score in each cell of the matrix. The magnitude of MAS gives the degree to which the observed residue interaction propensities deviate from random. In other words, MAS measures the 'information content' of the observed interaction propensities encoded in the scoring matrix. The value of MAS would be equal to zero if residue interactions occurred at random, i.e. without any particular interaction propensities. The mean absolute score for each different contact-type and contact-count scoring matrix is shown in Figure 2, and are described in detail below.

Figure 2

Magnitude of the contact-type/count scoring matrices. Each point gives the mean absolute score (MAS) of each cell in a particular residue interaction propensity scoring matrix. The different scoring matrices result from the different (given) residue-residue contact definitions used in matrix construction. The broken lines indicate the trend for the matrices with sequence separation filtering and the solid lines give the trend for the matrices without sequence separation filtering.


The sequence separation threshold has the biggest effect on the mean absolute score (MAS) of the contact-type scoring matrices (Figure 2). Without sequence separation filtering, the contact-type scoring matrices tend to have smaller values of MAS. This clearly shows the effect of including the inherently non-specific short-range contacts in the scoring matrix. The scoring matrices that include short-range contacts are 'more random', with respect to the observed contacts encoded in the matrix. A similar effect is seen with increasing contact distance threshold.


The values of the 'mean absolute scores' (MAS) of the contact-count matrices are consistently larger than those of the contact-type matrices (Figure 2). The number of elements in the contact-count scoring matrix may vary with the residue-residue contact definition used (Figure 3). However, the value of MAS is comparable between the different types of scoring matrix because MAS is the mean absolute score over all elements in the matrix. The comparison suggests that the 'number of contacts per residue type' is consistently more informative than the 'residue-residue contact-type', given any of the residue-residue contact definitions used here. Unlike the contact-type matrices, the contact-count matrices appear consistently more informative when short-range contacts are included.

Figure 3

Size of the contact-count scoring matrices. Each point gives the number of different 'number of contact-bins' in the different contact-count scoring matrices. The broken lines indicate the trend for the matrices with sequence separation filtering and the solid lines give the trend for the matrices without sequence separation filtering (as described in the text). The plot shows how the 'number of contact-bins' varies with the contact definition.


The 'Mean Absolute Score' shows that contact-count matrices are consistently more specific than the contact-type matrices. At all the distance thresholds, either with or without the sequence separation filtering, the values of MAS are largest for the C β -C β contacts, then the C α -C β contacts, then the C α -C α contacts. This shows that the C β atom captures the specific side chain interactions more accurately than either of the other two definitions.

Scoring native protein structures

Each of the scoring matrices was derived from a data set of 3, 070 monomers (see the Methods Section for details). As a simple test, each of the native proteins was scored using the contact-type and contact-count scoring matrices which had the highest value of MAS (as described above). Using either the contact-type or the contact-count scoring matrices, there were some proteins that scored significantly worse than average. Examination of the 130 worst cases showed that they were caused by a few anomalies and annotation errors.

There were 86 very small proteins and protein fragments. These included, for example, the structure of single C α -helices and extended, coiled-coil proteins. There were 25 membrane associated proteins, including alpha helical and beta-barrel lipoproteins. There were 12 proteins that adopted an extended conformation in complex with either DNA or several large ligand groups. Another four structures were found to be C α only models, containing only the backbone and no side-chain information.

Interestingly, in this group we found 3 structures of protein subunits from oligomeric proteins. These cases were incorrectly annotated monomers in the data set. These subunits appear 'non-native' because they would make many additional residue contacts in the native oligomer. For this reason the artificially isolated subunit is effectively 'non-native' and scored badly as a result.

An examination of 130 randomly selected proteins from the data set showed only a few protein fragments and DNA binding proteins. There were several proteins found binding large ligand groups, but the relative extent of the ligand was small compared to the cases found above. There were no trans-membrane structures found in the random sample.

For the above reasons, these 130 cases were removed from the data set giving a total of 2, 940 monomers. The matrices were re-calculated over this new data set for use in the following sections.

Scoring 'decoy' structures

In this section we describe a realistic benchmark of score performance [53] using several standard sets of 'near-native' protein decoys [54]. Here the scores are used to evaluate the decoys with reference to the C α RMSD of the decoy to its corresponding native structure. The C α RMSD is used as an independent measure of decoy quality in order to evaluate the various scores.

Description of the decoy sets used

Nine different sets of decoys were used in the current work. The structures of the decoys were taken from the Decoys-R-us database [54]. Each decoy set uses a particular method to generate several 'near-native' protein structures using a given native protein structure. Some additional information for the different decoy sets is given in Table 2. The different methods include: energy minimisation (lmds and vhp-mcmd), homology modelling (hg-structal, ig-structal, and ig-structal-hires), systematic randomisation with subsequent filtering (4-state-reduced and lattice-ssfit), ab-initio (semfold) and de-novo methods (fisa and fisa-casp3).

Table 2 Summary of the nine different decoys sets from decoys-R-us

The relationship between decoy 'quality' and score

When assessing the relationship between decoy 'quality' and the residue interaction propensity score, several different measures of score performance are important [53]. Here we apply two different measures of score performance, collected from the decoys as described below. The first measure is the Spearman rank correlation coefficient (S). The value of S shows whether the interaction propensity score can accurately discriminate between decoys of varying quality. The second measure is the Z-score of the native structure compared to the decoys with respect to interaction propensity score (Z). A large and positive Z indicates a clear discrimination of the native conformation from that of the decoys using the interaction propensity score.

The nine different decoy sets were analysed separately, and each decoy was scored using the different residue contact-type and contact-count scoring matrices as described above. The two different measures of score performance described above (S and Z) were calculated for each protein. For each decoy set we always report the mean value of S and Z over all the proteins in the set, given a particular residue-residue contact definition. In the following sections we refer to the 'best' score for a decoy set as the contact definition that had the best mean performance (on S or Z) over all the proteins in the set.

The best values of S per decoy set

Focusing only on the best performing scoring matrices, we saw considerable variation between decoy sets. The best values of S per decoy set varied between 0 and 0.7 (Table 3). Four of the nine decoy sets showed very little correlation (S below 0.30). Another four had some correlation (S between 0.3 and 0.6) and only two of the nine showed a reasonable correlation between score and quality (S above 0.6).

Table 3 The best results for each of the nine different decoy sets

The best scoring contact-type and contact-count scoring matrices have very similar performance over the nine different decoy sets. The nine different contact-type and contact-count S values in Table 3 have a Spearman rank correlation coefficient of 0.95. This clearly shows that the contact-type and contact-count scores have equivalent performance on the discrimination task. In all cases, a strong (or weak) correlation using the contact-type scores implies a strong (or weak) correlation using the contact-count scores.

The best values of Z per decoy set

The best values of Z for contact-type and contact-count are less strongly correlated, having a Spearman rank correlation coefficient of 0.7 (Table 3). In addition, the best Z do not correlate well with S. For one case in particular (lattice-ssfit) weak S is accompanied by a large Z (Table 3).

In general, the best contact-type scoring matrices have worse Z than the best contact-count scoring matrices (Table 3). In one case in particular (fisa-casp3), the best contact-type Z is very small (0) and the contact-count Z is moderate (1.4). However, the difference in the Z between the two different score types is not significant (p = 0.1, df = 8).

In the above two paragraphs we described the relative performance of the best contact-type and contact-count scoring matrices. The important question of which residue-residue contact definitions give the 'best' performance of S and Z is addressed in the following paragraph.

Choosing a specific residue-residue contact definition

The choice of a specific residue-residue contact definition can have a large and significant effect on the results of the scoring matrices. The performance can vary, not just between count and type scoring matrices, but also between different decoy sets. For example, the best S for the contact-count score occurs at C β -C β 8 Å without sequence separation for the 4-state-reduced decoy set, but at C β -C β 14 Å without sequence separation for the fisa decoy set. Using these contact definitions the values of S are 0.45 and 0.35 for the two decoy sets, respectively. Exchanging the contact definitions in these two cases, S falls to 0.40 and 0.25, respectively.

The ultimate aim of a scoring function is to rank near-native protein decoys according to their similarity to the native structure. The performance of the scoring function on this task should be independent of the method used to generate the decoys. For this reason, it is informative to look at the overall performance of each different scoring matrix across all the different decoy sets. The mean values of S and Z for each different scoring matrix over each decoy set are presented in Figure 4.

Figure 4

Mean performance. The mean value of S (upper) and Z (lower) for the contact-type (left) and contact-count (right) scoring matrices over the nine different decoy sets. The means were calculated as the mean of the mean value per decoy set, rather than the mean over the total set of proteins. The broken lines indicate the trend for the matrices with sequence separation filtering and the solid lines give the trend for the matrices without sequence separation filtering.

The best performance of the contact-type scoring matrices is obtained by defining residue-residue contact using C β -C β atoms with a distance threshold of 12 Å. This is obtained without sequence separation filtering, including short-range contacts. The best performance of the contact-count scoring matrices occurs at slightly longer distance threshold of 14 Å (Figure 4). Overall, the contact-count and contact-type matrices show a similar pattern of performance across different residue-residue contact definitions. The best Z are generally found using the contact-count scoring matrices. Both types of scoring matrix have a maximum in Z when using C β -C β atoms with a distance threshold of 8 Å. In addition, the Z of the contact-count scoring matrices is also high between 10 and 16 Å


Early work on single-body 'residue environment' potentials was very promising [2, 3, 5563]. However, the effectiveness of these potentials has never been directly compared to two-body 'residue pair' potentials in detail. Here we do not try to combine the two different types of residue interaction propensity score, but rather to assess their independent contribution to scoring protein structures. The objective is to examine how much information is stored in the two types of measure and to compare their performance on the realistic task of ranking a set of decoy structures.

The magnitude of the scoring matrices

To address the question of which type of residue contact propensity score contains the most specific information about protein structure, we assessed the mean absolute score (MAS) of the cells in the different scoring matrices. The score in each cell measures the strength of a certain residue contact propensity. In this sense, magnitude of the MAS gives the degree of 'non-randomness' or information content of the given residue contact propensity. The MAS suggests that, whatever the residue-residue contact definition used, the 'single-body' residue contact-count propensities were stronger or more informative than the 'two-body' residue contact-type propensities.

As the residue-residue contact definition was changed, we observed changes in MAS that were consistent with previous observations [64]. The most informative contact-type and contact-count matrices were obtained using C β -C β contacts at 6 Å without sequence separation filtering and using C β -C β at 7 Å with sequence separation filtering, respectively. However, the pattern of change in MAS that occurred as a consequence of changing residue-residue contact definition were not seen in the score performance on the task of scoring protein decoys.

Scoring native protein structures

Scoring the data set of 3, 070 native proteins highlighted some problematic structures. Some of the worst scoring proteins in this set when using either the contact-type or the contact-count scoring matrices were all found to be membrane proteins. It is not surprising that the residue contact propensities derived from a data set of mostly globular proteins are not generally the same as the propensities seen in membrane proteins.

Further down the list of the worst scoring native proteins, we find some protein subunits of oligomeric proteins that were incorrectly annotated monomers. These subunits appear 'non-native' because they would make many additional residue-residue contacts in the native oligomer. The artificially isolated subunit is effectively 'non-native'.

Ranking near-native protein decoys

Firstly, we observed that the decoys in some sets cannot be successfully ranked by either the contact-type or the contact-count scores. These sets of decoys are all considered equally 'native' (or equally 'non-native') by the residue contact propensity scores, despite having a range of different RMSD values to the native structure [54]. We observed that these decoy sets lacked decoys in the range of 1 to 5 Å RMSD, having less than 25% of the decoys below 5 Å. This observation suggests that the scores might perform better on decoys that are closer to native.

Secondly, and perhaps more importantly, we observe that the two different kinds of score perform equally well on the different decoy sets (Figure 5). The contact-type and contact-count performance in terms of both rank correlation coefficient (S) or Z-score (Z) are both highly correlated. The correlation of the best performance is 0.97 for S and 0.67 for Z.

Figure 5

Performance of the two different types of score. The best values of S and Z for the contact-type and contact-count scoring matrices. Each point represents one of the nine decoy sets studied in this work, encompassing the scores from 150 proteins and 244,794 decoys.

Finally, we observed that the specific residue-residue contact definition that gave the best performance varied between the different decoy sets studied. However, similar trends in performance were observed at any given residue-residue contact definition across all sets.

Several other groups have reported good performance on similar discrimination tasks using single-body residue burial terms. For example, in Godzik et. al. 1992 [60] it was reported that, in most cases, a burial term alone is a sufficient indicator of the native sequence compared to two- and three-body residue interaction terms. A Bayesian scoring function developed in Simons et. al. 1999 [9] suggested that residue burial scores have comparable performance to residue contact scores. Similarly, in Zhou et. al. 2004 [65] the authors concluded that the the residues solvent accessible surface area appears to be the most important among several different single-body terms tested. In addition, several groups have used a similar definition of residue contact-count as an approximation for burial [48, 50, 6567]).

The current work suggests that counting contacts between C β atoms using a distance threshold around 12 Å provides the most discriminative single body residue contact-count score (Figure 4). Similar observations have also been confirmed in the literature. For example, in Karchin et. al. 2004 [68] the best results were obtained with a 14 Å contact definition between C β atoms.

However, in a number of studies a distance threshold of 9 Å between C β atoms was used to count contacts [64, 69]. In one such study, it was stated that the 9 Å distance threshold used resulted in a slightly better performance than other cutoffs tested [69]. This difference may result from the specific count normalisation procedure applied in that work.

Only two different sequence separation filters were assessed in detail in this work, considering either all interactions or only the long-range interactions. Long range interactions were defined as interactions between residues that are more than 10 residues apart in the protein sequence. Results collected using alternative sequence separation thresholds of 5, 8 or 12 showed gave very little change in the scores collected. Using a sequence separation threshold of 2 showed scores that were roughly in between those of 0 and 10. It important to note that when scoring near-native decoys, sequence-separation filtering has very little effect on the performance of the score, as all decoys and the native protein have the same primary sequence.

Cooperativity in protein folding

It has long been suggested that pairwise potentials may not capture the inherent cooperativity of protein folding (for example see [14, 70]). Here we have presented results for the effectiveness of the contact-count score, suggesting that indeed higher order interactions are indeed important in protein structure. For example, it has been shown that contact-count can be estimated from a four-body residue-residue interaction potential [71]. However, the performance of such a four-body potential, assessed using the (SNAPP) score [72], is not significantly better than an equivalent two-body potential [73]. Despite this observation, four-body potentials are becoming much more commonly used as a way to better capture the cooperativity of protein interactions [74].

Future directions

The work presented here represents a basic comparison of contact-type and contact-count scores. There are several ways in which this basic work should be extended. However, it is important to note that the two scores compared in this work are far from optimal. It is known that distance dependent all-atom scores are more effective at discriminating between native and non-native protein structures [48, 69]. Developing the current work along these lines will be an important task for the future. In particular, it remains to be seen if the findings presented here at the residue level are consistent with observations at the atomic level.

We did not directly compare the statistical potentials derived in this work to similar potentials described by other authors in the literature. To extend the analysis presented here, our potentials should be compared directly with those in the literature (for a good example of this type of comparison see [75]) Additionally, a comparison of the important amino acid properties such as hydrophobicity and electrostatics should be performed.

In this work we did not address combinations of the two scores. The two types of potential studied perform equivalently, suggesting that they are based on a similar underlying principle. However, if a combination of scores improves the overall performance this would show that the scores carry different information. Although that an ideal scoring function should work in all possible cases, correlation between RMSD and score is usually only significant for RMSD below about 3 Å [76]. For this reason it would be useful to compare the scoring functions on decoys within specific ranges of RMSD from the native.


In this work we assessed the independent contribution of two different types of residue contact propensity to scoring protein structures. The main finding is that the contact-type and contact-count scores showed equivalent overall performance in the task of ranking protein decoys. Although the two different score types perform equivalently, the ability to automatically predict the number of contacts made by a residue [68, 77, 78] allows for a greater range of applications. In addition, a single-body term is amenable to an efficient dynamic programming method for alignment optimisation [3, 65, 79].

The work presented here represents our preliminary investigation of a multi-body potential for evaluation of protein structure. In future it should be possible to combine the contact-type and contact-count scores to better take into account the inherent cooperativity of protein folding.


The data set of native proteins and protein decoys

The non-redundant data set of monomers

The scoring matrices were derived from a non-redundant set of high-quality protein monomers. This set of 3, 070 proteins was selected using the following protocol. Only the monomers from the BioUnit section of the Protein Data Bank (PDB) [80, 81] were selected, excluding putative structures of dimers, trimers, and the other multi-subunit proteins. The resulting monomeric proteins were further filtered by size, having more than 20 amino acids, and by resolution, being better than 3 Å. Finally, the chains were made non-redundant at 30% sequence identity using BLASTClust [82]. The resulting set of 3, 070 monomers was used throughout this analysis.

The data set of 'near-native' protein decoys

Nine sets of protein decoys were taken from the Decoys R us database [54]. In total this data set included 244, 794 decoys derived from 150 native proteins.

Constructing the scoring matrices


Given the fraction of residues of type x and of type y (P x and P y ), the probability of randomly observing a contact of type xy is,P xy = P x ·P y

where P xy is the probability of a 'random' contact of type xy (for example see Table 1). This formula is obtained by assuming that contacts are made between randomly selected pairs of residues, assuming statistical independence. In this way, we make no assumptions about the distribution of contacts within the protein, such as the distribution of the number of contacts per residue.

The observed and expected probabilities of a contact of type xy can be combined into a score using the log-odds ratio;

S x y = log ( P x y o b s P x y e x p ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaem4uam1aaSbaaSqaaiabdIha4jabdMha5bqabaGccqGH9aqpcyGGSbaBcqGGVbWBcqGGNbWzdaqadaqcfayaamaalaaabaGaemiuaa1aa0baaeaacqWG4baEcqWG5bqEaeaacqWGVbWBcqWGIbGycqWGZbWCaaaabaGaemiuaa1aa0baaeaacqWG4baEcqWG5bqEaeaacqWGLbqzcqWG4baEcqWGWbaCaaaaaaGccaGLOaGaayzkaaaaaa@48A7@

The magnitude of the score S xy gives a measure of how 'non-randomly' the pair xy occurs. The score is positive when xy is observed more often than expected and negative when xy is observed less often than expected.


The contact-count scoring matrix is created in a similar way to the contact-type scoring matrix. However, instead of using the frequency P xy to denote the probability of a residue-residue contact between residue type x and type y, we use P xn to denote the probability of a residue of type x having exactly n residue-residue contacts. The probability P xn is defined as,P xn = N xn /N n

where N xn is the observed number of residues of type x having exactly n contacts and N n is the total number of residues with exactly n contacts (for example see Figure 1). The values of n were taken from those observed over all residues. The 'random' value of P xn is simply taken to be equal to P x , the fraction of residues of type x. Using this value assumes that there is no particular effect on the overall amino acid composition when filtered by a given number of contacts. Again the observed and expected probabilities can be combined as in Equation 2. Using this method, the magnitude of the score S xn can be easily interpreted as a measure of 'compositional bias' given a certain number of residue-residue contacts.


Certain residue-residue contact definitions could lead to low counts in the scoring matrices. For example, if the overall number of observed residue-residue contacts is low, certain contact-types may become rare. Similarly, if the number of different contact-counts spans a wide range, instances of a given residue type with a given contact-count may become rare. To address this issue of undersampling, if a cell of a scoring matrix was based on fewer than 5 observed or expected counts, that score was discarded. Those specific classes of contact were therefore ignored when scoring protein structures, being neither penalised nor rewarded.

Although a threshold of 5 observed or expected counts was used to filter undersampled classes in the results presented here, it should be noted that both the contact-type and contact-count scores appear very robust, showing only small changes in MAS when discarding cells with fewer than 50 observed or expected counts. The complete set of counts and scores for the contact-type and contact-count scoring matrices are included as additional files (see Additional file 1).


  1. 1.

    Narayana S, Argos P: Residue contacts in protein structures and implications for protein folding. Int J Pept Protein Res 1984, 24: 25–39.

    CAS  Article  Google Scholar 

  2. 2.

    Novotn'y J, Bruccoleri R, Karplus M: An analysis of incorrectly folded protein models. Implications for structure predictions. J Mol Biol 1984, 177: 787–818.

    Article  Google Scholar 

  3. 3.

    Bowie J, Lüthy R, Eisenberg D: A method to identify protein sequences that fold into a known three-dimensional structure. Science 1991, 253: 164–170.

    CAS  Article  Google Scholar 

  4. 4.

    Blundell T, Sibanda B, Sternberg M, Thornton J: Knowledge-based prediction of protein structures and the design of novel molecules. Nature 1987, 326: 347–352.

    CAS  Article  Google Scholar 

  5. 5.

    Sali A, Blundell T: Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol 1993, 234: 779–815.

    CAS  Article  Google Scholar 

  6. 6.

    Vajda S, Sippl M, Novotny J: Empirical potentials and functions for protein folding and binding. Curr Opin Struct Biol 1997, 7: 222–228.

    CAS  Article  Google Scholar 

  7. 7.

    Levitt M: A simplified representation of protein conformations for rapid simulation of protein folding. J Mol Biol 1976, 104: 59–107.

    CAS  Article  Google Scholar 

  8. 8.

    Tanaka S, Scheraga HA: Medium- and long-range interaction parameters between amino acids for predicting three-dimensional structures of proteins. Macromolecules 1976, 9(6):945–950.

    CAS  Article  Google Scholar 

  9. 9.

    Simons KT, Ruczinski I, Kooperberg C, Fox BA, Bystroff C, Baker D: Improved recognition of native-like protein structures using a combination of sequence-dependent and sequence-independent features of proteins. Proteins 1999, 34: 82–95.

    CAS  Article  Google Scholar 

  10. 10.

    Ponder J, Richards F: Tertiary templates for proteins. Use of packing criteria in the enumeration of allowed sequences for different structural classes. J Mol Biol 1987, 193: 775–791.

    CAS  Article  Google Scholar 

  11. 11.

    Wilson C, Doniach S: A computer model to dynamically simulate protein folding: studies with crambin. Proteins 1989, 6: 193–209.

    CAS  Article  Google Scholar 

  12. 12.

    Qin M, Wang J, Tang Y, Wang W: Folding behaviors of lattice model proteins with three kinds of contact potentials. Phys Rev E Stat Nonlin Soft Matter Phys 2003, 67: 061905.

    Article  Google Scholar 

  13. 13.

    Anfinsen C: Principles that govern the folding of protein chains. Science 1973, 181: 223–230.

    CAS  Article  Google Scholar 

  14. 14.

    Mayewski S: A multibody, whole-residue potential for protein structures, with testing by Monte Carlo simulated annealing. Proteins 2005, 59(2):152–169.

    CAS  Article  Google Scholar 

  15. 15.

    Zimmerman S, Pottle M, Némethy G, Scheraga H: Conformational analysis of the 20 naturally occurring amino acid residues using ECEPP. Macromolecules 1977, 10: 1–9.

    CAS  Article  Google Scholar 

  16. 16.

    Nemethy G, Gibson KD, Palmer KA, Yoon CN, Paterlini G, Zagari A, Rumsey S, Scheraga HA: Energy parameters in polypeptides. 10. Improved geometrical parameters and nonbonded interactions for use in the ECEPP/3 algorithm, with application to proline-containing peptides. Journal of Physical Chemistry 1992, 96(15):6472–6484.

    CAS  Article  Google Scholar 

  17. 17.

    Allinger NL: Conformational analysis. 130. MM2. A hydrocarbon force field utilizing V1 and V2 torsional terms. Journal of the American Chemical Society 1977, 99(25):8127–8134.

    CAS  Article  Google Scholar 

  18. 18.

    Allinger NL, Yuh YH, Lii JH: Molecular mechanics. The MM3 force field for hydrocarbons. 1. Journal of the American Chemical Society 1989, 111(23):8551–8566.

    CAS  Article  Google Scholar 

  19. 19.

    Weiner S, Kollman P, Case D, Singh U, Ghio C, Alagona G, Profeta S, Weiner P: A new force-field for molecular mechanical simulation of nucleic-acids and proteins. Journal of the American Chemical Society 1984, 106: 765–784.

    CAS  Article  Google Scholar 

  20. 20.

    Cornell WD, Cieplak P, Bayly CI, Gould IR, Merz KM, Ferguson DM, Spellmeyer DC, Fox T, Caldwell JW, Kollman PA: A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules. Journal of the American Chemical Society 1995, 117(19):5179–5197.

    CAS  Article  Google Scholar 

  21. 21.

    Brooks B, Bruccoleri R, Olafson B, States D, Swaminathan S, M K: A program for macromolecular energy, minimization, and dynamics calculations. J Comput Chem 1983, 4: 187–217.

    CAS  Article  Google Scholar 

  22. 22.

    MacKerell A, Bashford D, Bellott M, Dunbrack R, Evanseck J, Field M, Fischer S, Gao J, Guo H, Ha S, Joseph-McCarthy D, Kuchnir L, Kuczera K, Lau F, Mattos C, Michnick S, Ngo T, Nguyen D, Prodhom B, Reiher W, Roux B, Schlenkrich M, Smith J, Stote R, Straub J, Watanabe M, Wiorkiewicz-Kuczera J, Yin D, Karplus M: All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of Proteins. Journal of Physical Chemistry B 1998, 102(18):3586–3616.

    CAS  Article  Google Scholar 

  23. 23.

    Smith JC, Karplus M: Empirical force field study of geometries and conformational transitions of some organic molecules. Journal of the American Chemical Society 1992, 114(3):801–812.

    CAS  Article  Google Scholar 

  24. 24.

    Lindahl E, Hess B, Spoel D: GROMACS 3.0: a package for molecular simulation and trajectory analysis. Journal of Molecular Modeling 2001, 7(8):306–317.

    CAS  Google Scholar 

  25. 25.

    Miyazawa S, Jernigan RL: Estimation of effective interresidue contact energies from protein crystal structures: quasi-chemical approximation. Macromolecules 1985, 18(3):534–552.

    CAS  Article  Google Scholar 

  26. 26.

    Sippl MJ: Calculation of conformational ensembles from potentials of mean force. An approach to the knowledge-based prediction of local structures in globular proteins. J Mol Biol 1990, 213(4):859–883.

    CAS  Article  Google Scholar 

  27. 27.

    Godzik A, Kolinski A, Skolnick J: Are proteins ideal mixtures of amino acids? Analysis of energy parameter sets. Protein Sci 1995, 4(10):2107–2117.

    CAS  Article  Google Scholar 

  28. 28.

    Thomas P, Dill K: Statistical potentials extracted from protein structures: how accurate are they? J Mol Biol 1996, 257: 457–469.

    CAS  Article  Google Scholar 

  29. 29.

    Richardson J: The anatomy and taxonomy of protein structure. Adv Protein Chem 1981, 34: 167–339.

    CAS  Article  Google Scholar 

  30. 30.

    Efimov A: A novel super-secondary structure of proteins and the relation between the structure and the amino acid sequence. FEBS Lett 1984, 166: 33–38.

    CAS  Article  Google Scholar 

  31. 31.

    Finkelstein A, Ptitsyn O: Why do globular proteins fit the limited set of folding patterns? Prog Biophys Mol Biol 1987, 50: 171–190.

    CAS  Article  Google Scholar 

  32. 32.

    Chothia C, Finkelstein A: The classification and origins of protein folding patterns. Annu Rev Biochem 1990, 59: 1007–1039.

    CAS  Article  Google Scholar 

  33. 33.

    Harris N, Presnell S, Cohen F: Four helix bundle diversity in globular proteins. J Mol Biol 1994, 236: 1356–1368.

    CAS  Article  Google Scholar 

  34. 34.

    Walther D, Eisenhaber F, Argos P: Principles of helix-helix packing in proteins: the helical lattice superposition model. J Mol Biol 1996, 255: 536–553.

    CAS  Article  Google Scholar 

  35. 35.

    Hill T: Statistical mechanics: principles and selected applications. McGraw-Hill series in advanced chemistry, New York: McGraw-Hill; 1956.

    Google Scholar 

  36. 36.

    Xia Y, Levitt M: Extracting knowledge-based energy functions from protein structures by error rate minimization: Comparison of methods using lattice model. J Chem Phys 2000, 113(20):9318–9330.

    CAS  Article  Google Scholar 

  37. 37.

    Hendlich M, Lackner P, Weitckus S, Floeckner H, Froschauer R, Gottsbacher K, Casari G, Sippl MJ: Identification of native protein folds amongst a large number of incorrect models. The calculation of low energy conformations from potentials of mean force. J Mol Biol 1990, 216: 167–180.

    CAS  Article  Google Scholar 

  38. 38.

    Miyazawa S, Jernigan RL: Residue-residue potentials with a favorable contact pair term and an unfavorable high packing density term, for simulation and threading. J Mol Biol 1996, 256(3):623–644.

    CAS  Article  Google Scholar 

  39. 39.

    Betancourt MR, Thirumalai D: Pair potentials for protein folding: choice of reference states and sensitivity of predicted native states to variations in the interaction schemes. Protein Sci 1999, 8(2):361–369. [Comparative Study]

    CAS  Article  Google Scholar 

  40. 40.

    Simons KT, Kooperberg C, Huang E, Baker D: Assembly of protein tertiary structures from fragments with similar local sequences using simulated annealing and Bayesian scoring functions. J Mol Biol 1997, 268: 209–225.

    CAS  Article  Google Scholar 

  41. 41.

    Samudrala R, Moult J: An all-atom distance-dependent conditional probability discriminatory function for protein structure prediction. J Mol Biol 1998, 275(5):895–916.

    CAS  Article  Google Scholar 

  42. 42.

    Maiorov VN, Crippen GM: Contact potential that recognizes the correct folding of globular proteins. J Mol Biol 1992, 227(3):876–888.

    CAS  Article  Google Scholar 

  43. 43.

    Goldstein RA, Luthey-Schulten ZA, Wolynes PG: Optimal protein-folding codes from spin-glass theory. Proc Natl Acad Sci USA 1992, 89(11):4918–4922.

    CAS  Article  Google Scholar 

  44. 44.

    Hao MH, Scheraga H: Optimizing Potential Functions for Protein Folding. J Phys Chem 1996, 100(34):14540–14548.

    CAS  Article  Google Scholar 

  45. 45.

    Mirny LA, Shakhnovich EI: How to derive a protein folding potential? A new approach to an old problem. J Mol Biol 1996, 264(5):1164–1179. [Comparative Study]

    CAS  Article  Google Scholar 

  46. 46.

    Chiu TL, Goldstein RA: Optimizing energy potentials for success in protein tertiary structure prediction. Fold des 1998, 3(3):223–228.

    CAS  Article  Google Scholar 

  47. 47.

    Moult J: Comparison of database potentials and molecular mechanics force fields. Curr Opin Struct Biol 1997, 7(2):194–199.

    CAS  Article  Google Scholar 

  48. 48.

    Zhou H, Zhou Y: Distance-scaled, finite ideal-gas reference state improves structure-derived potentials of mean force for structure selection and stability prediction. Protein Sci 2002, 11: 2714–2726.

    CAS  Article  Google Scholar 

  49. 49.

    Jones D, Taylor W, Thornton J: A new approach to protein fold recognition. Nature 1992, 358: 86–89.

    CAS  Article  Google Scholar 

  50. 50.

    Sippl M: Boltzmann's principle, knowledge-based mean fields and protein folding. An approach to the computational determination of protein structures. J Comput Aided Mol Des 1993, 7: 473–501.

    CAS  Article  Google Scholar 

  51. 51.

    Tanaka S, Scheraga H: Model of protein folding: inclusion of short-, medium-, and long-range interactions. Proc Natl Acad Sci USA 1975, 72: 3802–3806.

    CAS  Article  Google Scholar 

  52. 52.

    Greene L, Higman V: Uncovering network systems within protein structures. J Mol Biol 2003, 334: 781–791.

    CAS  Article  Google Scholar 

  53. 53.

    Gilis D: Protein decoy sets for evaluating energy functions. J Biomol Struct Dyn 2004, 21: 725–736.

    CAS  Article  Google Scholar 

  54. 54.

    Samudrala R, Levitt M: Decoys 'R' Us: a database of incorrect conformations to improve protein structure prediction. Protein Sci 2000, 9: 1399–1401.

    CAS  Article  Google Scholar 

  55. 55.

    Eisenberg D, McLachlan A: Solvation energy in protein folding and binding. Nature 1986, 319: 199–203.

    CAS  Article  Google Scholar 

  56. 56.

    Novotný J, Rashin A, Bruccoleri R: Criteria that discriminate between native proteins and incorrectly folded models. Proteins 1988, 4: 19–30.

    Article  Google Scholar 

  57. 57.

    Baumann G, Frömmel C, Sander C: Polarity as a criterion in protein design. Protein Eng 1989, 2: 329–334.

    CAS  Article  Google Scholar 

  58. 58.

    Chiche L, Gregoret L, Cohen F, Kollman P: Protein model structure evaluation using the solvation free energy of folding. Proc Natl Acad Sci USA 1990, 87: 3240–3243.

    CAS  Article  Google Scholar 

  59. 59.

    Casari G, Sippl M: Structure-derived hydrophobic potential. Hydrophobic potential derived from X-ray structures of globular proteins is able to identify native folds. J Mol Biol 1992, 224: 725–732.

    CAS  Article  Google Scholar 

  60. 60.

    Godzik A, Skolnick J: Sequence-structure matching in globular proteins: application to supersecondary and tertiary structure determination. Proc Natl Acad Sci USA 1992, 89: 12098–12102.

    CAS  Article  Google Scholar 

  61. 61.

    Lüthy R, Bowie J, Eisenberg D: Assessment of protein models with three-dimensional profiles. Nature 1992, 356: 83–85.

    Article  Google Scholar 

  62. 62.

    Ouzounis C, Sander C, Scharf M, Schneider R: Prediction of protein structure by evaluation of sequence-structure fitness. Aligning sequences to contact profiles derived from three-dimensional structures. J Mol Biol 1993, 232: 805–825.

    CAS  Article  Google Scholar 

  63. 63.

    Huang E, Subbiah S, Levitt M: Recognizing native folds by the arrangement of hydrophobic and polar residues. J Mol Biol 1995, 252: 709–720.

    CAS  Article  Google Scholar 

  64. 64.

    Melo F, Sánchez R, Sali A: Statistical potentials for fold assessment. Protein Sci 2002, 11: 430–448.

    CAS  Article  Google Scholar 

  65. 65.

    Zhou H, Zhou Y: Single-body residue-level knowledge-based energy score combined with sequence-profile and secondary structure information for fold recognition. Proteins 2004, 55: 1005–1013.

    CAS  Article  Google Scholar 

  66. 66.

    Melo F, Feytmans E: Assessing protein structures with a non-local atomic interaction energy. J Mol Biol 1998, 277: 1141–1152.

    CAS  Article  Google Scholar 

  67. 67.

    Jones D: GenTHREADER: an efficient and reliable protein fold recognition method for genomic sequences. J Mol Biol 1999, 287: 797–815.

    CAS  Article  Google Scholar 

  68. 68.

    Karchin R, Cline M, Karplus K: Evaluation of local structure alphabets based on residue burial. Proteins 2004, 55(3):508–518.

    CAS  Article  Google Scholar 

  69. 69.

    Benkert P, Tosatto S, Schomburg D: QMEAN: A comprehensive scoring function for model quality assessment. Proteins 2008, 71(1):261–277.

    CAS  Article  Google Scholar 

  70. 70.

    Vendruscolo M, Domany E: Pairwise contact potentials are unsuitable for protein folding. J Chem Phys 1998, 109(24):11101–11108.

    CAS  Article  Google Scholar 

  71. 71.

    Munson P, Singh R: Statistical significance of hierarchical multi-body potentials based on Delaunay tessellation and their application in sequence-structure alignment. Protein Sci 1997, 6: 1467–1481.

    CAS  Article  Google Scholar 

  72. 72.

    Tropsha A, Carter C, Cammer S, Vaisman I: Simplicial neighborhood analysis of protein packing (SNAPP): a computational geometry approach to studying proteins. Meth Enzymol 2003, 374: 509–544.

    CAS  Article  Google Scholar 

  73. 73.

    Gan H, Tropsha A, Schlick T: Lattice protein folding with two and four-body statistical potentials. Proteins 2001, 43: 161–174.

    CAS  Article  Google Scholar 

  74. 74.

    Feng Y, Kloczkowski A, Jernigan R: Four-body contact potentials derived from two protein datasets to discriminate native structures from decoys. Proteins 2007, 68: 57–66.

    CAS  Article  Google Scholar 

  75. 75.

    Pokarowski P, Kloczkowski A, Jernigan R, Kothari N, Pokarowska M, Kolinski A: Inferring ideal amino acid interaction forms from statistical protein contact potentials. Proteins 2005, 59: 49–57.

    CAS  Article  Google Scholar 

  76. 76.

    Liu T, Samudrala R: The effect of experimental resolution on the performance of knowledge-based discriminatory functions for protein structure selection. Protein Eng Des Sel 2006, 19: 431–437.

    CAS  Article  Google Scholar 

  77. 77.

    Song J, Burrage K: Predicting residue-wise contact orders in proteins by support vector regression. BMC Bioinformatics 2006, 7: 425.

    Article  Google Scholar 

  78. 78.

    Ishida T, Nakamura S, Shimizu K: Potential for assessing quality of protein structure based on contact number prediction. Proteins 2006, 64: 940–947.

    CAS  Article  Google Scholar 

  79. 79.

    Gracy J, Chiche L, Sallantin J: Improved alignment of weakly homologous protein sequences using structural information. Protein Eng 1993, 6: 821–829.

    CAS  Article  Google Scholar 

  80. 80.

    Bernstein FC, Koetzle TF, Williams GJ, Meyer EFJ, Brice MD, Rodgers JR, Kennard O, Shimanouchi T, Tasumi M: The Protein Data Bank: a computer-based archival file for macromolecular structures. J Mol Biol 1977, 112(3):535–542.

    CAS  Article  Google Scholar 

  81. 81.

    Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucleic Acids Res 2000, 28: 235–242.

    CAS  Article  Google Scholar 

  82. 82.

    Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol 1990, 215(3):403–410.

    CAS  Article  Google Scholar 

Download references


DB wishes to thank: Ram Samudrala for making the 'Decoys-R-us' database available for public download, Linus Torvalds for the Linux OS, and all scientists who pursue their work with honesty and integrity.

Author information



Corresponding author

Correspondence to Dan M Bolser.

Additional information

Authors' contributions

DB conceived and designed the study, performed the statistical analysis and drafted the manuscript. IF helped design the study, helped write software for acquisition of data and helped draft the manuscript. HS helped write software for acquisition of data and helped draft the manuscript. JD wrote the software for acquisition of data. ML participated in design and coordination of the study and helped to draft the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

The complete set of counts and scores for the contact-type and contact-count scoring matrices

Additional file 1: . The zipped tar archive contains four directories containing all the residue contact data that was used in the analysis, along with the resulting scores. All data within the archive is provided as tab-delimited flat-files in text format. (GZ 6 MB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Bolser, D.M., Filippis, I., Stehr, H. et al. Residue contact-count potentials are as effective as residue-residue contact-type potentials for ranking protein decoys. BMC Struct Biol 8, 53 (2008).

Download citation


  • Distance Threshold
  • Potential Energy Function
  • Residue Interaction
  • Sequence Separation
  • Contact Definition