Structural characteristics of protein binding sites for calcium and lanthanide ions

Surveys of X-ray structures of Ca2+-containing and lanthanide ion-containing proteins and coordination complexes have been performed and structural features of the metal binding sites compared. A total of 515 structures of Ca2+-containing proteins were considered, although the final data set contained only 44 structures and 60 Ca2+ binding sites with a total of 323 ligands. Eighteen protein structures containing lanthanide ions were considered with a final data set containing eight structures and 11 metal binding sites. Structural features analysed include coordination numbers of the metal ions, the identity of their ligands, the denticity of carboxylate ligands, and the type of secondary structure from which the ligands are derived. Three general types of calcium binding site were identified in the final data set: class I sites supply the Ca2+ ligands from a continuous short sequence of amino acids; class II sites have one ligand supplied by a part of the amino acid sequence far removed from the main binding sequence; and class III sites are created by amino acids remote from one another in the sequence. The abundant EF-hand type of Ca2+ binding site was under-represented in the data set of structures analysed as far as its biological distribution is concerned, but was adequately represented for the chemical survey undertaken. A turn or loop structure was found to provide the bulk of the ligands to Ca2+, but helix and sheet secondary structures are slightly better providers of bidentate carboxylate ligation than turn or loop structures. The average coordination number for Ca2+ was 6.0, though for EF-hand sites it is 7. The average coordination number of a lanthanide ion in an intrinsic protein Ca2+ site was 7.2, but for the adventitious sites was only 4.4. A survey of the Cambridge Structural Database showed there are small-molecule lanthanide complexes with low coordination numbers but it is likely that water molecules, which do not appear in the electron density maps, are present for some lanthanide sites in proteins. A detailed comparison of the well-defined Ca2+ and lanthanide ion binding sites suggests that a reduction of hydrogen bonding associated with the ligating residues of the binding sites containing lanthanide ions may be a response to the additional positive charge of the lanthanide ion. Major structural differences between Ca2+ binding sites with weak and strong binding affinities were not obvious, a consequence of long-range electrostatic interactions and metal ion-induced protein conformational changes modulating affinities.     

Geometry of interaction of metal ions with histidine residues in protein structures

An analysis of the geometry and the orientation of metal ions bound to histidine residues in proteins is presented. Cations are found to lie in the imidazole plane along the lone pair on the nitrogen atom. Out of the two tautomeric forms of the imidazole ring, the NE2-protonated form is normally preferred. However, when bound to a metal ion the ND1-protonated form is predominant and NE2 is the ligand atom. When the metal coordination is through ND1, steric interactions shift the side chain torsional angle, chi 2 from its preferred value of 90 or 270 degrees. The orientation of histidine residues is usually stabilized through hydrogen bonding; ND1-protonated form of a helical residue can form a hydrogen bond with the carbonyl oxygen atom in the preceding turn of the helix. A considerable number of ligands are found in helices and beta-sheets. A helical residue bound to a heme group is usually found near the C-terminus of the helix. Two ligand groups four residues apart in a helix, or two residues apart in a beta-strand are used in many proteins to bind metal ions.     

Conformational studies of A23187 with mono-, di- and trivalent metal ions by circular dichroism spectroscopy

CD studies carried out on A23187 indicate a solvent-dependent conformation for the free acid. Alkali metal ions were found to bind to the ionophore weakly. Divalent metal ions such as Mg2+, Ca2+, Sr2+, Ba2+ and Co2+ and trivalent lanthanide metal ions like La3+ were found to form predominantly 2:1 (ionophore-metal ion) complexes at low concentrations of metal ions, but both 2:1 and 1:1 complexes were formed with increasing salt concentration. Mg2+ and Co2+ exhibit similar CD behaviour that differs from that observed for the other divalent and lanthanide metal ions. The structure of 2:1 complexes involves two ligand molecules coordinated to the metal ion through the carboxylate oxygen, benzoxazole nitrogen and keto-pyrrole oxygen from each ligand molecule along with one or more solvent molecules. Values of the binding constant were determined for 2:1 complexes of the ionophore with divalent and lanthanide metal ions.     

