The FunFOLD2 server for the prediction of protein-ligand interactions

The FunFOLD2 server is a new independent server that integrates our novel protein–ligand binding site and quality assessment protocols for the prediction of protein function (FN) from sequence via structure. Our guiding principles were, first, to provide a simple unified resource to make our function prediction software easily accessible to all via a simple web interface and, second, to produce integrated output for predictions that can be easily interpreted. The server provides a clean web interface so that results can be viewed on a single page and interpreted by non-experts at a glance. The output for the prediction is an image of the top predicted tertiary structure annotated to indicate putative ligand-binding site residues. The results page also includes a list of the most likely binding site residues and the types of predicted ligands and their frequencies in similar structures. The protein–ligand interactions can also be interactively visualized in 3D using the Jmol plug-in. The raw machine readable data are provided for developers, which comply with the Critical Assessment of Techniques for Protein Structure Prediction data standards for FN predictions. The FunFOLD2 webserver is freely available to all at the following web site: http://www.reading.ac.uk/bioinf/FunFOLD/FunFOLD_form_2_0.html.

FunFOLDQA: a quality assessment tool for protein-ligand binding site residue predictions

The estimation of prediction quality is important because without quality measures, it is difficult to determine the usefulness of a prediction. Currently, methods for ligand binding site residue predictions are assessed in the function prediction category of the biennial Critical Assessment of Techniques for Protein Structure Prediction (CASP) experiment, utilizing the Matthews Correlation Coefficient (MCC) and Binding-site Distance Test (BDT) metrics. However, the assessment of ligand binding site predictions using such metrics requires the availability of solved structures with bound ligands. Thus, we have developed a ligand binding site quality assessment tool, FunFOLDQA, which utilizes protein feature analysis to predict ligand binding site quality prior to the experimental solution of the protein structures and their ligand interactions. The FunFOLDQA feature scores were combined using: simple linear combinations, multiple linear regression and a neural network. The neural network produced significantly better results for correlations to both the MCC and BDT scores, according to Kendall’s τ, Spearman’s ρ and Pearson’s r correlation coefficients, when tested on both the CASP8 and CASP9 datasets. The neural network also produced the largest Area Under the Curve score (AUC) when Receiver Operator Characteristic (ROC) analysis was undertaken for the CASP8 dataset. Furthermore, the FunFOLDQA algorithm incorporating the neural network, is shown to add value to FunFOLD, when both methods are employed in combination. This results in a statistically significant improvement over all of the best server methods, the FunFOLD method (6.43%), and one of the top manual groups (FN293) tested on the CASP8 dataset. The FunFOLDQA method was also found to be competitive with the top server methods when tested on the CASP9 dataset. To the best of our knowledge, FunFOLDQA is the first attempt to develop a method that can be used to assess ligand binding site prediction quality, in the absence of experimental data.

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.

Probing protein−tannin interactions by isothermal titration microcalorimetry

Isothermal titration microcalorimetry (ITC) has been applied to investigate protein−tannin interactions. Two hydrolyzable tannins were studied, namely myrabolan and tara tannins, for their interaction with bovine serum albumin (BSA), a model globular protein, and gelatin, a model proline-rich random coil protein. Calorimetry data indicate that protein−tannin interaction mechanisms are dependent upon the nature of the protein involved. Tannins apparently interact nonspecifically with the globular BSA, leading to binding saturation at estimated tannin/BSA molar ratios of 48:1 for tara- and 178:1 for myrabolan tannins. Tannins bind to the random coil protein gelatin by a two-stage mechanism. The energetics of the first stage show evidence for cooperative binding of tannins to the protein, while the second stage indicates gradual saturation of binding sites as observed for interaction with BSA. The structure and flexibility of the tannins themselves alters the stoichiometry of the interaction, but does not appear to have any significant affect on the overall binding mechanism observed. This study demonstrates the potential of ITC for providing an insight into the nature of protein−tannin interactions.

Probing protein-tannin interactions by isothermal titration microcalorimetry

