A less‐biased analysis of metalloproteins reveals novel zinc coordination geometries

Zinc metalloproteins are involved in many biological processes and play crucial biochemical roles across all domains of life. Local structure around the zinc ion, especially the coordination geometry (CG), is dictated by the protein sequence and is often directly related to the function of the protein. Current methodologies in characterizing zinc metalloproteins' CG consider only previously reported CG models based mainly on nonbiological chemical context. Exceptions to these canonical CG models are either misclassified or discarded as “outliers.” Thus, we developed a less‐biased method that directly handles potential exceptions without pre‐assuming any CG model. Our study shows that numerous exceptions could actually be further classified and that new CG models are needed to characterize them. Also, these new CG models are cross‐validated by strong correlation between independent structural and functional annotation distance metrics, which is partially lost if these new CGs models are ignored. Furthermore, these new CG models exhibit functional propensities distinct from the canonical CG models. Proteins 2015; 83:1470–1487. © 2015 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

Insights into metalloproteins and metallodrugs from electron paramagnetic resonance spectroscopy

Metal ions play an important role in diverse biological processes, and much of the basic knowledge derived from studying native bioinorganic systems are applied in the synthesis of new molecules with the aim of diagnosing and treating diseases. At first glance, metalloproteins and metallodrugs are very different systems, but metal ion coordination, redox chemistry and substrate binding play essential roles in advancing both of these research fields. In this article, we discuss recent metalloprotein and metallodrug studies where electron paramagnetic resonance spectroscopy served as a major tool to gain a better understanding of metal-based structures and their function.     

Can the propensity of protein crystallization be increased by using systematic screening with metals?

Protein crystallization is one of the major bottlenecks in protein structure elucidation with new strategies being constantly developed to improve the chances of crystallization. Generally, well‐ordered epitopes possessing complementary surface and capable of producing stable inter‐protein interactions generate a regular three‐dimensional arrangement of protein molecules which eventually results in a crystal lattice. Metals, when used for crystallization, with their various coordination numbers and geometries, can generate such epitopes mediating protein oligomerization and/or establish crystal contacts. Some examples of metal‐mediated oligomerization and crystallization together with our experience on metal‐mediated crystallization of a putative rRNA methyltransferase from Sinorhizobium meliloti are presented. Analysis of crystal structures from protein data bank (PDB) using a non‐redundant data set with a 90% identity cutoff, reveals that around 67% of proteins contain at least one metal ion, with ～14% containing combination of metal ions. Interestingly, metal containing conditions in most commercially available and popular crystallization kits generally contain only a single metal ion, with combinations of metals only in a very few conditions. Based on the results presented in this review, it appears that the crystallization screens need expansion with systematic screening of metal ions that could be crucial for stabilizing the protein structure or for establishing crystal contact and thereby aiding protein crystallization.     

Bidentate and tridentate metal-ion coordination states within ternary complexes of RB69 DNA polymerase

Two divalent metal ions are required for primer‐extension catalyzed by DNA polymerases. One metal ion brings the 3′‐hydroxyl of the primer terminus and the α‐phosphorus atom of incoming dNTP together for bond formation so that the catalytically relevant conformation of the triphosphate tail of the dNTP is in an α,β,γ‐tridentate coordination complex with the second metal ion required for proper substrate alignment. A probable base selectivity mechanism derived from structural studies on Dpo4 suggests that the inability of mispaired dNTPs to form a substrate‐aligned, tridentate coordination complex could effectively cause the mispaired dNTPs to be rejected before catalysis. Nevertheless, we found that mispaired dNTPs can actually form a properly aligned tridentate coordination complex. However, complementary dNTPs occasionally form misaligned complexes with mutant RB69 DNA polymerases (RB69pols) that are not in a tridentate coordination state. Here, we report finding a β,γ‐bidentate coordination complex that contained the complementary dUpNpp opposite dA in the structure of a ternary complex formed by the wild type RB69pol at 1.88 Å resolution. Our observations suggest that several distinct metal‐ion coordination states can exist at the ground state in the polymerase active site and that base selectivity is unlikely to be based on metal‐ion coordination alone.     

Phosphorogenic and spontaneous formation of tris(bipyridine)ruthenium in peptide scaffolds

