The binding site distance test score: a robust method for the assessment of predicted protein binding sites

We propose a novel method for scoring the accuracy of protein binding site predictions – the Binding-site Distance Test (BDT) score. Recently, the Matthews Correlation Coefficient (MCC) has been used to evaluate binding site predictions, both by developers of new methods and by the assessors for the community wide prediction experiment – CASP8. Whilst being a rigorous scoring method, the MCC does not take into account the actual 3D location of the predicted residues from the observed binding site. Thus, an incorrectly predicted site that is nevertheless close to the observed binding site will obtain an identical score to the same number of nonbinding residues predicted at random. The MCC is somewhat affected by the subjectivity of determining observed binding residues and the ambiguity of choosing distance cutoffs. By contrast the BDT method produces continuous scores ranging between 0 and 1, relating to the distance between the predicted and observed residues. Residues predicted close to the binding site will score higher than those more distant, providing a better reflection of the true accuracy of predictions. The CASP8 function predictions were evaluated using both the MCC and BDT methods and the scores were compared. The BDT was found to strongly correlate with the MCC scores whilst also being less susceptible to the subjectivity of defining binding residues. We therefore suggest that this new simple score is a potentially more robust method for future evaluations of protein-ligand binding site predictions.

Interactions of decay-accelerating factor (DAF) with haemagglutinating human enteroviruses: utilizing variation in primate DAF to map virus binding sites

A cellular receptor for the haemagglutinating enteroviruses (HEV), and the protein that mediates haemagglutination, is the membrane complement regulatory protein decay accelerating factor (DAF; CD55). Although primate DAF is highly conserved, significant differences exist to enable cell lines derived from primates to be utilized for the characterization of the DAF binding phenotype of human enteroviruses. Thus, several distinct DAF-binding phenotypes of a selection of HEVs (viz. coxsackievirus A21 and echoviruses 6, 7, 11-13, 29) were identified from binding and infection assays using a panel of primate cells derived from human, orang-utan, African Green monkey and baboon tissues. These studies complement our recent determination of the crystal structure of SCR(34) of human DAF [Williams, P., Chaudhry, Y., Goodfellow, I. G., Billington, J., Powell, R., Spiller, O. B., Evans, D. J. & Lea, S. (2003). J Biol Chem 278, 10691-10696] and have enabled us to better map the regions of DAF with which enteroviruses interact and, in certain cases, predict specific virus-receptor contacts.

In silico identification and characterization of protein-ligand binding sites

Protein–ligand binding site prediction methods aim to predict, from amino acid sequence, protein–ligand interactions, putative ligands, and ligand binding site residues using either sequence information, structural information, or a combination of both. In silico characterization of protein–ligand interactions has become extremely important to help determine a protein’s functionality, as in vivo-based functional elucidation is unable to keep pace with the current growth of sequence databases. Additionally, in vitro biochemical functional elucidation is time-consuming, costly, and may not be feasible for large-scale analysis, such as drug discovery. Thus, in silico prediction of protein–ligand interactions must be utilized to aid in functional elucidation. Here, we briefly discuss protein function prediction, prediction of protein–ligand interactions, the Critical Assessment of Techniques for Protein Structure Prediction (CASP) and the Continuous Automated EvaluatiOn (CAMEO) competitions, along with their role in shaping the field. We also discuss, in detail, our cutting-edge web-server method, FunFOLD for the structurally informed prediction of protein–ligand interactions. Furthermore, we provide a step-by-step guide on using the FunFOLD web server and FunFOLD3 downloadable application, along with some real world examples, where the FunFOLD methods have been used to aid functional elucidation.

Binding of pentagalloyl glucose to two globular proteins occurs via multiple surface sites

The interaction between pentagalloyl glucose (PGG) and two globular proteins, bovine serum albumin (BSA) and ribulose-1,5-bisphosphate carboxylase oxygenase (rubisco), was investigated by isothermal titration calorimetry (ITC). ITC data fit to a binding model consisting of two sets of multiple binding sites, which reveal similarities in the mode of binding of PGG to BSA and rubisco. In both cases, the interaction is characterized by a high number of binding sites, which suggests that binding occurs by a surface adsorption mechanism that leads to coating of the protein surface, which promotes aggregation and precipitation of the PGG-protein complex. This model was confirmed by turbidimetry analysis of the PGG-BSA interaction. Analysis of tryptophan fluorescence quenching during the interaction of PGG with BSA suggests that binding of PGG leads to some conformational changes that are energetically closer to the unfolded state of the BSA structure, because small red shifts in the resulting emission spectra were observed.

