Email updates

Keep up to date with the latest news and content from BMC Bioinformatics and BioMed Central.

This article is part of the supplement: 22nd International Conference on Genome Informatics: Bioinformatics

Open Access Proceedings

A new protein-ligand binding sites prediction method based on the integration of protein sequence conservation information

Tianli Dai1, Qi Liu1, Jun Gao12, Zhiwei Cao13* and Ruixin Zhu14*

Author Affiliations

1 College of Life Science and Biotechnology, Tongji University, 200092, Shanghai, China

2 College of Information Engineering, Shanghai Maritime University, 201306, Shanghai, China

3 Shanghai Center for Bioinformation and Technology, 100 Qinzhou Road, Shanghai, 200235, China

4 Department of Chinese Material Medica, Liaoning University of Traditional Chinese Medicine, Shenyang, Liaoning 110032, China

For all author emails, please log on.

BMC Bioinformatics 2011, 12(Suppl 14):S9  doi:10.1186/1471-2105-12-S14-S9


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2105/12/S14/S9


Published:14 December 2011

© 2011 Dai et al; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Prediction of protein-ligand binding sites is an important issue for protein function annotation and structure-based drug design. Nowadays, although many computational methods for ligand-binding prediction have been developed, there is still a demanding to improve the prediction accuracy and efficiency. In addition, most of these methods are purely geometry-based, if the prediction methods improvement could be succeeded by integrating physicochemical or sequence properties of protein-ligand binding, it may also be more helpful to address the biological question in such studies.

Results

In our study, in order to investigate the contribution of sequence conservation in binding sites prediction and to make up the insufficiencies in purely geometry based methods, a simple yet efficient protein-binding sites prediction algorithm is presented, based on the geometry-based cavity identification integrated with sequence conservation information. Our method was compared with the other three classical tools: PocketPicker, SURFNET, and PASS, and evaluated on an existing comprehensive dataset of 210 non-redundant protein-ligand complexes. The results demonstrate that our approach correctly predicted the binding sites in 59% and 75% of cases among the TOP1 candidates and TOP3 candidates in the ranking list, respectively, which performs better than those of SURFNET and PASS, and achieves generally a slight better performance with PocketPicker.

Conclusions

Our work has successfully indicated the importance of the sequence conservation information in binding sites prediction as well as provided a more accurate way for binding sites identification.

Background

Proteins are the material basis of all life, the key components of body cells, and play important roles in the process of life activity. Since in most cellular processes, proteins interact with other molecules to perform their biological functions, the successful identification of ligand-binding sites on protein surfaces becomes vital and necessary to explore the proteins comprehensively [1]. In addition, as a result of various structural genomics projects performed, structural information of proteins with little or no functional annotations has been explosively increasing. Such increasingly accumulated data have become to attract much more interests in exploring the relationship between protein structure and function as well as elucidating the functions from their structures rather than merely from sequences.

In recent decades, many computational methods have been developed for candidate binding sites identification. Briefly, these algorithms can be divided into three categories, i.e.(1) purely geometry based methods, which follow the assumption that the protein-ligand binding sites are generally located at crevices on the protein surface or cavities in the protein. When the shapes of protein surface were calculated, it can be easily to predict the candidate protein-ligand binding sites without any ligands information. Methods following in this category include POCKET [2], LIGSITE [3], PASS [4], SURFNET [5] and, PocketPicker [6] etc. It is worth noting that this kind of methods focuses only on the shapes of protein surface without considering the physicochemical properties of amino acids. What’s more, a major number among these algorithms are based on the cubic grid representation, which means that their following results are often protein orientation dependent; (2) energetic based methods, which coat the protein surface with a layer of probes to calculate van der Waals interaction energies between the protein and probes. The energetically favorable probe sites are clustered according to their spatial proximity. Then the identified clusters are ranked according to the sum of the interaction energies within each cluster. As an example, Q-SiteFinder [7] is a classical tool following in this category; (3) knowledge based methods, which including various statistical methods [8], machine learning methods [9] and similarity comparison method set. Besides, a part of them predict protein-ligand binding sites by searching for clusters or patterns of conserved residues [10,11]. These method stake the assumption that the residues located in protein-ligand binding site usually being more important and more highly conserved than those located in other parts through evolution. Although the results for certain methods with only sequence conservation information are not satisfactory [12,13], it is still expected to be helpful in re-ranking the pockets in the process of prediction [14].

