The pKa Values of Acidic and Basic Residues Buried at the Same Internal Location in a Protein Are Governed by Different Factors

The pK_a values of internal ionizable groups are usually very different from the normal pK_a values of ionizable groups in water. To examine the molecular determinants of pK_a values of internal groups, we compared the properties of Lys, Asp, and Glu at internal position 38 in staphylococcal nuclease. Lys38 titrates with a normal or elevated pK_a, whereas Asp38 and Glu38 titrate with elevated pK_a values of 7.0 and 7.2, respectively. In the structure of the L38K variant, the buried amino group of the Lys38 side chain makes an ion pair with Glu122, whereas in the structure of the L38E variant, the buried carboxyl group of Glu38 interacts with two backbone amides and has several nearby carboxyl oxygen atoms. Previously, we showed that the pK_a of Lys38 is normal owing to structural reorganization and water penetration concomitant with ionization of the Lys side chain. In contrast, the pK_a values of Asp38 and Glu38 are perturbed significantly owing to an imbalance between favorable polar interactions and unfavorable contributions from dehydration and from Coulomb interactions with surface carboxylic groups. Their ionization is also coupled to subtle structural reorganization. These results illustrate the complex interplay between local polarity, Coulomb interactions, and structural reorganization as determinants of pK_a values of internal groups in proteins. This study suggests that improvements to computational methods for pK_a calculations will require explicit treatment of the conformational reorganization that can occur when internal groups ionize.     

Electrostatic effects in a network of polar and ionizable groups in staphylococcal nuclease

His121 and His124 are embedded in a network of polar and ionizable groups on the surface of staphylococcal nuclease. To examine how membership in a network affects the electrostatic properties of ionizable groups, the tautomeric state and the pK_a values of these histidines were measured with NMR spectroscopy in the wild-type nuclease and in 13 variants designed to disrupt the network. In the background protein, His121 and His124 titrate with pK_a values of 5.2 and 5.6, respectively. In the variants, where the network was disrupted, the pK_a values range from 4.03 to 6.46 for His121, and 5.04 to 5.99 for His124. The largest decrease in a pK_a was observed when the favorable Coulomb interaction between His121 and Glu75 was eliminated; the largest increase was observed when Tyr91 or Tyr93 was substituted with Ala or Phe. In all variants, the dominant tautomeric state at neutral pH was the N^ε2 state. At one level the network behaves as a rigid unit that does not readily reorganize when disrupted: crystal structures of the E75A or E75Q variants show that even when the pivotal Glu75 is removed, the overall configuration of the network was unaffected. On the other hand, a few key hydrogen bonds appear to govern the conformation of the network, and when these bonds are disrupted the network reorganizes. Coulomb interactions within the network report an effective dielectric constant of 20, whereas a dielectric constant of 80 is more consistent with the magnitude of medium to long-range Coulomb interactions in this protein. The data demonstrate that when structures are treated as static, rigid bodies, structure-based pK_a calculations with continuum electrostatics method are not useful to treat ionizable groups in cases where pK_a values are governed by short-range polar and Coulomb interactions.     

Carboxyl pKa Values and Acid Denaturation of BBL

The protein BBL undergoes structural transitions and acid denaturation between pH 1.2 and 8.0. Using NMR spectroscopy, we measured the pK_a values of all the carboxylic residues in this pH range. We employed ¹³C direct-detection two-dimensional IPAP (in-phase antiphase) CACO NMR spectroscopy to monitor the ionization state of different carboxylic groups and demonstrated its advantages over other NMR techniques in measuring pK_a values of carboxylic residues. The two residues Glu161 and Asp162 had significantly lowered pK_a values, showing that these residues are involved in a network of stabilizing electrostatic interactions, as is His166. The other carboxylates had unperturbed values. The pH dependence of the free energy of denaturation was described quantitatively by the ionizations of those three residues of perturbed pK_a, and, using thermodynamic cycles, we could calculate their pK_as in the native and denatured states as well as the equilibrium constants for denaturation of the different protonation states. We also measured ¹³C^α chemical shifts of individual residues as a function of pH. These shifts sense structural transitions rather than ionizations, and they titrated with pH consistent with the change in equilibrium constant for denaturation. Kinetic measurements of the folding of BBL E161Q indicated that, at pH 7, the stabilizing interactions with Glu161 are formed mainly in the transition state. We also found that local interactions still exist in the acid-denatured state of BBL, which attenuate somewhat the flexibility of the acid-denatured state.     

