Automated identification of binding sites for phosphorylated ligands in protein structures

Phosphorylation is a crucial step in many cellular processes, ranging from metabolic reactions involved in energy transformation to signaling cascades. In many instances, protein domains specifically recognize the phosphogroup. Knowledge of the binding site provides insights into the interaction, and it can also be exploited for therapeutic purposes. Previous studies have shown that proteins interacting with phosphogroups are highly heterogeneous, and no single property can be used to reliably identify the binding site. Here we present an energy‐based computational procedure that exploits the protein three‐dimensional structure to identify binding sites involved in the recognition of phosphogroups. The procedure is validated on three datasets containing more than 200 proteins binding to ATP, phosphopeptides, and phosphosugars. A comparison against other three generic binding site identification approaches shows higher accuracy values for our method, with a correct identification rate in the 80–90% range for the top three predicted sites. Addition of conservation information further improves the performance. The method presented here can be used as a first step in functional annotation or to guide mutagenesis experiments and further studies such as molecular docking. Proteins 2012;. © 2012 Wiley Periodicals, Inc.     

Construction of protein binding sites in scaffold structures

We have developed a strategy for grafting a protein-protein interface based on the known crystal structure of a native ligand and receptor proteins in a complex. The key interaction residues at the ligand protein binding interface are transferred onto a scaffold protein so that the mutated scaffold protein will bind the receptor protein in the same manner as the ligand protein. First, our method identifies key residues and atoms in the ligand protein, which strongly interact with the receptor protein. Second, this method searches the scaffold protein for combinations of candidate residues, among which the distance between any two candidate residues is similar to that between relevant key interaction residues in the ligand protein. These candidate residues are mutated to key interaction residues in the ligand protein respectively. The scaffold protein is superposed onto the ligand protein based upon the coordinates of corresponding atoms, which are assumed to strongly interact with the receptor protein. Complementarity between scaffold and receptor proteins is evaluated. Scaffold proteins with a low superposing rms difference and high complementary score are accepted for further analysis. Then, the relative position of the scaffold protein is adjusted so that the interfaces between the scaffold and receptor proteins have a reasonable packing density. Other mutations are also considered to reduce the desolvation energy or bad steric contacts. Finally, the scaffold protein is cominimized with the receptor protein and evaluated. To test the method, the binding interface of barstar, the inhibitor of barnase, was grafted onto small proteins. Four scaffold proteins with high complementary scores are accepted.     

Analysis of zinc binding sites in protein crystal structures

The geometrical properties of zinc binding sites in a dataset of high quality protein crystal structures deposited in the Protein Data Bank have been examined to identify important differences between zinc sites that are directly involved in catalysis and those that play a structural role. Coordination angles in the zinc primary coordination sphere are compared with ideal values for each coordination geometry, and zinc coordination distances are compared with those in small zinc complexes from the Cambridge Structural Database as a guide of expected trends. We find that distances and angles in the primary coordination sphere are in general close to the expected (or ideal) values. Deviations occur primarily for oxygen coordinating atoms and are found to be mainly due to H-bonding of the oxygen coordinating ligand to protein residues, bidentate binding arrangements, and multi-zinc sites. We find that H-bonding of oxygen containing residues (or water) to zinc bound histidines is almost universal in our dataset and defines the elec-His-Zn motif. Analysis of the stereochemistry shows that carboxyl elec-His-Zn motifs are geometrically rigid, while water elec-His-Zn motifs show the most geometrical variation. As catalytic motifs have a higher proportion of carboxyl elec atoms than structural motifs, they provide a more rigid framework for zinc binding. This is understood biologically, as a small distortion in the zinc position in an enzyme can have serious consequences on the enzymatic reaction. We also analyze the sequence pattern of the zinc ligands and residues that provide elecs, and identify conserved hydrophobic residues in the endopeptidases that also appear to contribute to stabilizing the catalytic zinc site. A zinc binding template in protein crystal structures is derived from these observations.     

Beyond structural genomics: computational approaches for the identification of ligand binding sites in protein structures

Structural genomics projects have revealed structures for a large number of proteins of unknown function. Understanding the interactions between these proteins and their ligands would provide an initial step in their functional characterization. Binding site identification methods are a fast and cost-effective way to facilitate the characterization of functionally important protein regions. In this review we describe our recently developed methods for binding site identification in the context of existing methods. The advantage of energy-based approaches is emphasized, since they provide flexibility in the identification and characterization of different types of binding sites.     