In summary, in this study, in order to investigate the contribution of sequence conservation information in binding sites prediction and to make up the insufficiencies for purely geometry based methods, we aims at designing a simple, yet efficient and practical binding site prediction algorithm based on the integration of sequence conservation information with geometry-based cleft identification.

Methods

Algorithm workflow

An overview of our method is shown in Figure 1. It is composed of three steps: (1) Calculation of geometrical characteristics of protein (cleft identification); (2) Filtering with sequence conservation information, and (3) Clustering potential atoms which will form the prediction binding sites according to their spatial distance-based similarity.

thumbnailFigure 1. The flowchart of proposed algorithm. Overview of our method. The prediction is based on the geometry-based cavity identification integrated with sequence conservation information.

Step 1: for each protein structure, the solvent accessible surface area (SASA) of each atom is calculated first. The SASA values can be available from many tools such as NACCESS [15], ASC [16], Surface Racer [17] etc. Here Surface Racer is used because it can exactly calculate accessible surface area on most common computer platforms. An atom is considered as an interface alone if its SASA is over 1Å [18]. In our study, these atoms are denoted as Set S (Surface). In addition, Qhull [19] is applied to calculate the convex hull based on all the protein atoms [20,21]. Considering the convex hull obtained here is just a smallest convex set of atom points, it is expanded by adding atoms which are within certain distance from the origin convex atoms (the parameter is adjustable, 6.5Å is used here). Such new convex set is denoted as Set Q (Qhull). Finally, after calculating those atoms which included in Set S and excluded from Set Q, it is convenient to get a set of protein atoms which locate at protein’s crevice regions.

Step 2: since the first step just identifies the clefts on the protein surface without any biological significance, the sequence conservation information are further added as a filter [8,14,22] to curate our results. This is achieved by the ConSurf-DB [23] which provides the pre-calculated evolutionary conservation profiles for proteins with known structures in the PDB. In ConSurf-DB, every residue in every corresponding protein is evaluated with a normalized conservation score. And then the normalized scores are binned into the 1-9 color scales for representing the conservation grades and projected on the 3D model of the query protein, where 1 corresponds to maximal variability and 9 to maximal conservation [23]. It is important to note that although the same color scale in ConSurf-DB is used in all the protein families, the conservation scores are not absolute and hence, defining the conservation scores as a filter between different protein families might be misleading. Accordingly, in this study, ConSurf-DB results are interpreted using the color scales rather than the conservation scores, and only the residues greater than or equal to certain conservation grade cutoff (such as 7, 8, 9. 7 is used here) will be retained [13]. We denote all the atoms of those conservative residues as Set C (Conservative atoms). After such physicochemical property as a biological factor, those atoms which appear on protein’s clefts will accordingly hold the sequence conservation information. These atoms are denoted as Set P (Potential atoms).

Step 3: a simple hierarchical clustering algorithm is applied to cluster those potential atoms according to their spatial distances. Each output cluster stands for a presumed protein-ligand binding site, and the center of each cluster represents the geometric center of each binding site. Geometric centers within a certain distance threshold (8Å used here) are grouped together as a new cluster [24] whereas the corresponding geometric center should be recalculated.

Besides the candidate protein-ligand binding site identification itself, binding site ranking is also a very important tissue. For instance, since there are often several presumed binding sites that can be detected on a protein surface, in order to select the more relevant ones, it is necessary to derive an approach to characterize and rank them. It is often said that the largest pocket tends to frequently correspond to the observed ligand binding site [25]. Based on this assumption, a most number of prediction methods rank the candidate sites by pocket size. On the other hand, different studies have tried to solve this ranking problem from other perspectives [14,26,27]. Among them, evolutionary information such as sequence conservation has been shown to be successful for re-ranking the binding sites [14]. Therefore, in our study, the candidate binding sites are ranked according to the conservation score of all residues in the same cluster.