Structural origins of high apparent dielectric constants experienced by ionizable groups in the hydrophobic core of a protein

The side chains of Lys66, Asp66, and Glu66 in staphylococcal nuclease are fully buried and surrounded mainly by hydrophobic matter, except for internal water molecules associated with carboxylic oxygen atoms. These ionizable side chains titrate with pK_a values of 5.7, 8.8, and 8.9, respectively. To reproduce these pK_a values with continuum electrostatics calculations, we treated the protein with high dielectric constants. We have examined the structural origins of these high apparent dielectric constants by using NMR spectroscopy to characterize the structural response to the ionization of these internal side chains. Substitution of Val66 with Lys66 and Asp66 led to increased conformational fluctuations of the microenvironments surrounding these groups, even under pH conditions where Lys66 and Asp66 are neutral. When Lys66, Asp66, and Glu66 are charged, the proteins remain almost fully folded, but resonances for a few backbone amides adjacent to the internal ionizable residues are broadened. This suggests that the ionization of the internal groups promotes a local increase in dynamics on the intermediate timescale, consistent with either partial unfolding or increased backbone fluctuations of helix 1 near residue 66, or, less likely, with increased fluctuations of the charged side chains at position 66. These experiments confirm that the high apparent dielectric constants reported by internal Lys66, Asp66, and Glu66 reflect localized changes in conformational fluctuations without incurring detectable global structural reorganization. To improve structure-based pK_a calculations in proteins, we will need to learn how to treat this coupling between ionization of internal groups and local changes in conformational fluctuations explicitly.     

Catalytic mechanism of inulinase from Arthrobacter sp. S37

Detailed catalytic roles of the conserved Glu323, Asp460, and Glu519 of Arthrobacter sp. S37 inulinase (EnIA), a member of the glycoside hydrolase family 32, were investigated by site-directed mutagenesis and pH-dependence studies of the enzyme efficiency and homology modeling were carried out for EnIA and for D460E mutant. The enzyme efficiency (k_cat/K_m) of the E323A and E519A mutants was significantly lower than that of the wild-type due to a substantial decrease in k_cat, but not due to variations in K_m, consistent with their putative roles as nucleophile and acid/base catalyst, respectively. The D460A mutant was totally inactive, whereas the D460E and D460N mutants were active to some extent, revealing Asp460 as a catalytic residue and demonstrating that the presence of a carboxylate group in this position is a prerequisite for catalysis. The pH-dependence studies indicated that the pK_a of the acid/base catalyst decreased from 9.2 for the wild-type enzyme to 7.0 for the D460E mutant, implicating Asp460 as the residue that interacts with the acid/base catalyst Glu519 and elevates its pK_a. Homology modeling and molecular dynamics simulation of the wild-type enzyme and the D460E mutant shed light on the structural roles of Glu323, Asp460, and Glu519 in the catalytic activity of the enzyme.     

Electrostatic analysis of the hepatitis C virus NS3 helicase reveals both active and allosteric site locations

