Thermodynamics of metal ion binding. 2. Metal ion binding by carbonic anhydrase variants

Inhibitor binding to metal-substituted metalloenzyme: Sulfonamide affinity for carbonic anhydrase IX

Transition metal ions are structural and catalytic cofactors of many proteins including human carbonic anhydrase (CA), a Zn-dependent hydrolase. Sulfonamide inhibitors of CA recognize and form a coordination bond with the Zn ion located in the active site of the enzyme. The Zn ion may be removed or substituted with other metal ions. Such CA protein retains the structure and could serve as a tool to study metal ion role in the recognition and binding affinity of inhibitor molecules. We measured the affinities of selected divalent transition metal ions, including Mn, Fe, Co, Ni, Cu, Cd, Hg, and Zn to metal-free CA isozymes CA I, CA II, and CAIX by fluorescence-based thermal shift assay, prepared metal-substituted CAs, and determined binding of diverse sulfonamide compounds. Sulfonamide inhibitor binding to metal substituted CA followed a U-shape pH dependence. The binding was dissected to contributing binding-linked reactions and the intrinsic binding reaction affinity was calculated. This value is independent of pH and protonation reactions that occur simultaneously upon binding native CA and as demonstrated here, to metal substituted CA. Sulfonamide inhibitor binding to cancer-associated isozyme CAIX diminished in the order: Zn > Co > Hg > Cu > Cd > Mn > Ni. Energetic contribution of the inhibitor-metal coordination bond was determined for all above metals. The understanding of the principles of metal influence on ligand affinity and selectivity should help design new drugs targeting metalloenzymes.

Continuous Injection Isothermal Titration Calorimetry for In Situ Evaluation of Thermodynamic Binding Properties of Ligand-Receptor Binding Models

We utilize a continuous injection approach (CIA) rather than the traditional incremental injection approach (IIA) to deliver ligand (or receptor) to the calorimeter cell to evaluate thermodynamic binding parameters for three common ligand-receptor binding models-single independent, competitive, and two independent binding sites-using isothermal titration calorimetry (ITC). A general mathematical expression for the binding isotherm for any binding stoichiometry under continuous delivery of ligand (or receptor) resulting in an analytical solution for the thermodynamic binding parameters is presented. The advantages of CIA include reduction in experimental time, estimation of thermodynamic binding parameter values, and automation of the experiment since thermodynamic parameters are estimated in situ. We demonstrate the inherent advantages of CIA over IIA for the three binding models. For the single independent site model, we utilized the binding of Ba²⁺ ions to ethylenediaminetetraacetic acid (EDTA), while competitive binding was captured by titration of Ca²⁺ ions into a buffered solution of Ba²⁺ and EDTA. We experimentally simulated a two independent binding site system by injecting Ca²⁺ into a solution of EDTA and 1,3-diaminopropane-N,N,N',N'-tetraacetic acid (DPTA). The results demonstrate estimation of thermodynamic parameters with greater confidence and simultaneous reduction in the experimental time of 83% and titrating reagent of 50%, as compared to IIA.

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.

Thermal and chemical unfolding pathways of PaSdsA1 sulfatase, a homo-dimer with topologically interlinked chains

Understanding the mechanisms as to how interlinked proteins entangle and fold is a challenge. PaSdsA1 sulfatase is a homo-dimer containing two zinc atoms per monomer. The monomer chains are interlinked in a dimerization domain. To study the unfolding pathways denaturation experiments were performed. In the native protein three forms coexist in chemical equilibrium, each with a different number of zinc atoms. In the chemical unfolding of the holo-dimers the entanglement of the chains is preserved and acts as a 'folding seed', allowing the unfolding process to be reversible. Thermal irreversible unfolding of the holo-dimers favours dissociation, producing monomers that are SDS-stabilized. The thermal unfolding of these monomers is reversible. However, it is not possible to form dimers from unfolded monomers.

Dissecting ITC data of metal ions binding to ligands and proteins

BACKGROUND

