Exploring the Composition of Protein-Ligand Binding Sites on a Large Scale

The residue composition of a ligand binding site determines the interactions available for diffusion-mediated ligand binding, and understanding general composition of these sites is of great importance if we are to gain insight into the functional diversity of the proteome. Many structure-based drug design methods utilize such heuristic information for improving prediction or characterization of ligand-binding sites in proteins of unknown function. The Binding MOAD database if one of the largest curated sets of protein-ligand complexes, and provides a source of diverse, high-quality data for establishing general trends of residue composition from currently available protein structures. We present an analysis of 3,295 non-redundant proteins with 9,114 non-redundant binding sites to identify residues over-represented in binding regions versus the rest of the protein surface. The Binding MOAD database delineates biologically-relevant “valid” ligands from “invalid” small-molecule ligands bound to the protein. Invalids are present in the crystallization medium and serve no known biological function. Contacts are found to differ between these classes of ligands, indicating that residue composition of biologically relevant binding sites is distinct not only from the rest of the protein surface, but also from surface regions capable of opportunistic binding of non-functional small molecules. To confirm these trends, we perform a rigorous analysis of the variation of residue propensity with respect to the size of the dataset and the content bias inherent in structure sets obtained from a large protein structure database. The optimal size of the dataset for establishing general trends of residue propensities, as well as strategies for assessing the significance of such trends, are suggested for future studies of binding-site composition.     

Recovering the true targets of specific ligands by virtual screening of the protein data bank

The Protein Data Bank (PDB) has been processed to extract a screening protein library (sc-PDB) of 2148 entries. A knowledge-based detection algorithm has been applied to 18,000 PDB files to find regular expressions corresponding to either protein, ions, co-factors, solvent, or ligand atoms. The sc-PDB database comprises high-resolution X-ray structures of proteins for which (i) a well-defined active site exists, (ii) the bound-ligand is a small molecular weight molecule. The database has been screened by an inverse docking tool derived from the GOLD program to recover the known target of four unrelated ligands. Both the database and the inverse screening procedures are accurate enough to rank the true target of the four investigated ligands among the top 1% scorers, with 70-100 fold enrichment with respect to random screening. Applying the proposed screening procedure to a small-sized generic ligand was much less accurate suggesting that inverse screening shall be reserved to rather selective compounds.     

Improving accuracy and efficiency of blind protein-ligand docking by focusing on predicted binding sites

The use of predicted binding sites (binding sites calculated from the protein structure alone) is evaluated here as a tool to focus the docking of small molecule ligands into protein structures, simulating cases where the real binding sites are unknown. The resulting approach consists of a few independent docking runs carried out on small boxes, centered on the predicted binding sites, as opposed to one larger blind docking run that covers the complete protein structure. The focused and blind approaches were compared using a set of 77 known protein-ligand complexes and 19 ligand-free structures. The focused approach is shown to: (1) identify the correct binding site more frequently than blind docking; (2) produce more accurate docking poses for the ligand; (3) require less computational time. Additionally, the results show that very few real binding sites are missed in spite of focusing on only three predicted binding sites per target protein. Overall the results indicate that, by improving the sampling in regions that are likely to correspond to binding sites, the focused docking approach increases accuracy and efficiency of protein ligand docking for those cases where the ligand-binding site is unknown. This is especially relevant in applications such as reverse virtual screening and structure-based functional annotation of proteins.     

Mapping of active site of alcohol dehydrogenase with low-molecular ligands

In search of an active alcohol dehydrogenase inhibitor, the structure of which may serve as the basis for a potential drug design, the active site of alcohol dehydrogenase containing NAD and Zn2+ ions was mapped using the method of molecular mechanics. Molecular docking was performed using a number of ligands containing characteristic functional groups: formate ion, ammonia, ammonium ion, methanol, and methylamine. Sites of preferable binding were revealed for each ligand and arranged in order of decreasing energy of binding to the enzyme. A comparison of the predicted ligand-binding sites and the experimental data on the location of water and inhibitor binding sites in the known structures of corresponding alcohol dehydrogenase complexes indicated a coincidence of the complex formation sites, which confirms the validity of the method and provides the requirements for a highly effective inhibitor (the pharmacophore model).     

Mapping of the active site of alcohol dehydrogenase with low-molecular ligands

