Zinc finger protein binding to DNA: an energy perspective using molecular dynamics simulation and free energy calculations on mutants of both zinc finger domains and their specific DNA bases

Energy calculations based on MM-GBSA were employed to study various zinc finger protein (ZF) motifs binding to DNA. Mutants of both the DNA bound to their specific amino acids were studied. Calculated energies gave evidence for a relationship between binding energy and affinity of ZF motifs to their sites on DNA. ΔG values were ?15.82(12), ?3.66(12), and ?12.14(11.6) kcal/mol for finger one, finger two, and finger three, respectively. The mutations in the DNA bases reduced the value of the negative energies of binding (maximum value for ΔΔG = 42Kcal/mol for F1 when GCG mutated to GGG, and ΔΔG = 22 kcal/mol for F2, the loss in total energy of binding originated in the loss in electrostatic energies upon mutation (r = .98). The mutations in key amino acids in the ZF motif in positions-1, 2, 3, and 6 showed reduced binding energies to DNA with correlation coefficients between total free energy and electrostatic was .99 and with Van der Waal was .93. Results agree with experimentally found selectivity which showed that Arginine in position-1 is specific to G, while Aspartic acid (D) in position 2 plays a complicated role in binding. There is a correlation between the MD calculated free energies of binding and those obtained experimentally for prepared ZF motifs bound to triplet bases in other reports (), our results may help in the design of ZF motifs based on the established recognition codes based on energies and contributing energies to the total energy.     

Binding of NDGA and morin with phospholipase A2: experimental and computational evidences

The effects of morin and nordihydroguaiaretic acid (NDGA), two plant secondary metabolites, on porcine pancreatic phospholipase A₂ (PLA₂) were investigated by isothermal titration calorimetry (ITC) and in silico docking analyses. The binding energies obtained for NDGA and morin from the ITC studies are ? 6.36 and ? 5.91 kcal mol^? 1, respectively. Similarly, the glide scores obtained for NDGA and morin towards PLA₂ were ? 7.32 and ? 7.23 kcal mol^? 1, respectively. Further the docked complexes were subjected to MD simulation in the presence of explicit water molecules to check the binding stability of the ligands in the active site of PLA₂. The bound ligands make hydrogen bonds with the active site residues of the enzyme and coordinate bonds with catalytically important Ca²⁺ ion. The binding of ligands at the active site of PLA₂ may also contribute to the reported anti-inflammatory properties of NDGA and morin.     

Development of a protein–ligand-binding site prediction method based on interaction energy and sequence conservation

We present a new method for predicting protein–ligand-binding sites based on protein three-dimensional structure and amino acid conservation. This method involves calculation of the van der Waals interaction energy between a protein and many probes placed on the protein surface and subsequent clustering of the probes with low interaction energies to identify the most energetically favorable locus. In addition, it uses amino acid conservation among homologous proteins. Ligand-binding sites were predicted by combining the interaction energy and the amino acid conservation score. The performance of our prediction method was evaluated using a non-redundant dataset of 348 ligand-bound and ligand-unbound protein structure pairs, constructed by filtering entries in a ligand-binding site structure database, LigASite. Ligand-bound structure prediction (bound prediction) indicated that 74.0 % of predicted ligand-binding sites overlapped with real ligand-binding sites by over 25 % of their volume. Ligand-unbound structure prediction (unbound prediction) indicated that 73.9 % of predicted ligand-binding residues overlapped with real ligand-binding residues. The amino acid conservation score improved the average prediction accuracy by 17.0 and 17.6 points for the bound and unbound predictions, respectively. These results demonstrate the effectiveness of the combined use of the interaction energy and amino acid conservation in the ligand-binding site prediction.     

Modeling of inhibitor–metalloenzyme interactions and selectivity using molecular mechanics grounded in quantum chemistry

We investigated the binding properties of the metalloprotease inhibitors hydroxamate, methanethiolate, and methylphosphoramidate to a model coordination site occurring in several Zn2+ metalloproteases, including thermolysin. This was carried out using both the SIBFA (sum of interactions between fragments ab initio-computed) molecular mechanics and the SCF/MP2 procedures for the purpose of evaluating SIBFA as a metalloenzyme modeling tool. The energy-minimized structures were closely similar to the X-ray crystallographic structures of related thermolysin-inhibitor complexes. We found that selectivity between alternative geometries and between inhibitors usually stemmed from multiple interaction components included in SIBFA. The binding strength sequence is hydroxamate > methanethiolate ≥ methylphosphoramidate from multiple interaction components included in SIBFA. The trends in interaction energy components, rankings, and preferences for mono- or bidentate binding were consistent in both computational procedures. We also compared the Zn2+ vs. Mg2+ selectivities in several other polycoordinated sites having various “hard” and “soft” qualities. This included a hexahydrate, a model representing Mg2+/Ca2+ binding sites, a chlorophyll-like structure, and a zinc finger model. The latter three favor Zn2+ over Mg2+ by a greater degree than the hydrated state, but the selectivity varies widely according to the ligand “softness.” SIBFA was able to match the ab initio binding energies by <2%, with the SIBFA terms representing dispersion and charge-transfer contributing the most to Zn2+/Mg2+ selectivity. These results showed this procedure to be a very capable modeling tool for metalloenzyme problems, in this case giving valuable information about details and limitations of “hard” and “soft” selectivity trends. Proteins 31:42–60, 1998. © 1998 Wiley-Liss, Inc.     

DFT studies of hydrocarbon combustion on metal surfaces

Catalytic combustion of hydrocarbons is an important technology to produce energy. Compared to conventional flame combustion, the catalyst enables this process to operate at lower temperatures; hence, reducing the energy required for efficient combustion. The reaction and activation energies of direct combustion of hydrocarbons (CH?→?C?+?H) on a series of metal surfaces were investigated using density functional theory (DFT). The data obtained for the Ag, Au, Al, Cu, Rh, Pt, and Pd surfaces were used to investigate the validity of the Brønsted-Evans-Polanyi (BEP) and transition state scaling (TSS) relations for this reaction on these surfaces. These relations were found to be valid (R²?=?0.94 for the BEP correlation and R²?=?1.0 for the TSS correlation) and were therefore used to estimate the energetics of the combustion reaction on Ni, Co, and Fe surfaces. It was found that the estimated transition state and activation energies (E_TS?=??69.70 eV and E_a?=?1.20 eV for Ni, E_TS?=??87.93 eV and E_a?=?1.08 eV for Co and E_TS?=??92.45 eV and E_a?=?0.83 eV for Fe) are in agreement with those obtained by DFT calculations (E_TS?=??69.98 eV and E_a?=?1.23 eV for Ni, E_TS?=??87.88 eV and E_a?=?1.08 eV for Co and E_TS?=??92.57 eV and E_a?=?0.79 eV for Fe). Therefore, these relations can be used to predict energetics of this reaction on these surfaces without doing the time consuming transition state calculations. Also, the calculations show that the activation barrier for CH dissociation decreases in the order Ag ? Au ? Al ? Cu ? Pt ? Pd ? Ni?>?Co?>?Rh?>?Fe.     

Thermal-stimulated pressure and current studies of bound water in lysozyme.

The state of bound water in crystalized lysozyme was studied by four techniques: electret thermal depolarization currents, thermal-stimulated pressure, isothermal polization decay, and thermogravimetry. Hydration levels ranged from 0 to 40 mg water/g protein. Desorption of bound water dipoles was found to be the main process responsible for electrical depolarization. Two different binding sites for water were identified with long relaxation times at room temperature (order 10(2)s) and activation energies of 0.34 plus or minus 0.02 eV and 0.55 plus or minus 0.04 eV.     

Binding Features of BCL2‐Targeted Oligodeoxynucleotides with K562 Cells

Delivery of various oligodeoxynucleotides into cells is mediated by binding to certain surface proteins followed by receptor‐mediated endocytosis. Moreover, oligonucleotides are able to provoke perturbation of cell surface proteins and growth factor receptors among them. Here we described binding sense BCL2 oligodeoxynucleotide targeted to translation start of BCL2 mRNA (ODN) with K562 cells. At low concentration ODN bound efficiently with K562 and penetrated into the cells via binding cell surface with rather high affinity and priming new binding sites. The loose binding constant at 4°C was 1.8 × 10⁹ M^? 1 both for binding ODN in solution and ODN‐associated liposome. The number of loose binding sites under both treatments was rather high: 4.6 to 6.6 pmoles per 10⁶ cells. The extent of ODN penetration into the cells showed higher potential site numbers than initially seen and reached 8.6 pmoles per 10⁶ cells for four hours incubation at 37°C.     

Cooperativity of complex salt bridges

The energetic contribution of complex salt bridges, in which one charged residue (anchor residue) forms salt bridges with two or more residues simultaneously, has been suggested to have importance for protein stability. Detailed analysis of the net energetics of complex salt bridge formation using double- and triple-mutant cycle analysis revealed conflicting results. In two cases, it was shown that complex salt bridge formation is cooperative, i.e., the net strength of the complex salt bridge is more than the sum of the energies of individual pairs. In one case, it was reported that complex salt bridge formation is anti-cooperative. To resolve these different findings, we performed analysis of the geometries of salt bridges in a representative set of structures from the PDB and found that over 87% of all complex salt bridges anchored by Arg/Lys have a geometry such that the angle formed by their Calpha atoms, Theta, is <90 degrees . This preferred geometry is observed in the two reported instances when the energetics of complex salt bridge formation is cooperative, while in the reported anti-cooperative complex salt bridge, Theta is close to 160 degrees . Based on these observations, we hypothesized that complex salt bridges are cooperative for Theta < 90 degrees and anti-cooperative for 90 degrees < Theta < 180 degrees . To provide a further experimental test for this hypothesis, we engineered a complex salt bridge with Theta = 150 degrees into a model protein, the activation domain of human procarboxypeptidase A2 (ADA2h). Experimentally derived stabilities of the ADA2h variants allowed us to show that the complex salt bridge in ADA2h is anti-cooperative.     

Ionic interactions and anion binding in the gramicidin channel. An electrostatic calculation

An electrostatic calculation is performed in order to examine basic features of interactions between ions (including anions) within the gramicidin channel. The calculation focuses on the effect of image forces. The substitute charge method is used for the calculation of image-force energies. Good arrangements of fictitious charges and contour points are described. Errors of calculated image-force energies are estimated not to greatly exceed 0.1%. The cases assumed are (i) the effective radius r of the channel is between 2.5 and 3.5 A, (ii) the binding site with the highest affinity is between 1 and 3.5 A from the channel end (outer site), and (iii) the dimple at the channel mouth is 0-5 A in depth. The induced energy of an ion placed at the outer site increases (and hence the affinity of the outer site decreases) with the increase in the depth of dimples, whereas the barrier height for translocation between the outer sites decreases in the presence of deeper dimples. The interactional energy between two monovalent cations placed at the outer sites is relatively small in the absence of dimples. It is large, however, in the presence of deeper dimples if the outer sites are 2.5 or 3.5 A from the ends; but, it is still relatively small even in the presence of dimples if the sites are 1 A from the ends. The interactional energy is very unfavorable for simultaneous occupancy by three cations. It is suggested that an ion pair may be formed at the channel mouth. Deeper positions of the outer site, smaller values of r, and deeper dimples favor the formation of the ion pair. In the presence of 5 A dimples, the binding constant of an anion (Cl-) for a cation (Na+) which has already been bound at an outer site (with no second cation at the opposite mouth) is estimated to be 0.4-8 molal-1 if local interactions between the bound anion and the channel wall is negligibly small. The anion binding constant increases in the presence of a cation (or an ion pair) at the opposite mouth. It decreases markedly in the absence of dimples. The interactional energy is considerably unfavorable for the binding of a third cation even if the presence of an ion pair is postulated. It is still large in the presence of an ion pair at each of the mouths.     

Tuning hydrogen storage of carbon nanotubes by mechanical bending: theoretical study

The effects of mechanical bending on tuning the hydrogen storage of titanium functionalised (4,0) carbon nanotube have been assessed using density functional theory calculations with reference to the ultimate targets of the US Department of Energy (DOE). The assessment has been carried out in terms of physisorption, gravimetric capacity, projected densities of states, statistical thermodynamic stability and reaction kinetics. The Ti atom binds at the hollow site of the hexagonal ring. The average adsorption energies (?0.54 eV) per hydrogen molecule meet the DOE target for physisorption (?0.20 to ?0.60 eV). The curvature attributed to the bending angle has no effect on the average adsorption energies per H₂ molecule. With no metal clustering, the system gravimetric capacities are expected to be as large as 9.0 wt%. The reactions of the deformed (bent) carbon nanotube have higher probabilities of occurring than those of the un-deformed carbon nanotube. The Gibbs free energies, enthalpies and entropies meet the ultimate targets of the DOE for all temperatures and pressures. The closest reactions to zero free energy occur at (378.15 K/2.961 atm.) and reverse at (340 and 360 K/1 atm.). The translational component is found to exact a dominant effect on the total entropy change with temperature. Favourable kinetics of the reactions at the temperatures targeted by DOE are reported regardless of the applied pressure. The more preferable thermodynamic properties assigned to the bending nanotube imply that hydrogen storage can be improved compared to the nonbending nanotube.     

Oxygen doped SiC nanocrystals: first principles study of the optical properties

We have studied a typical spherical SiC nanocrystal with a diameter of 1.2 nm (Si₄₃C₄₄H₇₆) using linear combination of atomic orbitals in combination with pesudopotential density functional calculation. The role of fluorine and oxygen impurities was investigated on the electronic and optical properties of the Si₄₃C₄₄H₇₆ nanocrystal. Total energy calculations show that the fluorine doped Si₄₃C₄₄H₇₆ nanocrystals are unstable. Oxygen doped Si₄₃C₄₄H₇₆ have different binding energies in various substitutional and interstitial defects. The maximum binding energy of the oxygen at carbon substitutional defect is about ?0.5 eV and at interstitial defect is ?0.18 eV. The HOMO-LUMO gap of the pure Si₄₃C₄₄H₇₆ is about 6.71 eV and after doping with oxygen changes on the order of 0.1 eV. Our studies show that the refractive index of the doped Si₄₃C₄₄H₇₆ nanocrystal significantly dispersed in comparison with pure SiC nanocrystal especially at the range of 6 to 8 eV.     

Allosteric sites can be identified based on the residue–residue interaction energy difference

Allosteric drugs act at a distance to regulate protein functions. They have several advantages over conventional orthosteric drugs, including diverse regulation types and fewer side effects. However, the rational design of allosteric ligands remains a challenge, especially when it comes to the identification allosteric binding sites. As the binding of allosteric ligands may induce changes in the pattern of residue–residue interactions, we calculated the residue–residue interaction energies within the allosteric site based on the molecular mechanics generalized Born surface area energy decomposition scheme. Using a dataset of 17 allosteric proteins with structural data for both the apo and the ligand‐bound state available, we used conformational ensembles generated by molecular dynamics simulations to compute the differences in the residue–residue interaction energies in known allosteric sites from both states. For all the known sites, distinct interaction energy differences (>25%) were observed. We then used CAVITY, a binding site detection program to identify novel putative allosteric sites in the same proteins. This yielded a total of 31 “druggable binding sites,” of which 21 exhibited >25% difference in residue interaction energies, and were hence predicted as novel allosteric sites. Three of the predicted allosteric sites were supported by recent experimental studies. All the predicted sites may serve as novel allosteric sites for allosteric ligand design. Our study provides a computational method for identifying novel allosteric sites for allosteric drug design. Proteins 2015; 83:1375–1384. © 2014 Wiley Periodicals, Inc.     

Exploring the Ion Selectivity Properties of a Large Number of Simplified Binding Site Models

The ability to discriminate between different cations efficiently is essential for the proper physiological functioning of many membrane transport proteins. One obvious mechanism of ion selectivity is when a binding site is structurally constrained by the protein architecture and its geometry is precisely adapted to fit an ion of a given size. This mechanism is not effective in the case of flexible protein binding sites that are able to deform structurally or to adapt to a bound ion. In this study, the concept of nontrivial ion selectivity arising in a highly flexible protein binding site conceptually represented as a microdroplet of ligands confined to a small volume is explored. The environment imposed by the spatial confinement is a critical feature of the reduced models. A large number of reduced binding site models (1077) comprising typical ion-coordinating ligands (carbonyl, hydroxyl, carboxylate, water) are constructed and characterized for Na⁺/K⁺ and Ca²⁺/Ba²⁺ size selectivity using free energy perturbation molecular dynamics simulations. Free energies are highly correlated with the sum of ion-ligand and ligand-ligand mean interactions, but the relative balance of those two contributions is different for K⁺-selective and Na⁺-selective binding sites. The analysis indicates that both the number and the type of ligands are important factors in ion selectivity.     

Interaction of flavonols with human serum albumin: a biophysical study showing structure–activity relationship and enhancement when coated on silver nanoparticles

Binding affinities of flavonols namely quercetin, myricetin, and kaempferol to human serum albumin (HSA) were determined fluorimetrically and the order was observed to be myricetin > quercetin > kaempferol demonstrating structure–activity relationship. Quercetin-coated silver nanoparticles (AgNPs) show higher binding affinity to HSA compared to free quercetin with binding constants 6.04 × 10⁷ M^?1 and 4.2 × 10⁶ M^?1, respectively. Using site-specific markers it is concluded that free quercetin and that coated on AgNPs bind at different sites. Significant structural changes in circular dichroism (CD) spectra of HSA were recorded with quercetin-coated AgNPs compared to free quercetin. These results were further substantiated by time-resolved fluorescence spectroscopy where fluorescence life time of the tryptophan residue in HSA–quercetin-coated AgNPs complex decreased to 3.63 ns from 4.22 ns in HSA–quercetin complex. Isothermal calorimetric studies reveal two binding modes for quercetin-coated AgNPs and also higher binding constants compared to free quercetin. These higher binding affinities are attributed to altered properties of quercetin when coated on AgNPs enabling it to reach the binding sites other than site II where free quercetin mainly binds.     

Revised model of calcium and magnesium binding to the bacterial cell wall

Metals bind to the bacterial cell wall, yet the binding mechanisms and affinity constants are not fully understood. The cell wall of gram positive bacteria is characterized by a thick layer of peptidoglycan and anionic teichoic acids anchored in the cytoplasmic membrane as lipoteichoic acid or covalently bound to the cell wall as wall teichoic acid. The polyphosphate groups of teichoic acid provide one-half of the metal binding sites for calcium and magnesium, which contradicts previous reports that calcium binding is 100 % dependent on teichoic acid. The remaining binding sites are formed with the carboxyl units of peptidoglycan. In this work we report equilibrium association constants and total metal binding capacities for the interaction of calcium and magnesium ions with the bacterial cell wall. Metal binding is much stronger than previously reported. Curvature of Scatchard plots from the binding data and the resulting two regions of binding affinity suggest the presence of negative cooperative binding, which means that the binding affinity decreases as more ions become bound to the sample. For Ca²⁺, Region I has a K_A = (1.0 ± 0.2) × 10⁶ M^?1 and Region II has a K_A = (0.075 ± 0.058) × 10⁶ M^?1. For Mg²⁺, K_A1 = (1.5 ± 0.1) × 10⁶ and K_A2 = (0.17 ± 0.10) × 10⁶. A binding capacity (η) is reported for both regions. However, since binding is still occurring in Region II, the total binding capacity is denoted by η₂, which are 0.70 ± 0.04 and 0.67 ± 0.03 µmol/mg for Ca²⁺ and Mg²⁺ respectively. These data contradict the current paradigm of only a single metal affinity value that is constant over a range of concentrations. We also find that measurement of equilibrium binding constants is highly sample dependent. This suggests a role for diffusion of metals through heterogeneous cell wall fragments. As a result, we are able to reconcile many contradictory theories that describe binding affinity and the binding mode of divalent metal cations.     

Preferential binding of fisetin to the native state of bovine serum albumin: spectroscopic and docking studies

We have investigated the binding of the biologically important flavonoid fisetin with the carrier protein bovine serum albumin using multi-spectroscopic and molecular docking methods. The binding constants were found to be in the order of 10⁴ M^?1 and the number of binding sites was determined as one. MALDI-TOF analyses showed that one fisetin molecule binds to a single bovine serum albumin (BSA) molecule which is also supported by fluorescence quenching studies. The negative Gibbs free energy change (?G°) values point to a spontaneous binding process which occurs through the presence of electrostatic forces with hydrophobic association that results in a positive entropy change (+51.69 ± 1.18 J mol^?1 K^?1). The unfolding and refolding of BSA in urea have been studied in absence and presence of fisetin using steady-state fluorescence and lifetime measurements. Urea denaturation studies indicate that fisetin is gradually released from its binding site on the protein. In the absence of urea, an increase in temperature that causes denaturation of the protein results in the release of fisetin from its bound state indicating that fisetin binds only to the native state of the protein. The circular dichroism (CD) and Fourier transform infrared (FTIR) spectroscopic studies showed an increase in % α-helix content of BSA after binding with fisetin. Site marker displacement studies in accordance with the molecular docking results suggested that fisetin binds in close proximity of the hydrophobic cavity in site 1 (subdomain IIA) of the protein. The PEARLS (Program of Energetic Analysis of Receptor Ligand System) has been used to estimate the interaction energy of fisetin with BSA and the results are in good correlation with the experimental findings.     

Influence of cysteine 164 on active site structure in rat cysteine dioxygenase

Cysteine dioxygenase is a non-heme mononuclear iron enzyme with unique structural features, namely an intramolecular thioether cross-link between cysteine 93 and tyrosine 157, and a disulfide bond between substrate l-cysteine and cysteine 164 in the entrance channel to the active site. We investigated how these posttranslational modifications affect catalysis through a kinetic, crystallographic and computational study. The enzyme kinetics of a C164S variant are identical to WT, indicating that disulfide formation at C164 does not significantly impair access to the active site at physiological pH. However, at high pH, the cysteine–tyrosine cross-link formation is enhanced in C164S. This supports the view that disulfide formation at position 164 can limit access to the active site. The C164S variant yielded crystal structures of unusual clarity in both resting state and with cysteine bound. Both show that the iron in the cysteine-bound complex is a mixture of penta- and hexa-coordinate with a water molecule taking up the final site (60 % occupancy), which is where dioxygen is believed to coordinate during turnover. The serine also displays stronger hydrogen bond interactions to a water bound to the amine of the substrate cysteine. However, the interactions between cysteine and iron appear unchanged. DFT calculations support this and show that WT and C164S have similar binding energies for the water molecule in the final site. This variant therefore provides evidence that WT also exists in an equilibrium between penta- and hexa-coordinate forms and the presence of the sixth ligand does not strongly affect dioxygen binding.     

The study on interactions between levofloxacin and model proteins by using multi-spectroscopic and molecular docking methods

The interactions of levofloxacin (LEV) with lysozyme (LYZ), trypsin and bovine hemoglobin (BHb) were investigated, respectively, by using multi-spectral techniques and molecular docking in vitro. Fluorescence studies showed that LEV quenched LYZ/trypsin fluorescence in a combined quenching ways and BHb fluorescence in a static quenching with binding constants of .14, .51 and .20 × 10⁵ L mol^?1 at 298 K, respectively. The thermodynamic parameters demonstrated that hydrophobic forces, hydrogen bonds, and van der Waals forces played the major role in the binding process. The binding distances between LEV and the inner tryptophan residues of LYZ, trypsin, and BHb were calculated to be 4.04, 3.38, and 4.52 nm, respectively. Furthermore, the results of circular dichroism spectra (CD), UV–vis, and three-dimensional fluorescence spectra indicated that the secondary structures of LYZ, trypsin, and BHb were partially changed by LEV with the α-helix percentage of LYZ-LEV system increased while that of BHb-LEV system was decreased, the β-sheet percentage of trypsin-LEV system increased from 41.3 to 42.9%. UV–vis spectral results showed that the binding interactions could cause conformational and some micro-environmental changes of LYZ, trypsin, and BHb. The results of molecular docking revealed that in LYZ and trypsin systems, LEV bound to the active sites residues GLU 35 and ASP 52 of LYZ and trypsin at the active site SER 195, and in BHb system, LEV was located in the central cavity, which was consistent with the results of synchronous fluorescence experiment. Besides, LEV made the activity of LYZ decrease while the activity of trypsin increased.     

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.