High resolution nuclear magnetic resonance studies of the active site of chymotrypsin. II. Polarization of histidine 57 by substrate analogues and competitive inhibitors

Biomolecular NMR: Past and future

The editors of this special volume suggested this topic, presumably because of the perspective lent by our combined >90-year association with biomolecular NMR. What follows is our personal experience with the evolution of the field, which we hope will illustrate the trajectory of change over the years. As for the future, one can confidently predict that it will involve unexpected advances. Our narrative is colored by our experience in using the NMR Facility for Biomedical Studies at Carnegie-Mellon University (Pittsburgh) and in developing similar facilities at Purdue (1977-1984) and the University of Wisconsin-Madison (1984-). We have enjoyed developing NMR technology and making it available to collaborators and users of these facilities. Our group's association with the Biological Magnetic Resonance data Bank (BMRB) and with the Worldwide Protein Data Bank (wwPDB) has also been rewarding. Of course, many groups contributed to the early growth and development of biomolecular NMR, and our brief personal account certainly omits many important milestones.

Reconstitution of Vanadium Haloperoxidase's Catalytic Activity by Boric Acid-Towards a Potential Biocatalytic Role of Boron

Boron's unusual properties inspired major advances in chemistry. In nature, the existence and importance of boron has been fairly explored (e.g. bacterial signaling, plant development) but its role as biological catalyst was never reported. Here, we show that boric acid [B(OH)₃ ] can restore chloroperoxidase activity of Curvularia inaequalis recombinant apo-haloperoxidase's (HPO) in the presence of hydrogen peroxide and chloride ions. Molecular modeling and semi-empirical PM7 calculations support a thermodynamically highly favored (bio)catalytic mechanism similarly to vanadium haloperoxidases (V-HPO) in which [B(OH)₃ ] is assumedly located in apo-HPO's active site and a monoperoxyborate [B(OH)₃ (OOH)^- ] intermediate is formed and stabilized by interaction with specific active site amino acids leading ultimately to the formation of HOCl. Thus, B(OH)₃ -HPO provides the first evidence towards the future exploitation of boron's role in biological systems.

Proton bridging in the interactions of thrombin with hirudin and its mimics

Thrombin is the pivotal serine protease enzyme in the blood cascade system and thus a target of drug design for control of its activity. The most efficient nonphysiologic inhibitor of thrombin is hirudin, a naturally occurring small protein. Hirudin and its synthetic mimics employ a range of hydrogen bonding, salt bridging, and hydrophobic interactions with thrombin to achieve tight binding with K(i) values in the nano- to femtomolar range. The one-dimensional (1)H nuclear magnetic resonance spectrum recorded at 600 MHz reveals a resonance 15.33 ppm downfield from silanes in complexes between human α-thrombin and r-hirudin in pH 5.6-8.8 buffers and between 5 and 35 °C. There is also a resonance between 15.17 and 15.54 ppm seen in complexes of human α-thrombin with hirunorm IV, hirunorm V, an Nα(Me)Arg peptide, RGD-hirudin, and Nα-2-naphthylsulfonyl-glycyl-DL-4-amidinophenylalanyl-piperidide acetate salt (NAPAP), while there is no such low-field resonance observed in a complex of porcine trypsin and NAPAP. The chemical shifts suggest that these resonances represent H-bonded environments. H-Donor-acceptor distances in the corresponding H-bonds are estimated to be <2.7 Å. Addition of Phe-Pro-Arg-chloromethylketone (PPACK) to a complex of human α-thrombin with r-hirudin results in an additional signal at 18.03 ppm, which is 0.10 ppm upfield from the observed signal [Kovach, I. M., et al. (2009) Biochemistry 48, 7296-7304] for thrombin covalently modified with PPACK. In contrast, the peak at 15.33 ppm remains unchanged. The fractionation factors for the thrombin-hirudin complexes are near 1.0 within 20% error. The most likely site of the short H-bond in complexes of thrombin with the hirudin family of inhibitors is in the hydrophobic patch of the C-terminus of hirudin where Glu(57') and Glu(58') are embedded and interact with Arg(75) and Arg(77) and their solvate water (on thrombin). Glu(57') and Glu(58') present in the hirudin family of inhibitors make up a key binding epitope of fibrinogen, thrombin's prime substrate, which lends substantial interest to the short hydrogen bond as a binding element at the fibrinogen recognition site.

Using the water signal to detect invisible exchanging protons in the catalytic triad of a serine protease

Chemical Exchange Saturation Transfer (CEST) is an MRI approach that can indirectly detect exchange broadened protons that are invisible in traditional NMR spectra. We modified the CEST pulse sequence for use on high-resolution spectrometers and developed a quantitative approach for measuring exchange rates based upon CEST spectra. This new methodology was applied to the rapidly exchanging Hδ1 and Hε2 protons of His57 in the catalytic triad of bovine chymotrypsinogen-A (bCT-A). CEST enabled observation of Hε2 at neutral pH values, and also allowed measurement of solvent exchange rates for His57-Hδ1 and His57-Hε2 across a wide pH range (3-10). Hδ1 exchange was only dependent upon the charge state of the His57 (k (ex,Im+) = 470 s(-1), k (ex,Im) = 50 s(-1)), while Hε2 exchange was found to be catalyzed by hydroxide ion and phosphate base (k(OH)⁻ = 1.7 × 10(10) M(-1) s(-1), K(HPO)²⁻₄ = 1.7 × 10(6) M(-1) s(-1)), reflecting its greater exposure to solute catalysts. Concomitant with the disappearance of the Hε2 signal as the pH was increased above its pK (a), was the appearance of a novel signal (δ = 12 ppm), which we assigned to Hγ of the nearby Ser195 nucleophile, that is hydrogen bonded to Nε2 of neutral His57. The chemical shift of Hγ is about 7 ppm downfield from a typical hydroxyl proton, suggesting a highly polarized O-Hγ bond. The significant alkoxide character of Oγ indicates that Ser195 is preactivated for nucleophilic attack before substrate binding. CEST should be generally useful for mechanistic investigations of many enzymes with labile protons involved in active site chemistry.

Catalytic reaction mechanism of acetylcholinesterase determined by Born-Oppenheimer ab initio QM/MM molecular dynamics simulations

Acetylcholinesterase (AChE) is a remarkably efficient serine hydrolase responsible for the termination of impulse signaling at cholinergic synapses. By employing Born-Oppenheimer molecular dynamics simulations with a B3LYP/6-31G(d) QM/MM potential and the umbrella sampling method, we have characterized its complete catalytic reaction mechanism for hydrolyzing neurotransmitter acetylcholine (ACh) and determined its multistep free-energy reaction profiles for the first time. In both acylation and deacylation reaction stages, the first step involves the nucleophilic attack on the carbonyl carbon, with the triad His447 serving as the general base, and leads to a tetrahedral covalent intermediate stabilized by the oxyanion hole. From the intermediate to the product, the orientation of the His447 ring needs to be adjusted very slightly, and then, the proton transfers from His447 to the product, and the break of the scissile bond happens spontaneously. For the three-pronged oxyanion hole, it only makes two hydrogen bonds with the carbonyl oxygen at either the initial reactant or the final product state, but the third hydrogen bond is formed and stable at all transition and intermediate states during the catalytic process. At the intermediate state of the acylation reaction, a short and low-barrier hydrogen bond (LBHB) is found to be formed between two catalytic triad residues His447 and Glu334, and the spontaneous proton transfer between two residues has been observed. However, it is only about 1-2 kcal/mol stronger than the normal hydrogen bond. In comparison with previous theoretical investigations of the AChE catalytic mechanism, our current study clearly demonstrates the power and advantages of employing Born-Oppenheimer ab initio QM/MM MD simulations in characterizing enzyme reaction mechanisms.