In search of an active alcohol dehydrogenase inhibitor, the structure of which may serve as the basis for a potential drug design, the active site of alcohol dehydrogenase containing NAD and Zn²⁺ ions was mapped using the method of molecular mechanics. Molecular docking was performed using a number of ligands containing characteristic functional groups: formate ion, ammonia, ammonium ion, methanol, and methylamine. Sites of preferable binding were revealed for each ligand and arranged in order of decreasing energy of binding to the enzyme. A comparison of the predicted ligand-binding sites and the experimental data on the location of water and inhibitor binding sites in the known structures of corresponding alcohol dehydrogenase complexes indicated a coincidence of the complex formation sites, which confirms the validity of the method and provides the requirements for a highly effective inhibitor (the pharmacophore model).     

Structural stability of binding sites: Consequences for binding affinity and allosteric effects

During the course of biological function, proteins interact with other proteins, ligands, substrates, inhibitors, etc. These interactions occur at precisely defined locations within the protein but their effects are sometimes propagated to distal regions, triggering highly specific responses. These effects can be used as signals directed to activate or inhibit other sites, modulate interactions with other molecules, and/or establish inter‐molecular communication networks. During the past decade, it has become evident that the energy of stabilization of the protein structure is not evenly distributed throughout the molecule and that, under native conditions, proteins lack global cooperativity and are characterized by the occurrence of multiple independent local unfolding events. From a biological point of view, it is important to assess if this uneven distribution reflects specific functional requirements. For example, are binding sites more likely to be found in well structured regions, unstable regions, or mixed regions? In this article, we have addressed these questions by performing a structure‐based thermodynamic stability analysis of non‐structurally homologous proteins for which high resolution structures of their complexes with specific ligands are available. The results of these studies indicate that for all 16 proteins considered, the binding sites have a dual character and are characterized by the presence of regions with very low structural stability and regions with high stability. In many cases the low stability regions are loops that become stable and cover a significant portion of low molecular weight ligands upon binding. For enzymes, catalytic residues are usually, but not always, located in regions with high structural stability. It is shown that this arrangement provides significant advantages for the optimization of binding affinity of small ligands. In allosteric enzymes, low stability regions in the regulatory site are shown to play a crucial role in the transmission of information to the catalytic site. Proteins 2000;41:63–71. © 2000 Wiley‐Liss, Inc.     

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.     

Characterization and prediction of alternative splice sites

Human alternative isoform, cryptic, skipped, and constitutive splice sites from the ALTEXTRON database were analysed regarding splice site strength, composition, GC content, position and binding site strength of polypyrimidine tract and branch site. Several features were identified which distinguish alternative isoform and cryptic splice sites, but not skipped splice sites from constitutive ones. These include splice site strength, introns GC content, U2AF35 binding site score, and oligonucleotide frequencies. For the predictive classification of splice sites, pattern recognition models for different splicing factor binding sites and oligonucleotide frequency models (OFMs) were combined using backpropagation networks. 67.45% of acceptor sites and 71.23% of donor sites are correctly classified by networks trained for classification of constitutive and alternative isoform/cryptic splice sites. A web-application for the prediction of alternative splice sites is available at http://es.embnet.org/~mwang/assp.html .     

Structural characteristics of protein binding sites for calcium and lanthanide ions

Surveys of X-ray structures of Ca2+-containing and lanthanide ion-containing proteins and coordination complexes have been performed and structural features of the metal binding sites compared. A total of 515 structures of Ca2+-containing proteins were considered, although the final data set contained only 44 structures and 60 Ca2+ binding sites with a total of 323 ligands. Eighteen protein structures containing lanthanide ions were considered with a final data set containing eight structures and 11 metal binding sites. Structural features analysed include coordination numbers of the metal ions, the identity of their ligands, the denticity of carboxylate ligands, and the type of secondary structure from which the ligands are derived. Three general types of calcium binding site were identified in the final data set: class I sites supply the Ca2+ ligands from a continuous short sequence of amino acids; class II sites have one ligand supplied by a part of the amino acid sequence far removed from the main binding sequence; and class III sites are created by amino acids remote from one another in the sequence. The abundant EF-hand type of Ca2+ binding site was under-represented in the data set of structures analysed as far as its biological distribution is concerned, but was adequately represented for the chemical survey undertaken. A turn or loop structure was found to provide the bulk of the ligands to Ca2+, but helix and sheet secondary structures are slightly better providers of bidentate carboxylate ligation than turn or loop structures. The average coordination number for Ca2+ was 6.0, though for EF-hand sites it is 7. The average coordination number of a lanthanide ion in an intrinsic protein Ca2+ site was 7.2, but for the adventitious sites was only 4.4. A survey of the Cambridge Structural Database showed there are small-molecule lanthanide complexes with low coordination numbers but it is likely that water molecules, which do not appear in the electron density maps, are present for some lanthanide sites in proteins. A detailed comparison of the well-defined Ca2+ and lanthanide ion binding sites suggests that a reduction of hydrogen bonding associated with the ligating residues of the binding sites containing lanthanide ions may be a response to the additional positive charge of the lanthanide ion. Major structural differences between Ca2+ binding sites with weak and strong binding affinities were not obvious, a consequence of long-range electrostatic interactions and metal ion-induced protein conformational changes modulating affinities.     

