Predicting calcium-binding sites in proteins - a graph theory and geometry approach

Identifying calcium-binding sites in proteins is one of the first steps towards predicting and understanding the role of calcium in biological systems for protein structure and function studies. Due to the complexity and irregularity of calcium-binding sites, a fast and accurate method for predicting and identifying calcium-binding protein is needed. Here we report our development of a new fast algorithm (GG) to detect calcium-binding sites. The GG algorithm uses a graph theory algorithm to find oxygen clusters of the protein and a geometric algorithm to identify the center of these clusters. A cluster of four or more oxygen atoms has a high potential for calcium binding. High performance with about 90% site sensitivity and 80% site selectivity has been obtained for three datasets containing a total of 123 proteins. The results suggest that a sphere of a certain size with four or more oxygen atoms on the surface and without other atoms inside is necessary and sufficient for quickly identifying the majority of the calcium-binding sites with high accuracy. Our finding opens a new avenue to visualize and analyze calcium-binding sites in proteins facilitating the prediction of functions from structural genomic information.     

Structure and dynamics of de novo proteins from a designed superfamily of 4-helix bundles

Libraries of de novo proteins provide an opportunity to explore the structural and functional potential of biological molecules that have not been biased by billions of years of evolutionary selection. Given the enormity of sequence space, a rational approach to library design is likely to yield a higher fraction of folded and functional proteins than a stochastic sampling of random sequences. We previously investigated the potential of library design by binary patterning of hydrophobic and hydrophilic amino acids. The structure of the most stable protein from a binary patterned library of de novo 4-helix bundles was solved previously and shown to be consistent with the design. One structure, however, cannot fully assess the potential of the design strategy, nor can it account for differences in the stabilities of individual proteins. To more fully probe the quality of the library, we now report the NMR structure of a second protein, S-836. Protein S-836 proved to be a 4-helix bundle, consistent with design. The similarity between the two solved structures reinforces previous evidence that binary patterning can encode stable, 4-helix bundles. Despite their global similarities, the two proteins have cores that are packed at different degrees of tightness. The relationship between packing and dynamics was probed using the Modelfree approach, which showed that regions containing a high frequency of chemical exchange coincide with less well-packed side chains. These studies show (1) that binary patterning can drive folding into a particular topology without the explicit design of residue-by-residue packing, and (2) that within a superfamily of binary patterned proteins, the structures and dynamics of individual proteins are modulated by the identity and packing of residues in the hydrophobic core.     

Tertiary templates for the design of diiron proteins.

Diiron proteins represent a diverse class of structures involved in the binding and activation of oxygen. This review explores the simple structural features underlying the common metal-ion-binding and oxygen-binding properties of these proteins. The backbone geometries of their active sites are formed by four-helix bundles, which may be parameterized to within approximately 1 A root mean square deviation. Such parametric models are excellent starting points for investigating how asymmetric deviations from an idealized geometry influence the functional properties of the metal ion centers. These idealized models also provide attractive frameworks for de novo protein design.     

Hidden partners: Using cross‐docking calculations to predict binding sites for proteins with multiple interactions

Protein‐protein interactions control a large range of biological processes and their identification is essential to understand the underlying biological mechanisms. To complement experimental approaches, in silico methods are available to investigate protein‐protein interactions. Cross‐docking methods, in particular, can be used to predict protein binding sites. However, proteins can interact with numerous partners and can present multiple binding sites on their surface, which may alter the binding site prediction quality. We evaluate the binding site predictions obtained using complete cross‐docking simulations of 358 proteins with 2 different scoring schemes accounting for multiple binding sites. Despite overall good binding site prediction performances, 68 cases were still associated with very low prediction quality, presenting individual area under the specificity‐sensitivity ROC curve (AUC) values below the random AUC threshold of 0.5, since cross‐docking calculations can lead to the identification of alternate protein binding sites (that are different from the reference experimental sites). For the large majority of these proteins, we show that the predicted alternate binding sites correspond to interaction sites with hidden partners, that is, partners not included in the original cross‐docking dataset. Among those new partners, we find proteins, but also nucleic acid molecules. Finally, for proteins with multiple binding sites on their surface, we investigated the structural determinants associated with the binding sites the most targeted by the docking partners.     