Proton bridging in the interactions of thrombin with small inhibitors

Thrombin is the pivotal serine protease enzyme in the blood cascade system. Phe-Pro-Arg-chloromethylketone (PPACK), phosphate, and phosphonate ester inhibitors form a covalent bond with the active-site Ser of thrombin. PPACK, a mechanism-based inhibitor, and the phosphate/phosphonate esters form adducts that mimic intermediates formed in reactions catalyzed by thrombin. Therefore, the dependence of the inhibition of human alpha-thrombin on the concentration of these inhibitors, pH, and temperature was investigated. The second-order rate constant (ki/Ki) and the inhibition constant (Ki) for inhibition of human alpha-thrombin by PPACK are (1.1 +/- 0.2) x 10(7) M(-1) s(-1) and (2.4 +/- 1.3) x 10(-8) M, respectively, at pH 7.00 in 0.05 M phosphate buffer and 0.15 M NaCl at 25.0 +/- 0.1 degrees C, in good agreement with previous reports. The activation parameters at pH 7.00 in 0.05 M phosphate buffer and 0.15 M NaCl are as follows: DeltaH = 10.6 +/- 0.7 kcal/mol, and DeltaS = 9 +/- 2 cal mol(-1) degrees C(-1). The pH dependence of the second-order rate constants of inhibition is bell-shaped. Values of pKa1 and pKa2 are 7.3 +/- 0.2 and 8.8 +/- 0.3, respectively, at 25.0 +/- 0.1 degrees C. A phosphate and a phosphonate ester inhibitor gave higher values, 7.8 and 8.0 for pKa1 and 9.3 and 8.6 for pKa2, respectively. They inhibit thrombin more than 6 orders of magnitude less efficiently than PPACK does. The deuterium solvent isotope effect for the second-order rate constant at pH 7.0 and 8.3 at 25.0 +/- 0.1 degrees C is unity within experimental error in all three cases, indicating the absence of proton transfer in the rate-determining step for the association of thrombin with the inhibitors, but in a 600 MHz 1H NMR spectrum of the inhibition adduct at pH 6.7 and 30 degrees C, a peak at 18.10 ppm with respect to TSP appears with PPACK, which is absent in the 1H NMR spectrum of a solution of the enzyme between pH 5.3 and 8.5. The peak at low field is an indication of the presence of a short-strong hydrogen bond (SSHB) at the active site in the adduct. The deuterium isotope effect on this hydrogen bridge is 2.2 +/- 0.2 (phi = 0.45). The presence of an SSHB is also established with a signal at 17.34 ppm for a dealkylated phosphate adduct of thrombin.

Unconventional serine proteases: variations on the catalytic Ser/His/Asp triad configuration

Serine proteases comprise nearly one-third of all known proteases identified to date and play crucial roles in a wide variety of cellular as well as extracellular functions, including the process of blood clotting, protein digestion, cell signaling, inflammation, and protein processing. Their hallmark is that they contain the so-called "classical" catalytic Ser/His/Asp triad. Although the classical serine proteases are the most widespread in nature, there exist a variety of "nonclassical" serine proteases where variations to the catalytic triad are observed. Such variations include the triads Ser/His/Glu, Ser/His/His, and Ser/Glu/Asp, and include the dyads Ser/Lys and Ser/His. Other variations are seen with certain serine and threonine peptidases of the Ntn hydrolase superfamily that carry out catalysis with a single active site residue. This work discusses the structure and function of these novel serine proteases and threonine proteases and how their catalytic machinery differs from the prototypic serine protease class.

NMR Study of the Inhibition of Pepsin by Glyoxal Inhibitors: Mechanism of Tetrahedral Intermediate Stabilization by the Aspartyl Proteases

Z-Ala-Ala-Phe-glyoxal (where Z is benzyloxycarbonyl) has been shown to be a competitive inhibitor of pepsin with a Ki = 89 +/- 24 nM at pH 2.0 and 25 degrees C. Both the ketone carbon (R13COCHO) and the aldehyde carbon (RCO13CHO) of the glyoxal group of Z-Ala-Ala-Phe-glyoxal have been 13C-enriched. Using 13C NMR, it has been shown that when the inhibitor is bound to pepsin, the glyoxal keto and aldehyde carbons give signals at 98.8 and 90.9 ppm, respectively. This demonstrates that pepsin binds and preferentially stabilizes the fully hydrated form of the glyoxal inhibitor Z-Ala-Ala-Phe-glyoxal. From 13C NMR pH studies with glyoxal inhibitor, we obtain no evidence for its hemiketal or hemiacetal hydroxyl groups ionizing to give oxyanions. We conclude that if an oxyanion is formed its pKa must be >8.0. Using 1H NMR, we observe four hydrogen bonds in free pepsin and in pepsin/Z-Ala-Ala-Phe-glyoxal complexes. In the pepsin/pepstatin complex an additional hydrogen bond is formed. We examine the effect of pH on hydrogen bond formation, but we do not find any evidence for low-barrier hydrogen bond formation in the inhibitor complexes. We conclude that the primary role of hydrogen bonding to catalytic tetrahedral intermediates in the aspartyl proteases is to correctly orientate the tetrahedral intermediate for catalysis.

Subangstrom crystallography reveals that short ionic hydrogen bonds, and not a His-Asp low-barrier hydrogen bond, stabilize the transition state in serine protease catalysis

To address questions regarding the mechanism of serine protease catalysis, we have solved two X-ray crystal structures of alpha-lytic protease (alphaLP) that mimic aspects of the transition states: alphaLP at pH 5 (0.82 A resolution) and alphaLP bound to the peptidyl boronic acid inhibitor, MeOSuc-Ala-Ala-Pro-boroVal (0.90 A resolution). Based on these structures, there is no evidence of, or requirement for, histidine-flipping during the acylation step of the reaction. Rather, our data suggests that upon protonation of His57, Ser195 undergoes a conformational change that destabilizes the His57-Ser195 hydrogen bond, preventing the back-reaction. In both structures the His57-Asp102 hydrogen bond in the catalytic triad is a normal ionic hydrogen bond, and not a low-barrier hydrogen bond (LBHB) as previously hypothesized. We propose that the enzyme has evolved a network of relatively short hydrogen bonds that collectively stabilize the transition states. In particular, a short ionic hydrogen bond (SIHB) between His57 Nepsilon2 and the substrate's leaving group may promote forward progression of the TI1-to-acylenzyme reaction. We provide experimental evidence that refutes use of either a short donor-acceptor distance or a downfield 1H chemical shift as sole indicators of a LBHB.

Low-barrier hydrogen bond hypothesis in the catalytic triad residue of serine proteases: correlation between structural rearrangement and chemical shifts in the acylation process

To elucidate the catalytic advantage of the low-barrier hydrogen bond (LBHB), we analyze the hydrogen bonding network of the catalytic triad (His57-Asp102-Ser195) of serine protease trypsin, one of the best examples of the LBHB reaction mechanism. Especially, we focus on the correlation between the change of the chemical shifts and the structural rearrangement of the active site in the acylation process. To clarify LBHB, we evaluate the two complementary properties. First, we calculate the NMR chemical shifts of the imidazole ring of His57 by the gauge-including atomic orbital (GIAO) approach within the ab initio QM/MM framework. Second, the free energy profile of the proton transfer from His57 to Asp102 in the tetrahedral intermediate is obtained by ab initio QM/MM calculations combined with molecular dynamics free energy perturbation (MD-FEP) simulations. The present analyses reveal that the calculated shifts reasonably reproduce the observed values for (1)H chemical shift of H(epsilon)(1) and H(delta)(1) in His57. The (15)N and (13)C chemical shifts are also consistent with the experiments. It is also shown that the proton between His57 and Asp102 is localized at the His57 side. This largely downfield chemical shift is originated from the strong electrostatic interaction, not a covalent-like bonding character between His57 and Asp102. Also, it is proved that a slight downfield character of H(epsilon)(1) is originated from a electrostatic interaction between His57 and the backbone carbonyl group of Val213 and Ser214. These downfield chemical shifts are observed only when the tetrahedral intermediate is formed in the acylation process.

