How enzymes harness highly unfavorable proton transfer reactions

Acid–base reactions that are exceedingly unfavorable under standard conditions can be catalytically important at enzyme active sites. For example, in triose phosphate isomerase, a glutamate side chain (nominal pK _a ≈ 4 in solution) can in fact deprotonate a CH group that is vicinal to a carbonyl (pK _a ≈ 18 in solution). This is true because of three distinct interactions: (a) ground state pK _a shifts due to environment polarity and electrostatics; (b) dramatic increases in effective molarity due to optimization of proximity and orientation; and (c) transition state pK _a shifts due to binding interactions and the formation of strong low barrier hydrogen bonds. In this report, we review the literature showing that the sum of these three effects supplies more than enough free energy to push forward proton transfer reactions that under standard conditions are exceedingly nonspontaneous and slow.     

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.     

Self-association of glutamic acid-rich fusion peptide analogs of influenza hemagglutinin in the membrane-mimic environments: Effects of positional difference of glutamic acids on side chain ionization constant and intra- and inter-peptide interactions deduced from NMR and gel electrophoresis measurements

Two glutamic acid-rich fusion peptide analogs of influenza hemagglutinin were synthesized to study the organization of the charged peptides in the membranous media. Fluorescence and gel electrophoresis experiments suggested a loose association between the monomers in the vesicles. A model was built which showed that a positional difference of 3, 7 and 4, 8 results in the exposure of Glu3 and Glu7 side chains to the apolar lipidic core. Supportive results include: first, pK_a values of two pH units higher than reference value in aqueous medium for Glu3 and Glu7 CγH, whereas the deviation of pK_a from the reference value for Glu4 and Glu8 CγH is substantially smaller; second, Hill coefficients of titration shift of these protons indicate anti-cooperativity for Glu3 and Glu7 side chain protons but less so for Glu4 and Glu8, implying a strong electrostatic interaction between Glu3 and Glu7 possibly resulting from their localization in an apolar environment; third, positive and larger titration shift for NH of Glu3 is observed compared to that of Glu4, suggesting stronger hydrogen bond between the NH and the carboxylic group of Glu3 than that of Glu4, consistent with higher degree of exposure to hydrophobic medium for the side chain of Glu3.     

The role of protonation in protein fibrillation

Many proteins fibrillate at low pH despite a high population of charged side chains. Therefore exchange of protons between the fibrillating peptide and its surroundings may play an important role in fibrillation. Here, we use isothermal titration calorimetry to measure exchange of protons between buffer and the peptide hormone glucagon during fibrillation. Glucagon absorbs or releases protons to an extent which allows it to attain a net charge of zero in the fibrillar state, both at acidic and basic pH. Similar results are obtained for lysozyme. This suggests that side chain pK_a values change dramatically in the fibrillar state.     

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.     

Mutation of non-conserved amino acids surrounding catalytic site to shift pH optimum of Bacillus circulans xylanase

Improvement of enzyme function by engineering pH dependence of enzymatic activity is of importance for industrial application of Bacillus circulans xylanases. Target mutation sites were selected by structural alignment between B. circulans xylanase and other xylanases having different pH optima. We selected non-conserved mutant sites within 8 Å from the catalytic residues, to see whether these residues have some role in modulating pK_as of the catalytic residues. We hypothesized that the non-conserved residues which may not have any role in enzyme catalysis might perturb pK_as of the catalytic residues. Change in pK_a of a titratable group due to change in electrostatic potential of a mutation was calculated and the change in pH optimum was predicted from the change in pK_a of the catalytic residues. Our strategy is proved to be useful in selection of promising mutants to shift the pH optimum of the xylanases towards desired side.     

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.     

Comparison of solution and crystal properties of Co(II)-substituted human carbonic anhydrase II