Crystal structure of calcium dodecin (Rv0379), from Mycobacterium tuberculosis with a unique calcium-binding site

In eukaryotes, calcium-binding proteins play a pivotal role in diverse cellular processes, and recent findings suggest similar roles for bacterial proteins at different stages in their life cycle. Here, we report the crystal structure of calcium dodecin, Rv0379, from Mycobacterium tuberculosis with a dodecameric oligomeric assembly and a unique calcium-binding motif. Structure and sequence analysis were used to identify orthologs of Rv0379 with different ligand-binding specificity.     

De novo design of α,β-didehydrophenylalanine containing peptides: from models to applications

The de novo design of peptides and proteins has emerged as an approach for investigating protein structure and function. The success relies heavily on the ability to design relatively short peptides that can adopt stable secondary structures. To this end, substitution with α,β-dehydroamino acids, especially α,β-didehydrophenylalanine (ΔPhe or ΔF) has blossomed in manifold directions, providing a rich diversity of well-defined structural motifs. Introduction of α,β-didehydrophenylalanine induces β-bends in small and 3(10)-helices in longer peptide sequences. Most favorable conformation of ΔF residues are (φ,ψ) ～(60°, 30°), (-60°, -30°), (-60°, 150°), and (60°, -150°). These features have been exploited in designing helix-turn-helix, helical bundle arrangements, and glycine zipper type super secondary structural motifs. The unusual capability of α,β-didehydrophenylalanine ring to form a variety of multicentered interactions (N-H…O, C-H…O, C-H…π, and N-H…π) suggests its possible exploitation for future de novo design of supramolecular structures. This work has now been extended to the de novo design of peptides with antibiotic, antifibrillization activity, etc. More recently, self-assembling properties of small dehydropeptides have been explored. This review focuses primarily on the structural and functional behavior of α,β-didehydrophenylalanine containing peptides.     

Evolution of binding sites for zinc and calcium ions playing structural roles

The geometry of metal coordination by proteins is well understood, but the evolution of metal binding sites has been less studied. Here we present a study on a small number of well-documented structural calcium and zinc binding sites, concerning how the geometry diverges between relatives, how often nonrelatives converge towards the same structure, and how often these metal binding sites are lost in the course of evolution. Both calcium and zinc binding site structure is observed to be conserved; structural differences between those atoms directly involved in metal binding in related proteins are typically less than 0.5 A root mean square deviation, even in distant relatives. Structural templates representing these conserved calcium and zinc binding sites were used to search the Protein Data Bank for cases where unrelated proteins have converged upon the same residue selection and geometry for metal binding. This allowed us to identify six "archetypal" metal binding site structures: two archetypal zinc binding sites, both of which had independently evolved on a large number of occasions, and four diverse archetypal calcium binding sites, where each had evolved independently on only a handful of occasions. We found that it was common for distant relatives of metal-binding proteins to lack metal-binding capacity. This occurred for 13 of the 18 metal binding sites we studied, even though in some of these cases the original metal had been classified as "essential for protein folding." For most of the calcium binding sites studied (seven out of eleven cases), the lack of metal binding in relatives was due to point mutation of the metal-binding residues, whilst for zinc binding sites, lack of metal binding in relatives always involved more extensive changes, with loss of secondary structural elements or loops around the binding site.     

EF-hand calcium-binding proteins

The EF-hand motif is the most common calcium-binding motif found in proteins. Several high-resolution structures containing different metal ions bound to EF-hand sites have given new insight into the modulation of their binding affinities. Recently determined structures of members of several newly identified protein families that contain the EF-hand motif in some of their domains, as well as of their complexes with target molecules, are throwing light on the surprising variety of functions that can be served by this simple and ingenious structural motif.     

Design principles for chlorophyll-binding sites in helical proteins

