On the role of Glu-68 in alcohol dehydrogenase.

Theoretical computations (molecular dynamics and combined quantum chemical and molecular mechanical geometry optimizations) have been performed on horse liver alcohol dehydrogenase. The results provide evidence that Glu-68, a highly conserved residue located 0.47 nm from the catalytic zinc ion, may intermittently coordinate to the zinc ion. Structures with Glu-68 coordinated to the zinc ion are almost as stable as structures with Glu-68 at the crystal position and the barrier between the two configurations of Glu-68 is so low that it can readily be bypassed at room temperature. There is a cavity behind the zinc ion that seems to be tailored to allow such coordination of Glu-68 to the zinc ion. It is suggested that Glu-68 may facilitate the exchange of ligands in the substrate site by coordinating to the zinc ion when the old ligand dissociates.     

The active-site metal coordination geometry of cadmium-substituted alcohol dehydrogenase

 The structure of eleven complexes of cadmium-substituted alcohol dehydrogenase with or without coenzyme and with different non-protein cadmium ligands has been estimated by combined quantum chemical and molecular mechanical geometry optimisations. The geometry of the optimised complexes is similar to the crystal structure of cadmium-substituted alcohol dehydrogenase, indicating that the method behaves well. The optimised structures do not differ significantly (except for the metal bond lengths) from those of the corresponding zinc complexes, which shows that cadmium is a good probe of zinc coordination geometries. The electric field gradients at the cadmium nucleus have been calculated quantum chemically at the MP2 level with a large cadmium basis set, and they have been used to interpret experimental data obtained by perturbed angular correlation of γ-rays. The experimental and calculated field gradients (all three eigenvalues) differ by less than 0.35 a.u. (3.4·10²¹ Vm^–2), the average error is 0.11 a.u., and the average relative error in the two largest eigenvalues of the field gradients is 9%. Calculated field gradients of four-coordinate structures agree better with the experimental results than do those of any five-coordinate model. Thus, the results indicate that the catalytic metal ion remains four-coordinate in all examined complexes. Two measurements are best explained by a four-coordinate cadmium ion with Glu-68 as the fourth ligand, indicating that Glu-68 probably coordinates intermittently to the catalytic metal ion in horse liver alcohol dehydrogenase under physiological conditions. Received: 10 January 1997 / Accepted: 24 May 1997     

Protein strain in blue copper proteins studied by free energy perturbations.

Free energy perturbations have been performed on two blue copper proteins, plastocyanin and nitrite reductase. By changing the copper coordination geometry, force constants, and charges, we have estimated the maximum energy with which the proteins may distort the copper coordination sphere. By comparing this energy with the quantum chemical energy cost for the same perturbation on the isolated copper complex, various hypotheses about protein strain have been tested. The calculations show that the protein can only modify the copper-methionine bond length by a modest amount of energy-<5 kJ/mol-and they lend no support to the suggestion that the quite appreciable difference in the copper coordination geometry encountered in the two proteins is a result of the proteins enforcing different Cu-methionine bond lengths. On the contrary, this bond is very flexible, and neither the geometry nor the electronic structure change appreciably when the bond length is changed. Moreover, the proteins are rather indifferent to the length of this bond. Instead, the Cu(II) coordination geometries in the two proteins represent two distinct minima on the potential surface of the copper ligand sphere, characterized by different electronic structures, a tetragonal, mainly sigma-bonded, structure in nitrite reductase and a trigonal, pi-bonded, structure in plastocyanin. In vacuum, the structures have almost the same energy, and they are stabilized in the proteins by a combination of geometric and electrostatic interactions. Plastocyanin favors the bond lengths and electrostatics of the trigonal structure, whereas in nitrite reductase, the angles are the main discriminating factor. Proteins 1999;36:157-174.     

Mechanistic inferences from the binding of ligands to LpxC, a metal-dependent deacetylase