Hydrolyzable tannin structures influence relative globular and random coil protein binding strengths

Binding parameters for the interactions of pentagalloyl glucose (PGG) and four hydrolyzable tannins (representing gallotannins and ellagitannins) with gelatin and bovine serum albumin (BSA) have been determined from isothermal titration calorimetry data. Equilibrium binding constants determined for the interaction of PGG and isolated mixtures of tara gallotannins and of sumac gallotannins with gelatin and BSA were of the same order of magnitude for each tannin (in the range of 10(4)-10(5) M-1 for stronger binding sites when using a binding model consisting of two sets of multiple binding sites). In contrast, isolated mixtures of chestnut ellagitannins and of myrabolan ellagitannins exhibited 3-4 orders of magnitude greater equilibrium binding constants for the interaction with gelatin (similar to 2 x 10(6) M-1) than for that with BSA (similar to 8 x 10(2) M-1). Binding stoichiometries revealed that the stronger binding sites on gelatin outnumbered those on BSA by a ratio of at least similar to 2:1 for all of the hydrolyzable tannins studied. Overall, the data revealed that relative binding constants for the interactions with gelatin and BSA are dependent on the structural flexibility of the tannin molecule.

Thermodynamic analysis of ferrous ion binding to Escherichia coli ferritin EcFtnA

Iron oxidation in the bacterial ferritin EcFtnA from Escherichia coli shows marked differences from its homologue human H-chain ferritin (HuHF). While the amino acid residues that constitute the dinuclear center in these proteins are highly conserved, EcFtnA has a third iron-binding site (C site) in close proximity to the dinuclear center that is seemingly responsible for these differences. Here, we describe the first thermodynamic study of Fe2+ binding to EcFtnA and its variants to determine the location of the primary ferrous ion-binding sites on the protein and to better understand the role of the third C site in iron binding. Isothermal titration calorimetric analyses of the wild-type protein reveal the presence of two main classes of binding sites in the pH range of 6.5-7.5, ascribed to Fe2+ binding, first at the A and then the B sites. Site-directed mutagenesis of ligands in the A, B, or C sites affects the apparent Fe2+-binding stoichiometries at the unaltered sites. The data imply some degree of inter- and intrasubunit negative cooperative interaction between sites. Unlike HuHF where only the A site initially binds Fe2+, both A and B sites in EcFtnA bind Fe2+, implying a role for the C site in influencing the binding of Fe2+ at the B site of the di-iron center of EcFtnA. The ITC equations describing a binding model for three classes of independent binding sites are reported here for the first time.

Development of novel brush-type chiral stationary phases based on terpenoid selectors: HPLC evaluation and theoretical investigation of enantioselective binding interactions

The terpenoid chiral selectors dehydroabietic acid, 12,14-dinitrodehydroabietic acid and friedelin have been covalently linked to silica gel yielding three chiral stationary phases CSP 1, CSP 2 and CSP 3, respectively. The enantiodiscriminating capability of each one of these phases was evaluated by HPLC with four families of chiral aromatic compounds composed of alcohols, amines, phenylalanine and tryptophan amino acid derivatives and beta-lactams. The CSP 3 phase, containing a selector with a large friedelane backbone is particularly suitable for resolving free alcohols and their derivatives bearing fluorine substituents, while CSP 2 with a dehydroabietic architecture is the only phase that efficiently discriminates 1, 1'-binaphthol atropisomers. CSP 3 also gives efficient resolution of the free amines. All three phases resolve well the racemates of N-trifluoracetyl and N-3,5-dinitrobenzoyl phenylalanine amino acid ester derivatives. Good enantioseparation of beta-lactams and N-benzoyl tryptophan amino acid derivatives was achieved on CSP 1. In order to understand the structural factors that govern the chiral molecular recognition ability of these phases, molecular dynamics simulations were carried out in the gas phase with binary diastereomeric complexes formed by the selectors of CSP 1 and CSP 2 and several amino acid derivatives. Decomposition of molecular mechanics energies shows that van der Waals interactions dominate the formation of the diastereomeric transient complexes while the electrostatic binding interactions are primarily responsible for the enantioselective binding of the (R)- and (S)-analytes. Analysis of the hydrogen bonds shows that electrostatic interactions are mainly associated with the formation of N-(HO)-O-...=C enantio selective hydrogen bonds between the amide binding sites from the selectors and the carbonyl groups of the analytes. The role of mobile phase polarity, a mixture of n-hexane and propan-2-ol in different ratios, was also evaluated through molecular dynamics simulations in explicit solvent. (c) 2006 Elsevier Ltd. All rights reserved.