ITC is a powerful technique that can reliably assess the thermodynamic underpinnings of a wide range of binding events. When metal ions are involved, complications arise in evaluating the data due to unavoidable solution chemistry that includes metal speciation and a variety of linked equilibria.

SCOPE OF REVIEW

This paper identifies these concerns, provides recommendations to avoid common mistakes, and guides the reader through the mathematical treatment of ITC data to arrive at a set of thermodynamic state functions that describe identical chemical events and, ideally, are independent of solution conditions. Further, common metal chromophores used in biological metal sensing studies are proposed as a robust system to determine unknown solution competition.

MAJOR CONCLUSIONS

Metal ions present several complications in ITC experiments. This review presents strategies to avoid these pitfalls and proposes and experimentally validates mathematical approaches to deconvolute complex equilibria that exist in these systems.

GENERAL SIGNIFICANCE

This review discusses the wide range of complications that exists in metal-based ITC experiments. It provides a starting point for scientists new to this field and articulates concerns that will help experienced researchers troubleshoot experiments.

Cadmium-containing carbonic anhydrase CDCA1 in marine diatom Thalassiosira weissflogii

The Carbon Concentration Mechanism (CCM) allows phytoplakton species to accumulate the dissolved inorganic carbon (DIC) necessary for an efficient photosynthesis even under carbon dioxide limitation. In this mechanism of primary importance for diatoms, a key role is played by carbonic anhydrase (CA) enzymes which catalyze the reversible hydration of CO2, thus taking part in the acquisition of inorganic carbon for photosynthesis. A novel CA, named CDCA1, has been recently discovered in the marine diatom Thalassiosira weissflogii. CDCA1 is a cambialistic enzyme since it naturally uses Cd2+ as catalytic metal ion, but if necessary can spontaneously exchange Cd2+ to Zn2+. Here, the biochemical and structural features of CDCA1 enzyme will be presented together with its putative biotechnological applications for the detection of metal ions in seawaters.

Structure-based studies on the metal binding of two-metal-dependent sugar isomerases

Two-metal-dependent sugar isomerases are important in the synthesis of rare sugars. Many of their properties, specifically their metal dependency, have not been sufficiently explored. Here we used X-ray crystallography, site-directed mutagenesis, isothermal titration calorimetry and electron paramagnetic resonance spectroscopy to investigate the molecular determinants of the metal-binding affinity of l-rhamnose isomerase, a two-Mn(2+) -dependent isomerase from Bacillus halodurans (BHRI). The crystal structure of BHRI confirmed the presence of two metal ion-binding sites: a structural metal ion-binding site for substrate binding, and a catalytic metal ion-binding site that catalyzes a hydride shift. One conserved amino acid, W38, in wild-type BHRI was identified as a critical residue for structural Mn(2+) binding and thus the catalytic efficiency of BHRI. This function of W38 was explored by replacing it with other amino acids. Substitution by Phe, His, Lys, Ile or Ala caused complete loss of catalytic activity. The role of W38 was further examined by analyzing the crystal structure of wild-type BHRI and two inactive mutants of BHRI (W38F and W38A) in complex with Mn(2+) . A structural comparison of the mutants and the wild-type revealed differences in their coordination of Mn(2+) , including changes in metal-ligand bond length and affinity for Mn(2+) . The role of W38 was further confirmed in another two-metal-dependent enzyme: xylose isomerase from Bacillus licheniformis. These data suggest that W38 stabilizes protein-metal complexes and in turn assists ligand binding during catalysis in two-metal-dependent isomerases.

STRUCTURED DIGITAL ABSTRACT

BHRI and BHRI bind by x-ray crystallography (View interaction).

The mononuclear metal center of type-I dihydroorotase from Aquifex aeolicus

BACKGROUND

