A new protein-ligand binding sites prediction method based on the integration of protein sequence conservation information
© Dai et al; licensee BioMed Central Ltd. 2011
Published: 14 December 2011
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.
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.
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.
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 . 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 , LIGSITE , PASS , SURFNET  and, PocketPicker  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  is a classical tool following in this category; (3) knowledge based methods, which including various statistical methods , machine learning methods  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 .
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.
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 , ASC , Surface Racer  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Å . In our study, these atoms are denoted as Set S (Surface). In addition, Qhull  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  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 . 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 . 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  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 . 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 . Therefore, in our study, the candidate binding sites are ranked according to the conservation score of all residues in the same cluster.
Results and discussion
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 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.
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.
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.
- 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.View ArticlePubMedGoogle Scholar
- 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. 10.1016/0263-7855(92)80074-NView ArticlePubMedGoogle Scholar
- 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 389 10.1016/S1093-3263(98)00002-3View ArticlePubMedGoogle Scholar
- 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. 10.1023/A:1008124202956View ArticlePubMedGoogle Scholar
- Laskowski RA: SURFNET: a program for visualizing molecular surfaces, cavities, and intermolecular interactions. J Mol Graph 1995, 13(5):323–330. 307–328 307-328 10.1016/0263-7855(95)00073-9View ArticlePubMedGoogle Scholar
- Weisel M, Proschak E, Schneider G: PocketPicker: analysis of ligand binding-sites with shape descriptors. Chem Cent J 2007, 1: 7. 10.1186/1752-153X-1-7PubMed CentralView ArticlePubMedGoogle Scholar
- Laurie AT, Jackson RM: Q-SiteFinder: an energy-based method for the prediction of protein-ligand binding sites. Bioinformatics 2005, 21(9):1908–1916. 10.1093/bioinformatics/bti315View ArticlePubMedGoogle Scholar
- 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. 10.1093/nar/gkl454PubMed CentralView ArticlePubMedGoogle Scholar
- Sonavane S, Chakrabarti P: Prediction of active site cleft using support vector machines. Journal of Chemical Information and Modeling 2010, 50(12):2266–2273. 10.1021/ci1002922View ArticlePubMedGoogle Scholar
- Capra JA, Singh M: Predicting functionally important residues from sequence conservation. Bioinformatics 2007, 23(15):1875–1882. 10.1093/bioinformatics/btm270View ArticlePubMedGoogle Scholar
- Manning JR, Jefferson ER, Barton GJ: The contrasting properties of conservation and correlated phylogeny in protein functional residue prediction. BMC Bioinformatics 2008, 9: 51. 10.1186/1471-2105-9-51PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1110/ps.03323604PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1007/s10822-010-9402-0PubMed CentralView ArticlePubMedGoogle Scholar
- Huang B, Schroeder M: LIGSITEcsc: predicting ligand binding sites using the Connolly surface and degree of conservation. BMC Struct Biol 2006, 6: 19. 10.1186/1472-6807-6-19PubMed CentralView ArticlePubMedGoogle Scholar
- Naccess Home Page[http://www.bioinf.manchester.ac.uk/naccess/]
- 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.View ArticleGoogle Scholar
- 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. 10.1002/jcc.10061View ArticlePubMedGoogle Scholar
- 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. 10.1186/1745-7580-6-S1-S1PubMed CentralView ArticlePubMedGoogle Scholar
- Barber CB, Dobkin DP, Huhdanpaa H: The Quickhull algorithm for convex hulls. Acm T Math Software 1996, 22(4):469–483. 10.1145/235815.235821View ArticleGoogle Scholar
- 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. 10.1186/1471-2105-7-316PubMed CentralView ArticlePubMedGoogle Scholar
- Coleman RG, Sharp KA: Travel depth, a new shape descriptor for macromolecules: application to ligand binding. J Mol Biol 2006, 362(3):441–458. 10.1016/j.jmb.2006.07.022View ArticlePubMedGoogle Scholar
- Magliery TJ, Regan L: Sequence variation in ligand binding sites in proteins. BMC Bioinformatics 2005, 6: 240. 10.1186/1471-2105-6-240PubMed CentralView ArticlePubMedGoogle Scholar
- 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 CentralView ArticlePubMedGoogle Scholar
- Huang B: MetaPocket:a meta approach to improve protein ligand binding site prediction. OMICS 2009, 13(4):325–330. 10.1089/omi.2009.0045View ArticlePubMedGoogle Scholar
- 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. 10.1002/prot.20897View ArticlePubMedGoogle Scholar
- An J, Totrov M, Abagyan R: Comprehensive identification of "druggable" protein ligand binding sites. Genome Inform 2004, 15(2):31–41.PubMedGoogle Scholar
- 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. 10.1021/ci700149kView ArticlePubMedGoogle Scholar
- 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. 10.1093/bioinformatics/btg243View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Jmol: an open-source Java viewer for chemical structures in 3D[http://www.jmol.org/]
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.