Localization of protein-binding sites within families of proteins

We address the question of whether or not the positions of protein-binding sites on homologous protein structures are conserved irrespective of the identities of their binding partners. First, for each domain family in the Structural Classification of Proteins (SCOP), protein-binding sites are extracted from our comprehensive database of structurally defined binary domain interactions (PIBASE). Second, the binding sites within each family are superposed using a structural alignment of its members. Finally, the degree of localization of binding sites within each family is quantified by comparing it with localization expected by chance. We found that 72% of the 1847 SCOP domain families in PIBASE have binding sites with localization values greater than expected by chance. Moreover, 554 (30%) of these families have localizations that are statistically significant (i.e., more than four standard deviations away from the mean expected by chance). In contrast, only 144 (8%) families have significantly low localization. The absence of a significant correlation of the binding site localization with the average sequence and structural conservations in a family suggests that localization can be helpful for describing the functional diversity of protein-protein interactions, complementing measures of sequence and structural conservation. Consideration of the binding site localization may also result in spatial restraints for the modeling of protein assembly structures.     

Nonequivalence of alpha-bungarotoxin binding sites in the native nicotinic receptor molecule

In the native, membrane-bound form of the nicotinic acetylcholine receptor (M-AcChR) the two sites for the cholinergic antagonist alpha-bungarotoxin (alpha-BGT) have different binding properties. One site has high affinity, and the M-AcChR/alpha-BGT complexes thus formed dissociate very slowly, similar to the complexes formed with detergent-solubilized AcChR (S-AcChR). The second site has much lower affinity (KD approximately 59 +/- 35 nM) and forms quickly reversible complexes. The nondenaturing detergent Triton X-100 is known to solubilize the AcChR in a form unable, upon binding of cholinergic ligands, to open the ion channel and to become desensitized. Solubilization of the AcChR in Triton X-100 affects the binding properties of this second site and converts it to a high-affinity, slowly reversible site. Prolonged incubation of M-AcChR at 4 degrees C converts the low-affinity site to a high-affinity site similar to those observed in the presence of Triton X-100. Although the two sites have similar properties when the AcChR is solubilized in Triton X-100, their nonequivalence can be demonstrated by the effect on alpha-BGT binding of concanavalin A, which strongly reduces the association rate of one site only. The Bmax of alpha-BGT to either Triton-solubilized AcChR or M-AcChR is not affected by the presence of concanavalin A. Occupancy of the high-affinity, slowly reversible site in M-AcChR inhibits the Triton X-100 induced conversion to irreversibility of the second site. At difference with alpha-BGT, the long alpha-neurotoxin from Naja naja siamensis venom (alpha-NTX) binds with high affinity and in a very slowly reversible fashion to two sites in the M-AcChR (Conti-Tronconi & Raftery, 1986). We confirm here that Triton-solubilized AcChR or M-AcChR binds in a very slowly reversible fashion the same amount of alpha-NTX.(ABSTRACT TRUNCATED AT 250 WORDS)     

Prediction of functional sites by analysis of sequence and structure conservation

