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.     

Residue-level prediction of DNA-binding sites and its application on DNA-binding protein predictions

Protein-DNA interactions are crucial to many cellular activities such as expression-control and DNA-repair. These interactions between amino acids and nucleotides are highly specific and any aberrance at the binding site can render the interaction completely incompetent. In this study, we have three aims focusing on DNA-binding residues on the protein surface: to develop an automated approach for fast and reliable recognition of DNA-binding sites; to improve the prediction by distance-dependent refinement; use these predictions to identify DNA-binding proteins. We use a support vector machines (SVM)-based approach to harness the features of the DNA-binding residues to distinguish them from non-binding residues. Features used for distinction include the residue's identity, charge, solvent accessibility, average potential, the secondary structure it is embedded in, neighboring residues, and location in a cationic patch. These features collected from 50 proteins are used to train SVM. Testing is then performed on another set of 37 proteins, much larger than any testing set used in previous studies. The testing set has no more than 20% sequence identity not only among its pairs, but also with the proteins in the training set, thus removing any undesired redundancy due to homology. This set also has proteins with an unseen DNA-binding structural class not present in the training set. With the above features, an accuracy of 66% with balanced sensitivity and specificity is achieved without relying on homology or evolutionary information. We then develop a post-processing scheme to improve the prediction using the relative location of the predicted residues. Balanced success is then achieved with average sensitivity, specificity and accuracy pegged at 71.3%, 69.3% and 70.5%, respectively. Average net prediction is also around 70%. Finally, we show that the number of predicted DNA-binding residues can be used to differentiate DNA-binding proteins from non-DNA-binding proteins with an accuracy of 78%. Results presented here demonstrate that machine-learning can be applied to automated identification of DNA-binding residues and that the success rate can be ameliorated as more features are added. Such functional site prediction protocols can be useful in guiding consequent works such as site-directed mutagenesis and macromolecular docking.     

Three-Dimensional Models of Immune Cell Surface Proteins and Identification of Binding Sites

The extracellular regions of many cell surface proteins of the immune system contain distinct domains that may be linked in many different ways and are often only loosely tethered to the transmembrane segment. In efforts to identify regions critical for binding, molecular models of these domains are used to select residues for mutagenesis and to map binding sites. Many immune cell surface proteins belong to protein superfamilies and display only limited sequence identity compared to proteins of known three-dimensional (3D) structure, often 30% or less. Therefore, detailed 3D structures are difficult to predict, and structure-based sequence analysis and model assessment are particularly important components of the model building process. In some cases, experimentally determined structures have made it possible to assess the accuracy of predictions, which illustrates the opportunities and shortcomings of the approach. Herein the model-based identification of binding sites in cell surface proteins is described and representative examples are discussed.     

Nucleotide sequence of the gene for alkaline phosphatase of Thermus caldophilus GK24 and characteristics of the deduced primary structure of the enzyme

The gene encoding Thermus caldophilus GK24 (Tca) alkaline phosphatase was cloned into Escherichia coli. The primary structure of Tca alkaline phosphatase was deduced from its nucleotide sequence. The Tca alkaline phosphatase precursor, including the signal peptide sequence, was comprised of 501 amino acid residues. Its molecular mass was determined to be 54? omitted?760 Da. On the alignment of the amino acid sequence, Tca alkaline phosphatase showed sequence homology with the microbial alkaline phosphatases, 20% identity with E. coli alkaline phosphatase and 22% Bacillus subtilis (Bsu) alkaline phosphatases. High sequence identity was observed in the regions containing the Ser-102 residue of the active site, the zinc and magnesium binding sites of E. coli alkaline phosphatase. Comparison of Tca alkaline phosphatase and E. coli alkaline phosphatase structures suggests that the reduced activity of the Tca alkaline phosphatase, in the presence of zinc, is directly involved in some of the different metal binding sites. Heat-stable Tca alkaline phosphatase activity was detected in E. coli YK537, harboring pJRAP.     