The metal-dependent deacetylase LpxC catalyzes the first committed step of lipid A biosynthesis in Gram-negative bacteria. Accordingly, LpxC is an attractive target for the development of inhibitors that may serve as potential new antibiotics for the treatment of Gram-negative bacterial infections. Here, we report the 2.7 A resolution X-ray crystal structure of LpxC complexed with the substrate analogue inhibitor TU-514 and the 2.0 A resolution structure of LpxC complexed with imidazole. The X-ray crystal structure of LpxC complexed with TU-514 allows for a detailed examination of the coordination geometry of the catalytic zinc ion and other enzyme-inhibitor interactions in the active site. The hydroxamate group of TU-514 forms a bidentate chelate complex with the zinc ion and makes hydrogen bond interactions with conserved active site residues E78, H265, and T191. The inhibitor C-4 hydroxyl group makes direct hydrogen bond interactions with E197 and H58. Finally, the C-3 myristate moiety of the inhibitor binds in the hydrophobic tunnel of the active site. These intermolecular interactions provide a foundation for understanding structural aspects of enzyme-substrate and enzyme-inhibitor affinity. Comparison of the TU-514 complex with cacodylate and imidazole complexes suggests a possible substrate diphosphate binding site and highlights residues that may stabilize the tetrahedral intermediate and its flanking transition states in catalysis. Evidence of a catalytic zinc ion in the native zinc enzyme coordinated by H79, H238, D242, and two water molecules with square pyramidal geometry is also presented. These results suggest that the native state of this metallohydrolase may contain a pentacoordinate zinc ion, which contrasts with the native states of archetypical zinc hydrolases such as thermolysin and carboxypeptidase A.     

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.     

First principle calculations of <Superscript>113</Superscript>Cd chemical shifts for proteins and model systems

¹¹³Cd isotropic NMR shieldings are calculated for a number of metal ion binding sites in proteins, using the GIAO-B3LYP and GIAO-HF methods with the uncontracted (19s15p9d4f) polarized basis set of Kellö and Sadlej on cadmium and 6-31G(d) on the ligands. The results compare favorably with experimental data, indicating that first principle calculations are a useful tool for structural interpretation of ¹¹³Cd chemical shift data from metal ion containing proteins. The effect of different ligand types (thiolate, imidazole, water, and monodentate carboxylate), coordination number, and deviations of the coordination geometry from ideal structures is evaluated. In particular, the ligand type and coordination number are important factors, but also changes in cadmium–ligand bond lengths may cause significant changes of the chemical shift.     

Zinc-mediated amino acid discrimination in cysteinyl-tRNA synthetase

Escherichia coli cysteinyl-tRNA synthetase (CysRS) achieves a high level of amino acid specificity without an editing reaction. The crystal structure of CysRS bound to substrate cysteine suggested that direct thiol coordination to a tightly bound zinc ion at the base of the active site is the primary determinant of selectivity against non-cognate amino acids. This hypothesis has now been supported by spectroscopic studies of cobalt-substituted CysRS. Binding of cysteine, but not non-cognate amino acids, induces high absorption in the ligand-to-metal charge transfer region, providing evidence for formation of a metal-thiolate bond. In addition, mutations in the zinc ligands alter the absorption spectrum without reducing the discrimination against non-cognate amino acids. These results argue strongly for a major role for the zinc ion in amino acid discrimination by CysRS, where the tight zinc-thiolate interaction and the strict structural geometry of the metal ion are sufficient to reject serine by more than 20,000-fold at the binding step.     

Structural basis of DNA flexibility

It has been recently shown by us, on the basis of crystal structure database that the flexibility of B-DNA double helices depends significantly on their base sequence. Our model building studies further indicated that the existence of bifurcated cross-strand hydrogen bonds between successive base pairs is possibly the main factor behind the sequence directed DNA flexibility. These cross-strand hydrogen bonds are, of course, weaker than the usual Watson-Crick hydrogen bonds and their bond geometry is characterized by relatively larger bond lengths and smaller bond angles. We have tried to improve our model structures by incorporating non-planarity of the amino groups in DNA bases due to the presence of lone pair electrons at the nitrogen atoms. Energy minimization studies have been carried out by using different quantum chemical methods, whereby it is found that in all cases of N-H....O type cross-strand hydrogen bonds, the bond geometry improves significantly. In the cases of N-H....N type hydrogen bonds, however, no such consistent improvements can be noticed. Perhaps the true picture would emerge only if all the other interactions present in the DNA macromolecule could be appropriately taken into account.     

Analysis of zinc binding sites in protein crystal structures