Dihydroorotase (DHO) is a zinc metalloenzyme, although the number of active site zinc ions has been controversial. E. coli DHO was initially thought to have a mononuclear metal center, but the subsequent X-ray structure clearly showed two zinc ions, α and β, at the catalytic site. Aquifex aeolicus DHO, is a dodecamer comprised of six DHO and six aspartate transcarbamoylase (ATC) subunits. The isolated DHO monomer, which lacks catalytic activity, has an intact α-site and conserved β-site ligands, but the geometry of the second metal binding site is completely disrupted. However, the putative β-site is restored when the complex with ATC is formed and DHO activity is regained. Nevertheless, the X-ray structure of the complex revealed a single zinc ion at the active site. The structure of DHO from the pathogenic organism, S. aureus showed that it also has a single active site metal ion.

RESULTS

Zinc analysis showed that the enzyme has one zinc/DHO subunit and the addition of excess metal ion did not stimulate catalytic activity, nor alter the kinetic parameters. The metal free apoenzyme was inactive, but the full activity was restored upon the addition of one equivalent of Zn2+ or Co2+. Moreover, deletion of the β-site by replacing the His180 and His232 with alanine had no effect on catalysis in the presence or absence of excess zinc. The 2.2 Å structure of the double mutant confirmed that the β-site was eliminated but that the active site remained otherwise intact.

CONCLUSIONS

Thus, kinetically competent A. aeolicus DHO has a mononuclear metal center. In contrast, elimination of the putative second metal binding site in amidohydrolyases with a binuclear metal center, resulted in the abolition of catalytic activity. The number of active site metal ions may be a consideration in the design of inhibitors that selectively target either the mononuclear or binuclear enzymes.

Mutation of outer-shell residues modulates metal ion co-ordination strength in a metalloenzyme

The metal ion co-ordination sites of many metalloproteins have been characterized by a variety of spectroscopic techniques and small-molecule model systems, revealing many important insights into the structural determinants of metal ion co-ordination. However, our understanding of this fundamentally and practically important phenomenon remains frustratingly simplistic; in many proteins it is essentially impossible to predict metal ion specificity and the effects of remote ‘outer-shell’ residues on metal ion co-ordination strength are also poorly defined. This is exemplified by our inability to explain why metalloenzymes with identical metal ion co-ordination spheres, such as the closely related orthologues of bacterial PTE (phosphotriesterase) from Agrobacterium radiobacter and Pseudomonas diminuta, display different metal ion specificity and co-ordination strength. In the present study, we present a series of PTE variants that all possess identical metal ion co-ordination spheres, yet display large differences in their metal ion co-ordination strength. Using measurement of the rates of metal ion dissociation from the active site alongside analysis of structural data obtained through X-ray crystallography, we show that ‘outer-shell’ residues provide essential support for the metal ion ligands, in effect buttressing them in their optimal orientation. Remote mutations appear to modulate metal ion interactions by increasing or decreasing the stabilizing effects of these networks. The present study therefore provides a description of how the greater protein fold can be modified to ‘tune’ the strength of metal ion co-ordination and metal ion specificity, as well as reinforcing the concept of proteins as ensembles of conformational states with unique structures and biochemical properties.

Kinetic and CD/MCD spectroscopic studies of the atypical, three-His-ligated, non-heme Fe2+ center in diketone dioxygenase: the role of hydrophilic outer shell residues in catalysis

Diketone cleaving enzyme (Dke1) is a dioxygenase with an atypical, three-histidine-ligated, mononuclear non-heme Fe(2+) center. To assess the role in enzyme catalysis of the hydrophilic residues in the active site pocket, residues Glu98, Arg80, Tyr70, and Thr107 were subjected to mutational analysis. Steady state and pre-steady state kinetics indicated a role for Glu98 in promoting both substrate binding and O(2) reduction. Additionally, the Glu98 substitution eliminated the pH dependence of substrate binding (k(cat)(app)/K(M)(app)-pH profile) present in wild-type Dke1 (pK(a) = 6.3 +/- 0.4 and 8.4 +/- 0.4). MCD spectroscopy revealed that the Glu98 --> Gln mutation leads to the conversion of the six-coordinate (6C) resting Fe(2+) center present in the wild-type enzyme at pH 7.0 to a mixture of five-coordinate (5C) and 6C sites. The 6C geometry was restored with a pH shift to 9.5 which also resulted in ligand field (LF) energy splittings identical to that found for wild-type (WT) Dke1 at pH 9.5. In WT Dke1, these LF transitions are shifted up in energy by approximately 300 cm(-1) at pH 9.5 relative to pH 7.0. These data, combined with CD pH titrations which reveal a pK(a) of approximately 8.2 for resting WT Dke1 and the Glu98 --> Gln variant, indicate the deprotonation of a metal-ligated water. Together, the kinetic and spectroscopic data reveal a stabilizing effect of Glu98 on the 6C geometry of the metal center, priming it for substrate ligation. Arg80 and Tyr70 are shown to promote O(2) reduction, while Thr107 stabilizes the Fe(II) cofactor.

