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.     

Structural influence of hydrophobic core residues on metal binding and specificity in carbonic anhydrase II

Aromatic residues in the hydrophobic core of human carbonic anhydrase II (CAII) influence metal ion binding in the active site. Residues F93, F95, and W97 are contained in a beta-strand that also contains two zinc ligands, H94 and H96. The aromatic amino acids contribute to the high zinc affinity and slow zinc dissociation rate constant of CAII [Hunt, J. A., and Fierke, C. A. (1997) J. Biol. Chem. 272, 20364-20372]. Substitution of these aromatic amino acids with smaller side chains enhances Cu(2+) affinity while decreasing Co(2+) and Zn(2+) affinity [Hunt, J. A., Mahiuddin, A., & Fierke, C. A. (1999) Biochemistry 38, 9054-9062]. Here, X-ray crystal structures of zinc-bound F93I/F95M/W97V and F93S/F95L/W97M CAIIs reveal the introduction of new cavities in the hydrophobic core, compensatory movements of surrounding side chains, and the incorporation of buried water molecules; nevertheless, the enzyme maintains tetrahedral zinc coordination geometry. However, a conformational change of direct metal ligand H94 as well as indirect (i.e., "second-shell") ligand Q92 accompanies metal release in both F93I/F95M/W97V and F93S/F95L/W97M CAIIs, thereby eliminating preorientation of the histidine ligands with tetrahedral geometry in the apoenzyme. Only one cobalt-bound variant, F93I/F95M/W97V CAII, maintains tetrahedral metal coordination geometry; F93S/F95L/W97M CAII binds Co(2+) with trigonal bipyramidal coordination geometry due to the addition of azide anion to the metal coordination polyhedron. The copper-bound variants exhibit either square pyramidal or trigonal bipyramidal metal coordination geometry due to the addition of a second solvent molecule to the metal coordination polyhedron. The key finding of this work is that aromatic core residues serve as anchors that help to preorient direct and second-shell ligands to optimize zinc binding geometry and destabilize alternative geometries. These geometrical constraints are likely a main determinant of the enhanced zinc/copper specificity of CAII as compared to small molecule chelators.     

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.     

Deciphering metal ion preference and primary coordination sphere robustness of a designed zinc finger with high‐resolution mass spectrometry

Small zinc finger (ZnF) motifs are promising molecular scaffolds for protein design owing to their structural robustness and versatility. Moreover, their characterization provides important insights into protein folding in general. ZnF motifs usually possess an exceptional specificity and high affinity towards Zn(II) ion to drive folding. While the Zn(II) ion is canonically coordinated by two cysteine and two histidine residues, many other coordination spheres also exist in small ZnFs, all having four amino acid ligands. Here we used high‐resolution mass spectrometry to study metal ion binding specificity and primary coordination sphere robustness of a designed zinc finger, named MM1. Based on the results, MM1 possesses high specificity for zinc with sub‐micromolar binding affinity. Surprisingly, MM1 retains metal ion binding affinity even in the presence of selective alanine mutations of the primary zinc coordinating amino acid residues.     

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

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

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.     

Functional role of the distal valine (E11) residue of alpha subunits in human haemoglobin

We have expressed human alpha-globin to a high level in Escherichia coli as a fusion protein, purified it and removed the N-terminal leader sequence by site-specific proteolysis with blood coagulation factor Xa. The apo globin has been refolded and reconstituted with haem and native beta-globin to form fully functional haemoglobin (Hb) with properties identical to those of native human Hb. By site-directed mutagenesis we have altered the distal residues of the alpha subunits and compared the functional properties of these mutant proteins. The rates of various ligands binding to these proteins in the R-state have been reported by Mathews et al. Here, we present the oxygen equilibrium curves of three E11 alpha mutants and the crystal structures of two of these mutants in the deoxy form. Replacing the distal valine residue of alpha-globin with alanine, leucine or isoleucine has no effect on the oxygen affinity of the protein in either quaternary state, in contrast to the equivalent mutations of beta subunits. The crystal structure of the valine E11 alpha----isoleucine mutant shows that the larger E11 residue excludes water from the haem pocket, but causes no significant movement of other amino acid residues. We conclude that the distal valine residue of alpha-globin does not control the oxygen affinity of the protein by sterically hindering ligand binding.     

Crystallographic analysis of Thr-200 → His human carbonic anhydrase II and its complex with the substrate,HCO