Redox‐active ruthenium complexes have been widely used in various fields; however, the harsh conditions required for their synthesis are not always conducive to their subsequent use in biological applications. In this study, we demonstrate the spontaneous formation of a derivative of tris(bipyridine)ruthenium at 37°C through the coordination of three bipyridyl ligands incorporated into a peptide to a ruthenium ion. Specifically, we synthesized six bipyridyl‐functionalized peptides with randomly chosen sequences. The six peptides bound to ruthenium ions and exhibited similar spectroscopic and electrochemical features to tris(bipyridine)ruthenium, indicating the formation of ruthenium complexes as we anticipated. The photo‐excited triplet state of the ruthenium complex formed in the peptides exhibited an approximately 1.6‐fold longer lifetime than that of tris(bipyridine)ruthenium. We also found that the photo‐excited state of the ruthenium complexes was able to transfer an electron to methyl viologen, indicating that the ruthenium complexes formed in the peptides had the same ability to transfer charge as tris(bipyridine)ruthenium. We believe that this strategy of producing ruthenium complexes in peptides under mild conditions will pave the way for developing new metallopeptides and metalloproteins containing functional metal‐complexes.     

Metal Ion availability in mitochondria

Transition metal ions are required for many aspects of mitochondrial physiology. Copper, iron, manganese and zinc are cofactors in metalloenzymes and metalloproteins within the organelle. Little is known how cells maintain optimal pools of these metal ions for mitochondrial function. This review documents the available literature on mitochondrial metal ion pools and protein metallation reactions. Upon perturbation in metal pools, mis-metallation reactions do occur. Thus, regulation of metal ion accessibility and bioavailability must exist.     

Functional role of metalloproteins in genome stability

Cells contain a large number of metalloproteins that commonly harbor at least one metal ion cofactor. In metalloproteins, metal ions are usually coordinated by oxygen, sulfur, or nitrogen centers belonging to amino acid residues in the protein. The presence of the metal ion in metalloproteins allows them to take part in diverse biological processes, such as genome stability, metabolic catalysis, and cell cycle progression. Clinically, alteration of the function of metalloproteins in mammals is genetically associated with diseases characterized by DNA damage and repair defects. The present review focuses on the current perspectives of metal ion homeostasis in different organisms and summarizes the most recent understanding on magnesium, copper, iron, and manganese-containing proteins and their functional involvement in the maintenance of genome stability.     

Metals on the move: zinc ions in cellular regulation and in the coordination dynamics of zinc proteins

Homeostatic control maintains essential transition metal ions at characteristic cellular concentrations to support their physiological functions and to avoid adverse effects. Zinc is especially widely used as a catalytic or structural cofactor in about 3000 human zinc proteins. In addition, the homeostatic control of zinc in eukaryotic cells permits functions of zinc(II) ions in regulation and in paracrine and intracrine signaling. Zinc ions are released from proteins through ligand-centered reactions in zinc/thiolate coordination environments, and from stores in cellular organelles, where zinc transporters participate in zinc loading and release. Muffling reactions allow zinc ions to serve as signaling ions (second messengers) in the cytosol that is buffered to picomolar zinc ion concentrations at steady-state. Muffling includes zinc ion binding to metallothioneins, cellular translocations of metallothioneins, delivery of zinc ions to transporter proteins, and zinc ion fluxes through cellular membranes with the result of removing the additional zinc ions from the cytosol and restoring the steady-state. Targets of regulatory zinc ions are proteins with sites for transient zinc binding, such as membrane receptors, enzymes, protein–protein interactions, and sensor proteins that control gene expression. The generation, transmission, targets, and termination of zinc ion signals involve proteins that use coordination dynamics in the inner and outer ligand spheres to control metal ion association and dissociation. These new findings establish critically important functions of zinc ions and zinc metalloproteins in cellular control.     

The best-fit size and geometry of metal ions for coordination with nitrogen-donor macrocycles

