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.     

Metal ions bound at the active site of the junction-resolving enzyme T7 endonuclease I

T7 endonuclease I is a nuclease that is selective for the structure of the four-way DNA junction. The active site is similar to those of a number of restriction enzymes. We have solved the crystal structure of endonuclease I with a wild-type active site. Diffusion of manganese ions into the crystal revealed two peaks of electron density per active site, defining two metal ion-binding sites. Site 1 is fully occupied, and the manganese ion is coordinated by the carboxylate groups of Asp55 and Glu65, and the main chain carbonyl of Thr66. Site 2 is partially occupied, and the metal ion has a single protein ligand, the remaining carboxylate oxygen atom of Asp55. Isothermal titration calorimetry showed the sequential exothermic binding of two manganese ions in solution, with dissociation constants of 0.58 +/- 0.019 and 14 +/- 1.5 mM. These results are consistent with a two metal ion mechanism for the cleavage reaction, in which the hydrolytic water molecule is contained in the first coordination sphere of the site 1-bound metal ion.     

Structural and dynamical properties of manganese catalase and the synthetic protein DF1 and their implication for reactivity from classical molecular dynamics calculations

There is a pressing need for accurate force fields to assist in metalloprotein analysis and design. Recent work on the design of mimics of dimetal proteins highlights the requirements for activity. DF1 is a de novo designed protein, which mimics the overall fold and active site geometry of a series of diiron and dimanganese proteins. Specifically, the dimanganese form of DF1 is a mimic of the natural enzyme manganese catalase, which catalyzes the dismutation reaction of hydrogen peroxide into water and oxygen. During catalytic turnover, the active site has to accommodate both the reduced and the oxidized state of the dimanganese core. The biomimetic protein DF1 is only stable in the reduced form and thus not active. Furthermore, the synthetic protein features an additional bridging glutamate sidechain, which occupies the substrate binding site. The goal of this study is to develop classical force fields appropriate for design of such important dimanganese proteins. To this aim, we use a nonbonded model to represent the metal-ligand interactions, which implicitly takes into account charge transfer and local polarization effects between the metal and its ligands. To calibrate this approach, we compare and contrast geometric and dynamical properties of manganese catalase and DF1. Having demonstrated a good correspondence with experimental structural data, we examine the effect of mutating the bridging glutamate to aspartate (M1) and serine (M2). Classical MD based on the refined forcefield shows that these point mutations affect not only the immediate coordination sphere of the manganese ions, but also the relative position of the helices, improving the similarity to Mn-catalase, especially in case of M2. On the basis of these findings, classical molecular dynamics calculations with the active site parameterization scheme introduced herein seem to be a promising addition to the protein design toolbox.     

Ni-Zn-[Fe4-S4] and Ni-Ni-[Fe4-S4] clusters in closed and open subunits of acetyl-CoA synthase/carbon monoxide dehydrogenase

The crystal structure of the tetrameric alpha2beta2 acetyl-coenzyme A synthase/carbon monoxide dehydrogenase from Moorella thermoacetica has been solved at 1.9 A resolution. Surprisingly, the two alpha subunits display different (open and closed) conformations. Furthermore, X-ray data collected from crystals near the absorption edges of several metal ions indicate that the closed form contains one Zn and one Ni at its active site metal cluster (A-cluster) in the alpha subunit, whereas the open form has two Ni ions at the corresponding positions. Alternative metal contents at the active site have been observed in a recent structure of the same protein in which A-clusters contained one Cu and one Ni, and in reconstitution studies of a recombinant apo form of a related acetyl-CoA synthase. On the basis of our observations along with previously reported data, we postulate that only the A-clusters containing two Ni ions are catalytically active.     

X-ray structure of the Escherichia coli periplasmic 5'-nucleotidase containing a dimetal catalytic site.