A complex of carbonic anhydrase (CA) with one of its substrates, bicarbonate, has been studied crystallographically. Human isoenzyme II was mutated at position 200 from threonine to histidine, which results in higher affinity for bicarbonate. The HCO ion binds in the active site to the zinc ion as a pseudo-bidentate ligand which gives the metal a coordination geometry between tetrahedral and trigonal bipyramide. The water/hydroxide normally bound with tetrahedral coordination to the zinc is probably replaced by the OH group of the bicarbonate ion. The importance of residues Thr-199 and Glu-106 in controlling the binding orientation of HCO is discussed as well as the catalytic mechanism. Both the complex as well as the uncomplexed mutant were studied at 1.9 Å resolution. © 1993 Wiley-Liss, Inc.     

Computational protocol for predicting the binding affinities of zinc containing metalloprotein-ligand complexes

Zinc is one of the most important metal ions found in proteins performing specific functions associated with life processes. Coordination geometry of the zinc ion in the active site of the metalloprotein-ligand complexes poses a challenge in determining ligand binding affinities accurately in structure-based drug design. We report here an all atom force field based computational protocol for estimating rapidly the binding affinities of zinc containing metalloprotein-ligand complexes, considering electrostatics, van der Waals, hydrophobicity, and loss in conformational entropy of protein side chains upon ligand binding along with a nonbonded approach to model the interactions of the zinc ion with all the other atoms of the complex. We examined the sensitivity of the binding affinity predictions to the choice of Lennard-Jones parameters, partial atomic charges, and dielectric treatments adopted for system preparation and scoring. The highest correlation obtained was R2 = 0.77 (r = 0.88) for the predicted binding affinity against the experiment on a heterogenous dataset of 90 zinc containing metalloprotein-ligand complexes consisting of five unique protein targets. Model validation and parameter analysis studies underscore the robustness and predictive ability of the scoring function. The high correlation obtained suggests the potential applicability of the methodology in designing novel ligands for zinc-metalloproteins. The scoring function has been web enabled for free access at www.scfbio-iitd.res.in/software/drugdesign/bapplz.jsp as BAPPL-Z server (Binding Affinity Prediction of Protein-Ligand complexes containing Zinc metal ions).     

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.     

Short and long spacer sequences and other structural features of zinc binding sites in zinc enzymes

The crystal structures of eleven zinc enzymes have served to identify common features of their Zn binding sites. Two of them have non-catalytic Zn sites, both of which contain four cysteine ligands closely spaced in the linear sequence of the protein with no bound water. In contrast, all the catalytic Zn sites have three protein ligands and, in addition, one coordinated, 'activated' water. Histidine is the predominant ligand. The spacing between the first two ligands (1-3 amino acids), the short spacer, ensures a nucleus for Zn binding. The third ligand, separated by from approximately 20 to approximately 120 amino acids, the long spacer, not only completes the coordination but also aligns protein residues for interaction with the substrate. The short and long spacing observed for catalytic zinc sites may also pertain to Fe and Cu proteins.     

Binding of ligands to the catalytic zinc ion in horse liver alcohol dehydrogenase

The affinity of nitrogen and sulfur ligands for the catalytic zinc ion in horse liver alcohol dehydrogenase has been investigated by their influence on the affinity labeling reaction with iodoacetate. All the nitrogen compounds including ammonia, a primary and a secondary amine, and heterocycles containing a pyridine-type nitrogen with the exception of 2,2-dipyridyl were found to activate the affinity labeling reaction. Activation results from inner-sphere ligand coordination to the catalytic zinc ion. Closely related pyridine compounds gave a regular increase in affinity for the enzyme with increasing basicity, as expected for coordination to a metal ion. The sulfur compounds penicillamine and mercaptoethanol also activated the affinity labeling reaction, but dimercaptopropanol bound very tightly as a bidentate inhibited the reaction. The anions hydrosulfide, diethyldithiocarbamate, and cyanide coordinated to the catalytic zinc ion, whereas azide, thiocyanate, tetrazole, and iodide complexed the anion-binding site. The anionic metal ligands increased the rate of inactivation of the enzyme with iodoacetamide by binding to the catalytic zinc ion, while the binding of iodoacetate to the anion-binding site was prevented.     

Probing the interaction of l‐captopril with metallo‐β‐lactamase CcrA by fluorescence spectra and molecular dynamic simulation