The technique of calculating the strain energy of metal ion complexes as a function of metal to ligand bond length (Hancock and McDougall, J. Am. Chem. Soc., 102 (1980) 6553) is used to study best-fit sizes of metal ions for coordinating with tetraaza and triaza macrocyles. In addition to varying the metal to ligand bond length in the calculations, different coordination geometries of the metal ion are also examined. The metal to nitrogen (M-N) bond lengths, and coordination geometries, that give lowest energies for several N-donor macrocyles, are calculated by molecular mechanics, and 16-aneN₄ (1,5,9,13-tetraazacyclohexadecane) is found, contrary to popular belief, to coordinate best with very small metal ions, with lowest energy occurring for a slightly flattened tetrahedral metal ion of M-N length = 1.81 Å. The best-fit size and geometry for coordination in 12-aneN₄ (1,4,7,10-tetraazacyclotetradecane) is an M-N length of 2.15 Å and square pyramidal geometry, and with cyclam (1,4,8,11-tetraazacyclotetradecane) an M-N length of 2.06 Å and planar coordination that is approximately square.     

Protein solution structure calculations in solution: Solvated molecular dynamics refinement of calbindin D9k

The three-dimensional solution structures of proteins determinedwith NMR-derived constraints are almost always calculated in vacuo. Thesolution structure of (Ca²⁺)_{_2}-calbindinD_9k has been redetermined by new restrained molecular dynamics(MD) calculations that include Ca²⁺ ions and explicit solventmolecules. Four parallel sets of MD refinements were run to provide accuratecomparisons of structures produced in vacuo, in vacuo withCa²⁺ ions, and with two different protocols in a solvent bathwith Ca²⁺ ions. The structural ensembles were analyzed interms of structural definition, molecular energies, packing density,solvent-accessible surface, hydrogen bonds, and the coordination of calciumions in the two binding loops. Refinement including Ca²⁺ ionsand explicit solvent results in significant improvements in the precisionand accuracy of the structure, particularly in the binding loops. Theseresults are consistent with results previously obtained in free MDsimulations of proteins in solution and show that the rMD refinedNMR-derived solution structures of proteins, especially metalloproteins, canbe significantly improved by these strategies.     

Ion‐Exchange and Adsorption of Fe(III) by Sepia Melanin

Sepia eumelanin is associated with many metal ions, yet little is known about its metal binding capacity and the chemical nature of the binding site(s). Herein, the natural concentrations of metal ions are presented and the ability to remove metals by exposure of the melanin granules to EDTA is quantified. The results reveal that the binding constants of melanin at pH 5.8 for Mg(II), Ca(II), Sr(II) and Cu(II) are, respectively, 5, 4, 14 and 34 times greater than the corresponding binding constants of these ions with EDTA. By exposing Sepia eumelanin to aqueous solutions of FeCl₃, the content of bound Fe(III) can be increased from a natural concentration of ～180 ppm to a saturation limit of ～80 000 ppm or 1.43 mmol/g of melanin. Similar saturation limits are found for Mg(II) and Ca(II). Exposure of Sepia melanin granules to aqueous solutions containing Ca(II) results in the stoichiometric replacement of the initially bound Mg(II), arguing that these two ions occupy the same binding site(s) in the pigment. The pH‐dependent binding of Mg(II) and Ca(II) suggests coordination of these ions to carboxylic acid groups in the pigment. Mg(II) and Ca(II) can be added to a Fe(III)‐saturated melanin sample without affecting the amount of Fe(III) pre‐adsorbed, clearly establishing Fe(III) and Mg(II)/Ca(II) occupy different binding sites. Taking recent Raman spectroscopic data into account, the binding of Fe(III) is concluded to involve coordination to o‐dihydroxyl groups. The effects of metal ion content on the surface morphology were analyzed. No significant changes were found over the full range of Fe(III) concentration studied, which is supported by the Brunauer–Emmett–Teller surface area analysis. These observations imply the existence of channels within the melanin granules that can serve to transport metal ions.     

Nickel(II) complexes stabilized by bis[N-(6-pivalamido-2-pyridylmethyl)]benzylamine: Synthesis and characterization of complexes stabilized by a hydrogen bonding network

Intramolecular hydrogen bonds in metalloproteins are key in directing reactivity yet these effects have been difficult to achieved in synthetic systems. We have been developing a synthetic system that uses hydrogen-bonding interactions to modulate the secondary coordination around a transition metal ion. This was accomplished with the ligand bis[N-(6-pivalamido-2-pyridylmethyl)]benzylamine (H₂pmb), which contains two carboxyamido units appended from pyridine rings. Several nickel complexes were prepared and structurally characterized. In particular, we found that the appended carboxyamido groups either provide intramolecular H-bond donors or can be converted to bind directly to a metal center. We established that the complex Ni^IIH₂pmb(Cl)₂ can be sequentially deprotonated with potassium tert-butoxide, causing coordination of the carboxyamido oxygen atoms and concomitant loss of the chloro ligands. The chloro ligands were also removed with silver(I) salts in the presence of acetate ions and the complex Ni^IIH₂pmb(κ²-OAc)(κ¹-OAc) was isolated, in which an intramolecular H-bonding network occurs between the H₂pmb ligand and the coordinate acetato ligands.     

