DynaBiS: A hierarchical sampling algorithm to identify flexible binding sites for large ligands and peptides

Knowing the ligand or peptide binding site in proteins is highly important to guide drug discovery, but experimental elucidation of the binding site is difficult. Therefore, various computational approaches have been developed to identify potential binding sites in protein structures. However, protein and ligand flexibility are often neglected in these methods due to efficiency considerations despite the recognition that protein–ligand interactions can be strongly affected by mutual structural adaptations. This is particularly true if the binding site is unknown, as the screening will typically be performed based on an unbound protein structure. Herein we present DynaBiS, a hierarchical sampling algorithm to identify flexible binding sites for a target ligand with explicit consideration of protein and ligand flexibility, inspired by our previously presented flexible docking algorithm DynaDock. DynaBiS applies soft-core potentials between the ligand and the protein, thereby allowing a certain protein–ligand overlap resulting in efficient sampling of conformational adaptation effects. We evaluated DynaBiS and other commonly used binding site identification algorithms against a diverse evaluation set consisting of 26 proteins featuring peptide as well as small ligand binding sites. We show that DynaBiS outperforms the other evaluated methods for the identification of protein binding sites for large and highly flexible ligands such as peptides, both with a holo or apo structure used as input.     

Role of structural water for prediction of cation binding sites in apoproteins

Structures of many metal-binding proteins are often obtained without structural cations in their apoprotein forms. Missing cation coordinates are usually updated from structural templates constructed from many holoprotein structures. Such templates usually do not include structural water, the important contributor to the ion binding energy. Structural templates are also inconvenient for taking into account structural modifications around the binding site at apo-/holo- transitions. An approach based upon statistical potentials readily takes into account structural modifications associated with binding as well as contribution of structural water molecules. Here, we construct a set of statistical potentials for Mg²⁺, Ca²⁺, and Zn²⁺ contacting with protein atoms of a different type or structural water oxygens. Each type of the cations tends to form tight contacts with protein atoms of specific types. Structural water contributes relatively more into the binding pseudo-energy of Mg²⁺ and Ca²⁺ than of Zn²⁺. We have developed PIONCA (Protein-Ion Calculator), a fast CUDA GPGPU-based algorithm that predicts ion-binding sites in apoproteins. Comparative tests demonstrate that PIONCA outperforms most of the tools based on structural templates or docking. Our software can be also used for locating bound cations in holoprotein structures with missing cation heteroatoms. PIONCA is equipped with an interactive web interface based upon JSmol.     

Sequence and structural determinants of mannose recognition

Mannose, an abundant cell surface monosaccharide binds to a diverse set of receptors, which are involved in a variety of important cellular processes. Structural analysis has been carried out on all the proteins containing non-covalently bound mannose as a monosaccharide in the Protein Data Bank, to identify common recognition principles. Proteins, highly specific to mannose, belonging to the super family of bulb lectins, are found to contain a consensus sequence motif QXDXNXVXY, which has been identified to be essential for mannose binding. Analysis of this motif in the crystal structures of bulb lectins has led to the understanding of the contribution of individual residues in mannose recognition. Comparison with other mannose binding proteins, reveals common hydrogen bonding patterns in all of them, despite differences in sequence, overall fold and the substructures at the binding sites of individual proteins. A database analysis also suggests that although the topology of the backbone, as at the binding site in bulb lectins, can generate mannose binding capability in a few other proteins, sequence and disposition of not only the residues in the motif, but also the residues in the neighborhood play a crucial role in allowing that property to be retained.     

Deciphering common recognition principles of nucleoside mono/di and tri‐phosphates binding in diverse proteins via structural matching of their binding sites

Nucleoside triphosphate (NTP) ligands are of high biological importance and are essential for all life forms. A pre‐requisite for them to participate in diverse biochemical processes is their recognition by diverse proteins. It is thus of great interest to understand the basis for such recognition in different proteins. Towards this, we have used a structural bioinformatics approach and analyze structures of 4677 NTP complexes available in Protein Data Bank (PDB). Binding sites were extracted and compared exhaustively using PocketMatch, a sensitive in‐house site comparison algorithm, which resulted in grouping the entire dataset into 27 site‐types. Each of these site‐types represent a structural motif comprised of two or more residue conservations, derived using another in‐house tool for superposing binding sites, PocketAlign. The 27 site‐types could be grouped further into 9 super‐types by considering partial similarities in the sites, which indicated that the individual site‐types comprise different combinations of one or more site features. A scan across PDB using the 27 structural motifs determined the motifs to be specific to NTP binding sites, and a computational alanine mutagenesis indicated that residues identified to be highly conserved in the motifs are also most contributing to binding. Alternate orientations of the ligand in several site‐types were observed and rationalized, indicating the possibility of some residues serving as anchors for NTP recognition. The presence of multiple site‐types and the grouping of multiple folds into each site‐type is strongly suggestive of convergent evolution. Knowledge of determinants obtained from this study will be useful for detecting function in unknown proteins. Proteins 2017; 85:1699–1712. © 2017 Wiley Periodicals, Inc.     

