Ligand binding by antibody IgE Lb4: assessment of binding site preferences using microcalorimetry, docking, and free energy simulations.

Antibody IgE Lb4 interacts favorably with a large number of different compounds. To improve the current understanding of the structural basis of this vast cross-reactivity, the binding of three dinitrophenyl (DNP) amino acids (DNP-alanine, DNP-glycine, and DNP-serine) is investigated in detail by means of docking and molecular dynamics free energy simulations. Experimental binding energies obtained by isothermal titration microcalorimetry are used to judge the results of the computational studies. For all three ligands, the docking procedure proposes two plausible subsites within the binding region formed by the antibody CDR loops. By subsequent molecular dynamics simulations and calculations of relative free energies of binding, one of these subsites, a tyrosine-surrounded pocket, is revealed as the preferred point of complexation. For this subsite, results consistent with experimental observations are obtained; DNP-glycine is found to bind better than DNP-serine, and this, in turn, is found to bind better than DNP-alanine. The suggested binding mode makes it possible to explain both the moderate binding affinity and the differences in binding energy among the three ligands.     

Analysis of binding site similarity, small-molecule similarity and experimental binding profiles in the human cytosolic sulfotransferase family

MOTIVATION: In the present work we combine computational analysis and experimental data to explore the extent to which binding site similarities between members of the human cytosolic sulfotransferase family correlate with small-molecule binding profiles. Conversely, from a small-molecule point of view, we explore the extent to which structural similarities between small molecules correlate to protein binding profiles. RESULTS: The comparison of binding site structural similarities and small-molecule binding profiles shows that proteins with similar small-molecule binding profiles tend to have a higher degree of binding site similarity but the latter is not sufficient to predict small-molecule binding patterns, highlighting the difficulty of predicting small-molecule binding patterns from sequence or structure. Likewise, from a small-molecule perspective, small molecules with similar protein binding profiles tend to be topologically similar but topological similarity is not sufficient to predict their protein binding patterns. These observations have important consequences for function prediction and drug design.     

A heme binding site on myelin basic protein: characterization, location, and significance

Myelin basic protein (MBP), an extrinsic membrane protein from the myelin sheath, binds dicyanohemin. The binding generates absorption bands in the Soret region and quenches the fluorescence emitted by the sole tryptophan residue. The absorption titration curves in the Soret demonstrate that the binding is stoichiometric, one heme per protein, with a large value of the extinction coefficient (8 X 10(4) M-1 cm-1 at 420 nm). Fluorescence quenching titration curves indicate an identical stoichiometry and a low quenching efficiency of 20%. From the heme titration curve the association constant between dicyanohemin and MBP is estimated to be greater than or equal to 10 nM-1 in 50 mM 4-morpholinepropanesulfonic acid buffer, pH 7.0, at 20 degrees C. Digestion of MBP by Staphylococcus aureus V8 protease yields a peptide (38-118) whose heme binding properties are identical to those of MBP. In contrast, peptides obtained by digestion of MBP with cathepsin D do not exhibit any specific binding of dicyanohemin. The cleavage of the Phe-Phe (42-43) bond appears to be critical in this respect. A comparison of the sequence immediately preceding, including these residues with a probable heme binding site of a mitochondrial cytochrome b, reveals a high degree of homology. The possible significance of heme binding is discussed.     

Probing the location of the substrate binding site of ascorbate oxidase near type 1 copper: an investigation through spectroscopic, inhibition and docking studies

The present investigation addresses the problem of the binding mode of phenolic inhibitors and the substrate ascorbate to the active site of ascorbate oxidase. The results from both types of compounds indicate that the binding site is located in a pocket near the type 1 copper center. This information is of general interest for blue multicopper oxidases. Docking calculations performed on the ascorbate oxidase-ascorbate complex show that binding of the substrate occurs in a pocket near type 1 Cu, and is stabilized by at least five hydrogen bonding interactions with protein residues, one of which involves the His512 Cu ligand. Similar docking studies show that the isomeric fluorophenols, which act as competitive inhibitors toward ascorbate, bind to the enzyme in a manner similar to ascorbate. The docking calculations are supported by 19F NMR relaxation measurements performed on fluorophenols in the presence of the enzyme, which show that the bound inhibitors undergo enhanced relaxation by the paramagnetic effect of a nearby Cu center. Unambiguous support to the location of the inhibitor close to type 1 Cu was obtained by comparative relaxation measurements of the fluorophenols in the presence of the ascorbate oxidase derivative where a Zn atom selectively replaces the paramagnetic type 2 Cu. The latter experiments show that contribution to relaxation of the bound inhibitors by the type 2 Cu site is negligible.     