The geometrical properties of zinc binding sites in a dataset of high quality protein crystal structures deposited in the Protein Data Bank have been examined to identify important differences between zinc sites that are directly involved in catalysis and those that play a structural role. Coordination angles in the zinc primary coordination sphere are compared with ideal values for each coordination geometry, and zinc coordination distances are compared with those in small zinc complexes from the Cambridge Structural Database as a guide of expected trends. We find that distances and angles in the primary coordination sphere are in general close to the expected (or ideal) values. Deviations occur primarily for oxygen coordinating atoms and are found to be mainly due to H-bonding of the oxygen coordinating ligand to protein residues, bidentate binding arrangements, and multi-zinc sites. We find that H-bonding of oxygen containing residues (or water) to zinc bound histidines is almost universal in our dataset and defines the elec-His-Zn motif. Analysis of the stereochemistry shows that carboxyl elec-His-Zn motifs are geometrically rigid, while water elec-His-Zn motifs show the most geometrical variation. As catalytic motifs have a higher proportion of carboxyl elec atoms than structural motifs, they provide a more rigid framework for zinc binding. This is understood biologically, as a small distortion in the zinc position in an enzyme can have serious consequences on the enzymatic reaction. We also analyze the sequence pattern of the zinc ligands and residues that provide elecs, and identify conserved hydrophobic residues in the endopeptidases that also appear to contribute to stabilizing the catalytic zinc site. A zinc binding template in protein crystal structures is derived from these observations.     

Active site electronic structure and dynamics during metalloenzyme catalysis

Zinc-dependent enzymes play important roles in many cellular processes. Assignment of their reaction mechanisms is often a subject of debate because the zinc ion is silent in several spectroscopic techniques. We have combined time-resolved X-ray absorption spectroscopy, pre-steady state kinetics and computational quantum chemistry to study the active site zinc ion of bacterial alcohol dehydrogenase during single substrate turnover. We detect a series of alternations in the coordination number and structure of the catalytic zinc ion with concomitant changes in metal-ligand bond distances. These structural changes are reflected in the effective charge of the metal ion. The present work emphasizes the flexibility of catalytic zinc sites during catalysis and provides novel mechanistic insights into alcohol dehydrogenase catalysis.     

Structure of acidic phospholipase A2 for the venom of Agkistrodon halys blomhoffii at 2.8 A resolution.

The crystal structure of acidic phospholipase A2 from the venom of Agkistrodon halys blomhoffii has been determined by molecular replacement methods based on the known structure of Crotalus atrox PLA2, a same group II enzyme. The overall structures, except the calcium-binding regions, are very similar to each other. A calcium ion is pentagonally ligated to two carboxylate oxygen atoms of Asp-49 and each carbonyl oxygen atoms of Tyr-28, Gly-30 and Ala-31. A reason why the former enzyme functions as monomeric form, while the latter one does as dimer, could be presumed by the structural comparison of these calcium-binding regions. Although Gly-32 is usually participated as a ligand in the coordination with calcium ion in group I PLA2, it is characteristically replaced to Ala-31 in the present structure, and thus the coordination geometry of calcium ion is rather different from the usually observed one.     

Higher metal-ligand coordination in the catalytic site of cobalt-substituted Thermoanaerobacter brockii alcohol dehydrogenase lowers the barrier for enzyme catalysis

Thermoanaerobacter brockii alcohol dehydrogenase (TbADH) is a zinc-dependent NADP(+)/H-linked class enzyme that reversibly catalyzes the oxidation of secondary alcohols to their corresponding ketones. Cobalt substitution studies of other members of the alcohol dehydrogenase (ADH) family showed that the cobalt-containing ADHs have a similar active site structure but slightly decreased activity compared to wild-type zinc ADHs. In contrast, the cobalt-substituted TbADH (Co-TbADH) exhibits an increase in specific activity compared to the native enzyme [Bogin, O., Peretz, M., and Burstein, Y. (1997) Protein Sci. 6, 450-458]. However, the structural basis underlying this behavior is not yet clear. To shed more light on this issue, we studied the local structure and electronics at the catalytic metal site in Co-TbADH by combining X-ray absorption (XAS) and quantum chemical calculations. Importantly, we show that the first metal-ligand coordination shell of Co-TbADH is distorted compared to its native tetrahedral coordination shell and forms an octahedral structure. This is mediated presumably by the addition of two water molecules and results in more positively charged catalytic metal ions. Recently, we have shown that the metal-ligand coordination number of the zinc ion in TbADH changes dynamically during substrate turnover. These structural changes are associated with a higher coordination number of the native catalytic zinc ion and the consequent buildup of a positive charge. Here we propose that the accumulation of a higher coordination number and positive charge at the catalytic metal ion in TbADH stabilizes the structure of the catalytic transition state and hence lowers the barrier for enzyme catalysis.     