Theoretical studies on the deacylation step of serine protease catalysis in the gas phase, in solution, and in elastase

The deacylation step of serine protease catalysis is studied using DFT and ab initio QM/MM calculations combined with MD/umbrella sampling calculations. Free energies of the entire reaction are calculated in the gas phase, in a continuum solvent, and in the enzyme elastase. The calculations show that a concerted mechanism in the gas phase is replaced by a stepwise mechanism when solvent effects or an acetate ion are added to the reference system, with the tetrahedral intermediate being a shallow minimum on the free energy surface. In the enzyme, the tetrahedral intermediate is a relatively stable species ( approximately 7 kcal/mol lower in energy than the transition state), mainly due to the electrostatic effects of the oxyanion hole and Asp102. It is formed in the first step of the reaction, as a result of a proton transfer from the nucleophilic water to His57 and of an attack of the remaining hydroxyl on the ester carbonyl. This is the rate-determining step of the reaction, which requires approximately 22 kcal/mol for activation, approximately 5 kcal/mol less than the reference reaction in water. In the second stage of the reaction, only small energy barriers are detected to facilitate the proton transfer from His57 to Ser195 and the breakdown of the tetrahedral intermediate. Those are attributed mainly to a movement of Ser195 and to a rotation of the His57 side chain. During the rotation, the imidazolium ion is stabilized by a strong H-bond with Asp102, and the C(epsilon)(1)-H...O H-bond with Ser214 is replaced by one with Thr213, suggesting that a "ring-flip mechanism" is not necessary as a driving force for the reaction. The movements of His57 and Ser195 are highly correlated with rearrangements of the binding site, suggesting that product release may be implicated in the deacylation process.

Role of Asp102 in the catalytic relay system of serine proteases: a theoretical study

The role of Asp102 in the catalytic relay system of serine proteases is studied theoretically by calculating the free energy profiles of the single proton-transfer reaction by the Asn102 mutant trypsin and the concerted double proton-transfer reaction (so-called the charge-relay mechanism) of the wild-type trypsin. For each reaction, the reaction free energy profile of the rate-determining step (the tetrahedral intermediate formation step) is calculated by using ab initio QM/MM electronic structure calculations combined with molecular dynamics-free energy perturbation method. In the mutant reaction, the free energy monotonically increases along the reaction path. The rate-determining step of the mutant reaction is the formation of tetrahedral intermediate complex, not the base (His57) abstraction of the proton from Ser195. In contrast to the single proton-transfer reaction of the wild-type, MD simulations of the enzyme-substrate complex show that the catalytically favorable alignment of the relay system (the hydrogen bonding network between the mutant triad, His57, Asn102, and Ser195) is rarely observed even in the presence of a substrate at the active site. In the double proton-transfer reaction, the energy barrier is observed at the proton abstraction step, which corresponds to the rate-determining step of the single proton-transfer reaction of the wild-type. Although both reaction profiles show an increase of the activation barrier by several kcals/mol, these increases have different energetic origins: a large energetic loss of the electrostatic stabilization between His57 and Asn102 in the mutant reaction, while the lack of stabilization by the protein environment in the double proton-transfer reaction. Comparing the present results with the single proton transfer of the wild-type, Asp102 is proven to play two important roles in the catalytic process. One is to stabilize the protonated His57, or ionic intermediate, formed during the acylation, and the other is to fix the configuration around the active site, which is favorable to promote the catalytic process. These two factors are closely related to each other and are indispensable for the efficient catalysis. Also the present calculations suggest the importance of the remote site interaction between His57 and Val213-Ser214 at the catalytic transition state.

Theoretical perspectives on the reaction mechanism of serine proteases: the reaction free energy profiles of the acylation process

The reaction mechanism of serine proteases (trypsin), which catalyze peptide hydrolysis, is studied theoretically by ab initio QM/MM electronic structure calculations combined with Molecular Dynamics-Free Energy Perturbation calculations. We have calculated the entire reaction free energy profiles of the first reaction step of this enzyme (acylation process). The present calculations show that the rate-determining step of the acylation is the formation of the tetrahedral intermediate, and the breakdown of this intermediate has a small energy barrier. The calculated activation free energy for the acylation is approximately 17.8 kcal/mol at QM/MM MP2/(aug)-cc-pVDZ//HF/6-31(+)G/AMBER level, and this reaction is an exothermic process. MD simulations of the enzyme-substrate (ES) complex and the free enzyme in aqueous phase show that the substrate binding induces slight conformational changes around the active site, which favor the alignment of the reactive fragments (His57, Asp102, and Ser195) together in a reactive orientation. It is also shown that the proton transfer from Ser195 to His57 and the nucleophilic attack of Ser195 to the carbonyl carbon of the scissile bond of the substrate occur in a concerted manner. In this reaction, protein environment plays a crucial role to lowering the activation free energy by stabilizing the tetrahedral intermediate compared to the ES complex. The polarization energy calculations show that the enzyme active site is in a very polar environment because of the polar main chain contributions of protein. Also, the ground-state destabilization effect (steric strain) is not a major catalytic factor. The most important catalytic factor of stabilizing the tetrahedral intermediate is the electrostatic interaction between the active site and particular regions of protein: the main chain NH groups in Gly193 and Ser195 (so-called oxyanion hole region) stabilize negative charge generated on the carbonyl oxygen of the scissile bond, and the main chain carbonyl groups in Ile212 approximately Ser214 stabilize a positive charge generated on the imidazole ring of His57.

Investigation of the role of the histidine-aspartate pair in the human exonuclease III-like abasic endonuclease, Ape1

Hydrogen bonded histidine-aspartate (His-Asp) pairs are critical constituents in several key enzymatic reactions. To date, the role that these pairs play in catalysis is best understood in serine and trypsin-like proteases, where structural and biochemical NMR studies have revealed important pK(a) values and hydrogen bonding patterns within the catalytic pocket. However, the role of the His-Asp pair in metal-assisted catalysis is less clear. Here, we apply liquid-state NMR to investigate the role of a critical histidine residue of apurinic endonuclease 1 (Ape1), a human DNA repair enzyme that cleaves adjacent to abasic sites in DNA using one or more divalent cations and an active-site His-Asp pair. The results of these studies suggest that the Ape1 His-Asp pair does not function as either a general base catalyst or a metal ligand. Rather, the pair likely stabilizes the pentavalent transition state necessary for phospho-transfer.

Tautomerism, acid-base equilibria, and H-bonding of the six histidines in subtilisin BPN' by NMR

We have determined by (15)N, (1)H, and (13)C NMR, the chemical behavior of the six histidines in subtilisin BPN' and their PMSF and peptide boronic acid complexes in aqueous solution as a function of pH in the range of from 5 to 11, and have assigned every (15)N, (1)H, C(epsilon 1), and C(delta2) resonance of all His side chains in resting enzyme. Four of the six histidine residues (17, 39, 67, and 226) are neutrally charged and do not titrate. One histidine (238), located on the protein surface, titrates with pK(a) = 7.30 +/- 0.03 at 25 degrees C, having rapid proton exchange, but restricted mobility. The active site histidine (64) in mutant N155A titrates with a pK(a) value of 7.9 +/- 0.3 and sluggish proton exchange behavior, as shown by two-site exchange computer lineshape simulation. His 64 in resting enzyme contains an extremely high C(epsilon 1)-H proton chemical shift of 9.30 parts per million (ppm) owing to a conserved C(epsilon 1)-H(.)O=C H-bond from the active site imidazole to a backbone carbonyl group, which is found in all known serine proteases representing all four superfamilies. Only His 226, and His 64 at high pH, exist as the rare N(delta1)-H tautomer, exhibiting (13)C(delta1) chemical shifts approximately 9 ppm higher than those for N(epsilon 2)-H tautomers. His 64 in the PMSF complex, unlike that in the resting enzyme, is highly mobile in its low pH form, as shown by (15)N-(1)H NOE effects, and titrates with rapid proton exchange kinetics linked to a pK(a) value of 7.47 +/- 0.02.