Carboxylate binding modes in zinc proteins: A theoretical study

The relative energies of different coordination modes (bidentate, monodentate, syn, and anti) of a carboxylate group bound to a zinc ion have been studied by the density functional method B3LYP with large basis sets on realistic models of the active site of several zinc proteins. In positively charged four-coordinate complexes, the mono- and bidentate coordination modes have almost the same energy (within 10 kJ/mol). However, if there are negatively charged ligands other than the carboxylate group, the monodentate binding mode is favored. In general, the energy difference between monodentate and bidentate coordination is small, 4-24 kJ/mol, and it is determined more by hydrogen-bond interactions with other ligands or second-sphere groups than by the zinc-carboxylate interaction. Similarly, the activation energy for the conversion between the two coordination modes is small, approximately 6 kJ/mol, indicating a very flat Zn-O potential surface. The energy difference between syn and anti binding modes of the monodentate carboxylate group is larger, 70-100 kJ/mol, but this figure again strongly depends on interactions with second-sphere molecules. Our results also indicate that the pK(a) of the zinc-bound water ligand in carboxypeptidase and thermolysin is 8-9.     

Binding ability of a HHP-tagged protein towards Ni2+ studied by paramagnetic NMR relaxation: The possibility of obtaining long-range structure information

The binding ability of a protein with a metal binding tag towards Ni(2+) was investigated by longitudinal paramagnetic NMR relaxation, and the possibility of obtaining long-range structure information from the paramagnetic relaxation was explored. A protein with a well-defined solution structure (Escherichia coli thioredoxin) was used as the model system, and the peptide His-His-Pro (HHP) fused to the N-terminus of the protein was used as the metal binding tag. It was found that the tag forms a stable dimer complex with the paramagnetic Ni(2+) ion, where each metal ion binds two HHP-tagged protein molecules. However, it was also found that additional sites in the protein compete with the HHP-tag for the binding of the metal ion. These binding sites were identified as the side chain carboxylate groups of the aspartic and glutamic acid residues. Yet, the carboxylate groups bind the Ni(2+) ions considerably weaker than the HHP-tag, and only protons spatially close to the carboxylate sites are affected by the Ni(2+) ions bound to these groups. As for the protons that are unaffected by the carboxylate-bound Ni(2+) ions, it was found that the long-range distances derived from the paramagnetic relaxation enhancements are in good agreement with the solution structure of thioredoxin. Specifically, the obtained long-range paramagnetic distance constraints revealed that the dimer complex is asymmetric with different orientations of the two protein molecules relative to the Ni(2+) ion.     

On the coordination properties of Eu3+ bound to tRNA

The luminescence properties of Eu3+ have been used to investigate the binding and coordination properties of the ion with tRNA, as an attempt to resolve the discussion of whether metal ions bind to tRNA in solution only by Debye-Hückel screening, or whether direct coordination to specific sites may occur. Binding studies with Escherichia coli tRNAmet/f (taking advantage of 4-thiouracil-sensitized Eu3+ emission) distinguish three classes of binding affinities. Two of these are single sites with affinities approx. 10(4) and approx. 10(3) tighter than the nonspecific affinity of Eu3+ for native DNA. Mg2+ competes for binding at both these sites. Measurement of the lifetime and excitation spectrum of Eu3+ bound to the highest affinity site shows that the ion has two to five non-phosphate ligands in its inner coordination sphere. The existence of this coordinated site demonstrates that electrostatic screening is not the only mechanism for metal ion interaction with tRNA. The coordination properties of the high-affinity Eu3+ site do not agree with the properties of any of the metal ion sites found in the two tRNAphe crystal forms. Possible reasons for this discrepancy are discussed; it may be that ions bind differently to isolated molecules in solution than to molecules packed in a crystal lattice.     