Genome-wide identification of protein binding sites on RNAs in mammalian cells

RNA-binding proteins (RBPs) are proteins that bind to the RNA and participate in forming ribonucleoprotein complexes. They have crucial roles in various biological processes such as RNA splicing, editing, transport, maintenance, degradation, intracellular localization and translation. The RBPs bind RNA with different RNA-sequence specificities and affinities, thus, identification of protein binding sites on RNAs (R-PBSs) will deeper our understanding of RNA-protein interactions. Currently, high-throughput sequencing of RNA isolated by crosslinking immunoprecipitation (HITS-CLIP, also known as CLIP-Seq) is one of the most powerful methods to map RNA-protein binding sites or RNA modification sites. However, this method is only used for identification of single known RBPs and antibodies for RBPs are required. Here we developed a novel method, called capture of protein binding sites on RNAs (RPBS-Cap) to identify genome-wide protein binding sites on RNAs without using antibodies. Double click strategy is used for the RPBS-Cap assay. Proteins and RNAs are UV-crosslinked in vivo first, then the proteins are crosslinked to the magnetic beads. The RNA elements associated with proteins are captured, reverse transcribed and sequenced. Our approach has potential applications for studying genome-wide RNA-protein interactions.     

Prediction of protein binding sites in protein structures using hidden Markov support vector machine

Background  

Predicting the binding sites between two interacting proteins provides important clues to the function of a protein. Recent research on protein binding site prediction has been mainly based on widely known machine learning techniques, such as artificial neural networks, support vector machines, conditional random field, etc. However, the prediction performance is still too low to be used in practice. It is necessary to explore new algorithms, theories and features to further improve the performance.     

PocketMatch: A new algorithm to compare binding sites in protein structures

Background  

Recognizing similarities and deriving relationships among protein molecules is a fundamental requirement in present-day biology. Similarities can be present at various levels which can be detected through comparison of protein sequences or their structural folds. In some cases similarities obscure at these levels could be present merely in the substructures at their binding sites. Inferring functional similarities between protein molecules by comparing their binding sites is still largely exploratory and not as yet a routine protocol. One of the main reasons for this is the limitation in the choice of appropriate analytical tools that can compare binding sites with high sensitivity. To benefit from the enormous amount of structural data that is being rapidly accumulated, it is essential to have high throughput tools that enable large scale binding site comparison.     

Novel metal sites in protein structures

Recently, the three-dimensional structures of several novel metalloenzymes have been solved. Of special interest are those containing uncommon and/or not well characterized metals such as molybdenum, tungsten, nickel, vanadium and cobalt. Modulated by the protein environment, the specific properties of these metals and of special metal-binding cofactors such as siroheme and topa quinone are used to catalyze a vast array of fascinating reactions.     

Automatic identification and representation of protein binding sites for molecular docking.

Molecular docking is a popular way to screen for novel drug compounds. The method involves aligning small molecules to a protein structure and estimating their binding affinity. To do this rapidly for tens of thousands of molecules requires an effective representation of the binding region of the target protein. This paper presents an algorithm for representing a protein's binding site in a way that is specifically suited to molecular docking applications. Initially the protein's surface is coated with a collection of molecular fragments that could potentially interact with the protein. Each fragment, or probe, serves as a potential alignment point for atoms in a ligand, and is scored to represent that probe's affinity for the protein. Probes are then clustered by accumulating their affinities, where high affinity clusters are identified as being the "stickiest" portions of the protein surface. The stickiest cluster is used as a computational binding "pocket" for docking. This method of site identification was tested on a number of ligand-protein complexes; in each case the pocket constructed by the algorithm coincided with the known ligand binding site. Successful docking experiments demonstrated the effectiveness of the probe representation.     

Localization of binding sites in protein structures by optimization of a composite scoring function