Multi-conformation continuum electrostatics (MCCE) was used to analyze various structures of the NS3 RNA helicase from the hepatitis C virus in order to determine the ionization state of amino acid side chains and their pK_as. In MCCE analyses of HCV helicase structures that lacked ligands, several active site residues were identified to have perturbed pK_as in both the nucleic acid binding site and in the distant ATP-binding site, which regulates helicase movement. In all HCV helicase structures, Glu493 was unusually basic and His369 was abnormally acidic. Both these residues are part of the HCV helicase nucleic acid binding site, and their roles were analyzed by examining the pH profiles of site-directed mutants. Data support the accuracy of MCCE predicted pK_a values, and reveal that Glu493 is critical for low pH enzyme activation. Several key residues, which were previously shown to be involved in helicase-catalyzed ATP hydrolysis, were also identified to have perturbed pK_as including Lys210 in the Walker-A motif and the DExD/H-box motif residues Asp290 and His293. When DNA was present in the structure, the calculated pK_as shifted for both Lys210 and Asp290, demonstrating how DNA binding might lead to electrostatic changes that stimulate ATP hydrolysis.     

pH dependence of conformational fluctuations of the protein backbone

Proton binding equilibria (pK_a values) of ionizable groups in proteins are exquisitely sensitive to their microenvironments. Apparent pK_a values measured for individual ionizable residues with NMR spectroscopy are actually population‐weighted averages of the pK_a in different conformational microstates. NMR spectroscopy experiments with staphylococcal nuclease were used to test the hypothesis that pK_a values of surface Glu and Asp residues are affected by pH‐sensitive fluctuations of the backbone between folded and locally unfolded conformations. ¹⁵N spin relaxation studies showed that as the pH decreases from the neutral into the acidic range the amplitudes of backbone fluctuations in the ps‐ns timescale increase near carboxylic residues. Hydrogen exchange experiments suggested that backbone conformational fluctuations promoted by decreasing pH also reflect slower local or sub‐global unfolding near carboxylic groups. This study has implications for structure‐based pK_a calculations: (1) The timescale of the backbone's response to ionization events in proteins can range from ps to ms, and even longer; (2) pH‐sensitive fluctuations of the backbone can be localized to both the segment the ionizable residue is attached to or the one that occludes the ionizable group; (3) Structural perturbations are not necessarily propagated through Coulomb interactions; instead, local fluctuations appear to be coupled through the co‐operativity inherent to elements of secondary structure and to networks of hydrogen bonds. These results are consistent with the idea that local conformational fluctuations and stabilities are important determinants of apparent pK_a values of ionizable residues in proteins. Proteins 2014; 82:3132–3143. © 2014 Wiley Periodicals, Inc.     

A Comprehensive Alanine-Scanning Mutagenesis Study Reveals Roles for Salt Bridges in the Structure and Activity of Pseudomonas aeruginosa Elastase

The relationship between salt bridges and stability/enzymatic activity is unclear. We studied this relationship by systematic alanine-scanning mutation analysis using the typical M4 family metalloprotease Pseudomonas aeruginosa elastase (PAE, also known as pseudolysin) as a model. Structural analysis revealed seven salt bridges in the PAE structure. We constructed ten mutants for six salt bridges. Among these mutants, six (Asp189Ala, Arg179Ala, Asp201Ala, Arg205Ala, Arg245Ala and Glu249Ala) were active and four (Asp168Ala, Arg198Ala, Arg253Ala, and Arg279Ala) were inactive. Five mutants were purified, and their catalytic efficiencies (k _cat/K _m), half-lives (t _1/2) and thermal unfolding curves were compared with those of PAE. Mutants Asp189Ala and Arg179Ala both showed decreased thermal stabilities and increased activities, suggesting that the salt bridge Asp189-Arg179 stabilizes the protein at the expense of catalytic efficiency. In contrast, mutants Asp201Ala and Arg205Ala both showed slightly increased thermal stability and slightly decreased activity, suggesting that the salt bridge Asp201-Arg205 destabilizes the protein. Mutant Glu249Ala is related to a C-terminal salt bridge network and showed both decreased thermal stability and decreased activity. Furthermore, Glu249Ala showed a thermal unfolding curve with three discernable states [the native state (N), the partially unfolded state (I) and the unfolded state (U)]. In comparison, there were only two discernable states (N and U) in the thermal unfolding curve of PAE. These results suggest that Glu249 is important for catalytic efficiency, stability and unfolding cooperativity. This study represents a systematic mutational analyses of salt bridges in the model metalloprotease PAE and provides important insights into the structure-function relationship of enzymes.     