Probing molecular docking in a charged model binding site

A model binding site was used to investigate charge-charge interactions in molecular docking. This simple site, a small (180A(3)) engineered cavity in cyctochrome c peroxidase (CCP), is negatively charged and completely buried from solvent, allowing us to explore the balance between electrostatic energy and ligand desolvation energy in a system where many of the common approximations in docking do not apply. A database with about 5300 molecules was docked into this cavity. Retrospective testing with known ligands and decoys showed that overall the balance between electrostatic interaction and desolvation energy was captured. More interesting were prospective docking scre"ens that looked for novel ligands, especially those that might reveal problems with the docking and energy methods. Based on screens of the 5300 compound database, both high-scoring and low-scoring molecules were acquired and tested for binding. Out of 16 new, high-scoring compounds tested, 15 were observed to bind. All of these were small heterocyclic cations. Binding constants were measured for a few of these, they ranged between 20microM and 60microM. Crystal structures were determined for ten of these ligands in complex with the protein. The observed ligand geometry corresponded closely to that predicted by docking. Several low-scoring alkyl amino cations were also tested and found to bind. The low docking score of these molecules owed to the relatively high charge density of the charged amino group and the corresponding high desolvation penalty. When the complex structures of those ligands were determined, a bound water molecule was observed interacting with the amino group and a backbone carbonyl group of the cavity. This water molecule mitigates the desolvation penalty and improves the interaction energy relative to that of the "naked" site used in the docking screen. Finally, six low-scoring neutral molecules were also tested, with a view to looking for false negative predictions. Whereas most of these did not bind, two did (phenol and 3-fluorocatechol). Crystal structures for these two ligands in complex with the cavity site suggest reasons for their binding. That these neutral molecules do, in fact bind, contradicts previous results in this site and, along with the alkyl amines, provides instructive false negatives that help identify weaknesses in our scoring functions. Several improvements of these are considered.     

Using protein binding site prediction to improve protein docking

Predicting protein interaction interfaces and protein complexes are two important related problems. For interface prediction, there are a number of tools, such as PPI-Pred, PPISP, PINUP, Promate, and SPPIDER, which predict enzyme-inhibitor interfaces with success rates of 23% to 55% and other interfaces with 10% to 28% on a benchmark dataset of 62 complexes. Here, we develop, metaPPI, a meta server for interface prediction. It significantly improves prediction success rates to 70% for enzyme-inhibitor and 44% for other interfaces. As shown with Promate, predicted interfaces can be used to improve protein docking. Here, we follow this idea using the meta server instead of individual predictions. We confirm that filtering with predicted interfaces significantly improves candidate generation in rigid-body docking based on shape complementarity. Finally, we show that the initial ranking of candidate solutions in rigid-body docking can be further improved for the class of enzyme-inhibitor complexes by a geometrical scoring which rewards deep pockets. A web server of metaPPI is available at scoppi.tu-dresden.de/metappi. The source code of our docking algorithm BDOCK is also available at www.biotec.tu-dresden.de/~bhuang/bdock.     

Protein function annotation by local binding site surface similarity

Hundreds of protein crystal structures exist for proteins whose function cannot be confidently determined from sequence similarity. Surflex‐PSIM, a previously reported surface‐based protein similarity algorithm, provides an alternative method for hypothesizing function for such proteins. The method now supports fully automatic binding site detection and is fast enough to screen comprehensive databases of protein binding sites. The binding site detection methodology was validated on apo/holo cognate protein pairs, correctly identifying 91% of ligand binding sites in holo structures and 88% in apo structures where corresponding sites existed. For correctly detected apo binding sites, the cognate holo site was the most similar binding site 87% of the time. PSIM was used to screen a set of proteins that had poorly characterized functions at the time of crystallization, but were later biochemically annotated. Using a fully automated protocol, this set of 8 proteins was screened against ～60,000 ligand binding sites from the PDB. PSIM correctly identified functional matches that predated query protein biochemical annotation for five out of the eight query proteins. A panel of 12 currently unannotated proteins was also screened, resulting in a large number of statistically significant binding site matches, some of which suggest likely functions for the podorly characterized proteins. Proteins 2014; 82:679–694. © 2013 Wiley Periodicals, Inc.     

Testing a flexible-receptor docking algorithm in a model binding site