Quantum chemical calculations on structural models of the catalytic site of chymotrypsin: comparison of calculated results with experimental data from NMR spectroscopy

Hybrid density functional quantum mechanical calculations were used to study the strength of the hydrogen bond between His(57) N(delta)(1) and Asp(102) O(delta)(1) in chymotrypsin and how it changes along the reaction coordinate. Comparison of experimental shifts with the results of chemical shift calculations on a variety of small molecules, including species containing very strong hydrogen bonds, has validated the overall approach and provided the means for calibrating and correcting the calculated values. Models of the active site of chymotrypsin in its resting state and tetrahedral intermediate state were derived from high-resolution X-ray structures. The distance between His(57) N(delta)(1) and Asp(102) O(delta)(1) in each model was varied between 2.77 A (weak hydrogen bond) and 2.50 A (extremely strong hydrogen bond), and the one-dimensional potential energy surface of the hydrogen-bonded proton (or deuteron/triton) was determined. The zero-point energy, probability distribution, and chemical shift were determined for each distance. Calculated values for NMR chemical shifts, NMR chemical shift differences between (1)H and (3)H, and (2)H/(1)H fractionation factors were compared with published experimental values. Energies provided by the calculations indicated that the hydrogen bond between His(57) N(delta)(1) and Asp(102) O(delta)(1) in the chymotrypsin active site increases in strength by 11 kcal mol(-)(1) in going from the resting state of the enzyme to the tetrahedral intermediate state. This result confirms the hypothesis that the strengthened hydrogen bond plays an important role in lowering the energy of the transition state and, hence, in the catalytic efficiency of the enzyme. Models of the transition state that best fit the experimental data are consistent with a "strong" hydrogen bond between His(57) N(delta)(1) and Asp(102) O(delta)(1) but apparently not a "low-barrier" or "very strong" hydrogen bond.

Formation of a trypsin-borate-4-aminobutanol ternary complex

The formation of ternary complexes involving serine proteases, borate, and an alcohol has important implications for understanding the physiological actions of borate and for the development of tight binding inhibitors for this class of enzymes. Recent studies of a related enzyme, gamma-glutamyl transpeptidase, which is subject to inhibition by a labile serine/borate mixture, have demonstrated that construction of a non-labile boronate analogue results in an inhibitor with nearly 10(5)-fold greater potency. To evaluate the generalization of this biochemistry to serine proteases, we have observed the ternary complex formed from 4-aminobutanol, borate, and trypsin. A combination of (11)B and (1)H NMR and spectrophotometric assays using acetylarginine p-nitroanilide (Ac-Arg-pNA) as the chromogenic substrate all indicate a cooperative binding interaction in which the borate is esterified by the oxygen atoms of the 4-aminobutanol and trypsin residue Ser(195). Two downfield-shifted proton resonances at 15.5 and 16.6 ppm are proposed to arise from the labile imidazolium protons on His(57), indicating a salt bridge interaction with the negatively charged borate. A cooperativity parameter alpha of 0.2 is derived from the assays. These results provide the first direct evidence for formation of a ternary complex involving a serine protease, borate, and an alcohol, and suggest that this represents a general approach for the development of tight binding ligands.

A catalytic diad involved in substrate-assisted catalysis: NMR study of hydrogen bonding and dynamics at the active site of phosphatidylinositol-specific phospholipase C

Phosphatidylinositol-specific phospholipase Cs (PI-PLCs, EC 3.1.4.10) are ubiquitous enzymes that cleave phosphatidylinositol or phosphorylated derivatives, generating second messengers in eukaryotic cells. A catalytic diad at the active site of Bacillus cereus PI-PLC composed of aspartate-274 and histidine-32 was postulated from the crystal structure to form a catalytic triad with the 2-OH group of the substrate [Heinz, D. W., et al. (1995) EMBO J. 14, 3855-3863]. This catalytic diad has been observed directly by proton NMR. The single low-field line in the 1H NMR spectrum is assigned by site-directed mutagenesis: The peak is present in the wild type but absent in the mutants H32A and D274A, and arises from the histidine Hdelta1 forming the Asp274-His32 hydrogen bond. This hydrogen is solvent-accessible, and exchanges slowly with H2O on the NMR time scale. The position of the low-field peak shifts from 16.3 to 13.8 ppm as the pH is varied from 4 to 9, reflecting a pKa of 8.0 at 6 degrees C, which is identified with the pKa of His32. The Hdelta1 signal is modulated by rapid exchange of the Hepsilon2 with the solvent. Estimates of the exchange rate as a function of pH and protection factors are derived from a line shape analysis. The NMR behavior is remarkably similar to that of the serine proteases. The postulated function of the Asp274-His32 diad is to hydrogen-bond with the 2-OH of phosphatidylinositol (PI) substrate to form a catalytic triad analogous to Asp-His-Ser of serine proteases. This is an example of substrate-assisted catalysis where the substrate provides the catalytic nucleophile of the triad. This hydrogen bond becomes shorter as the imidazole is protonated, suggesting it is stronger in the transition state, contributing further to the catalytic efficiency. The hydrogen bond fits the NMR criteria for a short, strong hydrogen bond, i.e., a highly deshielded proton resonance, bond length of 2.64 +/- 0.04 A at pH 6 measured by NMR, a D/H fractionation factor significantly lower than 1.0, and a protection factor > or = 100.

Short, strong hydrogen bonds at the active site of human acetylcholinesterase: proton NMR studies

Cholinesterases use a Glu-His-Ser catalytic triad to enhance the nucleophilicity of the catalytic serine. We have previously shown by proton NMR that horse serum butyryl cholinesterase, like serine proteases, forms a short, strong hydrogen bond (SSHB) between the Glu-His pair upon binding mechanism-based inhibitors, which form tetrahedral adducts, analogous to the tetrahedral intermediates in catalysis [Viragh, C., et al. (2000) Biochemistry 39, 16200-16205]. We now extend these studies to human acetylcholinesterase, a 136 kDa homodimer. The free enzyme at pH 7.5 shows a proton resonance at 14.4 ppm assigned to an imidazole NH of the active-site histidine, but no deshielded proton resonances between 15 and 21 ppm. Addition of a 3-fold excess of the mechanism-based inhibitor m-(N,N,N-trimethylammonio)trifluoroacetophenone (TMTFA) induced the complete loss of the 14.4 ppm signal and the appearance of a broad, deshielded resonance of equal intensity with a chemical shift delta of 17.8 ppm and a D/H fractionation factor phi of 0.76 +/- 0.10, consistent with a SSHB between Glu and His of the catalytic triad. From an empirical correlation of delta with hydrogen bond lengths in small crystalline compounds, the length of this SSHB is 2.62 +/- 0.02 A, in agreement with the length of 2.63 +/- 0.03 A, independently obtained from phi. Upon addition of a 3-fold excess of the mechanism-based inhibitor 4-nitrophenyl diethyl phosphate (paraoxon) to the free enzyme at pH 7.5, and subsequent deethylation, two deshielded resonances of unequal intensity appeared at 16.6 and 15.5 ppm, consistent with SSHBs with lengths of 2.63 +/- 0.02 and 2.65 +/- 0.02 A, respectively, suggesting conformational heterogeneity of the active-site histidine as a hydrogen bond donor to either Glu-327 of the catalytic triad or to Glu-199, also in the active site. Conformational heterogeneity was confirmed with the methylphosphonate ester anion adduct of the active-site serine, which showed two deshielded resonances of equal intensity at 16.5 and 15.8 ppm with phi values of 0.47 +/- 0.10 and 0.49 +/- 0.10 corresponding to average hydrogen bond lengths of 2.59 +/- 0.04 and 2.61 +/- 0.04 A, respectively. Similarly, lowering the pH of the free enzyme to 5.1 to protonate the active-site histidine (pK(a) = 6.0 +/- 0.4) resulted in the appearance of two deshielded resonances, at 17.7 and 16.4 ppm, consistent with SSHBs with lengths of 2.62 +/- 0.02 and 2.63 +/- 0.02 A, respectively. The NMR-derived distances agree with those found in the X-ray structures of the homologous acetylcholinesterase from Torpedo californica complexed with TMTFA (2.66 +/- 0.28 A) and sarin (2.53 +/- 0.26 A) and at low pH (2.52 +/- 0.25 A). However, the order of magnitude greater precision of the NMR-derived distances establishes the presence of SSHBs at the active site of acetylcholinesterase, and detect conformational heterogeneity of the active-site histidine. We suggest that the high catalytic power of cholinesterases results in part from the formation of a SSHB between Glu and His of the catalytic triad.