High-affinity small molecule-phospholipid complex formation: Binding of siramesine to phosphatidic acid

Siramesine (SRM) is a sigma-2 receptor agonist which has been recently shown to inhibit growth of cancer cells. Fluorescence spectroscopy experiments revealed two distinct binding sites for this drug in phospholipid membranes. More specifically, acidic phospholipids retain siramesine on the bilayer surface due to a high-affinity interaction, reaching saturation at an apparent 1:1 drug-acidic phospholipid stoichiometry, where after the drug penetrates into the hydrocarbon core of the membrane. This behavior was confirmed using Langmuir films. Of the anionic phospholipids, the highest affinity, comparable to the affinities for the binding of small molecule ligands to proteins, was measured for phosphatidic acid (PA, mole fraction Of X-PA = 0.2 in phosphatidylcholine vesicles), yielding a molecular partition coefficient of 240 +/- 80 x 10(6). An MD simulation on the siramesine:PA interaction was in agreement with the above data. Taking into account the key role of PA as a signaling molecule promoting cell growth our results suggest a new paradigm for the development of anticancer drugs, viz. design of small molecules specifically scavenging phospholipids involved in the signaling cascades controlling cell behavior.

Isothermal titration calorimetry study of epicatechin binding to serum albumin

The interaction of epicatechin with bovine serum albumin (BSA) was studied by isothermal titration calorimetry. The binding constant (K) and associated thermodynamic binding parameters (n, Delta H) were determined for the interaction at three solution concentrations of BSA using a binding model assuming independent binding sites. These data show weak non-covalent binding of epicatechin to BSA. The interaction energetics varied with BSA concentration in the calorimeter cell, suggesting that the binding of epicatechin induced BSA aggregation. The free energy (Delta G) remained constant within a range of 2 kJ mol(-1) and negative entropy was observed, indicating an enthalpy driven exothermic interaction. It is concluded that the non-covalent epicatechin-BSA complex is formed by hydrogen bonding. (c) 2006 Elsevier B.V. All rights reserved.

SARS coronavirus unique domain: three-domain molecular architecture in solution and RNA binding

Nonstructural protein 3 of the severe acute respiratory syndrome (SARS) coronavirus includes a "SARS-unique domain" (SUD) consisting of three globular domains separated by short linker peptide segments. This work reports NMR structure determinations of the C-terminal domain (SUD-C) and a two-domain construct (SUD-MC) containing the middle domain (SUD-M) and the C-terminal domain, and NMR data on the conformational states of the N-terminal domain (SUD-N) and the SUD-NM two-domain construct. Both SUD-N and SUD-NM are monomeric and globular in solution; in SUD-NM, there is high mobility in the two-residue interdomain linking sequence, with no preferred relative orientation of the two domains. SUD-C adopts a frataxin like fold and has structural similarity to DNA-binding domains of DNA-modifying enzymes. The structures of both SUD-M (previously determined) and SUD-C (from the present study) are maintained in SUD-MC, where the two domains are flexibly linked. Gel-shift experiments showed that both SUD-C and SUD-MC bind to single-stranded RNA and recognize purine bases more strongly than pyrimidine bases, whereby SUD-MC binds to a more restricted set of purine-containing RNA sequences than SUD-M. NMR chemical shift perturbation experiments with observations of (15)N-labeled proteins further resulted in delineation of RNA binding sites (i.e., in SUD-M, a positively charged surface area with a pronounced cavity, and in SUD-C, several residues of an anti-parallel beta-sheet). Overall, the present data provide evidence for molecular mechanisms involving the concerted actions of SUD-M and SUD-C, which result in specific RNA binding that might be unique to the SUD and, thus, to the SARS coronavirus.

Titanocene anticancer complexes and their binding mode of action to human serum albumin: a computational study

Due to the pivotal role played by human serum albumin (HSA) in the transport and cytotoxicity of titanocene complexes, a docking study has been performed on a selected set of titanocene complexes to aid in the current understanding of the potential mode of action of these titanocenes upon binding HSA. Analysis of the docking results has revealed potential binding at the known drug binding sites in HSA and has provided some explanation for the specificity and subsequent cytotoxicity of these titanocenes. Additionally, a new alternative binding site for these titanocenes has been postulated.