The crystal structure of 5'-nucleotidase (5'-NT) from E. coli, also known as UDP-sugar hydrolase, has been determined at 1.7 A resolution. Two zinc ions are present in the active site, which is located in a cleft between two domains. The dimetal center and a catalytic Asp-His dyad are the main players in the catalytic mechanism. Structure-based sequence comparisons show that the structure also provides a model for animal 5'-NTs, which together with other ectonucleotidases terminate the action of nucleotides as extracellular signaling substances in the nervous system.     

Nuclear magnetic resonance studies of the internal dynamics in Apo, (Cd2+)1 and (Ca2+)2 calbindin D9k. The rates of amide proton exchange with solvent.

The backbone dynamics of the EF-hand Ca(2+)-binding protein, calbindin D9k, has been investigated in the apo, (Cd2+)1 and (Ca2+)2 states by measuring the rate constants for amide proton exchange with solvent. 15N-1H correlation spectroscopy was utilized to follow direct 1H-->2H exchange of the slowly exchanging amide protons and to follow indirect proton exchange via saturation transfer from water to the rapidly exchanging amide protons. Plots of experimental rate constants versus intrinsic rate constants have been analyzed to give qualitative insight into the opening modes of the protein that lead to exchange. These results have been interpreted within the context of a progressive unfolding model, wherein hydrophobic interactions and metal chelation serve to anchor portions of the protein, thereby damping fluctuations and retarding amide proton exchange. The addition of Ca2+ or Cd2+ was found to retard the exchange of many amide protons observed to be in hydrogen-bonding environments in the crystal structure of the (Ca2+)2 state, but not of those amide protons that were not involved in hydrogen bonds. The largest changes in rate constant occur for residues in the ion-binding loops, with substantial effects also found for the adjacent residues in helices I, II and III, but not helix IV. The results are consistent with a reorganization of the hydrogen-bonding networks in the metal ion-binding loops, accompanied by a change in the conformation of helix IV, as metal ions are chelated. Further analysis of the results obtained for the three states of metal occupancy provides insight into the nature of the changes in conformational fluctuations induced by ion binding.     

The crystal structure of the carbohydrate-recognition domain of the glycoprotein sorting receptor p58/ERGIC-53 reveals an unpredicted metal-binding site and conformational changes associated with calcium ion binding

p58/ERGIC-53 is a calcium-dependent animal lectin that acts as a cargo receptor, binding to a set of glycoproteins in the endoplasmic reticulum (ER) and transporting them to the Golgi complex. It is similar in structure to calcium-dependent leguminous lectins. We have determined the structure of the carbohydrate-recognition domain of p58/ERGIC-53 in its calcium-bound form. The structure reveals localized but large conformational changes in relation to the previously determined metal ion-free structure, mapping mostly to the ligand-binding site. It reveals the presence of two calcium ion-binding sites located 6A apart, one of which has no equivalent in the plant lectins. The second metal ion-binding site present in that class of lectins, binding Mn(2+), is absent from p58/ERGIC-53. The absence of a short loop in the ligand-binding site in this protein suggests that it has adapted to optimally bind the high-mannose Man(8)(GlcNAc)(2) glycan common to glycoproteins at the ER exit stage.     

Significant change in the structure of a ribozyme upon introduction of a phosphorothioate linkage at P9: NMR reveals a conformational fluctuation in the core region of a hammerhead ribozyme

A modified hammerhead ribozyme (R32S) with a phosphorothioate linkage between G(8) and A(9), a site that is considered to play a crucial role in catalysis, was examined by high-resolution 1H and (31)P nuclear magnetic resonance (NMR) spectroscopy. Signals due to imino protons that corresponded to stems were observed, but the anticipated signals due to imino protons adjacent to the phosphorothioate linkage were not detected and the (31)P signal due to the phosphorothioate linkage was also absent irrespective of the presence or absence of the substrate. (31)P NMR is known to reflect backbone mobility, and thus the absence of signals indicated that the introduction of sulfur at P9 had increased the mobility of the backbone near the phosphorothioate linkage. The addition of metal ions did not regenerate the signals that had disappeared, a result that implied that the structure of the core region of the hammerhead ribozyme had fluctuated even in the presence of metal ions. Furthermore, kinetic analysis suggested that most of the R32S-substrate complexes generated in the absence of Mg(2+) ions were still in an inactive form and that Mg(2+) ions induced a further conformational change that converted such complexes to an activated state. Finally, according to available NMR studies, signals due to the imino protons of the central core region that includes the P9 metal binding site were broadened or not observed, suggesting that this catalytically important region might be intrinsically flexible. Our present analysis revealed a significant change in the structure of the ribozyme upon the introduction of the single phosphorothioate linkage at P9 that is in general considered to be a conservative modification.     