Test dataset

In this study, a regularly used dataset [14] is chosen as the standard test data, which consists of 210 non-duplicated protein-ligand complexes derived from the Protein Ligand Database (PLD) [28]. A rough statistics on the protein dataset classification has been shown in Figure 2. Furthermore, in order to assess the binding site prediction performance of our method, the identified sites are needed to be compared to the real binding sites. For the 210 bound proteins, the ligands are taken away when making predictions and then put back when performing evaluation. The PocketPicker criterion (PPc) [6] is adopted as the prediction criterion in this study. It defines the prediction to be a hit if the geometric center of the presumed binding site is within 4Å from any atom of the ligand.

thumbnailFigure 2. Distribution of the protein dataset by molecular function. This is a rough statistics on the protein dataset classification.

Results and discussion

Our algorithm is tested on 210 protein-ligand complexes. The presumed protein-ligand binding sites are ranked by conservation score. A successful prediction example [PDB: 6RNT] [29] visualized with Jmol [30] is presented in Figure 3. In order to check the contribution of sequence conservation as well as the feasibility of our method, three purely geometry based methods i.e. PocketPicker, SURFNET, and PASS with their own ranking methods are also tested for comparison [6,24]. The accuracy of the first one (TOP 1) and first three (TOP 3) in the prediction ranking lists have been calculated. It is indicated from Table 1 that our method obtained a 59% success rate for the top one prediction which means almost 124 of the 210 proteins are correctly predicted. The top one result is much higher than that of SURFNET and PASS. Although the top three success rate seems to be a little bit worse than PASS, our method still performs better than others. On the whole, the success rate in our study is comparable to that of PocketPicker which is one of the most popular prediction tools presented in 2007 while our method with the quick-reading operating process and grid-presentation independent.

thumbnailFigure 3. One case study of our method. PDB ID: 6RNT. (Red points: water molecule; Light blue: the whole protein; Golden: ligand molecular; Aquamarine : binding site’s center and Purple: predicted binding site constituted by amino acids.)

Table 1. Prediction success rate presented by different binding-sites prediction methods

As small molecular ligands are tend to combine with proteins in larger cavities on protein surface, the volume can be used as a ranking method to choose the likely candidates. In this study, the candidate binding sites are also ranked according to the space volume. The success rates are listed in the Volume column in Table 1. It can be seen that this kind of ranking method doesn’t show any advantage to that by conservation score. What’s more, the top one success rate ranked by volume hardly achieves 50%. It indicates that such volume ranking rule can’t be generalized with its own limit.

In addition, two factors, i.e. expand distance from origin convex hull set and the conservation color scale (ConCS) are tested for their influence in our study. The top one and top three success rates under different combinations of these two factors are shown in Figure 4 and Figure 5, respectively. The success rates are derived according to the cavities’ space volume in a sequential manner. It can be seen that no matter how the conservation score is set, the expand distance at 6.5 Å always performs a better success rate. The data of ConCS>=9cannot be available because when under such special condition, the number of the candidate atoms in some protein structures will be too small to form predicted protein-ligand binding sites. Moreover, compared to those candidate atoms without evolutionary information (the No ConCS line), we find that when the sequence conservation is introduced as a filter, the success rates in top 1 and top 3 are all improved significantly. However, it also can be seen that the higher conservation color scale the candidate atoms obtain doesn’t often result the higher prediction accuracy (the red line tends to be above the green line). We explained this by that the atoms forming the ligand binding sites are not only the most conservative ones but also the relatively conservative ones, which should be validated further more.

thumbnailFigure 4. TOP1 Success rates achieved by setting different parameters. The accuracy of the first one pocket sites (TOP 1) in the prediction ranking lists was different under different parameter combinations.

thumbnailFigure 5. TOP3 Success rates achieved by setting different parameters. The accuracy of the first three pocket sites (TOP 3) in the prediction ranking lists was different under different parameter combinations.

In summary, our study has validated the insufficiency of purely geometric methods, and at the same time, reflected the significance of sequence conservation in ligand binding sites prediction.

Conclusions

