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.

Relationship Between Expression of Voltage-Dependent Anion Channel (VDAC) Isoforms and Type of Hexokinase Binding Sites on Brain Mitochondria

Voltage-dependent anion channels (VDAC) are pore-forming proteins found in the outer mitochondrial membrane of eukaryotes. VDACs are known to play an essential role in cellular metabolism and in early stages of apoptosis. In mammals, three VDAC isoforms have been identified. A proteomic approach was exploited to study the expression of VDAC isoforms in rat, bovine, and chicken brain mitochondria. Given the importance of mitochondrially bound hexokinase in regulation of aerobic glycolysis in brain, we studied the possibility that differences in the relative expression of VDAC isoforms may be a factor in determining the species-dependent ratio of type A/type B hexokinase binding sites on brain mitochondria. The spots were characterized, and the signal intensities among spots were compared. VDAC1 was the most abundantly expressed of the three isoforms. Moreover the expression of VDAC1 plus VDAC2 was significantly higher in bovine than in rat brain. Chicken brain mitochondria showed the highest VDAC1 expression and the lowest of VDAC2. Bovine brain mitochondria had the highest VDAC2 levels. We concluded that the nature of hexokinase binding site is not determined by the expression of a single VDAC isoform.

Regulation of Na(+)/H(+) Exchanger Isoform 1 (NHE1) by Calmodulin-binding Sites: Role of Angiotensin II

We examined the effect of Angiotensin II (Ang II) on the interaction between the Ca(2+)/CaM complex and hNHE1. Considering that calmodulin binds to NHE1 at two sites (A and B), amino acids at both sites were modified and two mutants were constructed: SA(1K3R/4E) and SB(1K3R/4E). Wild type and mutants were transfected into PS120 cells and their activity was examined by H(+) flux (J(H+)). The basal J(H+) of wild type was 4.71 +/- 0.57 (mM/min), and it was similar in both mutants. However, the mutations partially impaired the binding of CaM to hNHE1. Ang II (10(-12) and 10(-9) M) increased the J(H+) in wild type and SB. Ang II (10(-6) M) increased this parameter only in SA. Ang II (10(-9) M) maintained the expression of calmodulin in wild type or mutants, and Ang II (10(-6) M) decreased it in wild type or SA, but not in SB. Dimethyl-Bapta-AM (10(-7) M), a calcium chelator, suppressed the effect of Ang II (10(-9) M) in wild type. With Ang II (10(-6) M), Bapta failed to affect wild type or SA, but it increased the J(H+) in SB. W13 or calmidazolium chloride (10(-5) M), two distinct calmodulin inhibitors, decreased the effect of Ang II (10(-9) M) in wild type or SB. With Ang II (10(-6) M), W13 or calmidazolium chloride decreased the J(H+) in wild type or SA and increased it in SB. Thus, with Ang II (10(-12) and 10(-9) M), site A seems to be responsible for the stimulation of hNHE1 and with Ang II (10(-6) M), site B is important to maintain its basal activity. Copyright (C) 2010 S. Karger AG, Basel

Novel Zn(2+)-binding Sites in Human Transthyretin IMPLICATIONS FOR AMYLOIDOGENESIS AND RETINOL-BINDING PROTEIN RECOGNITION

Human transthyretin (TTR) is a homotetrameric protein involved in several amyloidoses. Zn(2+) enhances TTR aggregation in vitro, and is a component of ex vivo TTR amyloid fibrils. We report the first crystal structure of human TTR in complex with Zn(2+) at pH 4.6-7.5. All four structures reveal three tetra-coordinated Zn(2+)-binding sites (ZBS 1-3) per monomer, plus a fourth site (ZBS 4) involving amino acid residues from a symmetry-related tetramer that is not visible in solution by NMR.Zn(2+) binding perturbs loop E-alpha-helix-loop F, the region involved in holo-retinol-binding protein (holo-RBP) recognition, mainly at acidic pH; TTR affinity for holo-RBP decreases similar to 5-fold in the presence of Zn(2+). Interestingly, this same region is disrupted in the crystal structure of the amyloidogenic intermediate of TTR formed at acidic pH in the absence of Zn(2+). HNCO and HNCA experiments performed in solution at pH 7.5 revealed that upon Zn(2+) binding, although the alpha-helix persists, there are perturbations in the resonances of the residues that flank this region, suggesting an increase in structural flexibility. While stability of the monomer of TTR decreases in the presence of Zn(2+), which is consistent with the tertiary structural perturbation provoked by Zn(2+) binding, tetramer stability is only marginally affected by Zn(2+). These data highlight structural and functional roles of Zn(2+) in TTR-related amyloidoses, as well as in holo-RBP recognition and vitamin A homeostasis.

