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.

Proteins and their interacting partners: an introduction to protein–ligand binding site prediction methods

Elucidating the biological and biochemical roles of proteins, and subsequently determining their interacting partners, can be difficult and time consuming using in vitro and/or in vivo methods, and consequently the majority of newly sequenced proteins will have unknown structures and functions. However, in silico methods for predicting protein–ligand binding sites and protein biochemical functions offer an alternative practical solution. The characterisation of protein–ligand binding sites is essential for investigating new functional roles, which can impact the major biological research spheres of health, food, and energy security. In this review we discuss the role in silico methods play in 3D modelling of protein–ligand binding sites, along with their role in predicting biochemical functionality. In addition, we describe in detail some of the key alternative in silico prediction approaches that are available, as well as discussing the Critical Assessment of Techniques for Protein Structure Prediction (CASP) and the Continuous Automated Model EvaluatiOn (CAMEO) projects, and their impact on developments in the field. Furthermore, we discuss the importance of protein function prediction methods for tackling 21st century problems.

The leptospiral antigen Lp49 is a two-domain protein with putative protein binding function

Pathogenic Leptospira is the etiological agent of leptospirosis, a life-threatening disease that affects populations worldwide. Currently available vaccines have limited effectiveness and therapeutic interventions are complicated by the difficulty in making an early diagnosis of leptospirosis. The genome of Leptospira interrogans was recently sequenced and comparative genomic analysis contributed to the identification of surface antigens, potential candidates for development of new vaccines and serodiagnosis. Lp49 is a membrane-associated protein recognized by antibodies present in sera from early and convalescent phases of leptospirosis patients. Its crystal structure was determined by single-wavelength anomalous diffraction using selenomethionine-labelled crystals and refined at 2.0 angstrom resolution. Lp49 is composed of two domains and belongs to the all-beta-proteins class. The N-terminal domain folds in an immunoglobulin-like beta-sandwich structure, whereas the C-terminal domain presents a seven-bladed beta-propeller fold. Structural analysis of Lp49 indicates putative protein-protein binding sites, suggesting a role in Leptospira-host interaction. This is the first crystal structure of a leptospiral antigen described to date. (C) 2008 Elsevier Inc. All rights reserved.

The binding of synthetic triiodo L-thyronine analogs to human transthyretin: Molecular basis of cooperative and non-cooperative ligand recognition

Transthyretin (TTR) is a tetrameric beta-sheet-rich transporter protein directly involved in human amyloid diseases. Several classes of small molecules can bind to TTR delaying its amyloid fibril formation, thus being promising drug candidates to treat TTR amyloidoses. In the present study, we characterized the interactions of the synthetic triiodo L-thyronine analogs and thyroid hormone nuclear receptor TR beta-selecfive agonists GC-1 and GC-24 with the wild type and V30M variant of human transthyretin (TTR). To achieve this aim, we conducted in vitro TTR acid-mediated aggregation and isothermal titration calorimetry experiments and determined the TTR:GC-1 and TTR:GC-24 crystal structures. Our data indicate that both GC-1 and GC-24 bind to TTR in a non-cooperative manner and are good inhibitors of TTR aggregation, with dissociation constants for both hormone binding sites (HBS) in the low micromolar range. Analysis of the crystal structures of TTRwt:GC-1(24) complexes and their comparison with the TTRwt X-ray structure bound to its natural ligand thyroxine (T4) suggests, at the molecular level, the basis for the cooperative process displayed by T4 and the non-cooperative process provoked by both GC-1 and GC-24 during binding to TTR. (C) 2010 Elsevier Inc. All rights reserved.

Flexibility and inhibitor binding in Cdc25 phosphatases

Cdc25 phosphatases involved in cell cycle checkpoints are now active targets for the development of anti-cancer therapies. Rational drug design would certainly benefit from detailed structural information for Cdc25s. However, only apo- or sulfate-bound crystal structures of the Cdc25 catalytic domain have been described so far. Together with previously available crystalographic data, results from molecular dynamics simulations, bioinformatic analysis, and computer-generated conformational ensembles shown here indicate that the last 30-40 residues in the C-terminus of Cdc25B are partially unfolded or disordered in solution. The effect of C-terminal flexibility upon binding of two potent small molecule inhibitors to Cdc25B is then analyzed by using three structural models with variable levels of flexibility, including an equilibrium distributed ensemble of Cdc25B backbone conformations. The three Cdc25B structural models are used in combination with flexible docking, clustering, and calculation of binding free energies by the linear interaction energy approximation to construct and validate Cdc25B-inhibitor complexes. Two binding sites are identified on top and beside the Cdc25B active site. The diversity of interaction modes found increases with receptor flexibility. Backbone flexibility allows the formation of transient cavities or compact hydrophobic units on the surface of the stable, folded protein core that are unexposed or unavailable for ligand binding in rigid and densely packed crystal structures. The present results may help to speculate on the mechanisms of small molecule complexation to partially unfolded or locally disordered proteins.