Intrinsically disordered inhibitor of glutamine synthetase is a functional protein with random‐coil‐like pKa values

The sequential action of glutamine synthetase (GS) and glutamate synthase (GOGAT) in cyanobacteria allows the incorporation of ammonium into carbon skeletons. In the cyanobacterium Synechocystis sp. PCC 6803, the activity of GS is modulated by the interaction with proteins, which include a 65‐residue‐long intrinsically disordered protein (IDP), the inactivating factor IF7. This interaction is regulated by the presence of charged residues in both IF7 and GS. To understand how charged amino acids can affect the binding of an IDP with its target and to provide clues on electrostatic interactions in disordered states of proteins, we measured the pK_a values of all IF7 acidic groups (Glu32, Glu36, Glu38, Asp40, Asp58, and Ser65, the backbone C‐terminus) at 100 mM NaCl concentration, by using NMR spectroscopy. We also obtained solution structures of IF7 through molecular dynamics simulation, validated them on the basis of previous experiments, and used them to obtain theoretical estimates of the pK_a values. Titration values for the two Asp and three Glu residues of IF7 were similar to those reported for random‐coil models, suggesting the lack of electrostatic interactions around these residues. Furthermore, our results suggest the presence of helical structure at the N‐terminus of the protein and of conformational changes at acidic pH values. The overall experimental and in silico findings suggest that local interactions and conformational equilibria do not play a role in determining the electrostatic features of the acidic residues of IF7.     

Energy cost for the proper ionization of active site residues in 6-phosphogluconate dehydrogenase from T. brucei

The catalytic mechanism of 6-phosphogluconate dehydrogenase requires the inversion of a Lys/Glu couple from its natural ionization state. The pK_a of these residues in free and substrate bound enzymes has been determined measuring by ITC the proton release/uptake induced by substrate binding at different pH values. Wt 6-phosphogluconate dehydrogenase from Trypanosoma brucei and two active site enzyme mutants, K185H and E192Q were investigated. Substrate binding was accompanied by proton release and was dependent on the ionization of a group with pK_a 7.07 which was absent in the E192Q mutant. Kinetic data highlighted two pK_a, 7.17 and 9.64, in the enzyme–substrate complex, the latter being absent in the E192Q mutant, suggesting that the substrate binding shifts Glu192 pK_a from 7.07 to 9.64. A comparison of wt and E192Q mutant appears to show that the substrate binding shifts Lys185 pK_a from 9.9 to 7.17. By comparing differences in proton release and the binding enthalpy of wt and mutant enzymes, the enthalpic cost of the change in the protonation state of Lys185 and Glu192 was estimated at ≈ 6.1 kcal/mol. The change in protonation state of Lys185 and Glu192 has little effect on Gibbs free energy, 240–325 cal/mol. However proton balance evidences the dissociation of other group(s) that can be collectively described by a single pK_a shift from 9.1 to 7.54. This further change in ionization state of the enzyme causes an increase of free energy with a total cost of 1.2–2.3 kcal/mol to set the enzyme into a catalytically competent form.     

Self‐guided Langevin dynamics study of regulatory interactions in NtrC