NMR evidence for a short, strong hydrogen bond at the active site of a cholinesterase

Cholinesterases (ChE), use a Glu-His-Ser catalytic triad to enhance the nucleophilicity of the catalytic serine. It has been shown that serine proteases, which employ an Asp-His-Ser catalytic triad for optimal catalytic efficiency, decrease the hydrogen bonding distance between the Asp-His pair to form a short, strong hydrogen bond (SSHB) upon binding mechanism-based inhibitors, which form tetrahedral Ser-adducts, analogous to the tetrahedral intermediates in catalysis, or at low pH when the histidine is protonated [Cassidy, C. S., Lin, J., Frey, P. A. (1997) Biochemistry 36, 4576-4584]. Two types of mechanism-based inhibitors were bound to pure equine butyrylcholinesterase (BChE), a 364 kDa homotetramer, and the complexes were studied by (1)H NMR at 600 MHz and 25-37 degrees C. The downfield region of the (1)H NMR spectrum of free BChE at pH 7.5 showed a broad, weak, deshielded resonance with a chemical shift, delta = 16.1 ppm, ascribed to a small amount of the histidine-protonated form. Upon addition of a 3-fold excess of diethyl 4-nitrophenyl phosphate (paraoxon) and subsequent dealkylation, the broad 16.1 ppm resonance increased in intensity 4.7-fold, and yielded a D/H fractionation factor phi = 0.72+/-0.10 consistent with a SSHB between Glu and His of the catalytic triad. From an empirical correlation of delta with hydrogen-bond length in small crystalline compounds, the length of this SSBH is 2.64+/-0.04 A, in agreement with the length of 2.62+/-0.02 A independently obtained from phi. The addition of a 3-fold excess of m-(N,N, N-trimethylammonio)trifluoroacetophenone to BChE yielded no signal at 16.1 ppm, and a 640 Hz broad, highly deshielded proton resonance with a chemical shift delta = 18.1 ppm and a D/H fractionation factor phi = 0.63+/-0.10, also consistent with a SSHB. The length of this SSHB is calculated to be 2.62+/-0.04 A from delta and 2.59+/-0.03 A from phi. These NMR-derived distances agree with those found in the X-ray structures of the homologous acetylcholinesterase complexed with the same mechanism-based inhibitors, 2.60+/-0.22 and 2.66+/-0.28 A. However, the order of magnitude greater precision of the NMR-derived distances establish the presence of SSHBs. We suggest that ChEs achieve their remarkable catalytic power in ester hydrolysis, in part, due to the formation of a SSHB between Glu and His of the catalytic triad.

Fractionation factors and activation energies for exchange of the low barrier hydrogen bonding proton in peptidyl trifluoromethyl ketone complexes of chymotrypsin

NMR investigations have been carried out of complexes between bovine chymotrypsin Aalpha and a series of four peptidyl trifluoromethyl ketones, listed here in order of increasing affinity for chymotrypsin: N-Acetyl-L-Phe-CF3, N-Acetyl-Gly-L-Phe-CF3, N-Acetyl-L-Val-L-Phe-CF3, and N-Acetyl-L-Leu-L-Phe-CF3. The D/H fractionation factors (phi) for the hydrogen in the H-bond between His 57 and Asp 102 (His 57-Hdelta1) in these four complexes at 5 degreesC were in the range phi = 0.32-0.43, expected for a low-barrier hydrogen bond. For this series of complexes, measurements also were made of the chemical shifts of His 57-Hepsilon1 (delta2,2-dimethylsilapentane-5-sulfonic acid 8.97-9. 18), the exchange rate of the His 57-Hdelta1 proton with bulk water protons (284-12.4 s-1), and the activation enthalpies for this hydrogen exchange (14.7-19.4 kcal.mol-1). It was found that the previously noted correlations between the inhibition constants (Ki 170-1.2 microM) and the chemical shifts of His 57-Hdelta1 (delta2, 2-dimethylsilapentane-5-sulfonic acid 18.61-18.95) for this series of peptidyl trifluoromethyl ketones with chymotrypsin [Lin, J., Cassidy, C. S. & Frey, P. A. (1998) Biochemistry 37, 11940-11948] could be extended to include the fractionation factors, hydrogen exchange rates, and hydrogen exchange activation enthalpies. The results support the proposal of low barrier hydrogen bond-facilitated general base catalysis in the addition of Ser 195 to the peptidyl carbonyl group of substrates in the mechanism of chymotrypsin-catalyzed peptide hydrolysis. Trends in the enthalpies for hydrogen exchange and the fractionation factors are consistent with a strong, double-minimum or single-well potential hydrogen bond in the strongest complexes. The lifetimes of His 57-Hdelta1, which is solvent shielded in these complexes, track the strength of the hydrogen bond. Because these lifetimes are orders of magnitude shorter than those of the complexes themselves, the enzyme must have a pathway for hydrogen exchange at this site that is independent of dissociation of the complexes.

Correlations of the basicity of His 57 with transition state analogue binding, substrate reactivity, and the strength of the low-barrier hydrogen bond in chymotrypsin

The basicity of His 57-Nepsilon2 within the low-barrier hydrogen-bonded (LBHB) diad His 57-Asp 102 and the 1H NMR chemical shift of the LBHB proton in tetrahedral, hemiketal complexes of chymotrypsin with peptidyl trifluoromethyl ketones (peptidyl-TFKs) have been studied. The following results were obtained with various peptidyl-TFKs at 5 degrees C: N-Ac-Gly-DL-Phe-CF3, pKa = 11.1 and deltaLBHB = 18.7 ppm; N-Ac-L-Val-DL-Phe-CF3, pKa = 11.8 and deltaLBHB = 18.9 ppm; N-Ac-L-Leu-DL-Val-CF3, pKa = 10.3 and deltaLBHB = 18.9 ppm; and N-Ac-L-Leu-DL-naphthyl-CF3, pKa = 10.9 and deltaLBHB = 19.0 ppm. Results for peptidyl-TFKs with Phe in the P1 position and N-Ac, N-Ac-Gly, N-Ac-L-Val, and N-Ac-L-Leu in the P2 position were well correlated with literature values for inhibition constants Ki and kcat/Km for the corresponding peptidyl methyl esters. The plot of log Ki versus the apparent pKa of His 57-Nepsilon2 displayed a slope of -0.77, and that of log kcat/Km for peptidyl methyl esters versus the pKa of His 57-Nepsilon2 in corresponding peptidyl-TFK complexes gave a slope of 0.68. The slope of a plot of pKa versus deltaLBHB was 3.7, and that of log kcat/Km for peptidyl methyl ester substrates versus deltaLBHB for the corresponding peptidyl-TFK-chymotrypsin complexes was 2.7. A plot of log Ki versus deltaLBHB displayed a slope of -3.0. These plots indicated that the pKa of His 57 and substrate reactivity were correlated with increasing strength of the low-barrier hydrogen bond. The apparent pKa of His 57-Nepsilon2 for the chymotrypsin-N-Ac-L-Leu-DL-Phe-CF3 complex is 10.6 at 25 degrees C, whereas it is 12.0 at 5 degrees C [Cassidy, C. S., Lin, J. L., and Frey, P. A. (1997) Biochemistry 36, 4576-4584]. The apparent discrepancy is likely to be due to a temperature dependence in the cooperative ionization of His 57 in peptidyl-TFK complexes, which appears to be coupled to inhibitor dissociation, hydration and ionization of free peptidyl-TFK, ionization of Ile 16, and a conformational change.