Isothermal titration microcalorimetry (ITC) has been applied to investigate protein-tannin interactions. Two hydrolyzable tannins were studied, namely myrabolan and tara tannins, for their interaction with bovine serum albumin (BSA), a model globular protein, and gelatin, a model proline-rich random coil protein. Calorimetry data indicate that protein-tannin interaction mechanisms are dependent upon the nature of the protein involved. Tannins apparently interact nonspecifically with the globular BSA, leading to binding saturation at estimated tannin/BSA molar ratios of 48:1 for tara- and 178:1 for myrabolan tannins. Tannins bind to the random coil protein gelatin by a two-stage mechanism. The energetics of the first stage show evidence for cooperative binding of tannins to the protein, while the second stage indicates gradual saturation of binding sites as observed for interaction with BSA. The structure and flexibility of the tannins themselves alters the stoichiometry of the interaction, but does not appear to have any significant affect on the overall binding mechanism observed. This study demonstrates the potential of ITC for providing an insight into the nature of protein-tannin interactions.

Characterization of protein-polyphenol interactions

Dietary polyphenols have received attention for their biologically significant functions as antioxidants, anticarcinogens or antimutagens, which have led to their recognition as potential nutraceuticals. Polyphenols also characteristically possess a significant binding affinity for proteins, which can lead to the formation of soluble and insoluble protein-polyphenol complexes. Questions remain concerning whether and to what extent the protein-polyphenol interaction influences functionality. For example, is the formation of protein-polyphenol complexes an obstacle to the nutritional bioavailability of either species? This article discusses the development of suitable methodologies to investigate the physicochemical basis of protein-polyphenol interactions and the influence of structure-activity relationships on binding affinities. (C) 2003 Elsevier Ltd. All rights reserved.

Protein-lipid interactions at the air/water interface

Surface pressure measurements and external reflection FTIR spectroscopy have been used to probe protein-lipid interactions at the air/water interface. Spread monomolecular layers of stearic acid and phosphocholine were prepared and held at different compressed phase states prior to the introduction of protein to the buffered subphase. Contrasting interfacial behaviour of the proteins, albumin and lysozyme, was observed and revealed the role of both electrostatic and hydrophobic interactions in protein adsorption. The rate of adsorption of lysozyme to the air/water interface increased dramatically in the presence of stearic acid, due to strong electrostatic interactions between the negatively charged stearic acid head group and lysozyme, whose net charge at pH 7 is positive. Introduction of albumin to the subphase resulted in solubilisation of the stearic acid via the formation of an albumin-stearic acid complex and subsequent adsorption of albumin. This observation held for both human and bovine serum albumin. Protein adsorption to a PC layer held at low surface pressure revealed adsorption rates similar to adsorption to the bare air/water interface and suggested very little interaction between the protein and the lipid. For PC layers in their compressed phase state some adsorption of protein occurred after long adsorption times. Structural changes of both lysozyme and albumin were observed during adsorption, but these were dramatically reduced in the presence of a lipid layer compared to that of adsorption to the pure air/water interface.

Probing protein-tannin interactions by isothermal titration microcalorimetry

Isothermal titration microcalorimetry (ITC) has been applied to investigate protein-tannin interactions. Two hydrolyzable tannins were studied, namely myrabolan and tara tannins, for their interaction with bovine serum albumin (BSA), a model globular protein, and gelatin, a model proline-rich random coil protein. Calorimetry data indicate that protein-tannin interaction mechanisms are dependent upon the nature of the protein involved. Tannins apparently interact nonspecifically with the globular BSA, leading to binding saturation at estimated tannin/BSA molar ratios of 48:1 for tara- and 178:1 for myrabolan tannins. Tannins bind to the random coil protein gelatin by a two-stage mechanism. The energetics of the first stage show evidence for cooperative binding of tannins to the protein, while the second stage indicates gradual saturation of binding sites as observed for interaction with BSA. The structure and flexibility of the tannins themselves alters the stoichiometry of the interaction, but does not appear to have any significant affect on the overall binding mechanism observed. This study demonstrates the potential of ITC for providing an insight into the nature of protein-tannin interactions.

IntFOLD: an integrated server for modelling protein structures and functions from amino acid sequences