Conformation and intramolecular hydrogen-bonding in the crystal structure of potassium D-gluconate monohydrate

The D-gluconate ion is found to have the planar, extended carbon-chain conformation in the crystal structure of potassium D-gluconate monohydrate, with an intramolecular hydrogen-bond between 0-2 and 0-4. The D-gluconate ions and water molecules are linked in puckered sheets by a series of intermolecular hydrogen-bonds that involve the water molecules, the carboxylate groups, and pairs of hydroxyl groups. One hydroxyl group in the ion does not form a hydrogen bond. The potassium ions lie between the puckered sheets, with an eight-fold coordination of six D-gluconate groups and two water oxygen atoms. The crystal structure was determined from three-dimensional, CuKα, X-ray diffraction data taken on an automatic diffractometer.     

Identification of the six ligands to manganese(II) in transition-state-analogue complexes of creatine kinase: oxygen-17 superhyperfine coupling from selectively labeled ligands

The complete coordination scheme for Mn(II) in transition-state-analogue complexes with creatine kinase has been determined by electron paramagnetic resonance (EPR) spectroscopy. Perturbations in the EPR spectra for Mn(II) due to superhyperfine coupling to 17O of selectively labeled ligands have been used to identify oxygen ligands in the first coordination sphere of the metal ion. The results show that in the complex of enzyme-MnADP-formate-creatine, Mn(II) is bound to oxygen ligands from both the alpha- and beta-phosphate groups of ADP, to an oxygen from the carboxylate group of formate, and to three water molecules. In the complex with thiocyanate replacing formate as the stabilizing anion, previous infrared experiments [Reed, G. H., Barlow, C. H., & Burns, R. A., Jr. (1978) J. Biol. Chem. 253, 4153-4158] indicated that the nitrogen from thiocyanate was bound to the Mn(II). The magnitudes of the 17O superphyperfine coupling constants from the O- ligands of the ADP phosphate groups and from the formate carboxylate are approximately equal and are larger than that for the water ligands. The symmetry of the zero-field-splitting tensor for Mn(II) indicates that the oxygens from the alpha- and beta-phosphate groups of ADP and the ligand donor atom from the anion occupy mutually cis positions in the octahedral coordination geometry. Water proton relaxation time measurements show that the three water molecules which are bound to Mn(II) are not in free exchange with the bulk solvent. Hence, an enclosed structure at the active site is indicated. The results suggest that for creatine kinase the activating metal ion is bound to all three phosphate groups in the transition state of the reaction.     

Metal substitution in transferrins: specific binding of cerium(IV) revealed by the crystal structure of cerium-substituted human lactoferrin

Proteins of the transferrin family play a key role in iron homeostasis through their extremely strong binding of iron, as Fe3+. They are nevertheless able to bind a surprisingly wide variety of other metal ions. To investigate how metal ions of different size, charge and coordination characteristics are accommodated, we have determined the crystal structure of human lactoferrin (Lf) complexed with Ce4+. The structure, refined at 2.2 A resolution (R=20.2%, Rfree=25.7%) shows that the two Ce4+ ions occupy essentially the same positions as do Fe3+, and that the overall protein structure is unchanged; the same closed structure is formed for Ce2Lf as for Fe2Lf. The larger metal ion is accommodated by small shifts in the protein ligands, made possible by the presence of water molecules adjacent to each binding site. The two Ce4+ sites are equally occupied, indicating that the known difference in the pH-dependent release of Ce4+ arises from a specific protonation event, possibly of the His ligand in one of the binding sites. Comparing the effects of binding Ce4+ with those for the binding of other metal ions, we conclude that the ability of transferrins to accommodate metal ions other than Fe3+ depends on an interplay of charge, size, coordination and geometrical preferences of the bound metal ion. However, it is the ability to accept the six-coordinate, approximately octahedral, site provided by the protein that is of greatest importance.     