Characterization of a buried neutral histidine in Bacillus circulans xylanase: internal dynamics and interaction with a bound water molecule

NMR spectroscopy was used to characterize the dynamic behavior of His149 in Bacillus circulans xylanase (BCX) and its interaction with an internal water molecule. Rate constants for the specific acid- and base-catalyzed exchange following bimolecular kinetics (EX2) of the nitrogen-bonded H epsilon 2 of this buried, neutral histidine were determined. At pDmin 7.0 and 30 degrees C, the lifetime for this proton is 9.9 h, corresponding to a protection factor of approximately 10(7) relative to that predicted for an exposed histidine. The apparent activation energies measured for specific acid and base catalysis (7.0 and 17.4 kcal/mol) indicate that exchange occurs via local structural fluctuations. Consistent with its buried environment, the N epsilon 2-H bond vector of His149 shows restricted mobility, as evidenced by an order parameter S2 = 0.83 determined from 15N relaxation measurements. The crystal structure of BCX reveals that a conserved, buried water hydrogen-bonds to the H epsilon 2 of His149. Strong support for this interaction in solution is provided by the observation of a negative nuclear Overhauser effect (NOE) and positive rotating-frame Overhauser effect (ROE) between His149 H epsilon 2 and a water molecule with the same chemical shift as the bulk solvent. However, the chemical shift of H epsilon 2 (12.2 ppm) and a D/H fractionation factor close to unity (0.89 +/- 0.02) indicate that this is not a so-called low-barrier hydrogen bond. Lower and upper bounds on the lifetime of the internal water are estimated to be 10(-8) and 10(-3) s. Therefore the chemical exchange of solvent protons with those of His149 H epsilon 2 and the diffusion or physical exchange of the internal water to which the histidine is hydrogen-bonded differ in rate by over 7 orders of magnitude.

A 1H NMR study of structurally relevant inter-segmental hydrogen bond in cytochrome c

NMR signal arising from His 26 N(epsilon)H proton in horse and tuna ferrocytochromes c has been assigned. This His residue is highly conserved in most mitochondrial cytochromes c and X-ray crystallographic studies strongly suggested that its side-chain imidazole participates in an internal hydrogen bond network which is relevant to the stability of the non-helical protein folding near the heme active site. The shift and line width of the assigned signal indicated that this NH hydrogen is indeed involved in an internal hydrogen bond. On the basis of the X-ray crystal structures, the carbonyl oxygen of the residue at 44 is thought to act as a proton-acceptor for this hydrogen. The observation of nuclear Overhauser effect correlation between His 26 C(epsilon)H and Asn 31 main-chain amide NH proton signals in the present proteins also demonstrated the formation of the hydrogen bond between these residues. Consequently, the presence of a unique triad hydrogen bond network in these cytochromes c in solution has been confirmed. Taking advantage of the sensitivity of His 26 N(epsilon)H proton signal to the structural properties of this hydrogen bond network, influences of the presence of high concentration of salt or various concentrations of denaturant on the protein folding were inferred from the analysis of the NMR spectral parameters of the signal.

Low barrier hydrogen bond is absent in the catalytic triads in the ground state but Is present in a transition-state complex in the prolyl oligopeptidase family of serine proteases

High frequency proton NMR spectra for two members of the prolyl oligopeptidase class of serine proteases, prolyl oligopeptidase and oligopeptidase B, showed that resonances corresponding to the active center histidine Ndelta1H and Nepsilon2H generally observed in this region, are absent in these enzymes. However, for both enzymes, as well as with the H652A and H652Q active center variants of oligopeptidase B, there are two resonances observed in this region that could be assigned to two protonated histidines with a noncatalytic function. The results indicate that these two histidines participate in strong hydrogen bonds. The absence of resonances pertinent to the active center histidine resonances suggests the absence of a low barrier hydrogen bond between the Asp and His in these two enzymes in their ground states. Addition of the peptide boronic acid t-butoxycarbonyl-(D)Val-Leu-(L)boroArg to oligopeptidase B resulted in potent, slow binding inhibition of the enzyme and the appearance of a new resonance at 15.8 ppm, whose chemical shift is appropriate for a tetrahedral boronate complex and a low barrier hydrogen bond. The results demonstrate important dissimilarities between the active centers of the prolyl oligopeptidase class of serine proteases and the pancreatic and subtilisin classes both in the ground state and in the transition-state analog complexes.

Unique push-pull mechanism of dealkylation in soman-inhibited cholinesterases

The pH-dependence and solvent isotope effects of dealkylation in diastereomeric adducts of Electric eel (Ee) and fetal bovine serum (FBS) acetylcholinesterase (AChE) inactivated with P(-)C(+) and P(-)C(-) 2-(3,3-dimethylbutyl) methylphosphonofluoridate (soman) were studied at 4.0 +/- 0.2 degrees C. The rate constant versus pH profiles were fit to a bell-shaped curve for all adducts. Best fit parameters are pK1 4.4-4.6 and pK2 6.3-6.5 for Ee AChE and pK1 4.8-5. 0 and pK2 5.8 for FBS AChE. The pKs are consistent with catalytic participation of the Glu199 anion and HisH+440. Maximal rate constants (kmax) are 13-16 x 10(-3) s-1 for Ee AChE and 8 x 10(-3) s-1 for FBS AChE. The solvent isotope effects at the pH maxima are 1.1-1.3, indicating unlikely proton transfer at the enzymic transition states for the dealkylation reaction. Slopes of log rate constant versus pH plots are near 1 at 25.0 +/- 0.2 degrees C between pH 7.0 and 10.0. In stark contrast, the corresponding adducts of trypsin are very stable even at 37.0 +/- 0.2 degrees C. The rate constants for diastereomers of soman-inhibited trypsin at 37.0 +/- 0.2 degrees C are pH independent and approximately 10(4) times smaller than kmax for analogous adducts with AChE. Dealkylation in soman-inhibited AChEs is estimated to occur at >10(10) times faster than a plausible nonenzymic reaction. Up to 40% of the catalytic acceleration can be attributed to an electrostatic push, and an electrostatic pull provides much of the balance. The results of this work together with results of a product analysis by Michel et al. (1969) can be explained by an initial and rate-determining methyl migration from Cbeta to Calpha. This is driven by the high electron density of residues (Glu199 and Trp84) at a crowded active site and may be concerted with C-O bond breaking. The positive charge at the rate-determining transition state is distributed between Cbeta and His440. A tertiary carbocation may have a fleeting existence before it is trapped by water or neighboring electrons which is likely to be promoted by Glu199 as the proton acceptor.

Effects of serpin binding on the target proteinase: global stabilization, localized increased structural flexibility, and conserved hydrogen bonding at the active site