Sampling receptor flexibility is challenging for database docking. We consider a method that treats multiple flexible regions of the binding site independently, recombining them to generate different discrete conformations. This algorithm scales linearly rather than exponentially with the receptor's degrees of freedom. The method was first evaluated for its ability to identify known ligands of a hydrophobic cavity mutant of T4 lysozyme (L99A). Some 200000 molecules of the Available Chemical Directory (ACD) were docked against an ensemble of cavity conformations. Surprisingly, the enrichment of known ligands from among a much larger number of decoys in the ACD was worse than simply docking to the apo conformation alone. Large decoys, accommodated in the larger cavity conformations sampled in the ensemble, were ranked better than known small ligands. The calculation was redone with an energy correction term that considered the cost of forming the larger cavity conformations. Enrichment improved, as did the balance between high-ranking large and small ligands. In a second retrospective test, the ACD was docked against a conformational ensemble of thymidylate synthase. Compared to docking against individual enzyme conformations, the flexible receptor docking approach improved enrichment of known ligands. Including a receptor conformational energy weighting term improved enrichment further. To test the method prospectively, the ACD database was docked against another cavity mutant of lysozyme (L99A/M102Q). A total of 18 new compounds predicted to bind this polar cavity and to change its conformation were tested experimentally; 14 were found to bind. The bound structures for seven ligands were determined by X-ray crystallography. The predicted geometries of these ligands all corresponded to the observed geometries to within 0.7A RMSD or better. Significant conformational changes of the cavity were observed in all seven complexes. In five structures, part of the observed accommodations were correctly predicted; in two structures, the receptor conformational changes were unanticipated and thus never sampled. These results suggest that although sampling receptor flexibility can lead to novel ligands that would have been missed when docking a rigid structure, it is also important to consider receptor conformational energy.     

BSSF: a fingerprint based ultrafast binding site similarity search and function analysis server

Background  

Genome sequencing and post-genomics projects such as structural genomics are extending the frontier of the study of sequence-structure-function relationship of genes and their products. Although many sequence/structure-based methods have been devised with the aim of deciphering this delicate relationship, there still remain large gaps in this fundamental problem, which continuously drives researchers to develop novel methods to extract relevant information from sequences and structures and to infer the functions of newly identified genes by genomics technology.     

Discovery of new quinolines as potent colchicine binding site inhibitors: design,synthesis, docking studies,and anti-proliferative evaluation

Discovering of new anticancer agents with potential activity against tubulin polymerisation is still a promising approach. Colchicine binding site inhibitors are the most relevant anti-tubulin polymerisation agents. Thus, new quinoline derivatives have been designed and synthesised to possess the same essential pharmacophoric features of colchicine binding site inhibitors. The synthesised compounds were tested in vitro against a panel of three human cancer cell lines (HepG-2, HCT-116, and MCF-7) using colchicine as a positive control. Comparing to colchicine (IC₅₀ = 7.40, 9.32, and 10.41 µM against HepG-2, HCT-116, and MCF-7, respectively), compounds 20, 21, 22, 23, 24, 25, 26, and 28 exhibited superior cytotoxic activities with IC₅₀ values ranging from 1.78 to 9.19 µM. In order to sightsee the proposed mechanism of anti-proliferative activity, the most active members were further evaluated in vitro for their inhibitory activities against tubulin polymerisation. Compounds 21 and 32 exhibited the highest tubulin polymerisation inhibitory effect with IC₅₀ values of 9.11 and 10.5 nM, respectively. Such members showed activities higher than that of colchicine (IC₅₀ = 10.6 nM) and CA-4 (IC₅₀ = 13.2 nM). The impact of the most promising compound 25 on cell cycle distribution was assessed. The results revealed that compound 25 can arrest the cell cycle at G2/M phase. Annexin V and PI double staining assay was carried out to explore the apoptotic effect of the synthesised compounds. Compound 25 induced apoptotic effect on HepG-2 thirteen times more than the control cells. To examine the binding pattern of the target compounds against the tubulin heterodimers active site, molecular docking studies were carried out.     

Resolving protein structure‐function‐binding site relationships from a binding site similarity network perspective

Functional annotation is seldom straightforward with complexities arising due to functional divergence in protein families or functional convergence between non‐homologous protein families, leading to mis‐annotations. An enzyme may contain multiple domains and not all domains may be involved in a given function, adding to the complexity in function annotation. To address this, we use binding site information from bound cognate ligands and catalytic residues, since it can help in resolving fold‐function relationships at a finer level and with higher confidence. A comprehensive database of 2,020 fold‐function‐binding site relationships has been systematically generated. A network‐based approach is employed to capture the complexity in these relationships, from which different types of associations are deciphered, that identify versatile protein folds performing diverse functions, same function associated with multiple folds and one‐to‐one relationships. Binding site similarity networks integrated with fold, function, and ligand similarity information are generated to understand the depth of these relationships. Apart from the observed continuity in the functional site space, network properties of these revealed versatile families with topologically different or dissimilar binding sites and structural families that perform very similar functions. As a case study, subtle changes in the active site of a set of evolutionarily related superfamilies are studied using these networks. Tracing of such similarities in evolutionarily related proteins provide clues into the transition and evolution of protein functions. Insights from this study will be helpful in accurate and reliable functional annotations of uncharacterized proteins, poly‐pharmacology, and designing enzymes with new functional capabilities. Proteins 2017; 85:1319–1335. © 2017 Wiley Periodicals, Inc.     