Zinc coordination sphere in biochemical zinc sites

Zinc is known to be indispensable to growth and development and transmission of the genetic message. It does this through a remarkable mosaic of zinc binding motifs that orchestrate all aspects of metabolism. There are now nearly 200 three dimensional structures for zinc proteins, representing all six classes of enzymes and covering a wide range of phyla and species. These structures provide standards of reference for the identity and nature of zinc ligands in other proteins for which only the primary structure is known. Three primary types of zinc sites are apparent from examination of these structures: structural, catalytic and cocatalytic. The most common amino acids that supply ligands to these sites are His, Glu, Asp and Cys. In catalytic sites zinc generally forms complexes with water and any three nitrogen, oxygen and sulfur donors with His being the predominant amino acid chosen. Water is always a ligand to such sites. Structural zinc sites have four protein ligands and no bound water molecule. Cys is the preferred ligand in such sites. Cocatalytic sites contain two or three metals in close proximity with two of the metals bridged by a side chain moiety of a single amino acid residue, such as Asp, Glu or His and sometimes a water molecule. Asp and His are the preferred amino acids for these sites. No Cys ligands are found in such sites. The scaffolding of the zinc sites is also important to the function and reactivity of the bound metal. The influence of zinc on quaternary protein structure has led to the identification of a fourth type of zinc binding site, protein inteface. In this case zinc sites are formed from ligands supplied from amino acid residues residing in the binding surface of two proteins. The resulting zinc site usually has the coordination properties of a catalytic or structural zinc binding site.     

Efficient recognition of protein fold at low sequence identity by conservative application of Psi-BLAST: application

Based on a study involving structural comparisons of proteins sharing 25% or less sequence identity, three rounds of Psi-BLAST appear capable of identifying remote evolutionary homologs with greater than 95% confidence provided that more than 50% of the query sequence can be aligned with the target sequence. Since it seems that more than 80% of all homologous protein pairs may be characterized by a lack of significant sequence similarity, the experimental biologist is often confronted with a lack of guidance from conventional homology searches involving pair-wise sequence comparisons. The ability to disregard levels of sequence identity and expect value in Psi-BLAST if at least 50% of the query sequence has been aligned allows for generation of new hypotheses by consideration of matches that are conventionally disregarded. In one example, we suggest a possible evolutionary linkage between the cupredoxin and immunoglobulin fold families. A thermostable hypothetical protein of unknown function may be a circularly permuted homolog to phosphotriesterase, an enzyme capable of detoxifying organophosphate nerve agents. In a third example, the amino acid sequence of another hypothetical protein of unknown function reveals the ATP binding-site, metal binding site, and catalytic sidechain consistent with kinase activity of unknown specificity. This approach significantly expands the utility of existing sequence data to define the primary structure degeneracy of binding sites for substrates, cofactors and other proteins.     

Genome‐wide computational determination of the human metalloproteome

Accurate prediction of protein function in humans is important for understanding biological processes at the molecular level in biomedicine and drug design. Over a third of proteins are commonly held to bind metal, and ～10% of human proteins, to bind zinc. Therefore, an initial step in protein function prediction frequently involves predicting metal ion binding. In recent years, methods have been developed to predict a set of residues in 3D space forming the metal‐ion binding site, often with a high degree of accuracy. Here, using extensions of these methods, we provide an extensive list of human proteins and their putative metal ion binding site residues, using translated gene sequences derived from the complete, resolved human genome. Under conditions of ～90% selectivity, over 900 new human putative metal ion binding proteins are identified. A statistical analysis of resolved metal ion binding sites in the human metalloproteome is furnished and the importance of remote homology analysis is demonstrated. As an example, a novel metal‐ion binding site involving a complex of a botulinum substrate with its inhibitor is presented. On the basis of the location of the predicted site and the interactions of the contacting residues at the complex interface, we postulate that metal ion binding in this region could influence complex formation and, consequently, the functioning of the protein. Thus, this work provides testable hypotheses about novel functions of known proteins. Proteins 2015; 83:931–939. © 2015 Wiley Periodicals, Inc.     