The cyclic tetrapyrroles, viz. chlorophylls (Chl), their bacterial analogs bacteriochlorophylls, and hemes are ubiquitous cofactors of biological catalysis that are involved in a multitude of reactions. One systematic approach for understanding how Nature achieves functional diversity with only this handful of cofactors is by designing de novo simple and robust protein scaffolds with heme and/or (bacterio)chlorophyll [(B)Chls]-binding sites. This strategy is currently mostly implemented for heme-binding proteins. To gain more insight into the factors that determine heme-/(B)Chl-binding selectivity, we explored the geometric parameters of (B)Chl-binding sites in a nonredundant subset of natural (B)Chl protein structures. Comparing our analysis to the study of a nonredundant database of heme-binding helical histidines by Negron et al. (Proteins 2009;74:400-416), we found a preference for the m-rotamer in (B)Chl-binding helical histidines, in contrast to the preferred t-rotamer in heme-binding helical histidines. This may be used for the design of specific heme- or (B)Chl-binding sites in water-soluble helical bundles, because the rotamer type defines the positioning of the bound cofactor with respect to the helix interface and thus the protein-binding site. Consensus sequences for (B)Chl binding were identified by combining a computational and database-derived approach and shown to be significantly different from the consensus sequences recommended by Negron et al. (Proteins 2009;74:400-416) for heme-binding helical proteins. The insights gained in this work on helix- (B)Chls-binding pockets provide useful guidelines for the construction of reasonable (B)Chl-binding protein templates that can be optimized by computational tools.     

Minimal functional sites allow a classification of zinc sites in proteins

Zinc is indispensable to all forms of life as it is an essential component of many different proteins involved in a wide range of biological processes. Not differently from other metals, zinc in proteins can play different roles that depend on the features of the metal-binding site. In this work, we describe zinc sites in proteins with known structure by means of three-dimensional templates that can be automatically extracted from PDB files and consist of the protein structure around the metal, including the zinc ligands and the residues in close spatial proximity to the ligands. This definition is devised to intrinsically capture the features of the local protein environment that can affect metal function, and corresponds to what we call a minimal functional site (MFS). We used MFSs to classify all zinc sites whose structures are available in the PDB and combined this classification with functional annotation as available in the literature. We classified 77% of zinc sites into ten clusters, each grouping zinc sites with structures that are highly similar, and an additional 16% into seven pseudo-clusters, each grouping zinc sites with structures that are only broadly similar. Sites where zinc plays a structural role are predominant in eight clusters and in two pseudo-clusters, while sites where zinc plays a catalytic role are predominant in two clusters and in five pseudo-clusters. We also analyzed the amino acid composition of the coordination sphere of zinc as a function of its role in the protein, highlighting trends and exceptions. In a period when the number of known zinc proteins is expected to grow further with the increasing awareness of the cellular mechanisms of zinc homeostasis, this classification represents a valuable basis for structure-function studies of zinc proteins, with broad applications in biochemistry, molecular pharmacology and de novo protein design.     

De novo heme proteins from designed combinatorial libraries.

We previously reported the design of a library of de novo amino acid sequences targeted to fold into four-helix bundles. The design of these sequences was based on a "binary code" strategy, in which the patterning of polar and nonpolar amino acids is specified explicitly, but the exact identities of the side chains is varied extensively (Kamtekar S, Schiffer JM, Xiong H, Babik JM, Hecht MH, 1993, Science 262:1680-1685). Because of this variability, the resulting collection of amino acid sequences may include de novo proteins capable of binding biologically important cofactors. To probe for such binding, the de novo sequences were screened for their ability to bind the heme cofactor. Among an initial collection of 30 binary code sequences, 15 are shown to bind heme and form bright red complexes. Characterization of several of these de novo heme proteins demonstrated that their absorption spectra and resonance Raman spectra resemble those of natural cytochromes. Because the design of these sequences is based on global features of polar/ nonpolar patterning, the finding that half of them bind heme highlights the power of the binary code strategy, and demonstrates that isolating de novo heme proteins does not require explicit design of the cofactor binding site. Because bound heme plays a key role in the functions of many natural proteins, these results suggest that binary code sequences may serve as initial prototypes for the development of large collections of functionally active de novo proteins.     

Stably folded de novo proteins from a designed combinatorial library