Stereochemistry of guanidine-metal interactions: implications for L-arginine-metal interactions in protein structure and function

The geometries of 150 guanidine-metal ion interactions retrieved from crystal structures deposited in the Cambridge Structural Database have been analyzed. Metal ions exhibit a preference for anti coordination stereochemistry in the plane of the unprotonated guanidine group, usually in chelate complexes with a diguanidine moiety, but syn-oriented interactions are occasionally found for single guanidine-metal interactions. Three L-arginine-metal coordination interactions are found in metalloenzyme structures deposited in the Protein Data Bank: biotin synthase from E. coli, His-67 --> Arg human carbonic anhydrase I, and inactivated B. caldovelox arginase complexed with L-arginine. In these proteins, L-arginine-metal coordination adopts syn/out-of-plane and anti/in-plane coordination stereochemistry. The implications of these results for L-arginine-metal interactions in protein structure and function are discussed. Although such interactions are rare, this analysis serves as a useful reference point for the growing interest in enzymes containing L-arginine residues that function as general bases or metal ligands.     

Structure of satellite tobacco necrosis virus after crystallographic refinement at 2.5 Å resolution

The structure of the protein subunit of satellite tobacco necrosis virus has been solved at 3.7 Å resolution. We have now crystallographically refined the original model and extended the resolution to 2.5 Å in order to get a model accurate enough to explain the details of the subunit interactions. The refinement was done with a novel method utilizing the icosahedral symmetry of the virus particle.The final model shows a complicated network of interactions, involving salt linkages, hydrogen bonds and hydrophobic contacts. In addition, we have located three different metal ion sites in the protein shell, linking the protein subunits together. These sites are probably occupied by calcium ions. One site is found in a general position near the icosahedral 3-fold axis of the virus. The ligands form an octahedral arrangement, with two main chain carbonyl oxygens (O-61 and O-64), one carboxylate oxygen (OD1 from Asp194) of the same subunit and a second carboxylate oxygen (OE1 of Glu25) from a 3-fold related subunit. Two water molecules complete the octahedral arrangement. A second site is on the icosahedral 3-fold axis and is liganded by the carboxylate oxygens of the 3-fold related Asp55 residues. The third metal ion site is found on the 5-fold axis, liganded by the five carbonyl oxygens of Thr138 and two water molecules.We are unable to locate the first 11 N-terminal amino acid residues, which point into the virus interior. No interpretable density for RNA has been found, indicating that the nucleic acid of the virus does not have a unique orientation in the crystal.     

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.     

Inhibition of HIV-1 proteinase by metal ions.

1. Certain metal ions have been identified as inhibitors (IC50 1-20 microM) of the aspartic proteinase of Human Immunodeficiency Virus Type 1 (HIV-PR). 2. By contrast most simple metal ions do not inhibit this enzyme. 3. Those that did inhibit have in common a high charge/size ratio or "hard" acidic nature, preferring to combine covalently with oxygen donor ligands. 4. Some evidence from independent X-ray crystal structure determinations suggests that the metalloinhibitors identified here may bind in the active site of the enzyme via coordination to the carboxylate side chains of the essential active site residues Asp 25 and 125. 5. Although the measured inhibition is only microM, very few enzyme-inhibitor interactions can be taking place and so more complex metalloinhibitors with ligands that can also bind to peptide side chains of the enzyme might be significantly more potent inhibitors of HIV-PR and of viral replication.     

Structure-based mechanism of ligand binding for periplasmic solute-binding protein of the Bug family