Multiple self‐guided Langevin dynamics (SGLD) simulations were performed to examine structural and dynamical properties of the receiver domain of nitrogen regulatory protein C (NtrC^r). SGLD and MD simulations of the phosphorylated active form structure suggest a mostly stable but broad structural ensemble of this protein. The finite difference Poisson–Boltzmann calculations of the pK_a values of the active site residues suggest an increase in the pK_a of His‐84 on phosphorylation of Asp‐54. In SGLD simulations of the phosphorylated active form with charged His‐84, the average position of the regulatory helix α4 is found closer to the starting structure than in simulations with the neutral His‐84. To model the transition pathway, the phosphate group was removed from the simulations. After 7 ns of simulations, the regulatory helix α4 was found approximately halfway between positions in the NMR structures of the active and inactive forms. Removal of the phosphate group stimulated loss of helix α4, suggesting that the pathway of conformational transition may involve partial unfolding mechanism. The study illustrates the potential utility of the SGLD method in studies of the coupling between ligand binding and conformational transitions. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Molecular mechanisms of pH-driven conformational transitions of proteins: insights from continuum electrostatics calculations of acid unfolding

The acid unfolding of staphylococcal nuclease (SNase) is very cooperative (Whitten and García-Moreno, Biochemistry 2000;39:14292-14304). As many as seven hydrogen ions (H+) are bound preferentially by the acid-unfolded state relative to the native (N) state in the pH range 3.2-3.9. To investigate the mechanism of acid unfolding, structure-based pKa calculations were performed with a variety of continuum electrostatic methods. The calculations reproduced successfully the H+ binding properties of the N state between pH 5 and 9, but they systematically overestimated the number of H+ bound upon acid unfolding. The calculated pKa values of all carboxylic residues in the N state were more depressed than they should be. The discrepancy between the observed and the calculated H+ uptake upon acid unfolding was not improved by using high protein dielectric constants, structures relaxed with molecular dynamics, or other empirical modifications implemented previously by others to maximize agreement between measured and calculated pKa values. This suggests an important role for conformational fluctuations of the backbone as important determinants of pKa values of carboxylic groups. Because no global or subglobal conformational changes have been observed previously for SNase under acidic conditions above the acid-unfolding region, these fluctuations must be local. The acid unfolding of SNase does not seem to involve the disruption of the N state by accruement of intramolecular repulsive interactions, nor the protonation of key ion paired carboxylic residues. It is more consistent with modest contributions from many H+ binding groups, with an important role for local conformational fluctuations in the coupling between H+ binding and the global structural transition.     

Backbone relaxation coupled to the ionization of internal groups in proteins: a self-guided Langevin dynamics study

Pathways of structural relaxation triggered by ionization of internal groups in staphylococcal nuclease (SNase) were studied through multiple self-guided Langevin dynamics (SGLD) simulations. Circular dichroism, steady-state Trp fluorescence, and nuclear magnetic resonance spectroscopy have shown previously that variants of SNase with internal Glu, Asp, and Lys at positions 66 or 92, and Arg at position 66, exhibit local reorganization or global unfolding when the internal ionizable group is charged. Except for Arg-66, these internal ionizable groups have unusual pK_a values and are neutral at physiological pH. The structural trends observed in the simulations are in general agreement with experimental observations. The I92D variant, which unfolds globally upon ionization of Asp-92, in simulations often exhibits extensive hydration of the protein core, and sometimes also significant perturbations of the β-barrel. In the crystal structure of the V66R variant, the β1 strand from the β-barrel is domain-swapped; in the simulations, the β1 strand is sometimes partially released. The V66K variant, which in solutions shows reorganization of six residues at the C-terminus of helix α1 and perturbations in the β-barrel structure, exhibits fraying of three residues of helix α1 in one simulation, and perturbations and partial unfolding of three β-strands in a few other simulations. In sharp contrast, very small structural changes were observed in simulations of the wild-type protein. The simulations indicate that charging of internal groups frequently triggers penetration of water into the protein interior. The pK_a values of Asp-92 and Arg-66 calculated with continuum methods on SGLD-relaxed structures reached the normal values in most simulations. Detailed analysis of accuracy and performance of SGLD demonstrates that SGLD outperforms LD in sampling of alternative protein conformations without loss of the accuracy and level of detail characteristic of regular LD.     

Microscopic mechanisms that govern the titration response and pKa values of buried residues in staphylococcal nuclease mutants