The visible absorption of crystals of Co(II)-substituted human carbonic anhydrase II (Co(II)-HCA II) were measured over a pH range of 6.0-11.0 giving an estimate of pK_a 8.4 for the ionization of the metal-bound water in the crystal. This is higher by about 1.2 pK_a units than the pK_a near 7.2 for Co(II)-CA II in solution. This effect is attributed to a nonspecific ionic strength effect of 1.4 M citrate in the precipitant solution used in the crystal growth. A pK_a of 8.3 for the aqueous ligand of the cobalt was measured for Co(II)-HCA II in solution containing 0.8 M citrate. Citrate is not an inhibitor of the catalytic activity of Co(II)-HCA II and was not observed in crystal structures. The X-ray structures at 1.5-1.6 Å resolution of Co(II)-HCA II were determined for crystals prepared at pH 6.0, 8.5 and 11.0 and revealed no conformational changes of amino-acid side chains as a result of the use of citrate. However, the studies of Co(II)-HCA II did reveal a change in metal coordination from tetrahedral at pH 11 to a coordination consistent with a mixed population of both tetrahedral and penta-coordinate at pH 8.5 to an octahedral geometry characteristic of the oxidized enzyme Co(III)-HCA II at pH 6.0.     

Structure and specific binding of trypsin: comparison of inhibited derivatives and a model for substrate binding

The high-resolution structure of bovine trypsin inhibited with DFP² was determined by Stroud et al. (1971 and R. M. Stroud, L. M. Kay, A. Cooper &; R. E. Dickerson, Abstr. 8th Int. Congr. Biochem. 1970). The experiments reported here were designed to study the specific side-chain binding pocket of trypsin using benzamidine, which is a competitive, specific inhibitor of trypsin. High-resolution electron density syntheses and difference syntheses unambiguously identify the side-chain binding pocket, which normally recognizes and binds the side chains of arginine or lysine during proteolysis. Several important conformational differences in the protein structure are apparent between DIP- and BA-trypsins, and these are discussed with particular reference to inhibition, the binding of lysine and arginine, subsequent orientation of the target at the active site, and the enhancement of tryptic activity towards non-specific substrates seen on binding small alkyl amines or guanidines in the specific binding pocket.The BA-trypsin structure provides a good model for the binding of real substrate side chains to trypsin during catalysis, explaining the sharp trypsin specificity for lysine or arginine side chains (Weinstein &; Doolittle, 1972) and the lack of specificity for stereochemically different basic side chains. Benzamidine is shown to inhibit trypsin by steric interference with the inferred position of good substrates, even when they do not carry any side chain.Apart from the substitution of benzamidine and DIP, the most significant differences between DIP-trypsin and BA-trypsin involve complete repositioning of the side chain of Gln192, alterations in the side chains of Asp102, His57 and Ser195 at the active site, and changes in the solvent structure around this region. The carboxyl group of Asp189, which is responsible for trypsin specificity, shows no movement on binding benzamidine. The amidinium cation of benzamidine forms a salt bridge with Asp189 in BA-trypsin; a similar salt bridge can be constructed between the side chains of model substrates with lysyl or arginyl side chains and Aspl89. The γ-oxygen of Ser190 is displaced by a 120 ° rotation about its αβ bond on binding benzamidine and the binding pocket closes to sandwich the inhibitor ring between the peptide planes of 190–191 and 215–216. These contacts are presumably found in the enzyme-substrate complex with specific substrates.The active site structure at pH 8.0 is discussed with particular reference to the microscopic pK_a values of Asp102 and His57, the pK_a of the Asp-His system, and the mechanistic consequences of these assignments.     

Modulation of an Active-Site Cysteine pKa Allows PDI to Act as a Catalyst of both Disulfide Bond Formation and Isomerization

Protein disulfide isomerase (PDI) plays a central role in disulfide bond formation in the endoplasmic reticulum. It is implicated both in disulfide bond formation and in disulfide bond reduction and isomerization. To be an efficient catalyst of all three reactions requires complex mechanisms. These include mechanisms to modulate the pK_a values of the active-site cysteines of PDI. Here, we examined the role of arginine 120 in modulating the pK_a values of these cysteines. We find that arginine 120 plays a significant role in modulating the pK_a of the C-terminal active-site cysteine in the a domain of PDI and plays a role in determining the reactivity of the N-terminal active-site cysteine but not via direct modulation of its pK_a. Mutation of arginine 120 and the corresponding residue, arginine 461, in the a′ domain severely reduces the ability of PDI to catalyze disulfide bond formation and reduction but enhances the ability to catalyze disulfide bond isomerization due to the formation of more stable PDI-substrate mixed disulfides. These results suggest that the modulation of pK_a of the C-terminal active cysteine by the movement of the side chain of these arginine residues into the active-site locales has evolved to allow PDI to efficiently catalyze both oxidation and isomerization reactions.     

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.     