Docking and molecular dynamics simulation of quinone compounds with trypanocidal activity

In this work, two different docking programs were used, AutoDock and FlexX, which use different types of scoring functions and searching methods. The docking poses of all quinone compounds studied stayed in the same region in the trypanothione reductase. This region is a hydrophobic pocket near to Phe396, Pro398 and Leu399 amino acid residues. The compounds studied displays a higher affinity in trypanothione reductase (TR) than glutathione reductase (GR), since only two out of 28 quinone compounds presented more favorable docking energy in the site of human enzyme. The interaction of quinone compounds with the TR enzyme is in agreement with other studies, which showed different binding sites from the ones formed by cysteines 52 and 58. To verify the results obtained by docking, we carried out a molecular dynamics simulation with the compounds that presented the highest and lowest docking energies. The results showed that the root mean square deviation (RMSD) between the initial and final pose were very small. In addition, the hydrogen bond pattern was conserved along the simulation. In the parasite enzyme, the amino acid residues Leu399, Met400 and Lys402 are replaced in the human enzyme by Met406, Tyr407 and Ala409, respectively. In view of the fact that Leu399 is an amino acid of the Z site, this difference could be explored to design selective inhibitors of TR.

The bacterial chromosome segregation protein Spo0J spreads along DNA from parS nucleation sites

Regulation of chromosome inheritance is essential to ensure proper transmission of genetic information. To accomplish accurate genome segregation, cells organize their chromosomes and actively separate them prior to cytokinesis. In Bacillus subtilis the Spo0J protein is required for accurate chromosome segregation and it regulates the developmental switch from vegetative growth to sporulation. Spo0J is a DNA-binding protein that recognizes at least eight identified parS sites located near the origin of replication. As judged by fluorescence microscopy, Spo0J forms discrete foci associated with the oriC region of the chromosome throughout the cell cycle. In an attempt to determine the mechanisms utilized by Spo0J to facilitate productive chromosome segregation, we have investigated the DNA binding activity of Spo0J. In vivo we find Spo0J associates with several kilobases of DNA flanking its specific binding sites (parS) through a parS-dependent nucleation event that promotes lateral spreading of Spo0J along the chromosome. Using purified components we find that Spo0J has the ability to coat non-specific DNA substrates. These 'Spo0J domains' provide large structures near oriC that could potentially demark, organize or localize the origin region of the chromosome.

H19-DMR allele-specific methylation analysis reveals epigenetic heterogeneity of CTCF binding site 6 but not of site 5 in head-and-neck carcinomas: A pilot case-control analysis

Aberrant methylation of seven potential binding sites of the CTCF factor in the differentially methylated region upstream of the H19 gene (H19-DMR) has been suggested as critical for the regulation of IGF2 and H19 imprinted genes. In this study, we analyzed the allele-specific methylation pattern of CTCF binding sites 5 and 6 using methylationsensitive restriction enzyme PCR followed by RFLP analysis in matched tumoral and lymphocyte DNA from head-and-neck squamous cell carcinoma (HNSCC) patients, as well as in lymphocyte DNA from control individuals who were cancer-free. The monoallelic methylation pattern was maintained in CTCF binding site 5 in 22 heterozygous out of 91 samples analyzed. Nevertheless, a biallelic methylation pattern was detected in CTCF binding site 6 in a subgroup of HNSCC patients as a somatic acquired feature of tumor cells. An atypical biallelic methylation was also observed in both tumor and lymphocyte DNA from two patients, and at a high frequency in the control group (29 out of 64 informative controls). Additionally, we found that the C/T transition detected by HhaI RFLP suppressed one dinucleotide CpG in critical CTCF binding site 6, of a mutation showing polymorphic frequencies. Although a heterogeneous methylation pattern was observed after DNA sequencing modified by sodium bisulfite, the biallelic methylation pattern was confirmed in 9 out of 10 HNSCCs. These findings are likely to be relevant in the epigenetic regulation of the DMR, especially in pathological conditions in which the imprinting of IGF2 and H19 genes is disrupted.

DNA methylation in the CTCF-binding site I and the expression pattern of the H19 gene: Does positive expression predict poor prognosis in early stage head and neck carcinomas?