Bug proteins form a large family of periplasmic solute-binding proteins well represented in beta-proteobacteria. They adopt a characteristic Venus flytrap fold with two globular domains bisected by a ligand-binding cleft. The structures of two liganded Bug proteins have revealed that the family is specific for carboxylated solutes, with a characteristic mode of binding involving two highly conserved beta strand-beta turn-alpha helix motifs originating from each domain. These two motifs form hydrogen bonds with a carboxylate group of the ligand, both directly and via conserved water molecules, and have thus been termed the carboxylate pincers. In both crystallized Bug proteins, the ligands were found enclosed between the two domains and inaccessible to solvent, suggesting an inter-domain hinge-bending motion upon ligand binding. We report here the first structures of an open, unliganded Bug protein and of the same protein with a citrate ion bound in the open cavity. One of the ligand carboxylate groups is bound to one half of the carboxylate pincers by the beta strand-beta turn-alpha helix motif from domain 1, and the citrate ion forms several additional interactions with domain 1. The ligand is accessible to solvent and has very few contacts with domain 2. In this open, liganded structure, the second part of the carboxylate pincers originating from domain 2 is not stabilized by ligand binding, and a loop replaces the beta turn. In the unliganded structure, both motifs of the carboxylate pincers are highly mobile, and neither of the two beta turns is formed. Thus, ligand recognition is performed by domain 1, with the carboxylate group serving as an initial anchoring point. Stabilization of the closed conformation requires proper interactions to be established with domain 2, and thus domain 2 discriminates between productively and non-productively bound ligands.     

The Ca2+ ion and membrane binding structure of the Gla domain of Ca-prothrombin fragment 1.

The structure of Ca-prothrombin fragment 1 (residues 1-156 prothrombin) has been solved and refined at 2.2-A resolution by X-ray crystallographic methods. The first two-thirds of the Gla domain (residues 1-48) and two carbohydrate chains (approximately 5 kDa) are disordered in crystals of apo-fragment 1. When crystals are grown in the presence of Ca2+ ions, the Gla domain exhibits a well-defined structure binding seven Ca2+ ions, but the carbohydrate is still disordered. Even so, the crystallographic R factor reduced to 0.171. The folding of the Gla domain is dominated by 9-10 turns of three different alpha-helices. These turns produce two internal carboxylate surfaces composed of Gla side chains. A polymeric array of five Ca2+ ions separated by about 4.0 A intercalates between the carboxylate surfaces. The coordination of the Ca2+ ions with Gla carboxylate oxygen atoms and water molecules leads to distorted polyhedral arrangements with mu-oxo bridges in a highly complex array that most likely orchestrates the folding of the domain. The overall mode of interaction of the Ca2+ ions is new and different from any Ca2+ ion-protein interactions heretofore observed or described. The fluorescence quenching event observed upon Ca2+ ion binding is due to a disulfide-pi-electron interaction that causes a 100 degrees reorientation of Trp42 of the Gla domain. The Ca2+ ion interaction also affords the N-terminus protection from acetylation because the latter is buried in the folded structure and makes hydrogen-bonding salt bridges with Gla17, Gla21, and Gla27. The Gla domain and its trailing disulfide unit associate intimately and together give rise to a domain-like structure. Electrostatic potential calculations indicate that the Gla domain is very electronegative. Since most of the carboxylate oxygen atoms of Gla residues are involved in Ca2+ ion binding, leaving only a few for bridging Ca2+ ion-phospholipid interactions, the role of bridging Ca2+ ions might be generally unspecific, with Ca2+ ions simply intervening between the negative Gla domain and negative head groups of the membrane surface. The folding of the kringle structure in apo- and Ca-fragment 1 is essentially the same. However, the Ser36-Ala47 helix of the Gla domain pivots around Cys48, shifting by approximately 30 degrees, and the helix encroaches on the kringle producing some concomitant changes. These might be related to the protection of carbohydrate carrying Asn101 from acetylation in the Ca-fragment 1 structure.     

Atomic resolution analysis of a 2:1 complex of CpG and acridine orange.