Carbonic anhydrase II-based metal ion sensing: Advances and new perspectives

Carbonic anhydrases are archetypical zinc metalloenzymes and as such, they have been developed as the recognition element of a family of fluorescent indicators (sensors) to detect metal ions, particularly Zn(2+) and Cu(2+). Subtle modification of the structure of human carbonic anhydrase II isozyme (CAII) alters the selectivity, sensitivity, and response time for these sensors. Sensors using CAII variants coupled with zinc-dependent fluorescent ligands demonstrate picomolar sensitivity, unmatched selectivity, ratiometric fluorescence signal, and near diffusion-controlled response times. Recently, these sensors have been applied to measuring the readily exchangeable concentrations of zinc in the cytosol and nucleus of mammalian tissue culture cells and concentrations of free Cu(2+) in seawater.

Zinc coordination geometry and ligand binding affinity: the structural and kinetic analysis of the second-shell serine 228 residue and the methionine 180 residue of the aminopeptidase from Vibrio proteolyticus

The chemical properties of zinc make it an ideal metal to study the role of coordination strain in enzymatic rate enhancement. The zinc ion and the protein residues that are bound directly to the zinc ion represent a functional charge/dipole complex, and polarization of this complex, which translates to coordination distortion, may tune electrophilicity, and hence, reactivity. Conserved protein residues outside of the charge/dipole complex, such as second-shell residues, may play a role in supporting the electronic strain produced as a consequence of functional polarization. To test the correlation between charge/dipole polarity and ligand binding affinity, structure-function studies were carried out on the dizinc aminopeptidase from Vibrio proteolyticus. Alanine substitutions of S228 and M180 resulted in catalytically diminished enzymes whose crystal structures show very little change in the positions of the metal ions and the protein residues. However, more detailed inspections of the crystal structures show small positional changes that account for differences in the zinc ion coordination geometry. Measurements of the binding affinity of leucine phosphonic acid, a transition state analogue, and leucine, a product, show a correlation between coordination geometry and ligand binding affinity. These results suggest that the coordination number and polarity may tune the electrophilicity of zinc. This may have provided the evolving enzyme with the ability to discriminate between reaction coordinate species.

Determination of zinc using carbonic anhydrase-based fluorescence biosensors

This chapter summarizes the use of carbonic anhydrase (CA)-based fluorescent indicators to determine free zinc in solution, in cells, and in subcellular organelles. Expression (both in situ and in vitro) and preparation of CA-based indicators are described, together with techniques of their use, and procedures to minimize contamination. Recipes for zinc buffers are supplied.

Dissecting the General Physicochemical Properties of Noncovalent Interactions Involving Tyrosine Side Chain as a Second-Shell Ligand in Biomolecular Metal-Binding Site Mimetics: An Experimental Study Combining Fluorescence,13C NMR Spectroscopy and ESI Mass Spectrometry