Fluorescence lifetime imaging microscopy (FLIM) to demonstrate the nuclear binding of flavanols and (--epigallocatechin gallate

The use of light microscopy and DMACA staining strongly suggested that plant and animal cell nuclei act as sinks for flavanols [1, 2]. Detailed uv-vis spectroscopic titration experiments indicated that histone proteins are the likely binding sites in the nucleus [2]. Here we report the development of a multi-photon excitation microscopy technique combined with fluorescent lifetime measurements of flavanols. Using this technique, (+) catechin, (-) epicatechin and (-) epigallocatechin gallate (EGCG) showed strikingly different excited state lifetimes in solution. Interaction of histone proteins with flavanols was indicated by the appearance of a significant τ2-component of 1.7 to 4.0ns. Tryptophan interference could be circumvented in the in vivo fluorescence lifetime imaging microscopy (FLIM) experiments with 2-photon excitation at 630nm. This enabled visualisation and semi-quantitative measurements that demonstrated unequivocally the absorption of (+)catechin, (-)epicatechin and EGCG by nuclei of onion cells. 3D FLIM revealed for the first time that externally added EGCG penetrated the whole nucleus in onion cells. The relative proportions of EGCG in cytoplasm: nucleus: nucleoli were ca. 1:10:100. FLIM experiments may therefore facilitate probing the health effects of EGCG, which is the major constituent of green tea.

Predicting metal-binding site residues in low-resolution structural models

The accurate prediction of the biochemical function of a protein is becoming increasingly important, given the unprecedented growth of both structural and sequence databanks. Consequently, computational methods are required to analyse such data in an automated manner to ensure genomes are annotated accurately. Protein structure prediction methods, for example, are capable of generating approximate structural models on a genome-wide scale. However, the detection of functionally important regions in such crude models, as well as structural genomics targets, remains an extremely important problem. The method described in the current study, MetSite, represents a fully automatic approach for the detection of metal-binding residue clusters applicable to protein models of moderate quality. The method involves using sequence profile information in combination with approximate structural data. Several neural network classifiers are shown to be able to distinguish metal sites from non-sites with a mean accuracy of 94.5%. The method was demonstrated to identify metal-binding sites correctly in LiveBench targets where no obvious metal-binding sequence motifs were detectable using InterPro. Accurate detection of metal sites was shown to be feasible for low-resolution predicted structures generated using mGenTHREADER where no side-chain information was available. High-scoring predictions were observed for a recently solved hypothetical protein from Haemophilus influenzae, indicating a putative metal-binding site.

Size and molecular flexibility affect the binding of ellagitannins to bovine serum albumin

Binding to bovine serum albumin of monomeric (vescalagin and pedunculagin) and dimeric ellagitannins (roburin A, oenothein B, and gemin A) was investigated by isothermal titration calorimetry and fluorescence spectroscopy, which indicated two types of binding sites. Stronger and more specific sites exhibited affinity constants, K1, of 104–106 M–1 and stoichiometries, n1, of 2–13 and dominated at low tannin concentrations. Weaker and less-specific binding sites had K2 constants of 103–105 M–1 and stoichiometries, n2, of 16–30 and dominated at higher tannin concentrations. Binding to stronger sites appeared to be dependent on tannin flexibility and the presence of free galloyl groups. Positive entropies for all but gemin A indicated that hydrophobic interactions dominated during complexation. This was supported by an exponential relationship between the affinity, K1, and the modeled hydrophobic accessible surface area and by a linear relationship between K1 and the Stern–Volmer quenching constant, KSV.

Binding of an oligomeric ellagitannin series to bovine serum albumin (BSA): analysis by isothermal titration calorimetry (ITC)

A unique series of oligomeric ellagitannins was used to study their interactions with bovine serum albumin (BSA) by isothermal titration calorimetry. Oligomeric ellagitannins, ranging from monomer to heptamer and a mixture of octamer–undecamers, were isolated as individual pure compounds. This series allowed studying the effects of oligomer size and other structural features. The monomeric to trimeric ellagitannins deviated most from the overall trends. The interactions of ellagitannin oligomers from tetramers to octa–undecamers with BSA revealed strong similarities. In contrast to the equilibrium binding constant, enthalpy showed an increasing trend from the dimer to larger oligomers. It is likely that first the macrocyclic part of the ellagitannin binds to the defined binding sites on the protein surface and then the “flexible tail” of the ellagitannin coats the protein surface. The results highlight the importance of molecular flexibility to maximize binding between the ellagitannin and protein surfaces.