Aggregation of ferrihaems. Dimerization and protolytic equilibria of protoferrihaem and deuteroferrihaem in aqueous solution

1. The absorption spectra of deutero- and proto-ferrihaem in aqueous solution at 25°C show marked changes with concentration and pH in the Soret band region. Quantitative studies of these phenomena imply that they are associated with ferrihaem dimerization and with protolytic equilibria involving monomeric (M) and dimeric (D) ferrihaem species according to the scheme: [Formula: see text] 2. For deuteroferrihaem we obtain K=1.9×10⁻², pK_a(M)=7.1, pK_a(D)=7.4. Protoferrihaem has a much higher dimerization constant, K=4.5 and pK_a(D)=7.5 (pK_a(M) is not accessible). 3. Possible structural relationships between monomeric and dimeric ferrihaem species in solution are discussed in relation to recent work on the oxo-bridged nature of crystalline ferrihaem dimers.     

Equilibrium constants for association of guanidinium and ammonium ions with oxyanions: The effect of changing basicity of the oxyanion

Using guanidinium and n-butylammonium cations (C⁺) as models for the positively charged side chains in arginine and lysine, we have determined the association constants with various oxyanions by potentiometric titration. For a dibasic acid, H₂A, three association complexes may exist: $\text{K}{}_{\mathrm{<mtext>1M</mtext>}}\text{=}\text{[CHA]}\text{[C}{}^{\mathrm{+}}\text{]}\text{[HA}{}^{\mathrm{?}}\text{]}$; $\text{K}{}_{\mathrm{<mtext>1D</mtext>}}\text{=}\text{[CA}{}^{\mathrm{?}}\text{]}\text{[C}{}^{\mathrm{+}}\text{]}\text{[A}{}^{\mathrm{2?}}\text{]}$; $\text{K}{}_{\mathrm{<mtext>2D</mtext>}}\text{=}\text{[C}{}_2\text{A]}\text{[C}{}^{\mathrm{+}}\text{]}\text{[CA}{}^{\mathrm{?}}\text{]}$. For guanidinium ion and phosphate, K_1M = 1.4, K_1D = 2.6, and K_2D = 5.1. The data for carboxylates indicate that the basicity of the oxyanion does not affect the association constant: acetate, pK_a = 4.8, K_1M = 0.37; formate, pK_a = 3.8, K_1M = 0.32; and chloroacetate, pK_a = 2.9, K_1M = 0.43, all with guanidinium ion. Association constants are also reported for carbonate, dimethylphosphinate, benzylphosphonate, and adenylate anions.     

Structure–activity relationships for 4-anilinoquinoline derivatives as inhibitors of the DNA methyltransferase enzyme DNMT1

A series of 4-anilinoquinoline derivatives related to the known inhibitor SGI-1027, containing side chains of varying pK_a, were prepared by acid-catalysed coupling of the pre-formed side chains with 4-chloroquinolines. The compounds were evaluated for their ability to reduce the level of DNMT1 protein in HCT116 human colon carcinoma cells by Western blotting. With a very strongly basic N-methylpyridinium side chain, only NHCO-linked compounds were effective, whereas less strongly basic ((diaminomethylene)hydrazono)ethyl or 3-methylpyrimidine-2,4-diamine side chains allowed both NHCO- and CONH-linked compounds to show activity. In contrast, the pK_a of the quinoline unit had little apparent influence on activity.     

Flavin oxidation in flavin‐dependent N‐monooxygenases

Siderophore A (SidA) from Aspergillus fumigatus is a flavin‐containing monooxygenase that hydroxylates ornithine (Orn) at the amino group of the side chain. Lysine (Lys) also binds to the active site of SidA; however, hydroxylation is not efficient and H₂O₂ is the main product. The effect of pH on steady‐state kinetic parameters was measured and the results were consistent with Orn binding with the side chain amino group in the neutral form. From the pH dependence on flavin oxidation in the absence of Orn, a pK_a value >9 was determined and assigned to the FAD‐N5 atom. In the presence of Orn, the pH dependence displayed a pK_a value of 6.7 ±0.1 and of 7.70 ±0.10 in the presence of Lys. Q102 interacts with NADPH and, upon mutation to alanine, leads to destabilization of the C4a‐hydroperoxyflavin (FAD_OOH). Flavin oxidation with Q102A showed a pK_a value of ~8.0. The data are consistent with the pK_a of the FAD N5‐atom being modulated to a value >9 in the absence of Orn, which aids in the stabilization of FAD_OOH. Changes in the FAD‐N5 environment lead to a decrease in the pK_a value, which facilitates elimination of H₂O₂ or H₂O. These findings are supported by solvent kinetic isotope effect experiments, which show that proton transfer from the FAD N5‐atom is rate limiting in the absence of a substrate, however, is significantly less rate limiting in the presence of Orn and or Lys.     

