Skip to main content

Discriminating the native structure from decoys using scoring functions based on the residue packing in globular proteins



Setting the rules for the identification of a stable conformation of a protein is of utmost importance for the efficient generation of structures in computer simulation. For structure prediction, a considerable number of possible models are generated from which the best model has to be selected.


Two scoring functions, Rs and Rp, based on the consideration of packing of residues, which indicate if the conformation of an amino acid sequence is native-like, are presented. These are defined using the solvent accessible surface area (ASA) and the partner number (PN) (other residues that are within 4.5 Å) of a particular residue. The two functions evaluate the deviation from the average packing properties (ASA or PN) of all residues in a polypeptide chain corresponding to a model of its three-dimensional structure. While simple in concept and computationally less intensive, both the functions are at least as efficient as any other energy functions in discriminating the native structure from decoys in a large number of standard decoy sets, as well as on models submitted for the targets of CASP7. Rs appears to be slightly more effective than Rp, as determined by the number of times the native structure possesses the minimum value for the function and its separation from the average value for the decoys.


Two parameters, Rs and Rp, are discussed that can very efficiently recognize the native fold for a sequence from an ensemble of decoy structures. Unlike many other algorithms that rely on the use of composite scoring function, these are based on a single parameter, viz., the accessible surface area (or the number of residues in contact), but still able to capture the essential attribute of the native fold.


Predicting the native structure of proteins from their amino acid sequences has yet remained an elusive goal. In general this entails the development of effective methods for conformation sampling and the design of an accurate function for structure discrimination [1, 2]. The functions could be based on elaborate calculations and analyses of forces between atoms [3, 4], or be knowledge-based that extract relevant parameters from a database of experimentally determined protein structures [5, 6]. One important area of application of knowledge-based potential functions has been in "protein threading" for the prediction of protein tertiary structure in the absence of detectable sequence homology. The technique involves threading a protein sequence onto the frameworks of known protein folds and finding the most energetically favorable conformation [710]. In addition to fold recognition applications, where the best conformation of a protein is selected from a database of known protein conformations, the knowledge-based scoring functions are also used in protein folding simulations [6, 1116]. Many statistical scoring functions assume that frequencies of non-bonded pairs of amino acids follow a Boltzmann-like distribution and the minimum value of the score occurs in the vicinity of the lowest energy structure. Additionally, a set of probability distributions can also be used to construct a scoring function such that it can identify the maximum probability structure.

For testing of empirical energy functions challenging and diverse datasets of decoy structures that are native-like in properties have been generated [12, 1719]. Models submitted in the community-wide experiment, CASP (Critical Assessment of techniques for protein Structure Prediction) [20] make up diverse sets of structures resulting from various computational approaches [21]. The most native-like structure needs to be identified from among these models [22]. An effective potential should be able to distinguish the native structure from decoy structures with a high degree of accuracy. Energy functions based on residue contact or compactness alone do not have enough discriminating power [12], or can rank the native structure highly only when the competing conformations are more random-coil like [23]. However, here we present two knowledge-based scoring functions based on the analysis of residue packing in protein structures that are quite robust in discriminating the native conformation from a number of misfolded conformations for a given primary protein sequence. The functions were also tested on ~ 19000 models from server predictions for 71 targets of CASP7 [20]. As a descriptor for the residue packing we use the average values of the accessible surface area or the number of other residues in contact around a given residue, calculated from a database of globular proteins. Each of the function then evaluates the cumulative value for the deviation of the parameter for individual residues from the corresponding average value over the whole polypeptide chain. The experimental structure is found to have the minimum deviation and thus the minimum value of the function, when applied to a set of decoys from which the native structure has to be identified. The success of the function indicates that the burial of each residue and its contact to the surrounding residues is optimized during folding and the average values of these parameters can be used as constraint to simulate folding process. Additionally, a surface patch with residues having a large overall deviation of these parameters from the average values may be indicative of the binding region on a protein structure, an issue that would be addressed in future to provide a common perception to both the folding and the binding processes.


Scoring functions have been used to validate X-ray crystal structures, assess and rank three-dimensional models generated for a protein sequence, predict the effect of mutations, etc. Here, we are concerned with the identification of the native structure from decoys. The idea of the use of the discriminatory function originated from the formula of R-factor in crystallography [24]. An exact equivalent formula would have meant the use of the expression (1) instead of (3), given in Methods.