Prediction of structures of zinc‐binding proteins through explicit modeling of metal coordination geometry

Metal ions play an essential role in stabilizing protein structures and contributing to protein function. Ions such as zinc have well‐defined coordination geometries, but it has not been easy to take advantage of this knowledge in protein structure prediction efforts. Here, we present a computational method to predict structures of zinc‐binding proteins given knowledge of the positions of zinc‐coordinating residues in the amino acid sequence. The method takes advantage of the “atom‐tree” representation of molecular systems and modular architecture of the Rosetta3 software suite to incorporate explicit metal ion coordination geometry into previously developed de novo prediction and loop modeling protocols. Zinc cofactors are tethered to their interacting residues based on coordination geometries observed in natural zinc‐binding proteins. The incorporation of explicit zinc atoms and their coordination geometry in both de novo structure prediction and loop modeling significantly improves sampling near the native conformation. The method can be readily extended to predict protein structures bound to other metal and/or small chemical cofactors with well‐defined coordination or ligation geometry.     

Metalloproteomics, metalloproteomes, and the annotation of metalloproteins

Metalloproteomics includes approaches that address the expression of metalloproteins and their changes in biological time and space. Metalloproteomes are investigated by a combination of approaches. Experimental approaches include structural genomics, which provides insights into the architecture of metal-binding sites in metalloproteins and establishes ligand signatures from the types and spacings of the metal ligands in the protein sequence. Theoretical approaches employ these ligand signatures as templates for homology searches in sequence databases. In this way, the number of metalloproteins in the iron, copper, and zinc metalloproteomes in various phyla of life has been estimated. Yet, manganese metalloproteomes remain poorly defined. Metals have catalytic and structural functions in proteins. However, additional functions have evolved. Proteins that control metal homeostasis and proteins that are metal-regulated bind metal ions transiently and are generally not accounted for in estimates from bioinformatics. Thus, metalloproteomes are dynamic and likely to be larger than present estimates suggest. This account discusses the assignment of transition metals in metalloproteins and the ensuing issues facing analytical chemists and structural and computational biologists. Biological and chemical selectivities render metal selection by metalloproteins either more stringent or less stringent depending on the metal homeostatic system of the organism, the subcellular location of the protein, and environmental factors. Failure to recognize the principles of metal utilization has led to assigning the wrong metal in metalloproteins and has missed some of the regulatory functions of transition metal ions.     

Coordination chemistry of glutathione

The metal ion coordination abilities of reduced and oxidized glutathione are reviewed. Reduced glutathione (GSH) is a very versatile ligand, forming stable complexes with both hard and soft metal ions. Several general binding modes of GSH are described. Soft metal ions coordinate exclusively or primarily through thiol sulfur. Hard ones prefer the amino acid-like moiety of the glutamic acid residue. Several transition metal ions can additionally coordinate to the peptide nitrogen of the gamma-Glu-Cys bond. Oxidized glutathione lacks the thiol function. Nevertheless, it proves to be a surprisingly efficient ligand for a range of metal ions, coordinating them primarily through the donors of the glutamic acid residue.     

Probing determinants of the metal ion selectivity in carbonic anhydrase using mutagenesis

Few studies measuring thermodynamic metal ion selectivity of metalloproteins have been performed, and the major determinants of metal ion selectivity in proteins are not yet well understood. Several features of metal ion binding sites and metal coordination have been hypothesized to alter the transition metal selectivity of chelators, including (1) the polarizability of the coordinating atom, (2) the relative sizes of the binding site and the metal ion, and (3) the metal ion binding site geometry. To test these hypotheses, we have measured the metal ion affinity and selectivity of a prototypical zinc enzyme, human carbonic anhydrase II (CAII), and a number of active site variants where one of the coordinating ligands is substituted by another side chain capable of coordinating metal. CAII and almost all of the variants follow the inherent metal ion affinity trend suggested by the Irving-Williams series, demonstrating that this trend operates within proteins as well as within small molecule chelators and may be a dominant factor in metal ion selectivity in biology. Neither the polarizability of the liganding side chains nor the size of the metal ion binding site correlates strongly with metal ion specificity; instead, changes in metal ion specificity in the variants correlate with the preferred coordination number and geometry of the metal ion. This correlation suggests that a primary feature driving deviations from the inherent ligand affinity trend is the positioning of active site groups such that a given metal ion can adopt a preferred coordination number/geometry.     