Primary structure of cobra complement component C3.

Complement component C3 is a multifunctional protein known to interact specifically with more than 10 different plasma proteins or cell surface receptors. Cobra venom contains cobra venom factor, a structural analogue of C3 that shares some properties with C3 (e.g., formation of a C3/C5 convertase) but differs in others (e.g., susceptibility to regulation by factors H and I). The elucidation of structural differences between C3 and cobra venom factor can be expected to help identify functionally important regions of C3 molecules. To that end we have undertaken the molecular cloning of both cobra C3 and cobra venom factor to take advantage of the unique biologic system where both proteins are produced by the same species. We report the primary structure of cobra C3 mRNA and the derived protein structure. Cobra C3 mRNA is 5211 bp in length. It contains an open reading frame of 4953 bp coding for a single pre-pro-C3 molecule, consisting of a 22-amino acid signal sequence, a 633-amino acid beta-chain (70 kDa), and a 992-amino acid alpha-chain (112 kDa) which is separated from the beta-chain by four arginine residues. There are no N-glycosylation sites in cobra C3. Cobra C3 exhibits approximately 58% nucleotide sequence identity with C3 from mammalian species. At the protein level, sequence identity is approximately 52% and sequence similarity approximately 71%. All 27 cysteine residues are highly conserved as are the C3 convertase cleavage site, the thioester site, and the factor B binding site. Cobra C3 also seems to have homologous binding sites for factor H and properdin, as well as a conserved sequence in the functionally important region of the C3a anaphylatoxin. The sequence homology at the CR2 and CR3 binding sites does not exceed the overall sequence homology. Accordingly, the existence of CR2 and CR3 binding sites can neither be deduced nor excluded.     

Identification of common structural features of binding sites in galactose-specific proteins

Galactose-binding proteins characterize an important subgroup of sugar-binding proteins that are involved in a variety of biological processes. Structural studies have shown that the Gal-specific proteins encompass a diverse range of primary and tertiary structures. The binding sites for galactose also seem to vary in different protein-galactose complexes. No common binding site features that are shared by the Gal-specific proteins to achieve ligand specificity are so far known. With the assumption that common recognition principles will exist for common substrate recognition, the present study was undertaken to identify and characterize any unique galactose-binding site signature by analyzing the three-dimensional (3D) structures of 18 protein-galactose complexes. These proteins belong to 7 nonhomologous families; thus, there is no sequence or structural similarity across the families. Within each family, the binding site residues and their relative distances were well conserved, but there were no similarities across families. A novel, yet simple, approach was adopted to characterize the binding site residues by representing their relative spatial dispositions in polar coordinates. A combination of the deduced geometrical features with the structural characteristics, such as solvent accessibility and secondary structure type, furnished a potential galactose-binding site signature. The signature was evaluated by incorporation into the program COTRAN to search for potential galactose-binding sites in proteins that share the same fold as the known galactose-binding proteins. COTRAN is able to detect galactose-binding sites with a very high specificity and sensitivity. The deduced galactose-binding site signature is strongly validated and can be used to search for galactose-binding sites in proteins. PROSITE-type signature sequences have also been inferred for galectin and C-type animal lectin-like fold families of Gal-binding proteins.     

Prediction of zinc binding sites in proteins using sequence derived information

Zinc is one the most abundant catalytic cofactor and also an important structural component of a large number of metallo-proteins. Hence prediction of zinc metal binding sites in proteins can be a significant step in annotation of molecular function of a large number of proteins. Majority of existing methods for zinc-binding site predictions are based on a data-set of proteins, which has been compiled nearly a decade ago. Hence there is a need to develop zinc-binding site prediction system using the current updated data to include recently added proteins. Herein, we propose a support vector machine-based method, named as ZincBinder, for prediction of zinc metal-binding site in a protein using sequence profile information. The predictor was trained using fivefold cross validation approach and achieved 85.37% sensitivity with 86.20% specificity during training. Benchmarking on an independent non-redundant data-set, which was not used during training, showed better performance of ZincBinder vis-à-vis existing methods. Executable versions, source code, sample datasets, and usage instructions are available at http://proteininformatics.org/mkumar/znbinder/     