Detailed physicochemical features inherent in the dynamic cation-pi interactions of aromatic amino acid side chains in the secondary coordination spheres around metal ions were extracted and mapped by intrinsic tyrosine fluorescence titration experiments with two homologous, artificially engineered metal-binding scaffolds which mimic metal-binding sites in metalloproteins. A newly formulated method for the treatment of fluorescence titration data allows straightforward assessment of both the magnitudes and properties of metal-chelation-assisted cation-aromatic interactions (K2) underlying a proposed two-step metallosupramolecular association process. The unprecedented linear platform-motif correlations between the two contrasting scaffolds in their changes in tyrosine fluorescence on binding of 3d metal cations help to elucidate the properties of general cation-arene recognition corresponding to the metal-responsive characteristics of the second-shell Tyr residue surrounding the metal-binding sites in the supramolecular context, and thereby define a new noncovalent design principle for metal-ion recognition in aqueous solution. As supported by NMR spectroscopic and ESI-MS analyses and molecular mechanics force field calculations, the systematic study exemplifies the concept of using steady-state tyrosine fluorescence as a powerful tool for comprehensive descriptions of cation-pi interactions in the extended environment of a metal-binding site. We established that the physicochemical properties pertaining to indirect metal-arene interactions are highly dependent on the electronic properties of the metal ions. This work suggests that second-shell cation-pi interactions may play more diverse roles, including modulation of structure, reactivity, and function of metal-binding sites, than the previously well-established direct cation-pi interactions involving hard cations (e.g., alkali metal ions). Moreover, such a study will continue to complement theoretical predications and/or the early experimental investigations in organic solvents.

Of Folding and Function: Understanding Active-Site Context through Metalloenzyme Design

In the emerging field of biomolecular design, the introduction of metal-binding sites into loop or turn regions of known protein scaffolds has been utilized to create unique metalloprotein and metallopeptide systems for study. This Forum Article highlights examples of the modular-turn-substitution approach to design and the range of structural and mechanistic questions to which this tool can be applied. Examples from the authors' laboratory are given to show that lanthanide-binding metallopeptides, and now a full metallohomeodomain, can be generated by modular substitution of a Ca-binding EF-hand loop into the unrelated scaffold, the engrailed helix-turn-helix motif. We have previously shown that these peptides bind trivalent Ln(III) ions and promote DNA and phosphate hydrolysis, the targeted function. Here, a series of chimeric peptides are presented that differ only in the ninth loop position [given in parentheses; Peptides P3N (Asn), P3E (Glu), P3A (Ala), and P3W(D) (Asp]. This residue, a putative second-shell ligand stabilizing a coordinated water, was found to influence not only metal affinity but also peptide folding. The affinity for Tb(III) was determined by Trp-Tb fluorescence resonance energy transfer and followed the order Ka = P3W(D) > P3A approximately P3E > P3N. However, circular dichroism (CD) titrations with EuCl3 showed that only P3W(D) and P3N folded to any extent upon metal binding, indicating that the Asp/Asn side chains stabilize the central loop structure and thus propagate folding of the peripheral helices, whereas neither Ala nor Glu appears to be interacting with the metal to organize the loop. Finally, we investigated the longer range context of a given loop substitution by cloning and expressing a lanthanide-binding homeodomain (C2), whose loop insertion sequence is analogous to that of peptide P3W(D). We find by CD that apo-C2 has a significant helical structure (approximately 25% alphahelicity), which increases further upon the addition of Tb(III) (approximately 32% alpha helicity). The protein's Tb(III) affinity is similar to that of the chimeric peptides. However, unlike previously reported metallopeptides, we find that EuC2 does not appreciably promote phosphate or DNA cleavage, which suggests a difference in metal accessibility in the context of the full domain. We have demonstrated that substituting beta turns with metal-binding turns does not necessarily require homologous parental scaffolds or small flexible peptides but rather relies on the structural similarity of the motifs flanking the turn.

Identification of Mn2+-binding aspartates from alpha, beta, and gamma subunits of human NAD-dependent isocitrate dehydrogenase