Spectroscopic and metal-binding properties of DF3: an artificial protein able to accommodate different metal ions

The design, synthesis, and metal-binding properties of DF3, a new de novo designed di-iron protein model are described (“DF” represents due ferri, Italian for “two iron,” “di-iron”). DF3 is the latest member of the DF family of synthetic proteins. They consist of helix–loop–helix hairpins, designed to dimerize and form an antiparallel four-helix bundle that encompasses a metal-binding site similar to those of non-heme carboxylate-bridged di-iron proteins. Unlike previous DF proteins, DF3 is highly soluble in water (up to 3 mM) and forms stable complexes with several metal ions (Zn, Co, and Mn), with the desired secondary structure and the expected stoichiometry of two ions per protein. UV–vis studies of Co(II) and Fe(III) complexes confirm a metal-binding environment similar to previous di-Co(II)- and di-Fe(III)-DF proteins, including the presence of a μ-oxo-di-Fe(III) unit. Interestingly, UV–vis, EPR, and resonance Raman studies suggest the interaction of a tyrosine adjacent to the di-Fe(III) center. The design of DF3 was aimed at increasing the accessibility of small molecules to the active site of the four-helix bundle. Indeed, binding of azide to the di-Fe(III) site demonstrates a more accessible metal site compared with previous DFs. In fact, fitting of the binding curve to the Hill equation allows us to quantify a 150% accessibility enhancement, with respect to DF2. All these results represent a significant step towards the development of a functional synthetic DF metalloprotein.     

Binding of calcium in the EF-hand of Escherichia coli lytic transglycosylase Slt35 is important for stability

The Escherichia coli lytic transglycosylase Slt35 contains a single metal ion-binding site that resembles EF-hand calcium-binding sites. The Slt35 EF-hand is only the second observation of such a domain in a prokaryotic protein. Two crystal structures at 2.1 A resolution show that both Ca2+ ions and Na+ ions can bind to the EF-hand domain, but in subtly different configurations. Heat-induced unfolding studies demonstrate that Ca2+ ions are preferentially bound, and that only Ca2+ ions significantly increase the melting temperature of Slt35. This shows that the EF-hand calcium-binding domain is important for the stability of Slt35.     

Crystal structure and mechanism of catalysis of a pyrazinamidase from Pyrococcus horikoshii.

Bacterial pyrazinamidase (PZAase)/nicotinamidase converts pyrazinamide (PZA) to ammonia and pyrazinoic acid, which is active against Mycobacterium tuberculosis. Loss of PZAase activity is the major mechanism of pyrazinamide-resistance by M. tuberculosis. We have determined the crystal structure of the gene product of Pyrococcus horikoshii 999 (PH999), a PZAase, and its complex with zinc ion by X-ray crystallography. The overall fold of PH999 is similar to that of N-carbamoylsarcosine amidohydrolase (CSHase) of Arthrobacter sp. and YcaC of Escherichia coli, a protein with unknown physiological function. The active site of PH999 was identified by structural features that are also present in the active sites of CSHase and YcaC: a triad (D10, K94, and C133) and a cis-peptide (between V128 and A129). Surprisingly, a metal ion-binding site was revealed in the active site and subsequently confirmed by crystal structure of PH999 in complex with Zn(2+). The roles of the triad, cis-peptide, and metal ion in the catalysis are proposed. Because of extensive homology between PH999 and PZAase of M. tuberculosis (37% sequence identity), the structure of PH999 provides a structural basis for understanding PZA-resistance by M. tuberculosis harboring PZAase mutations.     