We present a method for prediction of functional sites in a set of aligned protein sequences. The method selects sites which are both well conserved and clustered together in space, as inferred from the 3D structures of proteins included in the alignment. We tested the method using 86 alignments from the NCBI CDD database, where the sites of experimentally determined ligand and/or macromolecular interactions are annotated. In agreement with earlier investigations, we found that functional site predictions are most successful when overall background sequence conservation is low, such that sites under evolutionary constraint become apparent. In addition, we found that averaging of conservation values across spatially clustered sites improves predictions under certain conditions: that is, when overall conservation is relatively high and when the site in question involves a large macromolecular binding interface. Under these conditions it is better to look for clusters of conserved sites than to look for particular conserved sites.     

Zinc coordination sphere in biochemical zinc sites

Zinc is known to be indispensable to growth and development and transmission of the genetic message. It does this through a remarkable mosaic of zinc binding motifs that orchestrate all aspects of metabolism. There are now nearly 200 three dimensional structures for zinc proteins, representing all six classes of enzymes and covering a wide range of phyla and species. These structures provide standards of reference for the identity and nature of zinc ligands in other proteins for which only the primary structure is known. Three primary types of zinc sites are apparent from examination of these structures: structural, catalytic and cocatalytic. The most common amino acids that supply ligands to these sites are His, Glu, Asp and Cys. In catalytic sites zinc generally forms complexes with water and any three nitrogen, oxygen and sulfur donors with His being the predominant amino acid chosen. Water is always a ligand to such sites. Structural zinc sites have four protein ligands and no bound water molecule. Cys is the preferred ligand in such sites. Cocatalytic sites contain two or three metals in close proximity with two of the metals bridged by a side chain moiety of a single amino acid residue, such as Asp, Glu or His and sometimes a water molecule. Asp and His are the preferred amino acids for these sites. No Cys ligands are found in such sites. The scaffolding of the zinc sites is also important to the function and reactivity of the bound metal. The influence of zinc on quaternary protein structure has led to the identification of a fourth type of zinc binding site, protein inteface. In this case zinc sites are formed from ligands supplied from amino acid residues residing in the binding surface of two proteins. The resulting zinc site usually has the coordination properties of a catalytic or structural zinc binding site.     

Docking molecules by families to increase the diversity of hits in database screens: computational strategy and experimental evaluation

Molecular docking programs screen chemical databases for novel ligands that fit protein binding sites. When one compound fits the site well, close analogs typically do the same. Therefore, many of the compounds that are found in such screens resemble one another. This reduces the variety and novelty of the compounds suggested. In an attempt to increase the diversity of docking hit lists, the Available Chemicals Directory was grouped into families of related structures. All members of every family were docked and scored, but only the best scoring molecule of a high-ranking family was allowed in the hit list. The identity and scores of the other members of these families were recorded as annotations to the best family member, but they were not independently ranked. This family-based docking method was compared with molecule-by-molecule docking in screens against the structures of thymidylate synthase, dihydrofolate reductase (DHFR), and the cavity site of the mutant T4 lysozyme Leu99 --> Ala (L99A). In each case, the diversity of the hit list increased, and more families of known ligands were found. To investigate whether the newly identified hits were sensible, we tested representative examples experimentally for binding to L99A and DHFR. Of the six compounds tested against L99A, five bound to the internal cavity. Of the seven compounds tested against DHFR, six inhibited the enzyme with apparent K(i) values between 0.26 and 100 microM. The segregation of potential ligands into families of related molecules is a simple technique to increase the diversity of candidates suggested by database screens. The general approach should be applicable to most docking methods. Proteins 2001;42:279-293.     

Site-site interaction on mitochondrial F1-ATPase. Functional symmetry of the high-affinity nucleotide binding sites

Interactions between the high affinity binding sites on mitochondrial F1 were analysed by combined use of the nucleotide analogues 3'-O-(1-naphthoyl)-ADP (N-ADP) and 2'-3'-O-(2,4,6-trinitrophenyl)-ADP (TNP-ADP). The binding behaviour of F1 with respect to these ligands was studied by measuring the fluorescence of F1 and of TNP-ADP and the fluorescence anisotropy of N-ADP. A total of 3 high affinity binding sites can be occupied by TNP-ADP. By exchange experiments, it could be shown that binding of TNP-ADP to such a site considerably accelerates the dissociation of a ligand bound to a neighbouring site. These results support the notion that the functional behaviour of F1 is symmetric: during the catalytic cycle any individual site can successively assume different affinity states as has been predicted by hypotheses such as the binding change model.