Structure and stability of an unusual zinc-binding protein from Bacteroides thetaiotaomicron

The crystal structure of a putative protease from Bacteroides thetaiotaomicron (ppBat) suggested the presence of a zinc ion in each protomer of the dimer as well as a flavin in the dimer interface. Since the chemical identity of the flavin and the exact mode of binding remained unclear, we have determined the crystal structure of ppBat in complex with riboflavin. The obtained structure revealed that the isoalloxazine ring is sandwiched between two tryptophan residues (Trp164) from both chains and adopts two alternate orientations with the N(10)-ribityl side chain protruding from the binding site in opposite directions. In order to characterize the zinc-binding site, we generated two single variants and one double variant in which the two coordinating cysteine residues (Cys74 and Cys111) were replaced by alanine. All three variants were unable to bind zinc demonstrating that both cysteine residues are essential for binding. Moreover, the lack of zinc binding also resulted in drastically reduced thermal stability (11–15 °C). A similar effect was obtained when wild-type protein was incubated with EDTA supporting the conclusion that the zinc-binding site plays an important structural role in ppBat. On the other hand, attempts to identify proteolytic activity failed suggesting that the zinc may not act as a catalytic center in ppBat. Structurally similar zinc binding motives in other proteins were also found to play a structural rather than catalytic role and hence it appears that neither the flavin nor the zinc binding sites possess a catalytic function in ppBat.     

Remote homologue identification of Drosophila GAGA factor in mouse

GAGA factor (GAF) is involved in both gene activation and gene repression and plays a role in the modulation of chromatin structure. In Drosophila, Trithroax like (Trl) gene encodes the DNA binding protein called GAGA factor (GAF). Trl-GAF binds to GAGA sites through its C2H2 zinc finger domain and has an N-terminal BTB/POZ domain. Identification of Trl-GAF homologue in mouse helps in deeper understanding of the mechanism and function. Conventional alignment tools such as BLAST and FASTA cannot identify homologues in mouse genome as their sequence identity is below 30%. In the present study, various sequence and structure analyses were followed for the detection of remote homologues of Drosophila GAGA FACTOR in mouse to identify as Zbtb3. Through homology modeling and docking approach, the zinc finger region of mouse Zbtb3 showed conserved residues and favorable DNA binding sites with GAGA sites similar to that of Drosophila GAGA FACTOR.     

Solution structure of a CCHHC domain of neural zinc finger factor-1 and its implications for DNA binding

The structure of a CCHHC zinc-binding domain from neural zinc finger factor-1 (NZF-1) has been determined in solution though the use of NMR methods. This domain is a member of a family of domains that have the Cys-X(4)-Cys-X(4)-His-X(7)-His-X(5)-Cys consensus sequence. The structure determination reveals a novel fold based around a zinc(II) ion coordinated to three Cys residues and the second of the two conserved His residues. The other His residue is stacked between the metal-coordinated His residue and a relatively conserved aromatic residue. Analysis of His to Gln sequence variants reveals that both His residues are required for the formation of a well-defined structure, but neither is required for high-affinity metal binding at a tetrahedral site. The structure suggests that a two-domain protein fragment and a double-stranded DNA binding site may interact with a common two-fold axis relating the two domains and the two half-sites of the DNA-inverted repeat.     

Crystal structure of the zinc-binding transport protein ZnuA from Escherichia coli reveals an unexpected variation in metal coordination