The human NAD-dependent isocitrate dehydrogenase (IDH), with three types of subunits present in the ratio of 2alpha:1beta:1gamma, requires a divalent metal ion to catalyze the oxidative decarboxylation of isocitrate. With the aim of identifying ligands of the enzyme-bound Mn(2+), we mutated aspartates on the alpha, beta, or gamma subunits. Mutagenesis target sites were based on crystal structures of metal-isocitrate complexes of Escherichia coli and pig mitochondrial NADP-IDH and sequence alignments. Aspartates replaced by asparagine or cysteine were 206, 230, and 234 of the alpha subunit and those corresponding to alpha-Asp-206: 217 of the beta subunit and 215 of the gamma subunit. Each expressed, purified mutant enzyme has two wild-type subunits and one subunit with a single mutation. Specific activities of WT, alpha-D206N, alpha-D230C, alpha-D234C, beta-D217N, and gamma-D215N enzymes are 22, 29, 1.4, 0.2, 7.3 and 3.7 micromol of NADH/min/mg, respectively, whereas alpha-D230N and alpha-D234N enzymes showed no activity. The K(m,Mn(2+)) for alpha-D230C and gamma-D215N are increased 32- and 100-fold, respectively, along with elevations in K(m,isocitrate). The K(m,NAD) of alpha-D230C is increased 16-fold, whereas that of beta-D217N is elevated 10-fold. For all the mutants K(m,isocitrate) is decreased by ADP, indicating that these aspartates are not needed for normal ADP activation. This study demonstrates that alpha-Asp-230 and alpha-Asp-234 are critical for catalytic activity, but alpha-Asp-206 is not needed; alpha-Asp-230 and gamma-Asp-215 may interact directly with the Mn(2+); and alpha-Asp-230 and beta-Asp-217 contribute to the affinity of the enzyme for NAD. These results suggest that the active sites of the human NAD-IDH are shared between alpha and gamma subunits and between alpha and beta subunits.

Solvation and the hidden thermodynamics of a zinc finger probed by nonstandard repair of a protein crevice

The classical Zn finger contains a phenylalanine at the crux of its three architectural elements: a beta-hairpin, an alpha-helix, and a Zn(2+)-binding site. Surprisingly, phenylalanine is not required for high-affinity Zn2+ binding, but instead contributes to the specification of a precise DNA-binding surface. Substitution of phenylalanine by leucine leads to a floppy but native-like structure whose Zn affinity is maintained by marked entropy-enthalpy compensation (DeltaDeltaH -8.3 kcal/mol and -TDeltaDeltaS 7.7 kcal/mol). Phenylalanine and leucine differ in shape, size, and aromaticity. To distinguish which features correlate with dynamic stability, we have investigated a nonstandard finger containing cyclohexanylalanine at this site. The structure of the nonstandard finger is similar to that of the native domain. The cyclohexanyl ring assumes a chair conformation, and conformational fluctuations characteristic of the leucine variant are damped. Although the nonstandard finger exhibits a lower affinity for Zn2+ than does the native domain (DeltaDeltaG -1.2 kcal/mol), leucine-associated perturbations in enthalpy and entropy are almost completely attenuated (DeltaDeltaH -0.7 kcal/mol and -TDeltaDeltaS -0.5 kcal/mol). Strikingly, global changes in entropy (as inferred from calorimetry) are in each case opposite in sign from changes in configurational entropy (as inferred from NMR). This seeming paradox suggests that enthalpy-entropy compensation is dominated by solvent reorganization rather than nominal molecular properties. Together, these results demonstrate that dynamic and thermodynamic perturbations correlate with formation or repair of a solvated packing defect rather than type of physical interaction (aromatic or aliphatic) within the core.

Investigation of metal binding and activation of Escherichia coli glyoxalase I: kinetic, thermodynamic and mutagenesis studies