Exploring functional roles of multibinding protein interfaces

Cellular processes are highly interconnected and many proteins are shared in different pathways. Some of these shared proteins or protein families may interact with diverse partners using the same interface regions; such multibinding proteins are the subject of our study. The main goal of our study is to attempt to decipher the mechanisms of specific molecular recognition of multiple diverse partners by promiscuous protein regions. To address this, we attempt to analyze the physicochemical properties of multibinding interfaces and highlight the major mechanisms of functional switches realized through multibinding. We find that only 5% of protein families in the structure database have multibinding interfaces, and multibinding interfaces do not show any higher sequence conservation compared with the background interface sites. We highlight several important functional mechanisms utilized by multibinding families. (a) Overlap between different functional pathways can be prevented by the switches involving nearby residues of the same interfacial region. (b) Interfaces can be reused in pathways where the substrate should be passed from one protein to another sequentially. (c) The same protein family can develop different specificities toward different binding partners reusing the same interface; and finally, (d) inhibitors can attach to substrate binding sites as substrate mimicry and thereby prevent substrate binding.     

Evolution of binding sites for zinc and calcium ions playing structural roles

The geometry of metal coordination by proteins is well understood, but the evolution of metal binding sites has been less studied. Here we present a study on a small number of well-documented structural calcium and zinc binding sites, concerning how the geometry diverges between relatives, how often nonrelatives converge towards the same structure, and how often these metal binding sites are lost in the course of evolution. Both calcium and zinc binding site structure is observed to be conserved; structural differences between those atoms directly involved in metal binding in related proteins are typically less than 0.5 A root mean square deviation, even in distant relatives. Structural templates representing these conserved calcium and zinc binding sites were used to search the Protein Data Bank for cases where unrelated proteins have converged upon the same residue selection and geometry for metal binding. This allowed us to identify six "archetypal" metal binding site structures: two archetypal zinc binding sites, both of which had independently evolved on a large number of occasions, and four diverse archetypal calcium binding sites, where each had evolved independently on only a handful of occasions. We found that it was common for distant relatives of metal-binding proteins to lack metal-binding capacity. This occurred for 13 of the 18 metal binding sites we studied, even though in some of these cases the original metal had been classified as "essential for protein folding." For most of the calcium binding sites studied (seven out of eleven cases), the lack of metal binding in relatives was due to point mutation of the metal-binding residues, whilst for zinc binding sites, lack of metal binding in relatives always involved more extensive changes, with loss of secondary structural elements or loops around the binding site.     

Concentration-dependent inactivation of superoxide dismutase

1. Vasoactive intestinal peptide (VIP) receptors were identified in crude rat hepatic membranes by 125I-labelled VIP binding and by the ability of VIP to stimulate adenylate cyclase activity. The specificity of these receptors was evaluated by the capacity of secretin, synthetic secretin analogues, and secretin fragments to inhibit 125I-labelled VIP binding and to stimulate adenylate cyclase. 2. The results were compatible with the existence of two classes of VIP binding sites that could be distinguished according to their affinity for VIP and their specificity. High-affinity sites were more specific for VIP as secretin was 175 times less potent than VIP for recognition of these sites while being only 33 times less potent than VIP for recognition of low-affinity sites. 3. Secretin analogues, monosubstituted in position 2, 3, 4 or 6 were less potent than secretin for adenylate cyclase stimulation as well as for the recognition of the two classes of receptors. [Val5]secretin was more potent than secretin and appeared definitely more VIP-like than secretin; [Ala4, Val5] and [D-Ala4,Val5]secretin were equipotent to secretin. 4. The fragment secretin (7-27) was unable to recognize VIP receptors and to stimulate adenylate cyclase. The substituted fragment [Gln9,Asn15]secretin (5-27) recognized these receptors with weak potency but could not activate the enzyme.     

Lipid-binding surfaces of membrane proteins: evidence from evolutionary and structural analysis