The rise in the number of functionally uncharacterized protein structures is increasing the demand for structure-based methods for functional annotation. Here, we describe a method for predicting the location of a binding site of a given type on a target protein structure. The method begins by constructing a scoring function, followed by a Monte Carlo optimization, to find a good scoring patch on the protein surface. The scoring function is a weighted linear combination of the z-scores of various properties of protein structure and sequence, including amino acid residue conservation, compactness, protrusion, convexity, rigidity, hydrophobicity, and charge density; the weights are calculated from a set of previously identified instances of the binding-site type on known protein structures. The scoring function can easily incorporate different types of information useful in localization, thus increasing the applicability and accuracy of the approach. To test the method, 1008 known protein structures were split into 20 different groups according to the type of the bound ligand. For nonsugar ligands, such as various nucleotides, binding sites were correctly identified in 55%-73% of the cases. The method is completely automated (http://salilab.org/patcher) and can be applied on a large scale in a structural genomics setting.     

FTSite: high accuracy detection of ligand binding sites on unbound protein structures

MOTIVATION: Binding site identification is a classical problem that is important for a range of applications, including the structure-based prediction of function, the elucidation of functional relationships among proteins, protein engineering and drug design. We describe an accurate method of binding site identification, namely FTSite. This method is based on experimental evidence that ligand binding sites also bind small organic molecules of various shapes and polarity. The FTSite algorithm does not rely on any evolutionary or statistical information, but achieves near experimental accuracy: it is capable of identifying the binding sites in over 94% of apo proteins from established test sets that have been used to evaluate many other binding site prediction methods. AVAILABILITY: FTSite is freely available as a web-based server at http://ftsite.bu.edu.     

Recognition of functional sites in protein structures

Recognition of regions on the surface of one protein, that are similar to a binding site of another is crucial for the prediction of molecular interactions and for functional classifications. We first describe a novel method, SiteEngine, that assumes no sequence or fold similarities and is able to recognize proteins that have similar binding sites and may perform similar functions. We achieve high efficiency and speed by introducing a low-resolution surface representation via chemically important surface points, by hashing triangles of physico-chemical properties and by application of hierarchical scoring schemes for a thorough exploration of global and local similarities. We proceed to rigorously apply this method to functional site recognition in three possible ways: first, we search a given functional site on a large set of complete protein structures. Second, a potential functional site on a protein of interest is compared with known binding sites, to recognize similar features. Third, a complete protein structure is searched for the presence of an a priori unknown functional site, similar to known sites. Our method is robust and efficient enough to allow computationally demanding applications such as the first and the third. From the biological standpoint, the first application may identify secondary binding sites of drugs that may lead to side-effects. The third application finds new potential sites on the protein that may provide targets for drug design. Each of the three applications may aid in assigning a function and in classification of binding patterns. We highlight the advantages and disadvantages of each type of search, provide examples of large-scale searches of the entire Protein Data Base and make functional predictions.     

High-throughput identification of interacting protein-protein binding sites

Background  

With the advent of increasing sequence and structural data, a number of methods have been proposed to locate putative protein binding sites from protein surfaces. Therefore, methods that are able to identify whether these binding sites interact are needed.     

Searching for functional sites in protein structures

An ability to assign protein function from protein structure is important for structural genomics consortia. The complex relationship between protein fold and function highlights the necessity of looking beyond the global fold of a protein to specific functional sites. Many computational methods have been developed that address this issue. These include evolutionary trace methods, methods that involve the calculation and assessment of maximal superpositions, methods based on graph theory, and methods that apply machine learning techniques. Such function prediction techniques have been applied to the identification of enzyme catalytic triads and DNA-binding motifs.     

Geometric constraints for porphyrin binding in helical protein binding sites

Helical bundles which bind heme and porphyrin cofactors have been popular targets for cofactor-containing de novo protein design. By analyzing a highly nonredundant subset of the protein databank we have determined a rotamer distribution for helical histidines bound to heme cofactors. Analysis of the entire nonredundant database for helical sequence preferences near the ligand histidine demonstrated little preference for amino acid side chain identity, size, or charge. Analysis of the database subdivided by ligand histidine rotamer, however, reveals strong preferences in each case, and computational modeling illuminates the structural basis for some of these findings. The majority of the rotamer distribution matches that predicted by molecular simulation of a single porphyrin-bound histidine residue placed in the center of an all-alanine helix, and the deviations explain two prominent features of natural heme protein binding sites: heme distortion in the case of the cytochromes C in the m166 histidine rotamer, and a highly prevalent glycine residue in the t73 histidine rotamer. These preferences permit derivation of helical consensus sequence templates which predict optimal side chain-cofactor packing interactions for each rotamer. These findings thus promise to guide future design endeavors not only in the creation of higher affinity heme and porphyrin binding sites, but also in the direction of bound cofactor geometry.     

Quantitative identification of protein nitration sites

Several labelling strategies have been developed targeting specific amino acid residues and/or PTMs. Methods specifically tailored for the qualitative and sometimes quantitative determination of PTMs have emerged. Many research groups have focused their attention towards o‐nitrotyrosine residues, developing various methodologies for their identification, while direct quantification has remained elusive. So far the iTRAQ chemistry has been limited to primary amines. Here, we report a new strategy based on the use of iTRAQ reagents coupled to MS analysis for the selective labelling of o‐nitrotyrosine residues. This method was proved to lead to the simultaneous localisation and quantification of nitration sites both in model proteins and in biological systems.     

Similarity networks of protein binding sites

An increasing attention has been dedicated to the characterization of complex networks within the protein world. This work is reporting how we uncovered networked structures that reflected the structural similarities among protein binding sites. First, a 211 binding sites dataset has been compiled by removing the redundant proteins in the Protein Ligand Database (PLD) (http://www-mitchell.ch.cam.ac.uk/pld/). Using a clique detection algorithm we have performed all-against-all binding site comparisons among the 211 available ones. Within the set of nodes representing each binding site an edge was added whenever a pair of binding sites had a similarity higher than a threshold value. The generated similarity networks revealed that many nodes had few links and only few were highly connected, but due to the limited data available it was not possible to definitively prove a scale-free architecture. Within the same dataset, the binding site similarity networks were compared with the networks of sequence and fold similarity networks. In the protein world, indications were found that structure is better conserved than sequence, but on its own, sequence was better conserved than the subset of functional residues forming the binding site. Because a binding site is strongly linked with protein function, the identification of protein binding site similarity networks could accelerate the functional annotation of newly identified genes. In view of this we have discussed several potential applications of binding site similarity networks, such as the construction of novel binding site classification databases, as well as the implications for protein molecular design in general and computational chemogenomics in particular.     

Progress in the identification of interaction sites on the periplasmic maltose binding protein from E coli

The periplasmic maltose binding protein (MBP) is required for the high affinity transport of maltose and maltodextrins and for chemotaxis towards these sugars. In these functions, MBP interacts with proteins of the cytoplasmic membrane: MalF and MalG for transport, Tar for chemotaxis. A large number of MBP mutations have been isolated by us and other laboratories. We grouped these mutations into classes depending on the interactions affected and we represented the corresponding residues on the 3-D model for MBP so as to further identify the sites of MBP interacting with the MalF-MalG complex and with the Tar protein. MBP (like the other binding proteins) is composed of 2 lobes enclosing a cleft where the substrate binds. The face of the protein opposite the cleft seems to interact neither with MalF-MalG nor with Tar. The other face, corresponding to the cleft, contains sites for interactions with MalF-MalG and Tar. These sites appear to cover both sides of the cleft and may overlap in part. The present definition of the interaction sites suggests further that MBP has different in vivo orientations when it interacts with MalF-MalG or with Tar. This work constitutes an additional step in combining the use of genetic and structural analysis to define the interaction sites on MBP. Because of the structural similarities between periplasmic binding proteins, the regions of interaction defined could be relevant for other members of this family.     

LigBase: a database of families of aligned ligand binding sites in known protein sequences and structures

A database comprising all ligand-binding sites of known structure aligned with all related protein sequences and structures is described. Currently, the database contains approximately 50000 ligand-binding sites for small molecules found in the Protein Data Bank (PDB). The structure-structure alignments are obtained by the Combinatorial Extension (CE) program (Shindyalov and Bourne, Protein Eng., 11, 739-747, 1998) and sequence-structure alignments are extracted from the ModBase database of comparative protein structure models for all known protein sequences (Sanchez et al., Nucleic Acids Res., 28, 250-253, 2000). It is possible to search for binding sites in LigBase by a variety of criteria. LigBase reports summarize ligand data including relevant structural information from the PDB file, such as ligand type and size, and contain links to all related protein sequences in the TrEMBL database. Residues in the binding sites are graphically depicted for comparison with other structurally defined family members. LigBase provides a resource for the analysis of families of related binding sites.