Binary patterning of polar and nonpolar amino acids has been used as the key design feature for constructing large combinatorial libraries of de novo proteins. Each position in a binary patterned sequence is designed explicitly to be either polar or nonpolar; however, the precise identities of these amino acids are varied extensively. The combinatorial underpinnings of the "binary code" strategy preclude explicit design of particular side chains at specified positions. Therefore, packing interactions cannot be specified a priori. To assess whether the binary code strategy can nonetheless produce well-folded de novo proteins, we constructed a second-generation library based upon a new structural scaffold designed to fold into 102-residue four-helix bundles. Characterization of five proteins chosen arbitrarily from this new library revealed that (1) all are alpha-helical and quite stable; (2) four of the five contain an abundance of tertiary interactions indicative of well-ordered structures; and (3) one protein forms a well-folded structure with native-like features. The proteins from this new 102-residue library are substantially more stable and dramatically more native-like than those from an earlier binary patterned library of 74-residue sequences. These findings demonstrate that chain length is a crucial determinant of structural order in libraries of de novo four-helix bundles. Moreover, these results show that the binary code strategy--if applied to an appropriately designed structural scaffold--can generate large collections of stably folded and/or native-like proteins.     

Energy minimization method using automata network for sequence and side-chain conformation prediction from given backbone geometry

Globular proteins have high packing densities as a result of residue side chains in the core achieving a tight, complementary packing. The internal packing is considered the main determinant of native protein structure. From that point of view, we present here a method of energy minimization using an automata network to predict a set of amino acid sequences and their side-chain conformations from a desired backbone geometry for de novo design of proteins. Using discrete side-chain conformations, that is, rotamers, the sequence generation problem from a given backbone geometry becomes one of combinatorial problems. We focused on the residues composing the interior core region and predicted a set of amino acid Sequences and their side-chain conformations only from a given backbone geometry. The kinds of residues were restricted to six hydrophobic amino acids (Ala, Ile, Met, Leu, Phe, and Val) because the core regions are almost always composed of hydrophobic residues. The obtained sequences were well packed as was the native sequence. The method can be used for automated sequence generation in the de novo design of proteins. © 1994 Wiley-Liss, Inc.     

Restrained least squares refinement of native (calcium) and cadmium-substituted carp parvalbumin using X-ray crystallographic data at 1.6-A resolution

Carp parvalbumin coordinates calcium through one carbonyl oxygen atom and the oxygen-containing side chains of 5 amino acid residues, or 4 residues and a water molecule, in a helix-loop-helix structural motif. Other calcium-binding proteins, including calmodulin and troponin C, also possess this unique calcium-binding design, which is designated EF-hand or calmodulin fold. Parvalbumin has two such sites, labeled CD and EF. Each of the calcium-binding sites of refined structures of proteins belonging to this group has a 7-oxygen coordination sphere except those of the structure of parvalbumin as it was reported in 1975. This structure had been refined at 1.9 A using difference Fourier techniques on film data. The CD site appeared to be 6-coordinate and the EF site 8-coordinate. Results of NMR experiments using 113Cd-substituted parvalbumin, however, indicate that the sites are similar to one another with coordination number greater than 6. To resolve the inconsistency between crystallographic and NMR results, 1.6 A area detector data was collected for native and cadmium-substituted parvalbumin; the structures have been refined to R factors of 18.7% and 16.4%, respectively, with acceptable geometry and low errors in atomic coordinates. Differences between the parvalbumin structure described in 1975 and the present structure are addressed, including the discovery of 7-coordination for both the CD and EF sites.     

Q-SiteFinder: an energy-based method for the prediction of protein-ligand binding sites

MOTIVATION: Identifying the location of ligand binding sites on a protein is of fundamental importance for a range of applications including molecular docking, de novo drug design and structural identification and comparison of functional sites. Here, we describe a new method of ligand binding site prediction called Q-SiteFinder. It uses the interaction energy between the protein and a simple van der Waals probe to locate energetically favourable binding sites. Energetically favourable probe sites are clustered according to their spatial proximity and clusters are then ranked according to the sum of interaction energies for sites within each cluster. RESULTS: There is at least one successful prediction in the top three predicted sites in 90% of proteins tested when using Q-SiteFinder. This success rate is higher than that of a commonly used pocket detection algorithm (Pocket-Finder) which uses geometric criteria. Additionally, Q-SiteFinder is twice as effective as Pocket-Finder in generating predicted sites that map accurately onto ligand coordinates. It also generates predicted sites with the lowest average volumes of the methods examined in this study. Unlike pocket detection, the volumes of the predicted sites appear to show relatively low dependence on protein volume and are similar in volume to the ligands they contain. Restricting the size of the pocket is important for reducing the search space required for docking and de novo drug design or site comparison. The method can be applied in structural genomics studies where protein binding sites remain uncharacterized since the 86% success rate for unbound proteins appears to be only slightly lower than that of ligand-bound proteins. AVAILABILITY: Both Q-SiteFinder and Pocket-Finder have been made available online at http://www.bioinformatics.leeds.ac.uk/qsitefinder and http://www.bioinformatics.leeds.ac.uk/pocketfinder     