Membrane proteins function in the diverse environment of the lipid bilayer. Experimental evidence suggests that some lipid molecules bind tightly to specific sites on the membrane protein surface. These lipid molecules often act as co-factors and play important functional roles. In this study, we have assessed the evolutionary selection pressure experienced at lipid-binding sites in a set of α-helical and β-barrel membrane proteins using posterior probability analysis of the ratio of synonymous vs. nonsynonymous substitutions (ω-ratio). We have also carried out a geometric analysis of the membrane protein structures to identify residues in close contact with co-crystallized lipids. We found that residues forming cholesterol-binding sites in both β(2)-adrenergic receptor and Na(+)-K(+)-ATPase exhibit strong conservation, which can be characterized by an expanded cholesterol consensus motif for GPCRs. Our results suggest the functional importance of aromatic stacking interactions and interhelical hydrogen bonds in facilitating protein-cholesterol interactions, which is now reflected in the expanded motif. We also find that residues forming the cardiolipin-binding site in formate dehydrogenase-N γ-subunit and the phosphatidylglycerol binding site in KcsA are under strong purifying selection pressure. Although the lipopolysaccharide (LPS)-binding site in ferric hydroxamate uptake receptor (FhuA) is only weakly conserved, we show using a statistical mechanical model that LPS binds to the least stable FhuA β-strand and protects it from the bulk lipid. Our results suggest that specific lipid binding may be a general mechanism employed by β-barrel membrane proteins to stabilize weakly stable regions. Overall, we find that the residues forming specific lipid binding sites on the surfaces of membrane proteins often experience strong purifying selection pressure.     

Binding of [3H]forskolin to solubilized preparations of adenylate cyclase

The binding of [3H]forskolin to proteins solubilized from bovine brain membranes was studied by precipitating proteins with polyethylene glycol and separating [3H]forskolin bound to protein from free [3H]forskolin by rapid filtration. The Kd for [3H]forskolin binding to solubilized proteins was 14 nM which was similar to that for [3H]forskolin binding sites in membranes from rat brain and human platelets. Forskolin analogs competed for [3H]forskolin binding sites with the same rank potency in both brain membranes and in proteins solubilized from brain membranes. [3H]forskolin bound to proteins solubilized from membranes with a Bmax of 38 fmol/mg protein which increased to 94 fmol/mg protein when GppNHp was included in the binding assay. In contrast, GppNHp had no effect on [3H]forskolin binding to proteins solubilized from membranes preactivated with GppNHp. Solubilized adenylate cyclase from non-preactivated membranes had a basal activity of 130 pmol/mg/min which was increased 7-fold by GppNHp. In contrast, adenylate cyclase from preactivated membranes had a basal activity of 850 pmol/mg/min which was not stimulated by GppNHp or forskolin. Thus, the number of high affinity binding sites for [3H]forskolin in solubilized preparations correlated with the activation of adenylate cyclase by GppNHp via the guanine nucleotide binding protein (GS).     

Characterization of Imido [8-3H] Guanosine 5′-Triphosphate Binding Sites to Rat Brain Membranes

Besides their well-defined intracellular roles in transmembrane signals transduction, guanine derivatives play important roles by acting from the outside of neural cell membranes. These roles are mediated by two different pool sites in cell membranes: G proteins, which bind to specific (GDP and GTP) intracellular guanine derivatives, and sites that bind to extracellular guanine derivatives. In this study we investigated some methodological characteristics of both guanine derivatives binding sites (intracellular and extracellular) in rat brain neural membranes. By investigating the binding of a poorly hydrolyzed GTP analogue and the adenylate cyclase activity in neural membranes, we observed some distinctiveness of guanine derivatives binding sites: stability to washing procedures (extracellular) and modulation of adenylate cyclase activity (intracellular). These results allow dealing with each site separately, which could be useful for discriminating the roles of extracellular and intracellular guanine derivatives in the central nervous system.     

Hydrophilicity Matching – A Potential Prerequisite for the Formation of Protein-Protein Complexes in the Cell

A binding event between two proteins typically consists of a diffusional search of binding partners for one another, followed by a specific recognition of the compatible binding sites resulting in the formation of the complex. However, it is unclear how binding partners find each other in the context of the crowded, constantly fluctuating, and interaction-rich cellular environment. Here we examine the non-specific component of protein-protein interactions, which refers to those physicochemical properties of the binding partners that are independent of the exact details of their binding sites, but which can affect their localization or diffusional search for one another. We show that, for a large set of high-resolution experimental 3D structures of binary, transient protein complexes taken from the DOCKGROUND database, the binding partners display a surprising, statistically significant similarity in terms of their total hydration free energies normalized by a size-dependent variable. We hypothesize that colocalization of binding partners, even within individual cellular compartments such as the cytoplasm, may be influenced by their relative hydrophilicity, potentially in response to local hydrophilic gradients.