A model binding site for testing scoring functions in molecular docking

Prediction of interaction energies between ligands and their receptors remains a major challenge for structure-based inhibitor discovery. Much effort has been devoted to developing scoring schemes that can successfully rank the affinities of a diverse set of possible ligands to a binding site for which the structure is known. To test these scoring functions, well-characterized experimental systems can be very useful. Here, mutation-created binding sites in T4 lysozyme were used to investigate how the quality of atomic charges and solvation energies affects molecular docking. Atomic charges and solvation energies were calculated for 172,118 molecules in the Available Chemicals Directory using a semi-empirical quantum mechanical approach by the program AMSOL. The database was first screened against the apolar cavity site created by the mutation Leu99Ala (L99A). Compared to the electronegativity-based charges that are widely used, the new charges and desolvation energies improved ranking of known apolar ligands, and better distinguished them from more polar isosteres that are not observed to bind. To investigate whether the new charges had predictive value, the non-polar residue Met102, which forms part of the binding site, was changed to the polar residue glutamine. The structure of the resulting Leu99Ala and Met102Gln double mutant of T4 lysozyme (L99A/M102Q) was determined and the docking calculation was repeated for the new site. Seven representative polar molecules that preferentially docked to the polar versus the apolar binding site were tested experimentally. All seven bind to the polar cavity (L99A/M102Q) but do not detectably bind to the apolar cavity (L99A). Five ligand-bound structures of L99A/M102Q were determined by X-ray crystallography. Docking predictions corresponded to the crystallographic results to within 0.4A RMSD. Improved treatment of partial atomic charges and desolvation energies in database docking appears feasible and leads to better distinction of true ligands. Simple model binding sites, such as L99A and its more polar variants, may find broad use in the development and testing of docking algorithms.     

Scoring optimisation of unbound protein-protein docking including protein binding site predictions

The prediction of the structure of the protein-protein complex is of great importance to better understand molecular recognition processes. During systematic protein-protein docking, the surface of a protein molecule is scanned for putative binding sites of a partner protein. The possibility to include external data based on either experiments or bioinformatic predictions on putative binding sites during docking has been systematically explored. The external data were included during docking with a coarse-grained protein model and on the basis of force field weights to bias the docking search towards a predicted or known binding region. The approach was tested on a large set of protein partners in unbound conformations. The significant improvement of the docking performance was found if reliable data on the native binding sites were available. This was possible even if data for single key amino acids at a binding interface are included. In case of binding site predictions with limited accuracy, only modest improvement compared with unbiased docking was found. The optimisation of the protocol to bias the search towards predicted binding sites was found to further improve the docking performance resulting in approximately 40% acceptable solutions within the top 10 docking predictions compared with 22% in case of unbiased docking of unbound protein structures.     

Substrate binding site of naphthalene 1,2-dioxygenase: functional implications of indole binding

The three-dimensional structure of the aromatic hydroxylating enzyme naphthalene dioxygenase (NDO) from Pseudomonas sp. NCIB 9816-4 was recently determined. The refinement of the structure together with cyclic averaging showed that in the active site of the enzyme there is electron density for a flat aromatic compound. This compound appears to be an indole adduct, which in Escherichia coli is derived from tryptophan present in the rich culture medium. An indole-dioxygen adduct has been built which fits the electron density convincingly. Support for this interpretation was obtained from crystals of the enzyme purified from cells grown in the absence of tryptophan which had an empty substrate pocket. These types of crystals were soaked in indole solutions and the position of indole in this complex was similar to the corresponding part in the modelled indole-oxygen adduct. This suggests that a peroxide bound to iron end-on attacks the substrate and forms this intermediate. The substrate position has implications for the substrate specificity of the enzyme. Docking studies with indole, naphthalene and biphenyl inside the substrate pocket of NDO suggest the presence of subpockets where the one close to the active site iron is reserved for the binding of the aromatic ring which is hydroxylated upon catalysis. The plausible location for the binding of dioxygen is between this pocket and the catalytic iron. This is in accordance with the enantiospecificity of the products.     