Syntheses, characterization, density functional theory calculations, and activity of tridentate SNS zinc pincer complexes

A series of tridentate SNS ligand precursors were metallated with ZnCl₂ to give new tridentate SNS pincer zinc complexes. The zinc complexes serve as models for the zinc active site in liver alcohol dehydrogenase (LADH) and were characterized with single crystal X-ray diffraction, ¹H, ¹³C, and HSQC NMR spectroscopies and electrospray mass spectrometry. The bond lengths and bond angles of the zinc complexes correlate well to those in horse LADH. The zinc complexes feature SNS donor atoms and pseudotetrahedral geometry about the zinc center, as is seen for liver alcohol dehydrogenase. The SNS ligand precursors were characterized with ¹H, ¹³C, and HSQC NMR spectroscopies and cyclic voltammetry, and were found to be redox active. Gaussian calculations were performed and agree quite well with the experimentally observed oxidation potential for the pincer ligand. The zinc complexes were screened for the reduction of electron poor aldehydes in the presence of a hydrogen donor, 1-benzyl-1,4-dihydronicotinamide (BNAH). The zinc complexes enhance the reduction of electron poor aldehydes. Density functional theory calculations were performed to better understand why the geometry about the zinc center is pseudo-tetrahedral rather than pseudo-square planar, which is seen for most pincer complexes. For the SNS tridentate pincer complexes, the data indicate that the pseudo-tetrahedral geometry was 43.8 kcal/mol more stable than the pseudo-square planar geometry. Density functional theory calculations were also performed on zinc complexes with monodentate ligands and the data indicate that the pseudo-tetrahedral geometry was 30.6 kcal/mol more stable than pseudo-square planar geometry. Overall, the relative stabilities of the pseudo-tetrahedral and pseudo-square planar systems are the same for this coordination environment whether the ligand set is a single tridentate SNS system or is broken into three separate units. The preference of a d¹⁰ Zn center to attain a tetrahedral local environment trumps any stabilization gained by removal of constraints within the ligand set.     

A strong 13C chemical shift signature provides the coordination mode of histidines in zinc-binding proteins

Zinc is the second most abundant metal ion incorporated in proteins, and is in many cases a crucial component of protein three-dimensional structures. Zinc ions are frequently coordinated by cysteine and histidine residues. Whereas cysteines bind to zinc via their unique S(γ) atom, histidines can coordinate zinc with two different coordination modes, either N(δ1) or N(ε2) is coordinating the zinc ion. The determination of this coordination mode is crucial for the accurate structure determination of a histidine-containing zinc-binding site by NMR. NMR chemical shifts contain a vast amount of information on local electronic and structural environments and surprisingly their utilization for the determination of the coordination mode of zinc-ligated histidines has been limited so far to (15)N nuclei. In the present report, we observed that the (13)C chemical shifts of aromatic carbons in zinc-ligated histidines represent a reliable signature of their coordination mode. Using a statistical analysis of (13)C chemical shifts, we show that (13)C(δ2) chemical shift is sensitive to the histidine coordination mode and that the chemical shift difference δ{(13)C(ε1)} - δ{(13)C(δ2)} provides a reference-independent marker of this coordination mode. The present approach allows the direct determination of the coordination mode of zinc-ligated histidines even with non-isotopically enriched protein samples and without any prior structural information.     

Refined structure of porcine cytosolic adenylate kinase at 2.1 A resolution