Raman dispersion spectroscopy probes heme distortions in deoxyHb-trout IV involved in its T-state Bohr effect

The depolarization ratios of heme protein Raman lines arising from vibrations of the heme group exhibit significant dependence on the excitation wavelength. From the analysis of this depolarization ratio dispersion, one obtains information about symmetry-lowering distortions δQ^Γ of the heme group that can be classified in terms of the symmetry races Γ = A_1g, B_1g, B_2g, and A_2g in D_4h symmetry. The heme-protein interaction can be changed by the protonation of distinct amino acid side chains (i.e., for instance the Bohr groups in hemoglobin derivates), which gives rise to specific static heme distortions for each protonation state. From the Raman dispersion data, it is possible to obtain parameters by fitting to a theoretical expression of the Raman tensor, which provide information on these static distortions and also about the pK values of the involved titrable side chains. We have applied this method to the ν₄ (1,355 cm^-1) and ν₁₀ (1,620 cm^-1) lines of deoxygenated hemoglobin of the fourth component of trout and have measured their depolarization ratio dispersion as a function of pH between 6 and 9. From the pH dependence of the thus derived parameters, we obtain pK values identical to those of the Bohr groups, which were earlier derived from the corresponding O₂-binding isotherms. These are pK_α1 = pK_α2 = 8.5 for the α and pK_β1 = 7.5, pK_β2 = 7.4 for the β chains. We also obtain the specific distortion parameters for each protonation state. As shown in earlier studies, the ν₄ mode mainly probes distortions from interactions between the proximal histidine and atoms of the heme core (i.e., the nitrogens and the C_α atoms of the pyrroles). Group theoretical argumentation allows us to relate specific changes of the imidazole geometry as determined by its tilt and azimuthal angle and the iron-out-of-plane displacement to distinct variations of the normal distortions δQ^Γ derived from the Raman dispersion data. Thus, we found that the pH dependence of the heme distortions δQ^A1g (totally symmetric) and δQ^B1g (asymmetric) is caused by variations of the azimuthal rather than the tilt angle of the Fe-His (F8) bond. In contrast to this, the ν₁₀ line mainly monitors changes resulting from the interaction between peripheral substituents of the porphyrin macrocycle (vinyl). From the pH dependence of the parameters, it is possible to separately identify distortions δQ^Γ affecting the hemes in the α and β chains, respectively. From this, we find that in the α subunit structural changes induced on protonation of the corresponding Bohr groups are mainly transferred via the Fe—N_ε bond and give rise to changes in the azimuthal angle. In the β subunit, however, in addition, structural changes of the heme pocket arise, which most probably result from protonation of the imidazole of the COOH-terminal His (HC3 β). This rearranges the net of H bonds between His HC3 β, Ser (F9 β), and Glu (F7 β).     

The ionic environment determines ribozyme cleavage rate by modulation of nucleobase pK a

Several small ribozymes employ general acid–base catalysis as a mechanism to enhance site-specific RNA cleavage, even though the functional groups on the ribonucleoside building blocks of RNA have pK_a values far removed from physiological pH. The rate of the cleavage reaction is strongly affected by the identity of the metal cation present in the reaction solution; however, the mechanism(s) by which different cations contribute to rate enhancement has not been determined. Using the Neurospora VS ribozyme, we provide evidence that different cations confer particular shifts in the apparent pK_a values of the catalytic nucleobases, which in turn determines the fraction of RNA in the protonation state competent for general acid–base catalysis at a given pH, which determines the observed rate of the cleavage reaction. Despite large differences in observed rates of cleavage in different cations, mathematical models of general acid–base catalysis indicate that k₁, the intrinsic rate of the bond-breaking step, is essentially constant irrespective of the identity of the cation(s) in the reaction solution. Thus, in contrast to models that invoke unique roles for metal ions in ribozyme chemical mechanisms, we find that most, and possibly all, of the ion-specific rate enhancement in the VS ribozyme can be explained solely by the effect of the ions on nucleobase pK_a. The inference that k₁ is essentially constant suggests a resolution of the problem of kinetic ambiguity in favor of a model in which the lower pK_a is that of the general acid and the higher pK_a is that of the general base.     