The binding of human alpha1-proteinase inhibitor to rat trypsin was shown by NMR spectroscopy to raise the pKa' of His57 in the active site but not to disrupt the hydrogen bond between His57 and Asp102. Similar NMR results were observed for the Asp189 to serine mutant of rat trypsin, which is much more stable than wild-type trypsin against autoproteolysis as the result of mutation of the residue at the base of the specificity pocket. This mutant was used in further studies aimed at determining the extent of the conformational transition in trypsin that accompanies serpin binding and leads to disruption of the catalytic activity of the proteinase such that the inhibitor complex is trapped at the acyl enzyme intermediate stage. The stability of rat trypsin toward thermal denaturation was found to be lower in the free enzyme than in the complex with alpha1-proteinase inhibitor. This suggests that the complex contains extensive protein-protein interactions that stabilize overall folding. On the other hand, previous investigations have shown that the proteinase in serpin-proteinase complexes becomes more susceptible to limited proteolysis, suggesting that the conformational change that accompanies binding leads to the exposure of susceptible loops in the enzyme. The existence of this type of conformational change upon complex formation has been confirmed here by investigation of the rate of cleavage of disulfide linkages by added dithiothreitol. This study revealed that, despite the increased stability of trypsin in the complex, one or more of its disulfide bridges becomes much more easily reduced. We suggest that the process of complex formation with alpha1-proteinase inhibitor converts trypsin D189S into an inactive, loose structure, which serves as a "conformational trap" of the enzyme that prevents catalytic deacylation. It is also proposed that plastic region(s) of the activation domain of trypsin may play a crucial role in this inhibitor-induced structural rearrangement.

A low-barrier hydrogen bond in subtilisin: 1H and 15N NMR studies with peptidyl trifluoromethyl ketones

The N delta 1 proton of His 64 forms a hydrogen bond with Asp 32, as part of the catalytic triad in serine proteases of the subtilisin family. His 64 in subtilisin has been studied by 1H and 15N NMR spectroscopy in the presence and absence of peptidyl trifluoromethyl ketones (TFMKs) that are transition state analog inhibitors. For subtilisin Carlsberg, the downfield resonance of the imidazolium N delta 1 proton is approximately 18.3 ppm and the D/H fractionation factor is 0.55 +/- 0.04 at pH 5.5 (11 degrees C), and 0.63 +/- 0.04 (5 degrees C) and 0.68 +/- 0.04 at pH 6 (11 degrees C). In the complex between subtilisin Carlsberg and Z-L-leucyl-L-leucyl-L-phenylalanyltrifluoromethyl ketone (Z-LLF-CF3) at pH values between 6.5 and 10.6, His 64 remains positively charged, and the D/H fractionation factor of its N delta 1 proton is 0.85 +/- 0.05. In the complex between a subtilisin variant from Bacillus lentus and Z-LLF-CF3, the proton resonance at 18.8 ppm is correlated with a 15N resonance at 197.6 ppm downfield from liquid NH3 with a 1JNH of 81 Hz. The chemical shifts of subtilisin complexes with peptidyl TFMKs are among the most downfield shifts reported for any protein. At pH 9.5, His 64 is neutral and the D/H fractionation factor increases to 1.2 with a chemical shift of 15.0. His 64 is positively charged in the free enzyme at low pH, the inhibitor hemiketal complex at neutral pH, and the transition state for amide bond hydrolysis. These data thus provide indirect evidence for the presence of a low-barrier hydrogen bond in the catalytic mechanism of subtilisin proteases.

Protonation-state dependence of hydrogen bond strengths and exchange rates in a serine protease catalytic triad: bovine chymotrypsinogen A

Hydrogen-1 nuclear magnetic resonance spectroscopy was used to measure D/H fractionation factors and the temperature dependence of the rate of hydrogen exchange at two sites in the catalytic triad of chymotrypsinogen (hydrogen bond between aspartate-102 and histidine-57 and hydrogen bond between histidine-57 and serine-195) as a function of the protonation states of the constituent residues. Connectivities in one-dimensional spectra used to assign NMR data were collected at three pH values: pH 9, at which His-57 is neutral and Asp-102 is negatively charged; pH 3.5, at which His-57 is positively charged and Asp-102 is negatively charged; and pH 1, at which His-57 is positively charged and Asp-102 is neutral. The signal from H epsilon 2 of histidine-57 was assigned by reference to 1H-1H NOE connectivities at pH 3.5 to the previously assigned signals from the H epsilon 1 and H delta 2 of the same residue. The D/H fractionation factor, phi, for the hydrogen bond between Asp-102 and His-57 changed from phi = 2 at pH 9 to phi = 0.4 at pH 3.5. From studies of model systems, it may be concluded that a change of phi of this magnitude corresponds to a large increase in hydrogen bond strength. A signal from the hydrogen bond between Ser-195 and His-57 was detected only at the lower pH values studied. The D/H fractionation factor for this hydrogen bond was phi = 0.7 at pH 3.5, indicative of a moderately strong interaction. Data obtained at pH 1 indicate that the hydrogen bond between Asp-102 and His-57 is weakened but that the hydrogen bond between His-57 and Ser-195 persists. The results are consistent with the hypothesis that changes in hydrogen bonding strength serve to lower barriers along the reaction coordinate in the catalytic mechanism. Large pH-dependent changes were found in the activation enthalpy (delta H ++) for exchange with protons from the solvent at the hydrogen bond between aspartate-102 and histidine-57: delta H ++ was approximately 10-12 kcal.mol-1 higher at pH 3.5 than at pH 1 or 9.

Distortion of the active site of chymotrypsin complexed with a serpin

There is no complete understanding of how serine protease inhibitors of the serpin family inhibit their target enzymes. Structural and biochemical studies have suggested that serpins utilize a mechanism that is distinct from the standard mechanism of inhibition proposed for most small protein protease inhibitors. Proton nuclear magnetic resonance spectroscopy was used in the present study to demonstrate a fundamental difference in the atomic environment of the catalytic triad of enzyme in complex with serpins when compared to uncomplexed enzyme and enzyme in complex with standard mechanism inhibitors. This work demonstrates that the active site of chymotrypsin is distorted when complexed to a serpin and makes tenable a mechanism of inhibition in which the serpin induces a conformational change in the enzyme that dramatically reduces or completely abrogates the catalytic activity of the protease.

15N and 1H NMR spectroscopy of the catalytic histidine in chloromethyl ketone-inhibited complexes of serine proteases

The hemiketal hydroxyl groups in chloromethyl ketone (cmk) complexes of trypsin and chymotrypsin have been reported to ionize to the oxyanion with pK(a) values 2-4 pK(a) units below expectations for such a functional group on the basis of the behavior of the hemiketal carbon atom in 13C NMR spectra [Finucane, M. D., & Malthouse, J. P. G. (1992) Biochem. J. 286, 889-900]. The low pK(a) indicates the enzymes selectively stabilize the oxyanion form of the bound inhibitor, and therefore that cmk complexes may be good models of enzyme-mediated transition-state stabilization. However, the 13C NMR studies could not rule out His57 as the titrating group. Here we report the behavior of the ring 15N atoms of His57 in the Ala-Ala-Pro-Val-cmk complex of alpha-lytic protease. Both N(delta 1) and N(epsilon 2) of His57 respond to an ionization with a pK(a) of approximately 7.5, but His57 itself does not titrate as N(epsilon 2) remains alkylated and N(delta 1) remains bonded to a proton over the entire pH range. The species titrating with a pK(a) of approximately 7.5 must therefore be the hemiketal hydroxyl. The results also show that the 1H NMR signal from the proton in the Asp-His hydrogen bond behaves in a characteristic manner in cmk complexes and can be used diagnostically to confirm that His57 does not titrate and to measure the pK(a) of the hemiketal hydroxyl in cmk-protease complexes without resorting to 15N-labeling. We have used the behavior of this signal to directly confirm that His57 does not titrate in the trypsin and chymotrypsin complexes that were the subjects of the original 13C NMR studies.