Structural Analysis of Metal Sites in Proteins: Non-heme Iron Sites as a Case Study

In metalloproteins, the protein environment modulates metal properties to achieve the required goal, which can be protein stabilization or function. The analysis of metal sites at the atomic level of detail provided by protein structures can thus be of benefit in functional and evolutionary studies of proteins. In this work, we propose a structural bioinformatics approach to the study of metalloproteins based on structural templates of metal sites that include the PDB coordinates of protein residues forming the first and the second coordination sphere of the metal. We have applied this approach to non-heme iron sites, which have been analyzed at various levels. Templates of sites located in different protein domains have been compared, showing that similar sites can be found in unrelated proteins as the result of convergent evolution. Templates of sites located in proteins of a large superfamily have been compared, showing possible mechanisms of divergent evolution of proteins to achieve different functions. Furthermore, template comparisons have been used to predict the function of uncharacterized proteins, showing that similarity searches focused on metal sites can be advantageously combined with typical whole-domain comparisons. Structural templates of metal sites, finally, may constitute the basis for a systematic classification of metalloproteins in databases.     

Shaping mechanisms of metal specificity in a family of metazoan metallothioneins: evolutionary differentiation of mollusc metallothioneins

Background  

The degree of metal binding specificity in metalloproteins such as metallothioneins (MTs) can be crucial for their functional accuracy. Unlike most other animal species, pulmonate molluscs possess homometallic MT isoforms loaded with Cu⁺ or Cd²⁺. They have, so far, been obtained as native metal-MT complexes from snail tissues, where they are involved in the metabolism of the metal ion species bound to the respective isoform. However, it has not as yet been discerned if their specific metal occupation is the result of a rigid control of metal availability, or isoform expression programming in the hosting tissues or of structural differences of the respective peptides determining the coordinative options for the different metal ions. In this study, the Roman snail (Helix pomatia) Cu-loaded and Cd-loaded isoforms (HpCuMT and HpCdMT) were used as model molecules in order to elucidate the biochemical and evolutionary mechanisms permitting pulmonate MTs to achieve specificity for their cognate metal ion.     

Control of redox reactivity of flavin and pterin coenzymes by metal ion coordination and hydrogen bonding

The electron-transfer activities of flavin and pterin coenzymes can be fine-tuned by coordination of metal ions, protonation and hydrogen bonding. Formation of hydrogen bonds with a hydrogen-bond receptor in metal–flavin complexes is made possible depending on the type of coordination bond that can leave the hydrogen-bonding sites. The electron-transfer catalytic functions of flavin and pterin coenzymes are described by showing a number of examples of both thermal and photochemical redox reactions, which proceed by controlling the electron-transfer reactivity of coenzymes with metal ion binding, protonation and hydrogen bonding.     

Development of METAL‐ACTIVE SITE and ZINCCLUSTER tool to predict active site pockets

The advent of whole genome sequencing leads to increasing number of proteins with known amino acid sequences. Despite many efforts, the number of proteins with resolved three dimensional structures is still low. One of the challenging tasks the structural biologists face is the prediction of the interaction of metal ion with any protein for which the structure is unknown. Based on the information available in Protein Data Bank, a site (METALACTIVE INTERACTION) has been generated which displays information for significant high preferential and low‐preferential combination of endogenous ligands for 49 metal ions. User can also gain information about the residues present in the first and second coordination sphere as it plays a major role in maintaining the structure and function of metalloproteins in biological system. In this paper, a novel computational tool (ZINCCLUSTER) is developed, which can predict the zinc metal binding sites of proteins even if only the primary sequence is known. The purpose of this tool is to predict the active site cluster of an uncharacterized protein based on its primary sequence or a 3D structure. The tool can predict amino acids interacting with a metal or vice versa. This tool is based on the occurrence of significant triplets and it is tested to have higher prediction accuracy when compared to that of other available techniques.