The transient accumulation of the signaling state of photoactive yellow protein is controlled by the external pH

The signaling state of the photoreceptor photoactive yellow protein is the long-lived intermediate I₂′. The pH dependence of the equilibrium between the transient photocycle intermediates I₂ and I₂′ was investigated. The formation of I₂′ from I₂ is accompanied by a major conformational change. The kinetics and intermediates of the photocycle and of the photoreversal were measured by transient absorption spectroscopy from pH 4.6 to 8.4. Singular value decomposition (SVD) analysis of the data at pH 7 showed the presence of three spectrally distinguishable species: I₁, I₂, and I₂′. Their spectra were determined using the extrapolated difference method. I₂ and I₂′ have electronic absorption spectra, with maxima at 370 ± 5 and 350 ± 5 nm, respectively. Formation of the signaling state is thus associated with a change in the environment of the protonated chromophore. The time courses of the I₁, I₂, and I₂′ intermediates were determined from the wavelength-dependent transient absorbance changes at each pH, assuming that their spectra are pH-independent. After the formation of I₂′ (~2 ms), these three intermediates are in equilibrium and decay together to the initial dark state. The equilibrium between I₂ and I₂′ is pH dependent with a pK_a of 6.4 and with I₂′ the main species above this pK_a. Measurements of the pH dependence of the photoreversal kinetics with a second flash of 355 nm at a delay of 20 ms confirm this pK_a value. I₂ and I₂′ are photoreversed with reversal times of ~55 μs and several hundred microseconds, respectively. The corresponding signal amplitudes are pH dependent with a pK_a of ~6.1. Photoreversal from I₂′ dominates above the pK_a. The transient accumulation of I₂′, the active state of photoactive yellow protein, is thus controlled by the proton concentration. The rate constant k₃ for the recovery to the initial dark state also has a pK_a of ~6.3. This equality of the equilibrium and kinetic pK_a values is not accidental and suggests that k₃ is proportional to [I₂′].     

Origin of the Proton-transfer Step in the Cofactor-free (1H)-3-Hydroxy-4-oxoquinaldine 2,4-Dioxygenase: EFFECT OF THE BASICITY OF AN ACTIVE SITE HIS RESIDUE*

Dioxygenases catalyze a diverse range of chemical reactions that involve the incorporation of oxygen into a substrate and typically use a transition metal or organic cofactor for reaction. Bacterial (1H)-3-hydroxy-4-oxoquinaldine 2,4-dioxygenase (HOD) belongs to a class of oxygenases able to catalyze this energetically unfavorable reaction without any cofactor. In the quinaldine metabolic pathway, HOD breaks down its natural N-heteroaromatic substrate using a mechanism that is still incompletely understood. Experimental and computational approaches were combined to study the initial step of the catalytic cycle. We have investigated the role of the active site His-251/Asp-126 dyad, proposed to be involved in substrate hydroxyl group deprotonation, a critical requirement for subsequent oxygen reaction. The pH profiles obtained under steady-state conditions for the H251A and D126A variants show a strong pH effect on their k_cat and k_cat/K_m constants, with a decrease in k_cat/K_m of 5500- and 9-fold at pH 10.5, respectively. Substrate deprotonation studies under transient-state conditions show that this step is not rate-limiting and yield a pK_a value of ∼7.2 for WT HOD. A large solvent isotope effect was found, and the pK_a value was shifted to ∼8.3 in D₂O. Crystallographic and computational studies reveal that the mutations have a minor effect on substrate positioning. Computational work shows that both His-251 and Asp-126 are essential for the proton transfer driving force of the initial reaction. This multidisciplinary study offers unambiguous support to the view that substrate deprotonation, driven by the His/Asp dyad, is an essential requirement for its activation.