The individual term in Eq. (3) involves the absolute difference between the observed and the average values of ASA for a given residue, normalized by the average value. These terms are summed over the whole sequence. In Eq. (1) the numerator and the denominator are summed separately. Some other modified formulae, including the use of the standard deviation on the average values <ASAx> in the denominator, were also tried, but (2) and (3) were found most efficient to identify the native structure from a set of decoys. Depending on the structural context larger residues may have a considerable variation in their ASA values in protein structures (as indicated by larger standard deviations, Table 1) - normalization of the difference in the numerator in Eq. (3) has the effect of damping the contribution of such residues in the summation.

Table 1 Average values of partner number (<PN>) and accessible surface area (<ASA>) of different amino acid residues

Quantification of the overall packing of residues in protein structures

The average number of partner residues and the average accessible surface area for all twenty amino acids are provided in Table 1[25]. While the <ASA> values are almost identical to those calculated earlier [26], the values for the partner number are different, as the calculation is residue-based here, while in the earlier study the individual atoms constituted the partners.

As Rp and Rs indicate the extent of deviation of PN and ASA of residues from their average values, taken over the whole structure, these parameters can be used to judge the optimization of packing of residues in a structure [27]. We also wanted to see if there is any variation depending on the class of protein. However, as Rp and Rs provide cumulative values over all the residues in a structure, it is sensible to divide them by the number of residues in a structure before comparison. Individual protein structures in the dataset were classified according to CATH (Class, Architecture, Topology, Homologous superfamily; into 157 all-α, 142 all-β and 133 αβ (including α+β and α/β) classes of proteins. The normalized values (Table 2) are rather similar, except slightly higher values in the all-β class, indicating somewhat higher deviations from the optimum values of PN and ASA in these structures. The observation of higher values in β-proteins is in tune with a relatively lesser packing efficiency in these proteins, as is also demonstrated by the higher occurrence of cavities involving residues in β-sheets [28].

Table 2 Average values of Rs and Rp in various protein structural classesa

Identification of the native structure from misfolded decoys

PROSTAR decoy sets

The objective of this work is to discriminate between the native structure and one or more misfolded or low-resolution structures. The utility of Rp and Rs was tested on the decoy sets in the PROSTAR website and the results are shown in Table 3. When compared with other atomic or residue-based potentials, the present parameters, Rs and Rp have similar or better performance, except for 'Ifu'. Of the two parameters, Rs based on residue accessibility performs better than the one derived on the basis of partner number (Rp).

Table 3 Identification of the native structure from decoys in PROSTAR decoy sets using different scoring functionsa

The 'Misfold' decoy set, generated by Holm and Sander [17], consists of 24 examples of pairs of proteins with the same number of residues in the chain, but different sequences and conformations. Sequences are swapped between members of a pair, resulting in rather inappropriate environments for most of the side chains. For this set, Rs selects 100% of the structures correctly, but Rp fails in four. Attempts were made to see if the use of other cut-off distances (4.0, 5.0, 6.0 and 7.0 Å) in the definition of Rp improved the situation, but the performance of the parameter derived at 4.5 Å was found to be the best.

The 'Ifu' decoy set is based on a set of 43 peptides, 10-20 residues long, which are proposed to be independent folding units as determined by local hydrophobic burial and experimental evidence [29]. In this test set, Rs and Rp were unsuccessful to pick 21 and 22, respectively, out of 43 native structures. While performing the best, even the knowledge-based potential [14] failed in 11 cases in this test set. This is probably because the targets in these subsets are protein pieces and it is difficult for residue packing parameters derived from larger proteins to evaluate these structures.

The 'Asilomar' decoy set resulted from the first experiment on the Critical Assessment of Protein Structure Prediction methods (CASP), which produced a set of 41 comparative models of six different proteins [30]. The models vary in Cα rmsd to the corresponding experimental conformation, ranging from 0.53 to 7.40 Å, depending on the difficulty of the model building process. In this test set, the parameter Rs selects 100% native structures correctly, by far the best result from any discriminatory function. For Rp, missing 5 out of 41 cases, the performance is at par with other functions.

The 'Pdberr' decoy set consists of structures determined using X-ray crystallography that were later found to contain errors, and the corresponding corrected experimental conformations [31]. The 'sgpa' decoy set consists of the experimental structure Streptomyces griseus Protease A (2sga) and two conformations generated by molecular dynamics simulations starting with the experimental structure [32]. In these test sets, where the decoys are low-resolution X-ray structures, both the scoring functions Rs and Rp correctly picked the high-resolution structures in all cases, as did all other potential functions, except the one based on the residue contact potential with a composition-corrected scale [33].

Park and Levitt decoy set

The Park and Levitt decoy test set, available on the web site, consists of 7 sequences, each with nearly 600-700 decoys that cover structures showing an rmsd ranging from 0 (the correct fold) to 10 Å from the native structure [12]. The protein structures were generated by using four-state models (four discrete ϕ,ψ angles) to define the conformation of each of ten selected residues in each protein using an off-lattice model. From the very large number of conformations generated, only those compact structures were retained that scored well using a variety of scoring functions, as well as having a reasonable rmsd from the native structure. The 4-state-reduced decoy data set given in Additional file 1: Table S1 includes a range of small proteins from 54-75 residues with varying topological folds, with the numbers of decoys ranging from 630 for 1ctf to 687 for 4pti. A positive Z-score (Equations (4) and (5)) indicates that the value of the parameter for a particular native fold is lower than the average of the distribution. While considering the Zs, the native structure is well separated from the average of the distribution for all the structures, but Zp shows an inferior result for 1r69 and 1sn3. Figure 1 plots Rs vs rmsd for a representative dataset corresponding to the PDB file, 1ctf. The value of Rs is the minimum for the native structure. There is a good linear correlation between the two variables (R2 is 0.78), better than that (0.6) obtained using the knowledge-based potential of Lu and Skolnick [14]. While the various energy functions based on empirical contact, surface area and van der Waals energy did not perform consistently well to distinguish between correct and incorrect conformations and had to be used in combination for the proper identification of the correct fold [12], the rather simple parameter, Rs has a remarkable discriminatory power.

Figure 1

Scatter plot of R s vs. rmsd for a representative protein structure, 1ctf, along with its decoys.

The Levitt low-minima decoy sets (LMDS) also contain structural decoys (the number ranging from 343 to 500) for 7 small proteins, 36 to 68 residues long [19]. From an initial ten thousand structures, generated by randomly modifying only the loop dihedral angles, which were subjected to minimization using a modified ENCAD force field involving united and soft atoms [34], up to five hundred of the lowest energy conformations were retained to make up the decoy sets. For all the 7 cases the native structure has the minimum Rs value and the corresponding Z-score indicates that it is well separated from the decoys (Additional file 1: Table S1). However, Zp gives an inferior result for 1bba and 1fc2. Other energy functions also failed to identify the native structure for these two proteins [15, 22] due to the fact that the native conformation is simply not very well defined for the former [35] and the latter is a fragment of a larger protein and additionally, a constituent of a complex, and in the unbound form may have a structure different from that in the complex [36]. Interestingly however, based on Rs both the native structures are separated by about two standard deviations from the average of the distribution.

ROSETTA decoy sets

The ROSETTA all-atom decoy sets are composed of five different proteins ranging in size from 92 to 116 residues, and the number of decoys ranging from 994 to 999 (Additional file 1: Table S1) [37]. Fragments, between 3 and 9 residues, from known structures matched to the targets through a multiple sequence alignment process, were assembled into the protein structures via the fragment insertion-simulated annealing strategy [37]. The scoring functions used to select the lowest energy decoys included hydrophobic burial, electrostatics, the formation of β-sheets and the packing of α-helices and β-strands. The Z-scores based on Rs and Rp indicate that both the scoring functions perform well over all the 5 structures. The large Z-scores seen here, as compared to those in others, should be due to the high rmsds in the decoys used in this test set.

The original ROSETTA decoy set has been improved by increasing the number of proteins and frequency of near native models, providing 1,400 model structures for 78 diverse, single domain proteins with varying degrees of secondary structure and length from 25 to 87 residues for the evaluation of scoring functions [16]. The discriminatory ability of our scoring functions can be seen from the results on 41 cases (a subset of the complete dataset, which is downloadable) presented in Additional file 1: Table S2. The native structure did not have the minimum Rs value in 3 cases, while Rp failed in two additional cases. For these, the Z-score is also quite small, Zp even registering a negative value in two. It may be noted that two structures (1res and 1uxd) among the failed cases were derived from NMR experiments and the Rosetta energy functions are also less efficient in identifying the NMR structures as compared to X-ray crystal structures, probably because the former structures have greater deviation of side chain conformations from the canonical rotamer conformations [16].

Identification of the native structure from the native-like conformation constructed by homology modeling

Samudrala and Levitt [19] have a decoy set (hg_structal) for 29 globins. Each globin has been built by comparative modeling using 29 other globins as templates with the program segmod [38]; the rmsd for the modeled structures range from 1.96 to 8.57 Å. A similar decoy set (ig_structal_hires) involving 20 immunoglobulins and at a relatively higher resolution (1.7-2.2 Å, compared to the range of 1.7-3.1 Å for the full set of 61 proteins) is also available. The application of our scoring function on these two sets yields results given in Table 4. As with the other decoy sets, Rs performs better than Rp in identifying the native structure. Even though the homology built models in the 'ig_structal_hires' set are very close to the native structure, the latter was identifiable in 90% of the cases.

Table 4 Identification of native structure from decoys constructed by homology modeling

Score of the experimental structure relative to the solutions submitted to CASP7

The ability of our scoring function to identify the native structure from the best near-native solutions has been tested on the CASP7 dataset [20]. This is the most difficult test as the decoys are the best predicted near-native structures submitted by different groups participating in the CASP experiment. CASP7 experiment consists of 95 accepted targets for which about 22000 models were submitted. We have excluded the NMR structures and retained 71 targets (with ~ 19000 models) to evaluate our scoring functions (Additional file 1: Table S3). The rmsd between the native structure and the best predicted solution varies in the range 0.4 - 2.6 Å in the whole dataset. Zs identifies the best solution in 51 cases and Zp in 38. Table 5 compares the results of our study vis-à-vis those from other algorithms [22]. As we have seen before, Rs performs better than Rp. But even Rp outperforms other existing functions in locating the native structure among the top ten solutions. Rs identifies the native structure as the top solution in 72% of cases, which is considerably better than the next best performer (DFIRE and QMEAN3) at 62%.

Table 5 Performance of the different scoring function for predicting the native structure among the best near-native structures submitted in CASP7


There are many energy functions (knowledge based statistical scoring function or physics-based or a combination of both) which find the correct native conformation from misfolded decoys [3, 6, 9, 1215, 22, 3942]. However, it is rather nontrivial to develop a function that works across different decoy sets and a combination of functions is normally used [12, 13]. R-factor is the gold-standard for expressing the accuracy of crystallographic analysis, and as knowledge-based functions are mostly "trained" on crystal structures it is rather gratifying to develop functions similar to R-factor that can also be used to characterize the native structure (Table 2).

The present study demonstrates the development of scoring functions from the properties of residue packing that can be useful for discriminating the native conformation from various misfolded conformations for a given protein sequence. The algorithm assumes that a protein tries to take up a fold that has the minimum deviation of ASA (or PN) of each residue from the average value observed over all protein structures. The function Rs, based on residue accessibility, performs better than the one derived from the partner number, Rp, on decoy sets. The test on various decoy sets from the PROSTAR website demonstrated that the knowledge based scoring function developed in this study performs better or even at least of the same order than those previously derived by many authors [12, 14, 15]. Not only the present knowledge-based scoring functions pick the correct native structure in most cases, but the discrimination ratio is also better than that of the other potentials. However, as Equations (2) and (3) use the average values derived from a database of globular proteins, it is not likely to be very discriminatory for small proteins or peptides (as seen for the 'Ifu' set in Table 3). As such it would not be useful for checking local model quality in protein structures, as done by packages such as PROSA [43]. Along the same line it may be mentioned that the Verify3D server [44] for the visual analysis of the quality of a crystal structure works best on proteins with at least 100 residues.

The Park and Levitt decoy set had been shown to be quite a challenging dataset where the lowest-energy structures typically were 6-10 Å rmsd away from native ones [12]. The improved residue-based potential [18] also cannot recognize the native and near-native structures in all cases. The knowledge based scoring functions derived in this study are quite efficient to identify the near-native fold in Park and Levitt decoy sets. The correlation between the scoring function and rmsd is good in all cases and most of the cases the scoring functions have minimum value for the native structure. The scoring functions perform well also in the PROSTAR decoy sets, Levit's Local-Minima Decoy Sets (LMDS) and also in ROSETTA All-atom Decoy Sets. Considering 222 independent cases considered in this analysis Rs and Rp can efficiently discriminate native structures from all their corresponding decoys with a success rate greater than 85% and 74%, respectively. If we do not consider the 'Ifu' dataset, which comprises of small fragments of polypeptide chains, the success rate increases to 94% and 80%, respectively. The most rigorous test of a scoring function is to evaluate its performance in identifying the native structure with reference to the models submitted in CASP7 experiment. Even here, both Rs and Rp, the former in particular, stand out from all other methods (Table 5).

As our scoring functions depend on ASA or PN, these should be closely related to potentials of mean force derived from solvation or packing considerations. The performance of these potentials, however, depend critically on how the standard state is specified [6, 12, 23]. As the core and surface regions in proteins constitute distinct environments, potentials are sometimes divided into two parts, for the buried and the solvent-accessible regions [40]. The use of the average values of ASA or PN in globular proteins seems to have eliminated the need of such division, or the debate on the proper choice of the standard state.

A discussion on the uniqueness of our parameters vis-à-vis other knowledge-based discrimination functions is in order. First, a residue in the sequence is normally represented in these functions with one or two positions in the three-dimensional space and one or more of its properties, such as the secondary structure or backbone dihedral angle preferences, features in distance or sequence separation from other residues, etc. are considered [7, 23]. With such a coarse representation the function may not be as efficient as an all-atom discriminatory function, which takes into account the environment of all the atoms in a residue [13, 4547]. An all-atom representation is implicit in our method, as all the atoms are needed for the calculation of ASA or the partner number. However, each residue in the sequence contributes singly to the derivation of Rs or Rp. This is also in contrast to residue-residue interaction energy for each residue pair that is normally employed in other functions [12, 48, 49]. Furthermore, residue triplets and four-body contact potentials have also been developed [50, 51]. Secondly, the energy functions are generally less discriminatory when used individually, and the use of the hybrid scoring function is the norm for an enhanced performance [12, 16, 22]. While conceptually simple, Rs or Rp can work as efficiently. Thirdly, most formulations use energy as the criterion (with the assumption that the native structure is at a global free-energy minimum), while our function seeks to find the conformation that has the minimum deviation from the average value of the partner number or ASA. This way the selection of the most compact state of the polypeptide chain corresponding to a given sequence is achieved. The parameters are less likely to be fooled by over-abundance (which is penalized to the same extent as lower-abundance in equations 2 and 3) of contacts, as is the case with some functions [12]. Lastly, as the functions can identify the correct structure from the erroneous ones modeled from X-ray data ('Pdberr' set in Table 3) and vary within a narrow range in different protein classes (Table 2), these can be used for the validation of the structure determined crystallographically [52].

The functions developed here can also be used to delineate the compatibility of the sequence to a fold For example, azurin [53] and plastocyanin [54] are two small proteins having the same fold (a sandwich of two β-sheets having seven strands), but sequence identity of only 17% over an aligned length of 86 residues (Table 6). Expectedly, they have very similar Rs and Rp values. More interestingly however, when the sequence of plastocyanin is considered over the structure of azurin one gets a value of 0.97 for (Rs)azu/pcy, quite close to 0.89 obtained for the reverse process ((Rs)pcy/azu), thereby indicating the compatibility of the two sequences to the same fold.

Table 6 Rs and Rp for two proteins having the same fold belonging to the β class


This work demonstrates the effectiveness of a simple knowledge-based scoring function derived from residue packing for discriminating the native structures from a large set of decoys constructed by several groups. This knowledge-based scoring scheme is simple to derive and less computationally intensive than other energy functions and the performance is better (or at least at par) compared to others. Used in conjunction with other chemically intuitive parameter that captures the essence of the protein structure, it should be possible to achieve complete discrimination between the native structure and decoys.


Atomic coordinates were obtained from the Protein Data Bank (PDB) [55]. The analysis was carried out using the dataset of 432 polypeptide chains in 418 PDB files (given in [26]) with an R-factor ≤ 20%, a resolution ≤ 2.0 Å and sequence identity < 25%. Also the polypeptide chains with >40% of atoms with temperature factor (B-factor) >30 Å2 were excluded. The calculation of the partner number was restricted only to the well-ordered residues by excluding those with >40% atoms with temperature factor >30 Å2. The solvent accessible surface area (ASA) was computed using the program NACCESS [56], which is an implementation of the Lee and Richards algorithm [57]. The partner number of a residue is the number of other residues within a distance of 4.5 Å from any atom of the residue under consideration; the flanking residues were not considered as partner if the interaction was only with the main-chain atoms. The reason for the selection of the particular threshold value for the distance has been discussed [26, 58]. To be identified as a partner it is enough if just a pair of atoms is in contact.

Two parameters Rp and Rs based on the observed partner number and the accessibility at a given position in the protein sequence, as compared to the average value of the parameters for the same residue type in the whole database, were developed as given in the following two equations


where PNxi and ASAxi are the observed partner number and the solvent accessible surface area, respectively, for a residues of type x occurring at a particular position, i, in a PDB file and <PNx> and <ASAx> are the average values of the residue type x in the analyzed dataset. Considering (3), the function sums up the absolute value of the deviation of ASA at each position in the sequence from the average ASA of the residue type, each term being normalized by the average ASA value. The magnitude of each of the two parameters derived using (2) and (3) is used to discriminate the near native fold from the misfolded decoys. For the correct fold the values of these two parameters should be minimum.

A number of decoy datasets have been used from literature, the details of which are provided in Results. The Z-score of a native structure and the misfolded decoys was also evaluated. The Z-scores using the residue accessibility (Zs) and residue partner number (Zp) of a particular protein conformation are defined by the following equations


where Rs-nat (or Rp-nat) is the value of the parameter for the native conformation, and <Rs> (<Rp>) and σ are the average and the standard deviation of the distribution of the parameter in the set. The magnitude of the Z-score is an indication of how far that native conformation is separated from the near native structures in the distribution.


  1. 1.

    Bradley P, Misura KMS, Baker D: Toward high-resolution de novo structure prediction for small proteins. Science 2005, 309: 1868–1871. 10.1126/science.1113801

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Das R, Baker D: Macromolecular modeling with Rosetta. Annu Rev Biochem 2008, 77: 363–382. 10.1146/annurev.biochem.77.062906.171838

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Lazaridis T, Karplus M: Effective energy functions for protein structure prediction. Curr Opin Struct Biol 2000, 10: 139–145. 10.1016/S0959-440X(00)00063-4

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Jagielska A, Wroblewska L, Skolnick J: Protein model refinement using an optimized physics-based all-atom force field. Proc Natl Acad Sci USA 2008, 105: 8268–8273. 10.1073/pnas.0800054105

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  5. 5.

    Wodak SJ, Rooman MJ: Generating and testing protein folds. Curr Opin Struct Biol 1993, 3: 247–259. 10.1016/S0959-440X(05)80160-5

    CAS  Article  Google Scholar 

  6. 6.

    Sippl M: Knowledge based potentials for proteins. Curr Opin Struct Biol 1995, 5: 229–235. 10.1016/0959-440X(95)80081-6

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Sippl M: 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: 859–883. 10.1016/S0022-2836(05)80269-4

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Bowie J, Lüthy R, Eisenberg D: Method to identify protein sequences that fold into known three-dimensional structures. Science 1991, 253: 164–170. 10.1126/science.1853201

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Jones D, Taylor W, Thornton J: A new approach to protein fold recognition. Nature 1992, 258: 86–89. 10.1038/358086a0

    Article  Google Scholar 

  10. 10.

    Bryant S, Lawrence C: An empirical energy function for threading protein sequence through folding motif. Proteins 1993, 16: 92–112. 10.1002/prot.340160110

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Mirny LA, Shakhnovich EI: How to derive a protein folding potential? A new approach to an old problem. J Mol Biol 1996, 264: 1164–1179. 10.1006/jmbi.1996.0704

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Park B, Levitt M: Energy functions that discriminate X-ray and near-native folds from well-constructed decoys. J Mol Biol 1996, 258: 367–392. 10.1006/jmbi.1996.0256

    CAS  Article  PubMed  Google Scholar 

  13. 13.

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

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Lu H, Skolnick J: A distance-dependent atomic knowledge-based potential for improved protein structure selection. Proteins 2001, 44: 223–232. 10.1002/prot.1087

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Felts AK, Gallicchio E, Wallqvist A, Levy RM: Distinguishing native conformations of proteins from decoys with an effective free energy estimator based on the OPLS all-atom force field and the surface generalized Born solvent model. Proteins 2002, 48: 404–422. 10.1002/prot.10171

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Tsai J, Bonneau R, Morozov AV, Kuhlman B, Rohl CA, Baker D: An improved protein decoy set for testing energy functions for protein structure prediction. Proteins 2003, 53: 76–87. 10.1002/prot.10454

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Holm L, Sander CJ: Evaluation of protein models by atomic solvation preference. J Mol Biol 1992, 225: 93–105. 10.1016/0022-2836(92)91028-N

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Simons KT, Bonneau R, Ruczinski I, Baker D: Ab initio protein structure prediction of CASP III targets using ROSETTA. Proteins 1999, S3: 171–176. Publisher Full Text 10.1002/(SICI)1097-0134(1999)37:3+<171::AID-PROT21>3.0.CO;2-Z

    Article  Google Scholar 

  19. 19.

    Samudrala R, Levitt M: Decoys 'R' Us; A database of incorrect conformations to improved protein structure prediction. Protein Sci 2000, 9: 1399–1401. 10.1110/ps.9.7.1399

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  20. 20.

    Moult J, Fidelis K, Kryshtafovych A, Rost B, Hubbard T, Tramontano A: Critical assessment of methods of protein structure prediction-Round VII. Proteins 2007, 69(Suppl 8):3–9. 10.1002/prot.21767

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  21. 21.

    Fischer D: Servers for protein structure prediction. Curr Opin Struct Biol 2006, 6: 178–182. 10.1016/

    Article  Google Scholar 

  22. 22.

    Benkert P, Tosatto SC, Schomburg D: QMEAN: A comprehensive scoring function for model quality assessment. Proteins 2008, 71: 261–77. 10.1002/prot.21715

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Jernigan RL, Bahar I: Structure-derived potentials and protein simulations. Curr Opin Struct Biol 1996, 6: 195–209. 10.1016/S0959-440X(96)80075-3

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Glusker JP, Trueblood KN: Crystal Structure Analysis. A Primer. Oxford University Press, New York; 1985.

    Google Scholar 

  25. 25.

    Pal A, Bahadur RP, Ray PS, Chakrabarti P: Accessibility and partner number of protein residues, their relationship and a webserver, ContPlot for their display. BMC Bioinformatics 2009, 10: 103. 10.1186/1471-2105-10-103

    PubMed Central  Article  PubMed  Google Scholar 

  26. 26.

    Samanta U, Bahadur RP, Chakrabarti P: Quantifying the accessible surface area of protein residues in their local environment. Protein Eng 2002, 15: 659–667. 10.1093/protein/15.8.659

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Samanta U, Chakrabarti P: Assessing the role of tryptophan residues in the binding site. Protein Eng 2001, 14: 7–15. 10.1093/protein/14.1.7

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Sonavane S, Chakrabarti P: Cavities and atomic packing in protein structures and interfaces. PLoS Comput Biol 2008, 4(9):e1000188. 10.1371/journal.pcbi.1000188

    PubMed Central  Article  PubMed  Google Scholar 

  29. 29.

    Moult J, Unger R: An analysis of protein folding pathways. Biochemistry 1991, 30: 3816–3824. 10.1021/bi00230a003

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Mosimann S, Meleshko R, James MN: A critical assessment of comparative molecular modeling of tertiary structures of proteins. Proteins 1995, 23: 301–317. 10.1002/prot.340230305

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Braxenthaler M, Samudrala R, Pedersen J, Luo R, Milash B, Moult J: PROSTAR: The protein potential test site.1997. []

    Google Scholar 

  32. 32.

    Avbelj F, Moult J, Kitson DH, James MN, Hagler AT: Molecular dynamics study of the structure and dynamics of a protein molecule in a crystalline ionic environment, Streptomyces griseus protease A. Biochemistry 1990, 29: 8658–8676. 10.1021/bi00489a023

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Skolnick J, Kolinski A, Ortiz A: Derivation of protein-specific pair potentials based on weak sequence fragment similarity. Proteins 2000, 38: 3–16. Publisher Full Text 10.1002/(SICI)1097-0134(20000101)38:1<3::AID-PROT2>3.0.CO;2-S

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Levitt M, Hirshberg M, Sharon R, Daggett V: Potential energy function and parameters for simulation of the molecular dynamics of proteins and nucleic acids in solutions. Comput Phys Commun 1995, 91: 215–231. 10.1016/0010-4655(95)00049-L

    CAS  Article  Google Scholar 

  35. 35.

    Li X, Sutcliffe MJ, Schwartz TW, Dobson CM: Sequence-specific 1H NMR assignments and solution structure of bovine pancreatic polypeptide. Biochemistry 1992, 31: 1245–1253. 10.1021/bi00119a038

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Deisenhofer J: Crystallographic refinement and atomic models of a human Fc fragment and its complex with fragment B of protein A from Staphylococcus aureus at 2.9Å and 2.8 Å resolution. Biochemistry 1981, 20: 2361–2370. 10.1021/bi00512a001

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    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. 10.1006/jmbi.1997.0959

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    Levitt M: Accurate modeling of protein conformation by automatic segment matching. J Mol Biol 1992, 226: 507–533. 10.1016/0022-2836(92)90964-L

    CAS  Article  PubMed  Google Scholar 

  39. 39.

    Casari G, Sippl MJ: 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. 10.1016/0022-2836(92)90556-Y

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Kocher J-PA, Rooman MJ, Wodak SJ: Factors influencing the ability of knowledge-based potentials to identify native sequence-structure matches. J Mol Biol 1994, 235: 1598–1613. 10.1006/jmbi.1994.1109

    CAS  Article  PubMed  Google Scholar 

  41. 41.

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

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Melo F, Sanchez R, Sali A: Statistical potentials for fold assessment. Protein Sci 2002, 11: 430–448. 10.1110/ps.25502

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  43. 43.

    Wiederstein M, Sippl MJ: ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res 2007, (35 Web Server):W407–410. 10.1093/nar/gkm290

    Google Scholar 

  44. 44.

    Eisenberg D, Lüthy R, Bowie JU: VERIFY3D: assessment of protein models with three-dimensional profiles. Methods Enzymol 1997, 277: 396–404. full_text

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Melo F, Feytmans E: Novel knowledge-based mean force potential at atomic level. J Mol Biol 1997, 267: 207–222. 10.1006/jmbi.1996.0868

    CAS  Article  PubMed  Google Scholar 

  46. 46.

    McConkey BJ, Sobolev V, Edelman M: Discrimination of native protein structures using atom-atom contact scoring. Proc Natl Acad Sci USA 2003, 100: 3215–3220. 10.1073/pnas.0535768100

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  47. 47.

    Summa CM, Levitt M, DeGrado WF: An atomic environment potential for use in protein structure prediction. J Mol Biol 2005, 352: 986–1001. 10.1016/j.jmb.2005.07.054

    CAS  Article  PubMed  Google Scholar 

  48. 48.

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

    CAS  Article  Google Scholar 

  49. 49.

    Rajgaria R, McAllister SR, Floudas CA: Distance dependent centroid to centroid force fields using high resolution decoys. Proteins 2008, 70: 950–970. 10.1002/prot.21561

    CAS  Article  PubMed  Google Scholar 

  50. 50.

    Ngan SC, Inouye MT, Samudrala R: A knowledge-based scoring function based on residue triplets for protein structure prediction. Protein Eng Des Sel 2006, 19: 187–193. 10.1093/protein/gzj018

    CAS  Article  PubMed  Google Scholar 

  51. 51.

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

    CAS  Article  PubMed  Google Scholar 

  52. 52.

    Wilson KS, Butterworth S, Dauter Z, Lamzin VS, Walsh M, Wodak S, Pontius J, Richelle J, Vaguine A, Sander C, Hooft RWW, Vriend G, Thornton JM, Laskowski RA, MacArthur MW, Dodson EJ, Murshudov G, Oldfield TJ, Kaptein R, Rullmann JAC: Who checks the checkers? Four validation tools applied to eight atomic resolution structures. J Mo Biol 1998, 276: 417–436. 10.1006/jmbi.1997.1526

    Article  Google Scholar 

  53. 53.

    Adman ET, Jensen LH: Structural features of azurin at 2.7 Å resolution. Isr J Chem 1981, 21: 8–12.

    CAS  Article  Google Scholar 

  54. 54.

    Guss JM, Harrowell PR, Murata M, Norris VA, Freeman HC: Crystal structure analyses of reduced (CuI) poplar plastocyanin at six pH values. J Mol Biol 1986, 192: 361–387. 10.1016/0022-2836(86)90371-2

    CAS  Article  PubMed  Google Scholar 

  55. 55.

    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. 10.1093/nar/28.1.235

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  56. 56.

    Hubbard SJ: NACCESS: program for calculating accessibilities. Department of Biochemistry and Molecular Biology University College of London; 1992. []

    Google Scholar 

  57. 57.

    Lee B, Richards FM: The interpretation of protein structures: estimation of static accessibility. J Mol Biol 1971, 55: 379–400. 10.1016/0022-2836(71)90324-X

    CAS  Article  PubMed  Google Scholar 

  58. 58.

    Chakrabarti P, Bhattacharyya R: Geometry of nonbonded interactions involving planar groups in proteins. Prog Biophys Mol Biol 2007, 95: 83–137. 10.1016/j.pbiomolbio.2007.03.016

    CAS  Article  PubMed  Google Scholar 

  59. 59.

    Orengo CA, Michie AD, Jones S, Jones DT, Swindells MB, Thornton JM: CATH-A hierarchic classification of protein domain structures. Structure 1997, 5: 1093–1108. 10.1016/S0969-2126(97)00260-8

    CAS  Article  PubMed  Google Scholar 

Download references


We are grateful to the anonymous reviewers for their comments on the manuscript. The work was supported by a grant from the Department of Biotechnology, India. RPB thanks SRIC of IIT, Kharagpur for a startup grant.

Author information



Corresponding author

Correspondence to Pinak Chakrabarti.

Additional information

Authors' contributions

PC conceptualized the work that was carried out by RPB. RPB and PC participated in interpretation of the data and writing the manuscript. Both the authors read and approved the final manuscript.

Electronic supplementary material

Identification of native structure from decoys in different decoy sets

Additional file 1: . The file contains three tables, numbered S1 to S3. (DOC 261 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

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

Bahadur, R.P., Chakrabarti, P. Discriminating the native structure from decoys using scoring functions based on the residue packing in globular proteins. BMC Struct Biol 9, 76 (2009).

Download citation


  • Native Structure
  • Accessible Surface Area
  • Native Conformation
  • Solvent Accessible Surface Area
  • Protein Data Bank File