Electrostatic properties of protein-protein complexes

Statistical electrostatic analysis of 37 protein-protein complexes extracted from the previously developed database of protein complexes (ProtCom, http://www.ces.clemson.edu/compbio/protcom) is presented. It is shown that small interfaces have a higher content of charged and polar groups compared to large interfaces. In a vast majority of the cases the average pKa shifts for acidic residues induced by the complex formation are negative, indicating that complex formation stabilizes their ionizable states, whereas the histidines are predicted to destabilize the complex. The individual pKa shifts show the same tendency since 80% of the interfacial acidic groups were found to lower their pKas, whereas only 25% of histidines raise their pKa upon the complex formation. The interfacial groups have been divided into three sets according to the mechanism of their pKa shift, and statistical analysis of each set was performed. It was shown that the optimum pH values (pH of maximal stability) of the complex tend to be the same as the optimum pH values of the complex components. This finding can be used in the homology-based prediction of the 3D structures of protein complexes, especially when one needs to evaluate and rank putative models. It is more likely for a model to be correct if both components of the model complex and the entire complex have the same or at least similar values of the optimum pH.     

Histidine residues at the N- and C-termini of alpha-helices: perturbed pKas and protein stability.

A single histidine residue has been placed at either the N-terminus or the C-terminus of each of the two alpha-helices of barnase. The pKa of that histidine residue in each of the four mutants has been determined by 1H NMR. The pKas of the two residues at the C-terminus are, on average, 0.5 unit higher, and those of the residues at the N-terminus are 0.8 unit lower, than the pKa of histidines in unfolded barnase at low ionic strength. The conformational stability of the mutant proteins at different values of pH has been measured by urea denaturation. C-Terminal histidine mutants are approximately 0.6 kcal mol-1 more stable when the introduced histidine is protonated, both at low and high ionic strength. N-Terminal mutants with a protonated histidine residue are approximately 1.1 kcal mol-1 less stable at low ionic strength and 0.5 kcal mol-1 less stable at high ionic strength (1 M NaCl). The low-field 1H NMR spectra of the mutant proteins at low pH suggest that the C-terminal histidines form hydrogen bonds with the protein while the N-terminal histidines do not form the same. The perturbations of pKa and stability result from a combination of different electrostatic environments and hydrogen-bonding patterns at either ends of helices. The value of 0.6 kcal mol-1 represents a lower limit to the favorable electrostatic interaction between the alpha-helix dipole and a protonated histidine residue at the C-terminal end of the helix.(ABSTRACT TRUNCATED AT 250 WORDS)     

Very fast empirical prediction and rationalization of protein pKa values

A very fast empirical method is presented for structure-based protein pKa prediction and rationalization. The desolvation effects and intra-protein interactions, which cause variations in pKa values of protein ionizable groups, are empirically related to the positions and chemical nature of the groups proximate to the pKa sites. A computer program is written to automatically predict pKa values based on these empirical relationships within a couple of seconds. Unusual pKa values at buried active sites, which are among the most interesting protein pKa values, are predicted very well with the empirical method. A test on 233 carboxyl, 12 cysteine, 45 histidine, and 24 lysine pKa values in various proteins shows a root-mean-square deviation (RMSD) of 0.89 from experimental values. Removal of the 29 pKa values that are upper or lower limits results in an RMSD = 0.79 for the remaining 285 pKa values.     

Determination of pKa values of the histidine side chains of phosphatidylinositol-specific phospholipase C from Bacillus cereus by NMR spectroscopy and site-directed mutagenesis.

Two active site histidine residues have been implicated in the catalysis of phosphatidylinositol-specific phospholipase C (PI-PLC). In this report, we present the first study of the pKa values of histidines of a PI-PLC. All six histidines of Bacillus cereus PI-PLC were studied by 2D NMR spectroscopy and site-directed mutagenesis. The protein was selectively labeled with 13C epsilon 1-histidine. A series of 1H-13C HSQC NMR spectra were acquired over a pH range of 4.0-9.0. Five of the six histidines have been individually substituted with alanine to aid the resonance assignments in the NMR spectra. Overall, the remaining histidines in the mutants show little chemical shift changes in the 1H-13C HSQC spectra, indicating that the alanine substitution has no effect on the tertiary structure of the protein. H32A and H82A mutants are inactive enzymes, while H92A and H61A are fully active, and H81A retains about 15% of the wild-type activity. The active site histidines, His32 and His82, display pKa values of 7.6 and 6.9, respectively. His92 and His227 exhibit pKa values of 5.4 and 6.9. His61 and His81 do not titrate over the pH range studied. These values are consistent with the crystal structure data, which shows that His92 and His227 are on the surface of the protein, whereas His61 and His81 are buried. The pKa value of 6.9 corroborates the hypothesis of His82 acting as a general acid in the catalysis. His32 is essential to enzyme activity, but its putative role as the general base is in question due to its relatively high pKa.     

Assignment of histidine resonances in the 1H NMR (500 MHz) spectrum of subtilisin BPN' using site-directed mutagenesis

A spin-echo pulse sequence has been used to resolve the six histidine C-2H protons in the 500-MHz NMR spectrum of subtilisin BPN'. Five of these residues have been substituted by site-directed mutagenesis, and this has enabled a complete assignment of these protons to be obtained. Analysis of the pH titration curves of these signals has provided microscopic pKas for the six histidines in this enzyme. The pKas of the histidine residues in subtilisin BPN' have been compared with the values obtained for the histidines in the homologous enzyme from Bacillus licheniformis (subtilisin Carlsberg). Four of the five conserved histidines titrate with essentially identical pKa's in the two enzymes. It therefore appears that the assignments made for these residues in subtilisin BPN' can be transferred to subtilisin Carlsberg. On the basis of these assignments, the one histidine that titrates with a substantially different pKa in the two enzymes can be assigned to histidine-238. This difference in pKa has been attributed to a Trp to Lys substitution at position 241 in subtilisin Carlsberg.     

Relationship between electrostatics and redox function in human thioredoxin: characterization of pH titration shifts using two-dimensional homo- and heteronuclear NMR.

The electrostatic behavior of potentially titrating groups in reduced human thioredoxin was investigated using two-dimensional (2D) 1H and 15N nuclear magnetic resonance (NMR) spectroscopy. A total of 241 chemical shift titration curves were measured over the pH range of 2.1-10.6 from homonuclear 1H-1H Hartmann-Hahn (HOHAHA) and heteronuclear 1H-15N Overbodenhausen correlation spectra. Nonlinear least-squares fits of the data to simple relationships derived from the Henderson-Hasselbalch equation led to the determination of pKas for certain isolated ionizable groups, including the single histidine residue at position 43 (pKa = 5.5 +/- 0.1) and a number of aspartic and glutamic acid carboxylate groups. Many of the titration curves demonstrate complex behavior due to the effects of interacting titrating groups, the long range of electrostatic interactions through the protein interior, and, perhaps, pH-induced conformational changes on the chemical shifts. Unambiguous assignment of the pKas for most of the 38 potentially ionizing groups of human thioredoxin could therefore not be made. In addition, there was no clear evidence that Asp-26 titrates in a manner corresponding to that observed in the Escherichia coli protein [Dyson, H. J., Tennant, L. L., & Holmgren, A. (1991) Biochemistry 30, 4262-4268]. The pKas of the active site cysteines were measured, however, with Cys-32 having an anomalously low value of 6.3 +/- 0.1 and that of Cys-35 between 7.5 and 8.6. These pKas are in agreement with proposed mechanisms for redox catalysis of thioredoxin and previously measured pKas within the active site of E. coli thioredoxin [Kallis, G. B., & Holmgren, A. (1980) J. Biol. Chem. 255, 10261-10265]. The stabilization of a thiolate anion at physiological pH can be explained by the interaction of the S gamma of Cys-32 with the amide of Cys-35 observed in the previously determined high-resolution solution structure of reduced human thioredoxin [Forman-Kay, J. D., Clore, G. M., Wingfield, P. T., & Gronenborn, A. M. (1991) Biochemistry 30, 2685-2698].     

Electrostatic contribution of surface charge residues to the stability of a thermophilic protein: benchmarking experimental and predicted pKa values

Optimization of the surface charges is a promising strategy for increasing thermostability of proteins. Electrostatic contribution of ionizable groups to the protein stability can be estimated from the differences between the pKa values in the folded and unfolded states of a protein. Using this pKa-shift approach, we experimentally measured the electrostatic contribution of all aspartate and glutamate residues to the stability of a thermophilic ribosomal protein L30e from Thermococcus celer. The pKa values in the unfolded state were found to be similar to model compound pKas. The pKa values in both the folded and unfolded states obtained at 298 and 333 K were similar, suggesting that electrostatic contribution of ionizable groups to the protein stability were insensitive to temperature changes. The experimental pKa values for the L30e protein in the folded state were used as a benchmark to test the robustness of pKa prediction by various computational methods such as H++, MCCE, MEAD, pKD, PropKa, and UHBD. Although the predicted pKa values were affected by crystal contacts that may alter the side-chain conformation of surface charged residues, most computational methods performed well, with correlation coefficients between experimental and calculated pKa values ranging from 0.49 to 0.91 (p<0.01). The changes in protein stability derived from the experimental pKa-shift approach correlate well (r = 0.81) with those obtained from stability measurements of charge-to-alanine substituted variants of the L30e protein. Our results demonstrate that the knowledge of the pKa values in the folded state provides sufficient rationale for the redesign of protein surface charges leading to improved protein stability.     

Electrostatic coupling between retinal isomerization and the ionization state of Glu-204: a general mechanism for proton release in bacteriorhodopsin.

The pKa values of ionizable groups that lie between the active site region of bacteriorhodopsin (bR) and the extracellular surface of the protein are reported. Glu-204 is found to have an elevated pKa in the resting state of bR, suggesting that it corresponds to the proton-releasing group in bR. Its elevated pKa is predicted to be due in part to strong repulsive interactions with Glu-9. Following trans-cis isomerization of the retinal chromophore and the transfer of a proton to Asp-85, polar groups on the protein are able to interact more strongly with the ionized state of Glu-204, leading to a substantial reduction of its pKa. This suggests a general mechanism for proton release in which isomerization and subsequent charge separation initially produce a new electrostatic balance in the active site of bR. Here it is proposed that those events in turn drives a conformational change in the protein in which the ionized state of Glu-204 can be stabilized through interactions with groups that were previously inaccessible. Whether these groups should be identified with polar moieties in the protein, bound waters, or Arg-82 is an important mechanistic question whose elucidation will require further study.     

Identification by proton nuclear magnetic resonance of the histidines in cytochrome b5 modified by diethyl pyrocarbonate

Diethyl pyrocarbonate (DEP) is an electrophilic reagent that is used to modify reversibly the histidine residues of proteins. Unfortunately, the lability of the acylated histidine adduct usually does not permit the isolation and identification of the modified histidine. By use of 500-MHz proton NMR spectroscopy, it has been possible to identify the C-H resonances of the nonaxial histidines of trypsin-solubilized bovine, rabbit, and porcine cytochrome b5 and therefore observe the interaction of DEP with specific histidine residues of cytochrome b5. In addition, the pKa of the peripheral histidines of bovine and rabbit cytochrome b5 have been measured in D2O. In the bovine protein it was found that the histidines are modified sequentially with increasing DEP concentration in the order His-26 greater than His-15 greater than His-80. This order is maintained in the rabbit protein with the following additions: His-26 approximately His-27 greater than His-15 greater than or equal to His-17 greater than His-80. The relative reactivity of the peripheral histidines with DEP was rationalized by considering three of their characteristics: (1) the pKa of the histidine, (2) the fraction of the side chain exposed to the solvent, and (3) the hydrogen-bond interactions of the imidazole ring.     

Why are polar residues within the membrane core evolutionary conserved?

Here, we present a study of polar residues within the membrane core of alpha-helical membrane proteins. As expected, polar residues are less frequent in the membrane than expected. Further, most of these residues are buried within the interior of the protein and are only rarely exposed to lipids. However, the polar groups often border internal water filled cavities, even if the rest of the sidechain is buried. A survey of their functional roles in known structures showed that the polar residues are often directly involved in binding of small compounds, especially in channels and transporters, but other functions including proton transfer, catalysis, and selectivity have also been attributed to these proteins. Among the polar residues histidines often interact with prosthetic groups in photosynthetic- and oxidoreductase-related proteins, whereas prolines often are required for conformational changes of the proteins. Indeed, the polar residues in the membrane core are more conserved than other residues in the core, as well as more conserved than polar residues outside the membrane. The reason is twofold; they are often (i) buried in the interior of the protein and (ii) directly involved in the function of the proteins. Finally, a method to identify which polar residues are present within the membrane core directly from protein sequences was developed. Applying the method to the set of all human membrane proteins the prediction indicates that polar residues were most frequent among active transporter proteins and GPCRs, whereas infrequent in families with few transmembrane regions, such as non-GPCR receptors.     

Model for the differential stabilities of rhinovirus and poliovirus to mild acidic pH, based on electrostatics calculations.

Previous calculations of electrostatic interactions in the rhinovirus capsid have identified a subset of histidine residues, paired with lysine or arginine, that may be involved in pH-induced conformational changes related to viral uncoating. Further calculations with the finite difference method, accounting for the dielectric environment of the ionizable groups, suggest that charge burial in the crystal conformation will prevent protonation of these histidine residues in the pentamer-pentamer interface. Calculations with a modelled pentamer-pentamer interface in which three beta-strands are removed recover mildly acidic pKa values for the histidines. These results are discussed in the context of the structural interactions of these three beta-strands, which form a beta-sheet extension from the rest of the capsid, and with regard to the conformation of the homologous beta-sheet extension in poliovirus, which also possesses homologous histidine-lysine/arginine pairs. A model is developed in which the structural stability of the beta-sheet extension is related to the difference in acid stability of rhinovirus and poliovirus. It is suggested that, for poliovirus prior to cell receptor binding, the beta-sheet extension is stable at pH 3, the pentamer-pentamer interface histidines remain buried, and the virus is acid-stable. Cell receptor binding of poliovirus destabilizes the beta-sheet extension and the acid lability that is proposed to result could be involved in viral uncoating. For rhinovirus it is suggested that the observed conformational change in the absence of cell receptor binding involves a further acidic pH-activated process or conformational fluctuations that rearrange the beta-sheet extension and expose the pentamer-pentamer interface histidine residues to the acidic medium. Sequence analysis and electrostatics calculations reveal an aspartic acid in the beta-sheet extension that may have different pKa values in rhinovirus and poliovirus.     

A denaturation-induced proton-uptake study of horse ferricytochrome c.

The observation that 6 M-urea denatures horse ferricytochrome c in the pH range 4-6, but not horse ferrocytochrome c, has been exploited to determine the denaturation-induced proton uptake of ferricytochrome c. This is related to the pKa values of ionizable groups buried within the native protein. The data indicate that one of the haem propionic acid substituents of ferricytochrome c has a pKa of less than 4.5, whereas the other has a pKa of greater than 9.     

Structure-based pKa prediction provides a thermodynamic basis for the role of histidines in pH-induced conformational transitions in dengue virus

pH-induced conformational changes in dengue virus (DENV) are critical to its ability to infect host cells. The envelope protein heterodimers that make up the viral envelope shift from a dimer to a trimer conformation at low-pH during membrane fusion. Previous studies have suggested that the ionization of histidine residues at low-pH is central to this pH-induced conformational change. We sought out to use molecular modeling with structure-based pKa prediction to provide a quantitative basis for the role of histidines in pH-induced conformational changes and identify which histidine residues were primarily responsible for this transition. We combined existing crystallographic and cryo-electron microscopy data to construct templates of the dimer and trimer conformations for the mature and immature virus. We then generated homology models for the four DENV serotypes and carried out structure-based pKa prediction using Rosetta. Our results showed that the pKa values of a subset of conserved histidines in DENV successfully capture the thermodynamics necessary to drive pH-induced conformational changes during fusion. Here, we identified the structural determinants underlying these pKa values and compare our findings with previous experimental results.     

Electrostatic contributions to residue-specific protonation equilibria and proton binding capacitance for a small protein

Charge-charge interactions in proteins are important in a host of biological processes. Here we use 13C NMR chemical shift data for individual aspartate and glutamate side chain carboxylate groups to accurately detect site-specific protonation equilibria in a variant of the B1 domain of protein G (PGB1-QDD). Carbon chemical shifts are dominated by changes in the electron distribution within the side chain and therefore excellent reporters of the charge state of individual groups, and the data are of high precision. We demonstrate that it is possible to detect local charge interactions within this small protein domain that stretch and skew the chemical shift titration curves away from "ideal" behavior and introduce a framework for the analysis of such convoluted data to study local charge-charge interactions and electrostatic coupling. It is found that, due to changes in electrostatic potential, the proton binding affinity, Ka, of each carboxyl group changes throughout the titration process and results in a linearly pH dependent pKa value. This result could be readily explained by calculations of direct charge-charge interactions based on Coulomb's law. In addition, the slope of pKa versus pH was dependent on screening by salt, and this dependence allowed the selective study of charge-charge interactions. For PGB1-QDD, it was established that mainly differences in self-energy, and not direct charge-charge interactions, are responsible for shifted pKa values within the protein environment.     

Characterization of a buried neutral histidine residue in Bacillus circulans xylanase: NMR assignments, pH titration, and hydrogen exchange.

Bacillus circulans xylanase contains two histidines, one of which (His 156) is solvent exposed, whereas the other (His 149) is buried within its hydrophobic core. His 149 is involved in a network of hydrogen bonds with an internal water and Ser 130, as well as a potential weak aromatic-aromatic interaction with Tyr 105. These three residues, and their network of interactions with the bound water, are conserved in four homologous xylanases. To probe the structural role played by His 149, NMR spectroscopy was used to characterize the histidines in BCX. Complete assignments of the 1H, 13C, and 15N resonances and tautomeric forms of the imidazole rings were obtained from two-dimensional heteronuclear correlation experiments. An unusual spectroscopic feature of BCX is a peak near 12 ppm arising from the nitrogen bonded 1H epsilon 2 of His 149. Due to its solvent inaccessibility and hydrogen bonding to an internal water molecule, the exchange rate of this proton (4.0 x 10(-5) s-1 at pH*7.04 and 30 degrees C) is retarded by > 10(6)-fold relative to an exposed histidine. The pKa of His 156 is unperturbed at approximately 6.5, as measured from the pH dependence of the 15N- and 1H-NMR spectra of BCX. In contrast, His 149 has a pKa < 2.3, existing in the neutral N epsilon 2H tautomeric state under all conditions examined. BCX unfolds at low pH and 30 degrees C, and thus His 149 is never protonated significantly in the context of the native enzyme. The structural importance of this buried histidine is confirmed by the destablizing effect of substituting a phenylalanine or glutamine at position 149 in BCX.     

Simplified methods for pKa and acid pH-dependent stability estimation in proteins: removing dielectric and counterion boundaries

Much computational research aimed at understanding ionizable group interactions in proteins has focused on numerical solutions of the Poisson-Boltzmann (PB) equation, incorporating protein exclusion zones for solvent and counterions in a continuum model. Poor agreement with measured pKas and pH-dependent stabilities for a (protein, solvent) relative dielectric boundary of (4,80) has lead to the adoption of an intermediate (20,80) boundary. It is now shown that a simple Debye-Huckel (DH) calculation, removing both the low dielectric and counterion exclusion regions associated with protein, is equally effective in general pKa calculations. However, a broad-based discrepancy to measured pH-dependent stabilities is maintained in the absence of ionizable group interactions in the unfolded state. A simple model is introduced for these interactions, with a significantly improved match to experiment that suggests a potential utility in predicting and analyzing the acid pH-dependence of protein stability. The methods are applied to the relative pH-dependent stabilities of the pore-forming domains of colicins A and N. The results relate generally to the well-known preponderance of surface ionizable groups with solvent-mediated interactions. Although numerical PB solutions do not currently have a significant advantage for overall pKa estimations, development based on consideration of microscopic solvation energetics in tandem with the continuum model could combine the large deltapKas of a subset of ionizable groups with the overall robustness of the DH model.     

Do salt bridges stabilize proteins? A continuum electrostatic analysis

The electrostatic contribution to the free energy of folding was calculated for 21 salt bridges in 9 protein X-ray crystal structures using a continuum electrostatic approach with the DELPHI computer-program package. The majority (17) were found to be electrostatically destabilizing; the average free energy change, which is analogous to mutation of salt bridging side chains to hydrophobic isosteres, was calculated to be 3.5 kcal/mol. This is fundamentally different from stability measurements using pKa shifts, which effectively measure the strength of a salt bridge relative to 1 or more charged hydrogen bonds. The calculated effect was due to a large, unfavorable desolvation contribution that was not fully compensated by favorable interactions within the salt bridge and between salt-bridge partners and other polar and charged groups in the folded protein. Some of the salt bridges were studied in further detail to determine the effect of the choice of values for atomic radii, internal protein dielectric constant, and ionic strength used in the calculations. Increased ionic strength resulted in little or no change in calculated stability for 3 of 4 salt bridges over a range of 0.1-0.9 M. The results suggest that mutation of salt bridges, particularly those that are buried, to "hydrophobic bridges" (that pack at least as well as wild type) can result in proteins with increased stability. Due to the large penalty for burying uncompensated ionizable groups, salt bridges could help to limit the number of low free energy conformations of a molecule or complex and thus play a role in determining specificity (i.e., the uniqueness of a protein fold or protein-ligand binding geometry).     

Ionization properties of titratable groups in ribonuclease T1. I. pKa values in the native state determined by two-dimensional heteronuclear NMR spectroscopy

pKa values of amino acid side chains of ribonuclease T1 have been determined from the pH dependence of 13C and 15N resonances. It was possible to derive pKa values of single protonation or deprotonation sites of carboxylate and imidazole groups. Deviations from pKa values of free amino acids could be interpreted with electrostatic interactions of corresponding side chains with the protein environment. In particular, the interaction between H27 and E82 led to an increase of the H27 pKa and a decrease of the E82 pKa. The pKa of E28 at the C-terminal end of the alpha-helix was increased because of the dipolar character of the alpha-helix. D76 did not titrate in the investigated pH range of about 2-9. From the chemical shift value this buried side chain seems to be protonated. The pKa values of side chains in the active site deviate from a normal behaviour. The lower pKa value of E58 may be interpreted with the close proximity of this side chain with positively charged H40 and R77. A novel two-dimensional 1H(13Cdelta)13Cgamma correlation experiment was developed to observe the pH dependence of the chemical shifts of the Cgamma resonances of histidine residues. From the inspection of the Cgamma chemical shift-pH profiles it was possible to determine the predominant tautomeric form for the histidine residues at higher pH values.     

Electrostatic Environment of Proteorhodopsin Affects the pKa of Its Buried Primary Proton Acceptor

The protonation state of embedded charged residues in transmembrane proteins (TMPs) can control the onset of protein function. It is understood that interactions between an embedded charged residue and other charged or polar residues in the moiety would influence its pKa, but how the surrounding environment in which the TMP resides affects the pKa of these residues is unclear. Proteorhodopsin (PR), a light-responsive proton pump from marine bacteria, was used as a model to examine externally accessible factors that tune the pKa of its embedded charged residue, specifically its primary proton acceptor D97. The pKa of D97 was compared between PR reconstituted in liposomes with different net headgroup charges and equilibrated in buffer with different ion concentrations. For PR reconstituted in net positively charged compared to net negatively charged liposomes in low-salt buffer solutions, a drop of the apparent pKa from 7.6 to 5.6 was observed, whereas intrinsic pKa modeled with surface pH calculated from Gouy-Chapman predictions found an opposite trend for the pKa change, suggesting that surface pH does not account for the main changes observed in the apparent pKa. This difference in the pKa of D97 observed from PR reconstituted in oppositely charged liposome environments disappeared when the NaCl concentration was increased to 150 mM. We suggest that protein-intrinsic structural properties must play a role in adjusting the local microenvironment around D97 to affect its pKa, as corroborated with observations of changes in protein side-chain and hydration dynamics around the E-F loop of PR. Understanding the effect of externally controllable factors in tuning the pKa of TMP-embedded charged residues is important for bioengineering and biomedical applications relying on TMP systems, in which the onset of functions can be controlled by the protonation state of embedded residues.     

Active-site serine phosphate and histidine residues of phosphoglucomutase: pH titration studies monitored by 1H and 31P NMR spectroscopy

1H and 31P NMR pH titrations were conducted to monitor changes in the environment and protonation state of the histidine residues and phosphoserine group of rabbit muscle phosphoglucomutase on binding of metal ions at the activating site and of substrate (glucose phosphate) at the catalytic site. Imidazole C epsilon-H signals from 8 of the 10 histidines present in the free enzyme were observed in 1H NMR spectra obtained by a spin-echo pulse sequence at 470 MHz; their pH (uncorrected pH meter reading of a 2H2O solution measured with a glass electrode standardized with H2O buffer) titration properties (in 99% 2H2O) were determined. Three of these histidine residues, which have pKa values ranging from 6.5 to 7.9, exhibited an atypical pH-dependent perturbation of their chemical shifts with a pHmid of 5.8 and a Hill coefficient of about 2. Since none of the observed histidines has a pKa near 5.8, it appears that these three histidines interact with a cluster consisting of two or more groups which become protonated cooperatively at this pH. Binding of Cd2+ at the activating site of the enzyme abolishes the pH-dependent transition of these histidines; hence, the putative anion cluster may constitute the metal ion binding site, or part of it. Two separate 31P NMR peaks from phosphoserine-116 of the phosphoenzyme were observed between pH 6 and 9. Apparently, the metal-free enzyme exists as a pH-dependent mixture of conformers that provide two different environments, I and II, for the enzymic phosphate group; the transition of the phosphate group between these two environments is slow on the NMR time scale.(ABSTRACT TRUNCATED AT 250 WORDS)