IntFOLD is an independent web server that integrates our leading methods for structure and function prediction. The server provides a simple unified interface that aims to make complex protein modelling data more accessible to life scientists. The server web interface is designed to be intuitive and integrates a complex set of quantitative data, so that 3D modelling results can be viewed on a single page and interpreted by non-expert modellers at a glance. The only required input to the server is an amino acid sequence for the target protein. Here we describe major performance and user interface updates to the server, which comprises an integrated pipeline of methods for: tertiary structure prediction, global and local 3D model quality assessment, disorder prediction, structural domain prediction, function prediction and modelling of protein-ligand interactions. The server has been independently validated during numerous CASP (Critical Assessment of Techniques for Protein Structure Prediction) experiments, as well as being continuously evaluated by the CAMEO (Continuous Automated Model Evaluation) project. The IntFOLD server is available at: http://www.reading.ac.uk/bioinf/IntFOLD/

Immersive Projection for biochemical pattern matching

In all biological processes, protein molecules and other small molecules interact to function and form transient macromolecular complexes. This interaction of two or more molecules can be described by a docking event. Docking is an important phase for structure-based drug design strategies, as it can be used as a method to simulate protein-ligand interactions. Various docking programs exist that allow automated docking, but most of them have limited visualization and user interaction. It would be advantageous if scientists could visualize the molecules participating in the docking process, manipulate their structures and manually dock them before submitting the new conformations to an automated docking process in an immersive environment, which can help stimulate the design/docking process. This also could greatly reduce docking time and resources. To achieve this, we propose a new virtual modelling/docking program, whereby the advantages of virtual modelling programs and the efficiency of the algorithms in existing docking programs will be merged.

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.

An insight into the mechanism of protein separation by colloidal gas aphrons (CGA) generated from ionic surfactants

Colloidal gas aphrons (CGA), which are surfactant stabilised microbubbles, have been previously applied for the recovery of proteins from model mixtures and a few studies have demonstrated the potential of these dispersions for the selective recovery of proteins from complex mixtures. However there is a lack of understanding of the mechanism of separation and forces governing the selectivity of the separation. In this paper a mechanistic study is carried out to determine the main factors and forces influencing the selectivity of separation of whey proteins with CGA generated from ionic surfactants. Two different separation strategies were followed: (i) separation of lactoferrin and lactoperoxidase by anionic CGA generated from a solution of sodium bis-(2-ethyl hexyl) sulfosuccinate (AOT); (ii) separation of beta-lactoglobulin by cationic CGA generated from a solution of cetyltrimethylammonium bromide (CTAB). Separation results indicate that electrostatic interactions are the main forces determining the selectivity however these could not completely explain the selectivities obtained following both strategies. Protein-surfactant interactions were studied by measuring the zeta potential changes on individual proteins upon addition of surfactant and at varying pH. Interestingly strongest electrostatic interactions were measured at those pH and surfactant to protein mass ratios which were optimum for protein separation. Effect of surfactant on protein conformation was determined by measuring the change in fluorescence intensity upon addition of surfactant at varying pH. Differences in the fluorescence patterns were detected among proteins which were correlated to differences in their conformational features which could in turn explain their different separation behaviour. The effect of conformation on selectivity was further proven by experiments in which conformational changes were induced by pre-treatment of whey (heating) and by storage at 4 degrees C. Overall it can be concluded that separation of proteins by ionic CGA is driven mainly by electrostatic interactions however conformational features will finally determine the selectivity of the separation with competitive adsorption having also an effect. (c) 2006 Elsevier B.V. All rights reserved.

Molecular dynamics simulations of proteins: can the explicit water model be varied?

In molecular mechanics simulations of biological systems, the solvation water is typically represented by a default water model which is an integral part of the force field. Indeed, protein nonbonding parameters are chosen in order to obtain a balance between water-water and protein-water interactions and hence a reliable description of protein solvation. However, less attention has been paid to the question of whether the water model provides a reliable description of the water properties under the chosen simulation conditions, for which more accurate water models often exist. Here we consider the case of the CHARMM protein force field, which was parametrized for use with a modified TIP3P model. Using quantum mechanical and molecular mechanical calculations, we investigate whether the CHARMM force field can be used with other water models: TIP4P and TIP5P. Solvation properties of N-methylacetamide (NMA), other small solute molecules, and a small protein are examined. The results indicate differences in binding energies and minimum energy geometries, especially for TIP5P, but the overall description of solvation is found to be similar for all models tested. The results provide an indication that molecular mechanics simulations with the CHARMM force field can be performed with water models other than TIP3P, thus enabling an improved description of the solvent water properties.