Bacterial ATP-binding cassette transport systems for high-affinity uptake of zinc and manganese use a cluster 9 solute-binding protein. Structures of four cluster 9 transport proteins have been determined previously. However, the structural determinants for discrimination between zinc and manganese remain under discussion. To further investigate the variability of metal binding sites in bacterial transporters, we have determined the structure of the zinc-bound transport protein ZnuA from Escherichia coli to 1.75 A resolution. The overall structure of ZnuA is similar to other solute-binding transporters. A scaffolding alpha-helix forms the backbone for two structurally related globular domains. The metal-binding site is located at the domain interface. The bound zinc ion is coordinated by three histidine residues (His78, His161 and His225) and one glutamate residue (Glu77). The functional role of Glu77 for metal binding is unexpected, because this residue is not conserved in previously determined structures of zinc and manganese-specific transport proteins. The observed metal coordination by four protein residues differs significantly from the zinc-binding site in the ZnuA transporter from Synechocystis 6803, which binds zinc via three histidine residues. In addition, the E. coli ZnuA structure reveals the presence of a disulfide bond in the C-terminal globular domain that is not present in previously determined cluster 9 transport protein structures.     

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.     

Identification of the Zn2+ Binding Site and Mode of Operation of a Mammalian Zn2+ Transporter

Vesicular zinc transporters (ZnTs) play a critical role in regulating Zn²⁺ homeostasis in various cellular compartments and are linked to major diseases ranging from Alzheimer disease to diabetes. Despite their importance, the intracellular localization of ZnTs poses a major challenge for establishing the mechanisms by which they function and the identity of their ion binding sites. Here, we combine fluorescence-based functional analysis and structural modeling aimed at elucidating these functional aspects. Expression of ZnT5 was followed by both accelerated removal of Zn²⁺ from the cytoplasm and its increased vesicular sequestration. Further, activity of this zinc transport was coupled to alkalinization of the trans-Golgi network. Finally, structural modeling of ZnT5, based on the x-ray structure of the bacterial metal transporter YiiP, identified four residues that can potentially form the zinc binding site on ZnT5. Consistent with this model, replacement of these residues, Asp⁵⁹⁹ and His⁴⁵¹, with alanine was sufficient to block Zn²⁺ transport. These findings indicate, for the first time, that Zn²⁺ transport mediated by a mammalian ZnT is catalyzed by H⁺/Zn²⁺ exchange and identify the zinc binding site of ZnT proteins essential for zinc transport.     

Prediction of metal ion-binding sites in proteins using the fragment transformation method

The structure of a protein determines its function and its interactions with other factors. Regions of proteins that interact with ligands, substrates, and/or other proteins, tend to be conserved both in sequence and structure, and the residues involved are usually in close spatial proximity. More than 70,000 protein structures are currently found in the Protein Data Bank, and approximately one-third contain metal ions essential for function. Identifying and characterizing metal ion-binding sites experimentally is time-consuming and costly. Many computational methods have been developed to identify metal ion-binding sites, and most use only sequence information. For the work reported herein, we developed a method that uses sequence and structural information to predict the residues in metal ion-binding sites. Six types of metal ion-binding templates- those involving Ca(2+), Cu(2+), Fe(3+), Mg(2+), Mn(2+), and Zn(2+)-were constructed using the residues within 3.5 ? of the center of the metal ion. Using the fragment transformation method, we then compared known metal ion-binding sites with the templates to assess the accuracy of our method. Our method achieved an overall 94.6 % accuracy with a true positive rate of 60.5 % at a 5 % false positive rate and therefore constitutes a significant improvement in metal-binding site prediction.     

Identification of metal binding residues for the binuclear zinc phosphodiesterase reveals identical coordination as glyoxalase II