The prediction of protein-ligand binding sites has great significance for the protein function annotation and computer-aided drug design. Though many different outstanding studies have been carried out to solve this problem, some of them just use complicated calculation methods based on protein shape descriptor rather than considering other physicochemical and sequence properties with biological characteristics. In this paper, a simple yet efficient binding site prediction algorithm is designed based on the integration of geometry and sequence conservation information. The algorithm is tested on a regularly used benchmark dataset, and shows an encouraging result with the success rates come to 59% and 75% for the top one and top three pocket sites, respectively. Our algorithm performs comparative to PocketPicker while with more convenient prediction procedure. Last but not least, our result also reflects the un-ignorable importance of sequence conservation information which can be an effective attribute in ligand binding site prediction.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Conceived and designed the experiments: TD, ZC, RZ. Performed the experiments: TD, ZC, RZ. Analyzed the data: TD, QL, JG, ZC, RZ. Wrote the first draft: TD. Revised the draft: TD, QL, JG, ZC, RZ.

Acknowledgements

We are grateful to members of the same research group for helpful discussions and thank the reviewers for their valuable suggestions to improve the work. This work was supported in part by grants from National Natural Science Foundation of China (30976611, 31100956, 61173117), Program for New Century Excellent Talents in University (NCET-08-0399), Research Fund for the Doctoral Program of Higher Education of China (20100072110008, 20100072120050), Shanghai Magnolia Funding (2010B127), Shanghai Pujiang Talents Funding (Grant No. 11PJ1407400) and Tongji Excellent Young Scientist Funding (Grant No. 2000219052).

This article has been published as part of BMC Bioinformatics Volume 12 Supplement 14, 2011: 22nd International Conference on Genome Informatics: Bioinformatics. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/12?issue=S14.