The crystal structure of porcine cytosolic adenylate kinase has been established at 2.1 A resolution using a restrained least-squares refinement method. Based on 11,251 independent reflections of better than 10 A resolution, a final R-factor of 19.3% was obtained with a model obeying standard geometry within 0.026 A in bond lengths and 3.3 degrees in bond angles. In comparison with the previous structure at 3 A resolution, there is a significant improvement. The high resolution structure has been used to rationalize the strictly conserved residues in the adenylate kinase family. Among these is the glycine-rich loop, which forms a giant anion hole accommodating a sulfate ion which mimics a phosphoryl group of a substrate. Such a structure seems to occur in a large group of mononucleotide binding proteins. Moreover, a conserved cis-proline has been detected in the active center. A structural comparison with the complex between adenylate kinase from yeast and a substrate-analog at medium resolution indicates that this kinase performs appreciable mechanical movements during a catalytic cycle. The reported structure presumably represents an open form of the enzyme, similar to that in solution in the absence of substrates. However, since there are large intermolecular contacts in the crystal, some deviation from the solution structure has to be expected.     

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.     

Mechanism of cyanamide hydration catalyzed by carbonic anhydrase II suggested by cryogenic X-ray diffraction

The three-dimensional structure of a possible intermediate in the hydration reaction of cyanamide to urea catalyzed by human carbonic anhydrase II (hCAII) has been determined by cryocrystallographic techniques. The crystal structure shows that two different adducts are formed under the experimental conditions and that they have different occupancy in the crystal. The high occupancy form consists of a binary hCAII-cyanamide complex where the substrate has replaced the zinc-bound hydroxide anion present in the native enzyme, maintaining the tetrahedral geometry around the metal ion. The second, low-occupancy form consists of a hCAII-cyanamide-water ternary complex where the catalytic zinc ion, still being bound to cyanamide, is approached by a water molecule in a five-coordinate adduct. While the first form can be considered a nonproductive complex, the second form may represent an intermediate state of the catalyzed reaction where the water molecule is about to perform a nucleophilic attack on the zinc-activated cyanamide substrate. The structural evidence is consistent with the kinetic data previously reported about this recently described hydrolytic reaction catalyzed by hCAII, and indicates that a different mechanism with respect to that generally accepted for the physiologic carbon dioxide hydration reaction may be adopted by the enzyme, depending on the substrate chemical properties.     

EXAFS investigation of the active site of iron superoxide dismutase of Escherichia coli and Propionibacterium shermanii

The local structure of the iron site in ferric superoxide dismutase from P. shermanii was analyzed by X-ray absorption spectroscopy. The metal-ligand cluster of the enzyme is found to be similar to the crystallographically investigated ferric superoxide dismutase from E. coli. At pH 6.4 the enzyme is five-fold coordinated with three histidines, an aspartate and a water molecule. The average bond lengths between the metal and the histidines are about 2.10 Å, between metal and aspartate they are about 1.86 Å and between metal and water 1.96 Å. With an increase in pH a change in the coordination number from five to six is observed both in pre-edge peak and EXAFS spectra analysis. However, the bond lengths of the ligands do not change dramatically, they are conserved for the aspartate and increase slightly to 2.13 Å for the average metal - histidine distance at pH 9.3. The observation of the increase in coordination number is correlated with a decrease in enzymatic activity which occurs in the high pH range. The zinc EXAFS spectra of P. shermanii superoxide dismutase have shown that zinc can be incorporated in the active center instead of the iron.     

The refined structure of the complex between adenylate kinase from beef heart mitochondrial matrix and its substrate AMP at 1.85 A resolution

The crystal structure of the complex between adenylate kinase from bovine mitochondrial matrix and its substrate AMP has been refined at 1.85 A resolution (1 A = 0.1 nm). Based on 42,519 independent reflections of better than 10 A resolution, a final R-factor of 18.9% was obtained with a model obeying standard geometry within 0.016 A in bond lengths and 3.2 degrees in bond angles. There are two enzyme: substrate complexes in the asymmetric unit, each consisting of 226 amino acid residues, one AMP and one sulfate ion. A superposition of the two full-length polypeptides revealed deviations that can be described as small relative movements of three domains. Best superpositions of individual domains yielded a residual overall root-mean-square deviation of 0.3 A for the backbone atoms and 0.5 A for the sidechains. The final model contains 381 solvent molecules in the asymmetric unit, 2 x 72 = 144 of which occupy corresponding positions in both complexes.