GlxI (glyoxalase I) isomerizes the hemithioacetal formed between glutathione and methylglyoxal. Unlike other GlxI enzymes, Escherichia coli GlxI exhibits no activity with Zn(2+) but maximal activation with Ni(2+). To elucidate further the metal site in E. coli GlxI, several approaches were undertaken. Kinetic studies indicate that the catalytic metal ion affects the k (cat) without significantly affecting the K (m) for the substrate. Inductively coupled plasma analysis and isothermal titration calorimetry confirmed one metal ion bound to the enzyme, including Zn(2+), which produces an inactive enzyme. Isothermal titration calorimetry was utilized to determine the relative binding affinity of GlxI for various bivalent metals. Each metal ion examined bound very tightly to GlxI with an association constant ( K (a))>10(7) M(-1), with the exception of Mn(2+) ( K (a) of the order of 10(6) M(-1)). One of the ligands to the catalytic metal, His(5), was altered to glutamine, a side chain found in the Zn(2+)-active Homo sapiens GlxI. The affinity of the mutant protein for all bivalent metals was drastically decreased. However, low levels of activity were now observed for Zn(2+)-bound GlxI. Although this residue has a marked effect on metal binding and activation, it is not the sole factor determining the differential metal activation between the human and E. coli GlxI enzymes.

First-second shell interactions in metal binding sites in proteins: a PDB survey and DFT/CDM calculations

The role of the second shell in the process of metal binding and selectivity in metalloproteins has been elucidated by combining Protein Data Bank (PDB) surveys of Mg, Mn, Ca, and Zn binding sites with density functional theory/continuum dielectric methods (DFT/CDM). Peptide backbone groups were found to be the most common second-shell ligand in Mg, Mn, Ca, and Zn binding sites, followed (in decreasing order) by Asp/Glu, Lys/Arg, Asn/Gln, and Ser/Thr side chains. Aromatic oxygen- or nitrogen-containing side chains (Tyr, His, and Trp) and sulfur-containing side chains (Cys and Met) are seldom found in the second coordination layer. The backbone and Asn/Gln side chain are ubiquitous in the metal second coordination layer as their carbonyl oxygen and amide hydrogen can act as a hydrogen-bond acceptor and donor, respectively, and can therefore partner practically every first-shell ligand. The second most common outer-shell ligand, Asp/Glu, predominantly hydrogen bonds to a metal-bound water or Zn-bound histidine and polarizes the H-O or H-N bond. In certain cases, a second-shell Asp/Glu could affect the protonation state of the metal ligand. It could also energetically stabilize a positively charged metal complex more than a neutral ligand such as the backbone and Asn/Gln side chain. As for the first shell, the second shell is predicted to contribute to the metal selectivity of the binding site by discriminating between metal cations of different ionic radii and coordination geometries. The first-shell-second-shell interaction energies decay rapidly with increasing solvent exposure of the metal binding site. They are less favorable but are of the same order of magnitude as compared to the respective metal-first-shell interaction energies. Altogether, the results indicate that the structure and properties of the second shell are dictated by those of the first layer. The outer shell is apparently designed to stabilize/protect the inner-shell and complement/enhance its properties.

The hidden thermodynamics of a zinc finger

The Zn finger provides a model for studies of protein structure and stability. Its core contains a conserved phenylalanine residue adjoining three architectural elements: a beta-hairpin, an alpha-helix and a tetrahedral Zn(2+)-binding site. Here, we demonstrate that the consensus Phe is not required for high-affinity Zn(2+) binding but contributes to the specification of a precise DNA-binding surface. Substitution of Phe by leucine in a ZFY peptide permits Zn(2+)-dependent folding. Although a native-like structure is retained, structural fluctuations lead to attenuation of selected nuclear Overhauser enhancements and accelerated amide proton exchange. Surprisingly, wild-type Zn affinity is maintained by entropy-enthalpy compensation (EEC): a hidden entropy penalty (TDeltaDeltaS 7kcal/mol) is balanced by enhanced enthalpy of association (DeltaDeltaH -7kcal/mol) at 25 degrees C. Because the variant is less well ordered than the Phe-anchored domain, the net change in entropy is opposite to the apparent change in configurational entropy. By analogy to the thermodynamics of organometallic complexation, we propose that EEC arises from differences in solvent reorganization. Exclusion of Leu among biological sequences suggests an evolutionary constraint on the dynamics of a Zn finger.