Structural and binding studies of the three-metal center in two mycobacterial PPM Ser/Thr protein phosphatases

Phospho-Ser/Thr protein phosphatases (PPs) are dinuclear metalloenzymes classed into two large families, PPP and PPM, on the basis of sequence similarity and metal ion dependence. The archetype of the PPM family is the α isoform of human PP2C (PP2Cα), which folds into an α/β domain similar to those of PPP enzymes. The recent structural studies of three bacterial PPM phosphatases, Mycobacterium tuberculosis MtPstP, Mycobacterium smegmatis MspP, and Streptococcus agalactiae STP, confirmed the conservation of the overall fold and dinuclear metal center in the family, but surprisingly revealed the presence of a third conserved metal-binding site in the active site. To gain insight into the roles of the three-metal center in bacterial enzymes, we report structural and metal-binding studies of MtPstP and MspP. The structure of MtPstP in a new trigonal crystal form revealed a fully active enzyme with the canonical dinuclear metal center but without the third metal ion bound to the catalytic site. The absence of metal correlates with a partially unstructured flap segment, indicating that the third manganese ion contributes to reposition the flap, but is dispensable for catalysis. Studies of metal binding to MspP using isothermal titration calorimetry revealed that the three Mn²⁺-binding sites display distinct affinities, with dissociation constants in the nano- and micromolar range for the two catalytic metal ions and a significantly lower affinity for the third metal-binding site. In agreement, the structure of inactive MspP at acidic pH was determined at atomic resolution and shown to lack the third metal ion in the active site. Structural comparisons of all bacterial phosphatases revealed positional variations in the third metal-binding site that are correlated with the presence of bound substrate and the conformation of the flap segment, supporting a role of this metal ion in assisting enzyme-substrate interactions.     

Artificial di-iron proteins: solution characterization of four helix bundles containing two distinct types of inter-helical loops

Peptide-based models have an enormous impact for the development of metalloprotein models, as they seem appropriate candidates to mimic both the structural characteristics and reactivity of the natural systems. Through the de novo design of four-helix bundles, we developed the DF (Due Ferri) family of artificial proteins, as models of di-iron and di-manganese metalloproteins. The goal of our research is to elucidate how the electrostatic environment, polarity and solvent accessibility of the metal-binding site, influence the functional properties of di-iron proteins. The first two subsets of the DF protein family, DF1 and DF2, consist of two non-covalently associated helix-loop-helix motifs, which bind the di-metal cofactor near the center of the structure. The DF2 subset was designed to improve the properties of DF1: DF2 and DF2t have several changes in their sequences to improve solubility and metal ion access, as well as a change in the loop connecting the two helices. In order to evaluate how these changes affect the overall structure of the model proteins, we solved the NMR structures of the di-Zn(II) complexes of DF2 and DF2t, and compared these structures with those recently obtained from X-ray crystallography. Further, we examined the thermodynamic consequences associated with the mutations, by measuring the stability of DF2t in the presence of different metal ions, and comparing the results with the data already obtained for DF2. Taken together, analysis of all the data showed the importance of the turn conformation in the design and stability of four-helix bundle.Electronic Supplementary Material Supplementary material is available for this article at     

Crystallographic and X-ray absorption spectroscopic characterization of Helicobacter pylori UreE bound to Ni²⁺ and Zn²⁺ reveals a role for the disordered C-terminal arm in metal trafficking