Geometric constraints for porphyrin binding in helical protein binding sites

Helical bundles which bind heme and porphyrin cofactors have been popular targets for cofactor-containing de novo protein design. By analyzing a highly nonredundant subset of the protein databank we have determined a rotamer distribution for helical histidines bound to heme cofactors. Analysis of the entire nonredundant database for helical sequence preferences near the ligand histidine demonstrated little preference for amino acid side chain identity, size, or charge. Analysis of the database subdivided by ligand histidine rotamer, however, reveals strong preferences in each case, and computational modeling illuminates the structural basis for some of these findings. The majority of the rotamer distribution matches that predicted by molecular simulation of a single porphyrin-bound histidine residue placed in the center of an all-alanine helix, and the deviations explain two prominent features of natural heme protein binding sites: heme distortion in the case of the cytochromes C in the m166 histidine rotamer, and a highly prevalent glycine residue in the t73 histidine rotamer. These preferences permit derivation of helical consensus sequence templates which predict optimal side chain-cofactor packing interactions for each rotamer. These findings thus promise to guide future design endeavors not only in the creation of higher affinity heme and porphyrin binding sites, but also in the direction of bound cofactor geometry.     

Surface grafting onto template-assembled synthetic protein scaffolds in molecular recognition

Creating functional biological molecules de novo requires a detailed understanding of the intimate relationship between primary sequence, folding mechanism, and packing topology, and remains up to now a most challenging goal in protein design and mimicry. As a consequence, the use of well-defined robust macromolecules as scaffolds for the introduction of function by grafting surface residues has become a major objective in protein engineering and de novo design. In this article, the concept of scaffolds is demonstrated on some selected examples, illustrating that novel types of functional molecules can be generated. Reengineered proteins and, most notably, de novo designed peptide scaffolds exhibiting molecular function, are ideal tools for structure-function studies and as leads in drug design.     

Heme proteins—Diversity in structural characteristics,function, and folding

The characteristics of heme prosthetic groups and their binding sites have been analyzed in detail in a data set of nonhomologous heme proteins. Variations in the shape, volume, and chemical composition of the binding site, in the mode of heme binding and in the number and nature of heme–protein interactions are found to result in significantly different heme environments in proteins with different functions in biology. Differences are also seen in the properties of the apo states of the proteins. The apo states of proteins that bind heme permanently in their functional form show some disorder, ranging from local unfolding in the heme binding pocket to complete unfolding to give a random coil. In contrast, proteins that bind heme transiently are fully folded in their apo and holo states, presumably allowing both apo and holo forms to remain biologically active resisting aggregation or proteolysis. The principles identified here provide a framework for the design of de novo proteins that will exhibit tight heme ligand binding and for the identification of the function of structural genomic target proteins with heme ligands. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

OptGraft: A computational procedure for transferring a binding site onto an existing protein scaffold

One of the many challenging tasks of protein design is the introduction of a completely new function into an existing protein scaffold. In this study, we introduce a new computational procedure OptGraft for placing a novel binding pocket onto a protein structure so as its geometry is minimally perturbed. This is accomplished by introducing a two‐level procedure where we first identify where are the most appropriate locations to graft the new binding pocket into the protein fold by minimizing the departure from a set of geometric restraints using mixed‐integer linear optimization. On identifying the suitable locations that can accommodate the new binding pocket, CHARMM energy calculations are employed to identify what mutations in the neighboring residues, if any, are needed to ensure that the minimum energy conformation of the binding pocket conserves the desired geometry. This computational framework is benchmarked against the results available in the literature for engineering a copper binding site into thioredoxin protein. Subsequently, OptGraft is used to guide the transfer of a calcium‐binding pocket from thermitase protein (PDB: 1thm) into the first domain of CD2 protein (PDB:1hng). Experimental characterization of three de novo redesigned proteins with grafted calcium‐binding centers demonstrated that they all exhibit high affinities for terbium (K_d ～ 22, 38, and 55 μM) and can selectively bind calcium over magnesium.