To probe the microscopic mechanisms that govern the titration behavior of buried ionizable groups, microsecond explicit solvent molecular dynamics simulations are carried out for several mutants of Staphylococcal nuclease using both fixed charge and polarizable force fields. While the ionization of Asp 66, Glu 66, and Lys 125 lead to enhanced structural fluctuations and partial unfolding of adjacent α‐helical regions, the ionization of Lys 25 causes local unfolding of adjacent β sheets. Using the sampled conformational ensembles, good agreement with experimental pK_a values is obtained with Poisson–Boltzmann calculations using a protein dielectric constant of 2–4 for V66D/E; slightly larger dielectric constants are needed for Lys mutants especially L25K, suggesting that structural responses beyond microseconds are involved in ionization of Lys 25. Overall, the set of unbiased simulations provides insights into the spatial and temporal scales of protein and solvent motions that dictate the diverse titration behaviors of buried protein residues. Proteins 2017; 85:268–281. © 2016 Wiley Periodicals, Inc.     

Origin of the pKa shift of the catalytic lysine in acetoacetate decarboxylase

The pK_a value of Lys115, the catalytic residue in acetoacetate decarboxylate, was calculated using atomic coordinates of the X-ray crystal structure with consideration of the protonation states of all titratable sites in the protein. The calculated pK_a value of Lys115 (pK_a(Lys115)) was unusually low (≈6) in agreement with the experimentally measured value. Although charged residues impact pK_a(Lys115) considerably in the native protein, the significant pK_a(Lys115) downshift in the protein with respect to aqueous solution was mainly due to loss of the solvation energy in the catalytic active site relative to bulk water.     

Electrostatic contributions to the stability of the GCN4 leucine zipper structure

Ion pairs are ubiquitous in X-ray structures of coiled coils, and mutagenesis of charged residues can result in large stability losses. By contrast, pK_a values determined by NMR in solution often predict only small contributions to stability from charge interactions. To help reconcile these results we used triple-resonance NMR to determine pK_a values for all groups that ionize between pH 1 and 13 in the 33 residue leucine zipper fragment, GCN4p. In addition to the native state we also determined comprehensive pK_a values for two models of the GCN4p denatured state: the protein in 6 M urea, and unfolded peptide fragments of the protein in water. Only residues that form ion pairs in multiple X-ray structures of GCN4p gave large pK_a differences between the native and denatured states. Moreover, electrostatic contributions to stability were not equivalent for oppositely charged partners in ion pairs, suggesting that the interactions between a charge and its environment are as important as those within the ion pair. The pH dependence of protein stability calculated from NMR-derived pK_a values agreed with the stability profile measured from equilibrium urea-unfolding experiments as a function of pH. The stability profile was also reproduced with structure-based continuum electrostatic calculations, although contributions to stability were overestimated at the extremes of pH. We consider potential sources of errors in the calculations, and how pK_a predictions could be improved. Our results show that although hydrophobic packing and hydrogen bonding have dominant roles, electrostatic interactions also make significant contributions to the stability of the coiled coil.     

Solvent interactions with the active site of the pig heart cytosolic aspartate aminotransferase.

Aqueous solvent interactions with the chromophoric pyridoxal phosphate prosthetic group of aspartate aminotransferase (EC 2.6.1.1) were analyzed quantitatively with ethylene glycol, glycerol, dimethylsulfoxide (DMSO), sucrose, and xylitol as cosolvents. The smaller cosolvents perturb the visible absorption and visible dichroic spectra of the free enzyme, but this solvent perturbation is not observed with the acidic enzymeglutarate complex. Addition of cosolvents caused an increase in the enzyme's affinity for glutarate. This increase in affinity resulted from an increase in the acidic dissociation constant (pK₂) of the enzyme-glutarate complex. The changes in the acidic dissociation constant of the enzyme-glutarate complex, upon addition of cosolvents, correlate well with the changes observed in the pK_a's of carboxylic acids in comparable solvents. Since these solvents have little effect on the pK_a of the enzyme itself, it is concluded that the increase in affinity is due to a specific solvation effect on a carboxyl group of the enzymebound glutarate, rather than resulting from a conformational change in the protein.     