Ligand binding sites in Escherichia coli inorganic pyrophosphatase: effects of active site mutations

Type I soluble inorganic pyrophosphatases (PPases) are well characterized both structurally and mechanistically. Earlier we measured the effects of active site substitutions on pH--rate profiles for the type I PPases from both Escherichia coli (E-PPase) and Saccharomyces cerevisae (Y-PPase). Here we extend these studies by measuring the effects of such substitutions on the more discrete steps of ligand binding to E-PPase, including (a) Mg(2+) and Mn(2+) binding in the absence of added ligand; (b) Mg(2+) binding in the presence of either P(i) or hydroxymethylbisphosphonate (HMBP), a competitive inhibitor of E-PPase; and (c) P(i) binding in the presence of Mn(2+). The active site of a type I PPase has well-defined subsites for the binding of four divalent metal ions (M1--M4) and two phosphates (P1, P2). Our results, considered in light of pertinent results from crystallographic studies on both E-PPase and Y-PPase and parallel functional studies on Y-PPase, allow us to conclude the following: (a) residues E20, D65, D70, and K142 play key roles in the functional organization of the active site; (b) the major structural differences between the product and substrate complexes of E-PPase are concentrated in the lower half of the active site; (c) the M1 subsite is functionally isolated from the rest of the active site; and (d) the M4 subsite is an especially unconstrained part of the active site.     

Molecular dynamics simulation of the P2Y14 receptor. Ligand docking and identification of a putative binding site of the distal hexose moiety

A rhodopsin-based homology model of the P2Y14 receptor was inserted into a phospholipid bilayer and refined by molecular dynamics (MD) simulation. The binding modes of several known agonists, namely UDP-glucose and its analogues, were proposed using automatic molecular docking combined with Monte Carlo Multiple Minimum calculations. Compared to other P2Y receptors, the P2Y14 receptor has an atypical binding mode of the nucleobase, ribose, and phosphate moieties. The diphosphate moiety interacts with only one cationic residue, namely Lys171 of EL2, while in other P2Y receptor subtypes three Arg or Lys residues interact with the phosphate chain. Two other conserved cationic residues, namely Arg253 (6.55) and Lys277 (7.35) of the P2Y14 receptor together with two anionic residues (Glu166 and Glu174, located in EL2), are likely involved in interactions with the distal hexose moiety.     

The Poisson Index: a new probabilistic model for protein ligand binding site similarity

MOTIVATION: The large-scale comparison of protein-ligand binding sites is problematic, in that measures of structural similarity are difficult to quantify and are not easily understood in terms of statistical similarity that can ultimately be related to structure and function. We present a binding site matching score the Poisson Index (PI) based upon a well-defined statistical model. PI requires only the number of matching atoms between two sites and the size of the two sites-the same information used by the Tanimoto Index (TI), a comparable and widely used measure for molecular similarity. We apply PI and TI to a previously automatically extracted set of binding sites to determine the robustness and usefulness of both scores. RESULTS: We found that PI outperforms TI; moreover, site similarity is poorly defined for TI at values around the 99.5% confidence level for which PI is well defined. A difference map at this confidence level shows that PI gives much more meaningful information than TI. We show individual examples where TI fails to distinguish either a false or a true site paring in contrast to PI, which performs much better. TI cannot handle large or small sites very well, or the comparison of large and small sites, in contrast to PI that is shown to be much more robust. Despite the difficulty of determining a biological 'ground truth' for binding site similarity we conclude that PI is a suitable measure of binding site similarity and could form the basis for a binding site classification scheme comparable to existing protein domain classification schema.     

Interaction of c(60)-fullerene and carboxyfullerene with proteins: docking and binding site alignment

The unique properties of fullerenes have raised the interest of using them for biomedical applications. Within this framework, the interactions of fullerenes with proteins have been an exciting research target, yet little is known about how native proteins can bind fullerenes, and what is the nature of these interactions. Moreover, though some proteins have been shown to interact with fullerenes, up to date, no crystal structure of such complexes was obtained. Here we report docking studies aimed at examining the interactions of fullerene in two forms (C60 nonsubstituted fullerene and carboxyfullerene) with four proteins that are known to bind fullerene derivatives: HIV protease, fullerene-specific antibody, human serum albumin, and bovine serum albumin. Our work provides docking models with detailed binding pockets information, which closely match available experimental data. We further compare the predicted binding sites using a novel multiple binding site alignment method. A high similarity between the physicochemical properties and surface geometry was found for fullerene's binding sites of HIV protease and the human and bovine serum albumins.