The survival and growth of the pathogen Helicobacter pylori in the gastric acidic environment is ensured by the activity of urease, an enzyme containing two essential Ni2? ions in the active site. The metallo-chaperone UreE facilitates in vivo Ni2? insertion into the apoenzyme. Crystals of apo-HpUreE (H. pylori UreE) and its Ni?- and Zn?-bound forms were obtained from protein solutions in the absence and presence of the metal ions. The crystal structures of the homodimeric protein, determined at 2.00 ? (apo), 1.59 ? (Ni2?) and 2.52 ? (Zn2?) resolution, show the conserved proximal and solvent-exposed His1?2 residues from two adjacent monomers invariably involved in metal binding. The C-terminal regions of the apoprotein are disordered in the crystal, but acquire significant ordering in the presence of the metal ions due to the binding of His1?2. The analysis of X-ray absorption spectral data obtained using solutions of Ni2?- and Zn2?-bound HpUreE provided accurate information of the metal-ion environment in the absence of solid-state effects. These results reveal the role of the histidine residues at the protein C-terminus in metal-ion binding, and the mutual influence of protein framework and metal-ion stereo-electronic properties in establishing co-ordination number and geometry leading to metal selectivity.     

Structure of the bacteriophage lambda Ser/Thr protein phosphatase with sulfate ion bound in two coordination modes

The protein phosphatase encoded by bacteriophage lambda (lambda PP) belongs to a family of Ser/Thr phosphatases (Ser/Thr PPases) that includes the eukaryotic protein phosphatases 1 (PP1), 2A (PP2A), and 2B (calcineurin). These Ser/Thr PPases and the related purple acid phosphatases (PAPs) contain a conserved phosphoesterase sequence motif that binds a dinuclear metal center. The mechanisms of phosphoester hydrolysis by these enzymes are beginning to be unraveled. To utilize lambda PP more effectively as a model for probing the catalytic mechanism of the Ser/Thr PPases, we have determined its crystal structure to 2.15 A resolution. The overall fold resembles that of PP1 and calcineurin, including a conserved beta alpha beta alpha beta structure that comprises the phosphoesterase motif. Substrates and inhibitors probably bind in a narrow surface groove that houses the active site dinuclear Mn(II) center. The arrangement of metal ligands is similar to that in PP1, calcineurin, and PAP, and a bound sulfate ion is present in two novel coordination modes. In two of the three molecules in the crystallographic asymmetric unit, sulfate is coordinated to Mn2 in a monodentate, terminal fashion, and the two Mn(II) ions are bridged by a solvent molecule. Two additional solvent molecules are coordinated to Mn1. In the third molecule, the sulfate ion is triply coordinated to the metal center with one oxygen coordinated to both Mn(II) ions, one oxygen coordinated to Mn1, and one oxygen coordinated to Mn2. The sulfate in this coordination mode displaces the bridging ligand and one of the terminal solvent ligands. In both sulfate coordination modes, the sulfate ion is stabilized by hydrogen bonding interactions with conserved arginine residues, Arg 53 and Arg 162. The two different active site structures provide models for intermediates in phosphoester hydrolysis and suggest specific mechanistic roles for conserved residues.     

Solution structure and thermodynamics of a divalent metal ion binding site in an RNA pseudoknot.