pH-dependent random coil 1H, 13C,and 15N chemical shifts of the ionizable amino acids: a guide for protein pK a measurements

The pK _a values and charge states of ionizable residues in polypeptides and proteins are frequently determined via NMR-monitored pH titrations. To aid the interpretation of the resulting titration data, we have measured the pH-dependent chemical shifts of nearly all the ¹H, ¹³C, and ¹⁵N nuclei in the seven common ionizable amino acids (X = Asp, Glu, His, Cys, Tyr, Lys, and Arg) within the context of a blocked tripeptide, acetyl-Gly-X-Gly-amide. Alanine amide and N-acetyl alanine were used as models of the N- and C-termini, respectively. Together, this study provides an essentially complete set of pH-dependent intra-residue and nearest-neighbor reference chemical shifts to help guide protein pK _a measurements. These data should also facilitate pH-dependent corrections in algorithms used to predict the chemical shifts of random coil polypeptides. In parallel, deuterium isotope shifts for the side chain ¹⁵N nuclei of His, Lys, and Arg in their positively-charged and neutral states were also measured. Along with previously published results for Asp, Glu, Cys, and Tyr, these deuterium isotope shifts can provide complementary experimental evidence for defining the ionization states of protein residues.     

pH-rate profiles support a general base mechanism for galactokinase (Lactococcus lactis)

Galactokinase (GALK), a member the Leloir pathway for normal galactose metabolism, catalyzes the conversion of α-d-galactose to galactose-1-phosphate. For this investigation, we studied the kinetic mechanism and pH profiles of the enzyme from Lactococcus lactis. Our results show that the mechanism for its reaction is sequential in both directions. Mutant proteins D183A and D183N are inactive (<10 000 fold), supporting the role of Asp183 as a catalytic base that deprotonates the C-1 hydroxyl group of galactose. The pH-k_cat profile of the forward reaction has a pK_a of 6.9 ± 0.2 that likely is due to Asp183. The pH-k_cat/K_Gal profile of the reverse reaction further substantiates this role as it is lacking a key pK_a required for a direct proton transfer mechanism. The R36A and R36N mutant proteins show over 100-fold lower activity than that for the wild-type enzyme, thus suggesting that Arg36 lowers the pK_a of the C-1 hydroxyl to facilitate deprotonation.     

Structure of the D142N mutant of the family 18 chitinase ChiB from Serratia marcescens and its complex with allosamidin

Catalysis by ChiB, a family 18 chitinase from Serratia marcescens, involves a conformational change of Asp142 which is part of a characteristic D₁₄₀XD₁₄₂XE₁₄₄ sequence motif. In the free enzyme Asp142 points towards Asp140, whereas it rotates towards the catalytic acid, Glu144, upon ligand binding. Mutation of Asp142 to Asn reduced k_cat and affinity for allosamidin, a competitive inhibitor. The X-ray structure of the D142N mutant showed that Asn142 points towards Glu144 in the absence of a ligand. The active site also showed other structural adjustments (Tyr10, Ser93) that had previously been observed in the wild-type enzyme upon substrate binding. The X-ray structure of a complex of D142N with allosamidin, a pseudotrisaccharide competitive inhibitor, was essentially identical to that of the wild-type enzyme in complex with the same compound. Thus, the reduced allosamidin affinity in the mutant is not caused by structural changes but solely by the loss of electrostatic interactions with Asp142. The importance of electrostatics was further confirmed by the pH dependence of catalysis and allosamidin inhibition. The pH-dependent apparent affinities for allosamidin were not correlated with k_cat, indicating that it is probably better to view the inhibitor as a mimic of the oxazolinium ion reaction intermediate than as a transition state analogue.