Zinc phosphodiesterase (ZiPD) is a member of the metallo-beta-lactamase family with a binuclear zinc binding site. As an experimental attempt to identify the metal ligands of Escherichia coli ZiPD and to investigate their function in catalysis, we mutationally exchanged candidate metal coordinating residues and performed kinetic and X-ray absorption spectroscopic analysis of the mutant proteins. All mutants (H66E, H69A, H141A, D212A, D212C, H231A, H248A, and H270A) show significantly lower catalytic rates toward the substrate bis(p-nitrophenyl)phosphate. Substrate binding, represented by the kinetic value K', remains unchanged for six mutants, whereas it is increased 3-4-fold for H231A and H270A. Accordingly, these two residues are supposed to be involved in substrate binding, whereas the others are more important for catalytic turnover and thus are assumed to be involved in zinc ligation. Structural insight into the metal binding of D212 was gained by zinc K-edge extended X-ray absorption fine structure (EXAFS). The sulfur coordination number of the cysteine mutant was found to be 1, demonstrating binding to both zinc metals in a bridging mode. Taken together with two residues from a strictly conserved sequence region within the metallo-beta-lactamase family, the metal site of ZiPD is proposed with H64, H66, and H141 coordinating ZnA, D68, H69, and H248 coordinating ZnB, and D212 bridging both metals. Surprisingly, the same coordination sphere is found in glyoxalase II. This is further substantiated by comparable EXAFS spectra of both native enzymes. This is the first example of the same metal site in two members of the metallo-beta-lactamase domain proteins catalyzing different reactions. The kinetic analysis of mutants provides unexpected insights into the reaction mechanism of ZiPD.     

A Structure-Based Approach for Detection of Thiol Oxidoreductases and Their Catalytic Redox-Active Cysteine Residues

Cysteine (Cys) residues often play critical roles in proteins, for example, in the formation of structural disulfide bonds, metal binding, targeting proteins to the membranes, and various catalytic functions. However, the structural determinants for various Cys functions are not clear. Thiol oxidoreductases, which are enzymes containing catalytic redox-active Cys residues, have been extensively studied, but even for these proteins there is little understanding of what distinguishes their catalytic redox Cys from other Cys functions. Herein, we characterized thiol oxidoreductases at a structural level and developed an algorithm that can recognize these enzymes by (i) analyzing amino acid and secondary structure composition of the active site and its similarity to known active sites containing redox Cys and (ii) calculating accessibility, active site location, and reactivity of Cys. For proteins with known or modeled structures, this method can identify proteins with catalytic Cys residues and distinguish thiol oxidoreductases from the enzymes containing other catalytic Cys types. Furthermore, by applying this procedure to Saccharomyces cerevisiae proteins containing conserved Cys, we could identify the majority of known yeast thiol oxidoreductases. This study provides insights into the structural properties of catalytic redox-active Cys and should further help to recognize thiol oxidoreductases in protein sequence and structure databases.     

Metal binding specificity in carbonic anhydrase is influenced by conserved hydrophobic core residues.

The role of highly conserved aromatic residues surrounding the zinc binding site of human carbonic anhydrase II (CAII) in determining the metal ion binding specificity of this enzyme has been examined by mutagenesis. Residues F93, F95, and W97 are located along a beta-strand containing two residues that coordinate zinc, H94 and H96, and these aromatic amino acids contribute to the high zinc affinity and slow zinc dissociation rate constant of CAII [Hunt, J. A., and Fierke, C. A. (1997) J. Biol. Chem. 272, 20364-20372]. Substitutions of these aromatic amino acids with smaller side chains enhance the copper affinity (up to 100-fold) while decreasing the affinity of both cobalt and zinc, thereby altering the metal binding specificity up to 10(4)-fold. Furthermore, the free energy of the stability of native CAII, determined by solvent-induced denaturation, correlates positively with increased hydrophobicity of the amino acids at positions 93, 95, and 97 as well as with cobalt and zinc affinity. Conversely, increased copper affinity correlates with decreased protein stability. Zinc specificity is therefore enhanced by formation of the native enzyme structure. These data suggest that the hydrophobic cluster in CAII is important for orienting the histidine residues to stabilize metals bound with a distorted tetrahedral geometry and to destabilize the trigonal bipyramidal geometry of bound copper. Knowledge of the structural factors that lead to high metal ion specificity will aid in the design of metal ion biosensors and de novo catalytic sites.