Identification and characterization of a metal ion binding site in an RNA pseudoknot was accomplished using cobalt (III) hexammine, Co(NH3)63+, as a probe for magnesium (II) hexahydrate, Mg(H2O)62+, in nuclear magnetic resonance (NMR) structural studies. The pseudoknot causes efficient -1 ribosomal frameshifting in mouse mammary tumor virus. Divalent metal ions, such as Mg2+, are critical for RNA structure and function; Mg2+preferentially stabilizes the pseudoknot relative to its constituent hairpins. The use of Co(NH3)63+as a substitute for Mg2+was investigated by ultraviolet absorbance melting curves, NMR titrations of the imino protons, and analysis of NMR spectra in the presence of Mg2+or Co (NH3)63+. The structure of the pseudoknot-Co(NH3)63+complex reveals an ion-binding pocket formed by a short, two-nucleotide loop and the major groove of a stem. Co(NH3)63+stabilizes the sharp loop-to-stem turn and reduces the electrostatic repulsion of the phosphates in three proximal strands. Hydrogen bonds are identified between the Co(NH3)63+protons and non-bridging phosphate oxygen atoms, 2' hydroxyl groups, and nitrogen and oxygen acceptors on the bases. The binding site is significantly different from that previously characterized in the major groove surface of tandem G.U base-pairs, but is similar to those observed in crystal structures of a fragment of the 5 S rRNA and the P5c helix of the Tetrahymena thermophila group I intron. Changes in chemical shifts occurred at the same pseudoknot protons on addition of Mg2+as on addition of Co(NH3)63+, indicating that both ions bind at the same site. Ion binding dissociation constants of approximately 0.6 mM and 5 mM (in 200 mM Na+and a temperature of 15 degrees C) were obtained for Co(NH3)63+and Mg2+, respectively, from the change in chemical shift as a function of metal ion concentration. An extensive array of non-sequence-specific hydrogen bond acceptors coupled with conserved structural elements within the binding pocket suggest a general mode of divalent metal ion stabilization of this type of frameshifter pseudoknot. These results provide new thermodynamic and structural insights into the role divalent metal ions play in stabilizing RNA tertiary structural motifs such as pseudoknots.     

Structural and functional characterization of a novel phosphodiesterase from Methanococcus jannaschii

Methanococcus jannaschii MJ0936 is a hypothetical protein of unknown function with over 50 homologs found in many bacteria and Archaea. To help define the molecular (biochemical and biophysical) function of MJ0936, we determined its crystal structure at 2.4-A resolution and performed a series of biochemical screens for catalytic activity. The overall fold of this single domain protein consists of a four-layered structure formed by two beta-sheets flanked by alpha-helices on both sides. The crystal structure suggested its biochemical function to be a nuclease, phosphatase, or nucleotidase, with a requirement for some metal ions. Crystallization in the presence of Ni(2+) or Mn(2+) produced a protein containing a binuclear metal center in the putative active site formed by a cluster of conserved residues. Analysis of MJ0936 against a panel of general enzymatic assays revealed catalytic activity toward bis-p-nitrophenyl phosphate, an indicator substrate for phosphodiesterases and nucleases. Significant activity was also found with two other phosphodiesterase substrates, thymidine 5'-monophosphate p-nitrophenyl ester and p-nitrophenylphosphorylcholine, but no activity was found for cAMP or cGMP. Phosphodiesterase activity of MJ0936 had an absolute requirement for divalent metal ions with Ni(2+) and Mn(2+) being most effective. Thus, our structural and enzymatic studies have identified the biochemical function of MJ0936 as that of a novel phosphodiesterase.     

Structural studies of the Alzheimer's amyloid precursor protein copper-binding domain reveal how it binds copper ions

Alzheimer's disease (AD) is the major cause of dementia. Amyloid beta peptide (Abeta), generated by proteolytic cleavage of the amyloid precursor protein (APP), is central to AD pathogenesis. APP can function as a metalloprotein and modulate copper (Cu) transport, presumably via its extracellular Cu-binding domain (CuBD). Cu binding to the CuBD reduces Abeta levels, suggesting that a Cu mimetic may have therapeutic potential. We describe here the atomic structures of apo CuBD from three crystal forms and found they have identical Cu-binding sites despite the different crystal lattices. The structure of Cu(2+)-bound CuBD reveals that the metal ligands are His147, His151, Tyr168 and two water molecules, which are arranged in a square pyramidal geometry. The site resembles a Type 2 non-blue Cu center and is supported by electron paramagnetic resonance and extended X-ray absorption fine structure studies. A previous study suggested that Met170 might be a ligand but we suggest that this residue plays a critical role as an electron donor in CuBDs ability to reduce Cu ions. The structure of Cu(+)-bound CuBD is almost identical to the Cu(2+)-bound structure except for the loss of one of the water ligands. The geometry of the site is unfavorable for Cu(+), thus providing a mechanism by which CuBD could readily transfer Cu ions to other proteins.