NMR studies of the alpha-chymotrypsin-(R)-1-acetamido-2-(4-fluorophenyl)ethane-1-boronic acid complex at pH 7

The interaction of (R)-1-acetamido-2-(4-fluorophenyl)ethane-1-boronic acid with alpha-chymotrypsin at pH 7 was studied by a variety of fluorine and proton NMR experiments and the results compared to observations made at pH 4. It was demonstrated that this compound forms a complex with a 1:1 stoichiometry at pH 7; proton NMR indicates that the boronic acid likely is coordinated to the serine-195 residue at the active site. Analysis of fluorine T1 relaxation behavior and 19F(1H) NOE data shows that the rate constant for dissociation of the complex is 4.7 s-1, somewhat faster than is observed at pH 4. The data analysis and the results of two-dimensional 19F(1H) NOE experiments show that interactions between the fluoroaromatic ring of the inhibitor and the enzyme are weaker at the higher pH value, although the motion of the fluoroaromatic ring within the complex appears to be just as restricted as is the case at pH 4.

A low-barrier hydrogen bond in the catalytic triad of serine proteases

Spectroscopic properties of chymotrypsin and model compounds indicate that a low-barrier hydrogen bond participates in the mechanism of serine protease action. A low-barrier hydrogen bond between N delta 1 of His57 and the beta-carboxyl group of Asp102 in chymotrypsin can facilitate the formation of the tetrahedral adduct, and the nuclear magnetic resonance properties of this proton indicate that it is a low-barrier hydrogen bond. These conclusions are supported by the chemical shift of this proton, the deuterium isotope effect on the chemical shift, and the properties of hydrogen-bonded model compounds in organic solvents, including the hydrogen bond in cis-urocanic acid, in which the imidazole ring is internally hydrogen-bonded to the carboxyl group.

1H NMR spectroscopic studies of selenosubtilisin

Anomalously low-field signals in 1H NMR spectra of serine proteases provide valuable information on the protonation state of the catalytic histidine residue. We have examined the pH dependence of the deshielded protons of three different oxidation states of selenosubtilisin, a semisynthetic selenoenzyme with significant peroxidase activity, in order to evaluate the influence of the selenium prosthetic group on the hydrogen-bonding network in the modified active site. In the spectra of the anionic seleninate and selenolate derivatives, two resonances were observed at 18.0 and 15.5/14.0 ppm, assigned respectively to the N delta 1 and N epsilon 2 protons of protonated His64. These signals were apparent from pH 4 to above pH 10, indicating that the negatively charged prosthetic group increases the stability of the imidazolium dramatically, raising its pKa by at least 3-4 pH units. In contrast, a neutral selenenyl sulfide species exhibits no deshielded proton signals at 18 ppm at any pH but has a weak signal at 14.1 ppm above pH 7 which was assigned to the N delta 1 imidazole proton of neutral His64. While the pKa of His64 appears normal (approximately 7) in this derivative, the selenenyl sulfide substitution may alter the orientation of the imidazole ring within the active site for steric reasons. Together with data on the influence of pH on peroxidase activity, these results suggest that selenosubtilisin's His64 acts as a general acid facilitating the reduction of the selenenyl sulfide to selenolate by thiols.

Identification of serine and histidine adducts in complexes of trypsin and trypsinogen with peptide and nonpeptide boronic acid inhibitors by 1H NMR spectroscopy

We have previously shown, in 15N NMR studies of the enzyme's active site histidine residue, that boronic acid inhibitors can form two distinct types of complexes with alpha-lytic protease. Inhibitors that are structural analogs of good alpha-lytic protease substrates form transition-state-like tetrahedral complexes with the active site serine whereas those that are not form complexes in which N epsilon 2 of the active site histidine is covalently bonded to the boron of the inhibitor. This study also demonstrated that the serine and histidine adduct complexes exhibit quite distinctive and characteristic low-field 1H NMR spectra [Bachovchin, W. W., Wong, W. Y. L., Farr-Jones, S., Shenvi, A. B., & Kettner, C. A. (1988) Biochemistry 27, 7689-7697]. Here we have used low-field 1H NMR diagnostically for a series of boronic acid inhibitor complexes of trypsin and trypsinogen. The results show that H-D-Val-Leu-boroArg and Ac-Gly-boroArg, analogs of good trypsin substrates, form transition-state-like serine adducts with trypsin, whereas the nonsubstrate analog inhibitors boric acid, methane boronic acid, butane boronic acid, and triethanolamine borate all form histidine adducts, thereby paralleling the previous results obtained with alpha-lytic protease. However, with trypsinogen, Ac-Gly-boroArg forms predominantly a histidine adduct while H-D-Val-Leu-boroArg forms both histidine and serine adducts, with the histidine adduct predominating below pH 8.0 and the serine adduct predominating above pH 8.0.(ABSTRACT TRUNCATED AT 250 WORDS)

A study of the stabilization of tetrahedral adducts by trypsin and delta-chymotrypsin

delta-Chymotrypsin has been alkylated by 1-13C- and 2-13C-enriched tosylphenylalanylchloromethane. In the intact inhibitor derivative, signals due to the 1-13C- and 2-13C-enriched carbon atoms have chemical shifts which titrate from 55.10 to 59.50 p.p.m. and from 99.10 to 103.66 p.p.m. respectively with similar pKa values of 8.99 and 8.85 respectively. These signals are assigned to a tetrahedral adduct formed between the hydroxy group of serine-195 and the inhibitor. An additional signal at 58.09 p.p.m. and at 204.85 p.p.m. in the 1-13C- and 2-13C-enzyme-inhibitor derivatives respectively does not titrate when the pH is changed and it is assigned to alkylated methionine-192. On denaturation/autolysis of the 1-13C-enriched enzyme-inhibitor derivative these signals associated with the intact inhibitor derivative are no longer detected, and a new signal, which titrates from 56.28 to 54.84 p.p.m. with a pKa of 5.26, is detected. The titration shift of this signal is assigned to the deprotonation of the imidazolium cation of alkylated histidine-57 in the denatured/autolysed enzyme-inhibitor derivative. Model compounds which form stable hydrates and hemiketals in aqueous solutions have been synthesized. By comparing the 13C titration shifts of these model compounds with those of the 13C enriched trypsin- and delta-chymotrypsin-inhibitor derivatives, we deduce that, in both of the intact enzyme-inhibitor derivatives, the zwitterionic tetrahedral adduct containing the imidazolium cation of histidine-57 and the hemiketal oxyanion predominates at alkaline pH values. It is estimated that in both the trypsin and delta-chymotrypsin-inhibitor derivatives the concentration of this zwitterionic tetrahedral adduct is 10,000-fold greater than it would be in water. We conclude that the pKa of the oxyanion of the hemiketal in the presence of the imidazolium cation of histidine-57 is 7.9 and 8.9 in the trypsin and delta-chymotrypsin-inhibitor derivatives respectively and that the pKa of the imidazolium cation of histidine-57 is greater than 7.9 and greater than 8.9 when the oxyanion is present as its conjugate acid, whereas, when the oxyanion is present, the pKa of the imidazolium cation is greater than 11 in both enzyme-inhibitor derivatives. We discuss how these enzymes preferentially stabilize zwitterionic tetrahedral adducts in the intact enzyme-inhibitor derivatives and how they could stabilize similar tetrahedral intermediates during catalysis. It is suggested that substrate binding could raise the pKa of the imidazolium cation of histidine-57 before tetrahedral-intermediate formation which would explain the enhanced nucleophilicity of the hydroxy group of serine-195.(ABSTRACT TRUNCATED AT 400 WORDS)