Both the molecular recognition and interaction of metallo‐β‐lactamase CcrA with l ‐captopril were studied by the combined use of fluorescence spectra and molecular dynamic simulation. The results showed that the binding constant was 8.89 × 10⁴ L mol^?1 at 296 K. Both Zn1 and Zn2 displayed tetrahedral coordination geometries in the CcrA–Lcap complex, the S atom in l ‐captopril displaced the nucleophilic hydroxide in apo CcrA and occupied the fourth coordination site for each ion, resulting in a competitively inhibited CcrA enzyme. Strong electrostatic interaction between the two zinc ions in CcrA and negatively charged l ‐captopril provided the main driving force for the binding affinity. Through a partly structural transformation from β‐sheet to random coil, loop 1 (residues 24–34) completely opened the binding pocket of CcrA to allow an induced fit of the newly introduced ligand. This study may provide some valuable information for designing and developing a more tightly binding inhibitor to resist superbugs.     

Metal binding specificity in carbonic anhydrase is influenced by conserved hydrophobic core residues.

The role of highly conserved aromatic residues surrounding the zinc binding site of human carbonic anhydrase II (CAII) in determining the metal ion binding specificity of this enzyme has been examined by mutagenesis. Residues F93, F95, and W97 are located along a beta-strand containing two residues that coordinate zinc, H94 and H96, and these aromatic amino acids contribute to the high zinc affinity and slow zinc dissociation rate constant of CAII [Hunt, J. A., and Fierke, C. A. (1997) J. Biol. Chem. 272, 20364-20372]. Substitutions of these aromatic amino acids with smaller side chains enhance the copper affinity (up to 100-fold) while decreasing the affinity of both cobalt and zinc, thereby altering the metal binding specificity up to 10(4)-fold. Furthermore, the free energy of the stability of native CAII, determined by solvent-induced denaturation, correlates positively with increased hydrophobicity of the amino acids at positions 93, 95, and 97 as well as with cobalt and zinc affinity. Conversely, increased copper affinity correlates with decreased protein stability. Zinc specificity is therefore enhanced by formation of the native enzyme structure. These data suggest that the hydrophobic cluster in CAII is important for orienting the histidine residues to stabilize metals bound with a distorted tetrahedral geometry and to destabilize the trigonal bipyramidal geometry of bound copper. Knowledge of the structural factors that lead to high metal ion specificity will aid in the design of metal ion biosensors and de novo catalytic sites.     

Human glutathione-dependent formaldehyde dehydrogenase. Structural changes associated with ternary complex formation

Human glutathione-dependent formaldehyde dehydrogenase plays an important role in the metabolism of glutathione adducts such as S-(hydroxymethyl)glutathione and S-nitrosoglutathione. The role of specific active site residues in binding these physiologically important substrates and the structural changes during the catalytic cycle of glutathione-dependent formaldehyde dehydrogenase was examined by determining the crystal structure of a ternary complex with S-(hydroxymethyl)glutathione and the reduced coenzyme to 2.6 A resolution. The formation of the ternary complex caused the movement of the catalytic domain toward the coenzyme-binding domain. This represents the first observation of domain closure in glutathione-dependent formaldehyde dehydrogenase in response to substrate binding. A water molecule adjacent to the 2'-ribose hydroxyl of NADH suggests that the alcohol proton is relayed to solvent directly from the coenzyme, rather than through the action of the terminal histidine residue as observed in the proton relay system for class I alcohol dehydrogenases. S-(Hydroxymethyl)glutathione is directly coordinated to the active site zinc and forms interactions with the highly conserved residues Arg114, Asp55, Glu57, and Thr46. The active site zinc has a tetrahedral coordination environment with Cys44, His66, and Cys173 as the three protein ligands in addition to S-(hydroxymethyl)glutathione. This is in contrast to zinc coordination in the binary coenzyme complex where all of the ligands were contributed by the enzyme and included Glu67 as the fourth protein ligand. This change in zinc coordination is accomplished by an approximately 2.3 A movement of the catalytic zinc.     

Individual metal ligands play distinct functional roles in the zinc sensor Staphylococcus aureus CzrA