Cytidylyl-3', 5'-guanosine and acridine orange crystallize in a highly-ordered triclinic lattice which diffracts X-rays to 0.85 angstrom resolution. The crystal structure has been solved and refined to a residual factor of 9.5%. The two dinucleoside phosphate molecules form an antiparallel double helix with the acridine orange intercalated between them. The two base pairs of the double helical fragment have a twist angle of 10 degrees and it is found to have a C3' endo-(3', 5')-C2' endo mixed sugar puckering along the nucleotide backbone as has been observed for other simple intercalator complexes. Twenty-five water molecules have been located in the lattice together with a sodium ion. The intercalator double helical fragments form sheets which are held together by van der Waals interactions in one direction and hydrogen bonding interactions in the other. The crystal lattice contains aqueous channels in which sixteen water molecules are hydrogen bonded to the nucleotide, none to the intercalator, five water molecules are coordinated about the sodium ion and four water molecules bind solely to other water molecules. The bases in the base pairs have a dihedral angle of 7 to 8 degrees between them.     

Systematics in the interaction of metal ions with the main-chain carbonyl group in protein structures

An analysis of the geometry of metal binding by peptide carbonyl groups in proteins is presented. Such metal ions are predominantly calcium in known protein structures. Cations tend to be located in the peptide plane, near the C = O bond direction. This distribution differs from that observed for water molecules bound to carbonyl oxygens. Most metal ions are bound to carbonyl oxygens of peptides in turns or in regions with no regular secondary structure. More infrequent binding interactions occur at the C-terminal end of alpha-helices or at the edges and sides of beta-sheets, where the geometrical preferences of the metal-carbonyl interaction may be satisfied. In many proteins carbonyl groups that are one, two, or three residues apart along the polypeptide chain bind to the same cation; these structures show a limited number of main-chain conformations around the metal center.     

Metal-ion environment in solid Mn(II), Co(II) and Ni(II) hyaluronates

Amorphous powders and films of some metal hyaluronate complexes of general composition (C14H20O11N)2 x xH2O (M = Mn2+, Ni2+ and Co2+) have been prepared at pH 5.5-6.0. The coordination geometry around the metal ions has been analyzed by EXAFS (extended X-ray absorption fine structure) and FTIR spectroscopy. Mn2+, Ni2+, and Co2+ ions are coordinated to carboxylate oxygen atoms and water molecules. The process of local geometry formation round the metal ions is sensitive to sample preparation.     

Differential inhibition of cytosolic PEPCK by substrate analogues. Kinetic and structural characterization of inhibitor recognition

The mechanisms of molecular recognition of phosphoenolpyruvate (PEP) and oxaloacetate (OAA) by cytosolic phosphoenolpyruvate carboxykinase (cPEPCK) were investigated by the systematic evaluation of a variety of PEP and OAA analogues as potential reversible inhibitors of the enzyme against PEP. The molecules that inhibit the enzyme in a competitive fashion were found to fall into two general classes. Those molecules that mimic the binding geometry of PEP, namely phosphoglycolate and 3-phosphonopropionate, are found to bind weakly (millimolar Ki values). In contrast, those competitive inhibitors that mimic the binding of OAA (oxalate and phosphonoformate) coordinate directly to the active site manganese ion and bind an order of magnitude more tightly (micromolar Ki values). The competitive inhibitor sulfoacetate is found to be an outlier of these two classes, binding in a hybrid fashion utilizing modes of recognition of both PEP and OAA in order to achieve a micromolar inhibition constant in the absence of direct coordination to the active site metal. The kinetic studies in combination with the structural characterization of the five aforementioned competitive inhibitors demonstrate the molecular requirements for high affinity binding of molecules to the active site of the enzyme. These features include cis-planar carbonyl groups that are required for coordination to the active site metal, a bridging electron rich atom at the position corresponding to the C2 methylene group of OAA to facilitate interactions with R405, a carboxylate or sulfonate moiety at a position corresponding to the C1 carboxylate of OAA, and the edge-on aromatic interaction between a carboxylate and Y235.