Kinetic Analysis of Gill (Na+,K+)-ATPase Activity in Selected Ontogenetic Stages of the Amazon River Shrimp, (Decapoda, Palaemonidae): Interactions at ATP- and Cation-Binding Sites

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

CONCERTED ACTION OF THE HIGH-AFFINITY CALCIUM-BINDING SITES IN SKELETAL-MUSCLE TROPONIN-C

Mutants of each of the four divalent cation binding sites of chicken skeletal muscle troponin C (TnC) were constructed using site directed mutagenesis to convert Asp to Ala at the first coordinating position in each site. With a view to evaluating the importance of site-site interactions both within and between the N- and C-terminal domains, in this study the mutants are examined for their ability to associate with other components of the troponin-tropomyosin regulatory complex and to regulate thin filaments. The functional effects of each mutation in reconstitution assays are largely confined to the domain in which it occurs, where the unmutated site is unable to compensate for the defect, Thus the mutants of sites I and II bind to the regulatory complex but are impaired in ability to regulate tension and actomyosin ATPase activity, whereas the mutants of sites III and IV regulate activity but are unable to remain bound to thin filaments unless Ca2+ is present. When all four sites are intact, free Mg2+ causes a 50-60-fold increase in TnC's affinity for the other components of the regulatory complex, allowing it to attach firmly to thin filaments. Calcium can replace Mg2+ at a concentration ratio of 1:5000, and at this ratio the Ca2 . TnC complex is more tightly bound to the filaments than the Mg2 . TnC form, In the C-terminal mutants, higher concentrations of Ca2+ (above tension threshold) are required to effect this transformation than in the recombinant wild-type protein, suggesting that the mutants reveal an attachment mediated by Ca2+ in the N-domain sites.

Conformational analyses and docking studies of a series of 5-nitrofuran- and 5-nitrothiophen-semicarbazone derivatives in three possible binding sites of trypanothione and glutathione reductases

To explore three possible binding sites of trypanothione and glutathione reductase, namely, the active, the dimer interface and the coenzyme NADPH binding site, a series of eight compounds, nitrofurans and nitrothiophenes derivatives, were docked, using their crystallographic and modeled conformations. Docking results showed that, for both families and both enzymes, compounds are more likely to bind in the interface site, even though there is some probability of binding in the active site. These studies are in agreement with experimental data, which suggest that these class of compounds can act either as uncompetitive or mixed type inhibitors, and also with the finding that there is an alpha-helix which connects the active with the interface site, thus allowing charge transference between them. (c) 2005 Elsevier B.V. All rights reserved.

Mapping eIF5A binding sites for Dys1 and Lia1: In vivo evidence for regulation of eIF5A hypusination

The evolutionarily conserved factor eIF5A is the only protein known to undergo hypusination, a unique posttranslational modification triggered by deoxyhypusine synthase (Dys1). Although eIF5A is essential for cell viability, the function of this putative translation initiation factor is still obscure. To identify eIF5A-binding proteins that could clarify its function, we screened a two-hybrid library and identified two eIF-5A partners in S. cerevisiae: Dys1 and the protein encoded by the gene YJR070C, named Lia1 (Ligand of eIF5A). The interactions were confirmed by GST pulldown. Mapping binding sites for these proteins revealed that both eIF5A domains can bind to Dys1, whereas the C-terminal domain is sufficient to bind Lia1. We demonstrate for the first time in vivo that the N-terminal α-helix of Dys1 can modulate enzyme activity by inhibiting eIF5A interaction. We suggest that this inhibition be abrogated in the cell when hypusinated and functional eIF5A is required. © 2003 Published by Elsevier B.V. on behalf of the Federation of European Biochemical Societies.