Recent studies on metalloregulatory proteins suggest that coordination number/geometry and metal ion availability in a host cytosol are key determinants for biological specificity. Here, we investigate the contribution that individual metal ligands of the alpha5 sensing site of Staphylococcus aureus CzrA (Asp84, His86, His97', and His100') make to in vitro metal ion binding affinity, coordination geometry, and allosteric negative regulation of DNA operator/promoter region binding. All ligand substitution mutants exhibit significantly reduced metal ion binding affinity (K(Me)) by > or =10(3) M(-1). Substitutions of Asp84 and His97 give rise to non-native coordination geometries upon metal binding and are non-functional in allosteric coupling of metal and DNA binding (DeltaG(coupling) approximately 0 kcal mol(-1)). In contrast, His86 and His100 could be readily substituted with potentially liganding (Asp, Glu) and poorly liganding (Asn, Gln) residues with significant native-like tetrahedral metal coordination geometry retained in these mutants, leading to strong functional coupling (DeltaG(coupling) > or = +3.0 kcal mol(-1)). 1H-(15)N heteronuclear single quantum coherence (HSQC) spectra of wild-type and mutant CzrAs reveal that all H86 and H100 substitution mutants undergo 4 degrees structural switching on binding Zn(II), while D84N, H97N and H97D CzrAs do not. Thus, only those variant CzrAs that retain some tetrahedral coordination geometry characteristic of wild-type CzrA upon metal binding are capable of driving 4 degrees structural conformational changes linked to allosteric regulation of DNA binding in vitro, irrespective of the magnitude of K(Me).     

Using protein design to dissect the effect of charged residues on metal binding and protein stability

Ca2+ controls biological processes by interacting with proteins with different affinities, which are largely influenced by the electrostatic interaction from the local negatively charged ligand residues in the coordination sphere. We have developed a general strategy for rationally designing stable Ca2+- and Ln3+-binding proteins that retain the native folding of the host protein. Domain 1 of cluster differentiation 2 (CD2) is the host for the two designed proteins in this study. We investigate the effect of local charge on Ca2+-binding affinity based on the folding properties and metal-binding affinities of the two proteins that have similarly located Ca2+-binding sites with two shared ligand positions. While mutation and Ca2+ binding do not alter the native structure of the protein, Ca2+ binding specifically induced changes around the designed Ca2+-binding site. The designed protein with a -5 charge at the binding sphere displays a 14-, 20-, and 12-fold increase in the binding affinity for Ca2+, Tb3+, and La3+, respectively, compared to the designed protein with a -3 charge, which suggests that higher local charges are preferred for both Ca2+ and Ln3+ binding. The localized charged residues significantly decrease the thermal stability of the designed protein with a -5 charge, which has a T(m) of 41 degrees C. Wild-type CD2 has a T(m) of 61 degrees C, which is similar to the designed protein with a -3 charge. This decrease is partially restored by Ca2+ binding. The effect on the protein stability is modulated by the environment and the secondary structure locations of the charged mutations. Our study demonstrates the capability and power of protein design in unveiling key determinants to Ca2+-binding affinity without the complexities of the global conformational changes, cooperativity, and multibinding process found in most natural Ca2+-binding proteins.     

Generation of novel copper sites by mutation of the axial ligand of amicyanin. Atomic resolution structures and spectroscopic properties

Amicyanin from Paracoccus denitrificans is a type 1 copper protein with three strong equatorial copper ligands provided by nitrogens of His53 and His95 and the sulfur of Cys92, with an additional weak axial ligand provided by the sulfur of Met98. Met98 was replaced with either Gln or Ala. As isolated, the M98A and M98Q mutant proteins contain zinc in the active site. The zinc is then removed and replaced with copper so that the copper-containing proteins may be studied. Each of the mutant amicyanins exhibits a marked decrease in thermal stability relative to that of native amicyanin, consistent with the weaker affinity for copper. Crystal structures were obtained for the oxidized and reduced forms of M98A and M98Q amicyanins at atomic resolution (相似文献   

Mutagenesis and molecular modeling reveal the importance of the 5-HT3 receptor F-loop

The 5-HT(3) receptor is a member of the Cys-loop family of ligand-gated ion channels. The extracellular domains of these proteins contain six amino acid loops (A-F) that converge to form the ligand binding site. In this study we have mutated 21 residues in or close to the 5-HT(3) receptor F-loop (Ile(192) to Gly(212)) to Ala or to a residue with similar chemical properties. Mutant receptors were expressed in HEK293 cells, and binding affinity was measured using [(3)H]granisetron. Two regions displayed decreases in binding affinity when mutated to Ala (Ile(192)-Arg(196) and Asp(204)-Ser(206)), but only one region was sensitive when mutated to chemically similar residues (Ile(192)-Val(201)). Homology modeling using acetylcholine-binding protein crystal structures with a variety of different bound ligands suggests there may be distinct movements of Trp(195) and Asp(204) upon ligand binding, indicating that these residues and their immediate neighbors have the ability to interact differently with different ligands. The models suggest predominantly lateral movement around Asp(204) and rotational movement around Trp(195), indicating the former is in a more flexible region. Overall our results are consistent with a flexible 5-HT(3) receptor F-loop with two regions that have specific but distinct roles in ligand binding.