Loss of allele-specific expression by the imprinted genes IGF2 and H19 has been correlated with a differentially methylated region (DMR) upstream to the H19 gene. The H19-DMR contains seven potential CCCTC-binding factor (CTCF) binding sites. CTCF is a chromatin insulator and a multifunctional transcription factor whose binding to the H19-DMR is suppressed by DNA methylation. Our study included a group of 41 head and neck squamous cell carcinoma (HNSCC) samples. The imprinting status of the H19 gene was analyzed in 11 out of 35 positive cases for H19 gene expression, and only 1 of them showed loss of imprinting. We detected a significant correlation (P=0.041, Fisher's exact test) between H19 expression and tumor recurrence. Among H19 positive cases, six were T2, in which five developed recurrence and/or metastasis. Inversely, in the group of tumors that showed no H19 gene expression, 5 out of 24 were T2 and only I presented regional recurrence. These data support the hypothesis that H19 expression could be used as a prognostic marker to indicate recurrence in early stage tumors. We also examined the methylation of the CTCF binding site 1 in a subgroup of these samples. The H19 gene silencing and loss of imprinting were not correlated with the methylation pattern of the CTCF binding site 1. However, the significant correlation between H19 expression and tumor recurrence suggest that this transcript could be a marker for the progression of HNSCC. (c) 2005 Wiley-Liss, Inc.

The Influence of Solutes on the Enthalpy/Entropy Change of the Actinomycin D Binding to DNA: Hydration, Energy Compensation and Long-Range Deformation on DNA

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

CHARACTERISTICS OF ZINC-BINDING TO HUMAN RED-BLOOD-CELL MEMBRANES

The objective of the present study was to standardize the analysis of zinc binding on human red blood cell (RBC) membranes in 20 normal adults. The displacement studies revealed that at the maximal stable zinc concentration tested (600 muM), 57% (mean) of the bound Zn-65 was displaced and to displace half maximal Zn-65, the stable zinc concentration was 300 muM. Scatchard plots revealed two classes of binding sites for zinc on RBC membranes: one with higher affinity, Kd = 1.20 x 10(-5) M (site I), and the other with lower affinity, Kd = 2.77 x 10(-4) M (site II). Binding sites occupancy was 97% means and 58.5% means for sites I and 11, respectively. The displacement was affected by temperature, membrane protein concentration, freezing, thawing, and dialysis. Other metal cations, including Co++, Fe++, and Mn++, had very little effect on Zn-65 displacement, in contrast copper displaced Zn-65 from its binding sites on RBC membranes. Zinc binding to RBC membranes was rapid and readily reversible in a dynamic equilibrium with its binding sites. It is anticipated that this method will be applicable to studies of a wide variety of diseases specifically related to zinc metabolism in humans as well as in animals. (C) 1994 Wiley-Liss, Inc.

Water regulation of actinomycin-D binding to DNA: the interplay among drug affinity, DNA long-range conformation, and hydration

Actiaomycin-D (actD) binds to natural DNA at two different classes of binding sites, weak and strong. The affinity for these sites is highly dependent on DNA se(sequence and solution conditions, and the interaction appears to be purely entropic driven Although the entropic character of this reaction has been attributed to the release of water molecules upon drug to DNA complex formation, the mechanism by which hydration regulates actD binding and discrimination between different classes of binding sites on natural DNA is still unknown. In this work, we investigate the role of hydration on this reaction using the osmotic stress method. We skew that the decrease of solution water activity, due to the addition of sucrose, glycerol ethylene glycol, and betaine, favors drug binding to the strong binding sites on DNA by increasing both the apparent binding affinity Delta G, and the number of DNA base pairs apparently occupied by the bound drug n(bp/actD). These binding parameters vary linearly with the logarithm of the molar fraction of water in solution log(X-w), which indicates the contribution of water binding to the energetic of the reaction. It is demonstrated that the hydration change measured upon binding increases proportionally to the apparent size of the binding site n(bp/uctD). This indicates that n(bp/actD) measured from the Scatchard plod is a measure of the size of the DNA molecule changing conformation due to ligand binding. We also find that the contribution of DNA deformation, gauged by n(bp/act) to the total free energy of binding Delta G, is given by Delta G = Delta G(local) + n(bp/actD) x delta G(DNA), where Delta G(local), = -8020 +/- 51 cal/mol of actD bound and delta G(DNa) = -24.1 +/- 1.7cal/mol of base pair at 25 degrees C. We interpret Delta G(local), as the energetic contribution due to the direct interactions of actD with the actual tetranucleotide binding site, and it n(bp/actB) X delta G(DNA) as that due to change inconformation, induced by binding, of it n(bp/actD) DNA base pairs flanking the local site. This interpretation is supported by the agreement found between the value of delta G(DNA) and the torsional free energy change measured independently. We conclude suggesting an allosteric model for ligand binding to DNA, such that the increase in binding affinity is achieved by increasing the relaxation of the unfavorable free energy of binding storage at the local site through a larger number of DNA base pairs. The new aspect on this model is that the size of the complex is not fixed but determined by solutions conditions, such as water activity, which modulate the energetic barrier to change helix conformation. These results may suggest that long-range allosteric transitions of duplex DNA are involved in the inhibition of RNA synthesis by actD, and more generally, in the regulation of transcription. (C) 2000 John Wiley & Sons, Inc.

Exploring the binding mechanism of Guaijaverin to human serum albumin: Fluorescence spectroscopy and computational approach

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)