References

  1. Thornton JM, Todd AE, Milburn D, Borkakoti N, Orengo CA: From structure to function: approaches and limitations.

    Nat Struct Biol 2000, 7(Suppl):991-994. PubMed Abstract | Publisher Full Text OpenURL

  2. Levitt DG, Banaszak LJ: POCKET: a computer graphics method for identifying and displaying protein cavities and their surrounding amino acids.

    J Mol Graph 1992, 10(4):229-234. PubMed Abstract | Publisher Full Text OpenURL

  3. Hendlich M, Rippmann F, Barnickel G: LIGSITE: automatic and efficient detection of potential small molecule-binding sites in proteins.

    J Mol Graph Model 1997, 15(6):359-363.

    389

    PubMed Abstract | Publisher Full Text OpenURL

  4. Brady GP Jr., Stouten PF: Fast prediction and visualization of protein binding pockets with PASS.

    J Comput Aided Mol Des 2000, 14(4):383-401. PubMed Abstract | Publisher Full Text OpenURL

  5. Laskowski RA: SURFNET: a program for visualizing molecular surfaces, cavities, and intermolecular interactions.

    J Mol Graph 1995, 13(5):323-330.

    307-328

    PubMed Abstract | Publisher Full Text OpenURL

  6. Weisel M, Proschak E, Schneider G: PocketPicker: analysis of ligand binding-sites with shape descriptors.

    Chem Cent J 2007, 1:7. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  7. Laurie AT, Jackson RM: Q-SiteFinder: an energy-based method for the prediction of protein-ligand binding sites.

    Bioinformatics 2005, 21(9):1908-1916. PubMed Abstract | Publisher Full Text OpenURL

  8. Zhou YQ, Liang SD, Zhang C, Liu S: Protein binding site prediction using an empirical scoring function.

    Nucleic Acids Res 2006, 34(13):3698-3707. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Sonavane S, Chakrabarti P: Prediction of active site cleft using support vector machines.

    Journal of Chemical Information and Modeling 2010, 50(12):2266-2273. PubMed Abstract | Publisher Full Text OpenURL

  10. Capra JA, Singh M: Predicting functionally important residues from sequence conservation.

    Bioinformatics 2007, 23(15):1875-1882. PubMed Abstract | Publisher Full Text OpenURL

  11. Manning JR, Jefferson ER, Barton GJ: The contrasting properties of conservation and correlated phylogeny in protein functional residue prediction.

    BMC Bioinformatics 2008, 9:51. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  12. Caffrey DR, Somaroo S, Hughes JD, Mintseris J, Huang ES: Are protein-protein interfaces more conserved in sequence than the rest of the protein surface?

    Protein Sci 2004, 13(1):190-202. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Prymula K, Jadczyk T, Roterman I: Catalytic residues in hydrolases: analysis of methods designed for ligand-binding site prediction.

    J Comput Aided Mol Des 2011, 25(2):117-133. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Huang B, Schroeder M: LIGSITEcsc: predicting ligand binding sites using the Connolly surface and degree of conservation.

    BMC Struct Biol 2006, 6:19. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  15. Naccess Home Page [http://www.bioinf.manchester.ac.uk/naccess/] webcite

  16. Eisenhaber F, P Argos: Improved strategy in analytic surface calculation for molecular systems: handling of singularities and computational efficiency.

    Journal of Computational Chemistry 2004, 14(11):1272-1280. OpenURL

  17. Tsodikov OV, Record MT Jr., Sergeev YV: Novel computer program for fast exact calculation of accessible and molecular surface areas and average surface curvature.

    J Comput Chem 2002, 23(6):600-609. PubMed Abstract | Publisher Full Text OpenURL

  18. Wu D, Sun J, Xu T, Wang S, Li G, Li Y, Cao Z: Stacking and energetic contribution of aromatic islands at the binding interface of antibody proteins.

    Immunome Res 2010, 6(Suppl 1):S1. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  19. Barber CB, Dobkin DP, Huhdanpaa H: The Quickhull algorithm for convex hulls.

    Acm T Math Software 1996, 22(4):469-483. Publisher Full Text OpenURL

  20. Petrek M, Otyepka M, Banas P, Kosinova P, Koca J, Damborsky J: CAVER: a new tool to explore routes from protein clefts, pockets and cavities.

    BMC Bioinformatics 2006, 7:316. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  21. Coleman RG, Sharp KA: Travel depth, a new shape descriptor for macromolecules: application to ligand binding.

    J Mol Biol 2006, 362(3):441-458. PubMed Abstract | Publisher Full Text OpenURL

  22. Magliery TJ, Regan L: Sequence variation in ligand binding sites in proteins.

    BMC Bioinformatics 2005, 6:240. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  23. Goldenberg O, Erez E, Nimrod G, Ben-Tal N: The ConSurf-DB: pre-calculated evolutionary conservation profiles of protein structures.

    Nucleic Acids Res 2009, 37(Database issue):D323-327. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Huang B: MetaPocket:a meta approach to improve protein ligand binding site prediction.

    OMICS 2009, 13(4):325-330. PubMed Abstract | Publisher Full Text OpenURL

  25. Nayal M, Honig B: On the nature of cavities on protein surfaces: application to the identification of drug-binding sites.

    Proteins 2006, 63(4):892-906. PubMed Abstract | Publisher Full Text OpenURL

  26. An J, Totrov M, Abagyan R: Comprehensive identification of "druggable" protein ligand binding sites.

    Genome Inform 2004, 15(2):31-41. PubMed Abstract OpenURL

  27. Zhong S, MacKerell AD Jr.: Binding response: a descriptor for selecting ligand binding site on protein surfaces.

    J Chem Inf Model 2007, 47(6):2303-2315. PubMed Abstract | Publisher Full Text OpenURL

  28. Puvanendrampillai D, Mitchell JB: L/D Protein Ligand Database (PLD): additional understanding of the nature and specificity of protein-ligand complexes.

    Bioinformatics 2003, 19(14):1856-1857. PubMed Abstract | Publisher Full Text OpenURL

  29. Ding J, Koellner G, Grunert HP, Saenger W: Crystal structure of ribonuclease T1 complexed with adenosine 2'-monophosphate at 1.8-A resolution.

    J Biol Chem 1991, 266(23):15128-15134. PubMed Abstract | Publisher Full Text OpenURL

  30. Jmol: an open-source Java viewer for chemical structures in 3D [http://www.jmol.org/] webcite