Proton Bridging in Catalysis by and Inhibition of Serine Proteases of the Blood Cascade System.

Inquiries into the participation of short hydrogen bonds in stabilizing transition states and intermediate states in the thrombin, factor Xa, plasmin and activated protein C-catalyzed reactions revealed that specific binding of effectors at Sn, n = 1-4 and S'n, n = 1-3 and at remote exosites elicit complex patterns of hydrogen bonding and involve water networks. The methods employed that yielded these discoveries include; (1) kinetics, especially partial or full kinetic deuterium solvent isotope effects with short cognate substrates and also with the natural substrates, (2) kinetic and structural probes, particularly low-field high-resolution nuclear magnetic resonance (1H NMR), of mechanism-based inhibitors and substrate-mimic peptide inhibitors. Short hydrogen bonds form at the transition states of the catalytic reactions at the active site of the enzymes as they do with mechanism-based covalent inhibitors of thrombin. The emergence of short hydrogen bonds at the binding interface of effectors and thrombin at remote exosites has recently gained recognition. Herein, I describe our contribution, a confirmation of this discovery, by low-field 1H NMR. The principal conclusion of this review is that proton sharing at distances below the sum of van der Waals radii of the hydrogen and both donor and acceptor atoms contribute to the remarkable catalytic prowess of serine proteases of the blood clotting system and other enzymes that employ acid-base catalysis. Proton bridges also play a role in tight binding in proteins and at exosites, i.e., allosteric sites, of enzymes.

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.

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.

Locating the rate-determining step(s) for three-step hydrolase-catalyzed reactions with DYNAFIT.

Hydrolytic reactions of oligopeptide 4-nitroanilides catalyzed by human-alpha-thrombin, human activated protein C and human factor Xa were studied at pH 8.0-8.4 and 25.0+/-0.1 degrees C by the progress curve method and individual rate constants were calculated mostly within 10% internal error using DYNAFITV. A systematic strategy has been developed for fitting a three-step consecutive mechanism to eighteen hundred to six thousand time-course data points polled from two to four independent kinetic experiments. Enzyme and substrate concentrations were also calculated. Individual rate constants well reproduce published values obtained under comparable conditions and the Michaelis-Menten kinetic parameters calculated from these elementary rate constants are also within reasonable limits of published values. For comparison, the integrated Michaelis-Menten equation was also fitted to data from twelve sets. Both the k(cat) and k(cat)/K(m) values are within 15% agreement with those calculated using the elementary rate constants obtained with DYNAFITV. Rate constants for the second and third consecutive steps are within 3-4 fold indicating that both determine the overall rate. The Factor Xa-catalyzed hydrolysis of N-alpha-Z-D-Arg-Gly-Arg-pNA.2HCl at pH 8.4 in a series of buffers containing increasing fractions of deuterium at 25.0+/-0.1 degrees C shows a very strong dependence of k(3) and a moderate dependence of k(2) on D content in the buffer: the fractionation factors are: 0.49+/-0.03 for K(1,) 0.70+/-0.05 for k(2), and (0.32+/-0.03)(2) for k(3).

Deuterium solvent isotope effect and proton-inventory studies of factor Xa-catalyzed reactions.

Kinetic solvent isotope effects (KSIEs) for the factor Xa (FXa)-catalyzed activation of prothrombin in the presence and absence of factor Va (FVa) and 5.0 x 10(-5) M phospholipid vesicles are slightly inverse, 0.82-0.93, when substrate concentrations are at 0.2 Km. This is consistent with the rate-determining association of the enzyme-prothrombin assembly, rather than the rate-limiting chemical transformation. FVa is known to effect a major conformational change to expose the first scissile bond in prothrombin, which is the likely event triggering a major solvent rearrangement. At prothrombin concentrations > 5 Km, the KSIE is 1.6 +/- 0.3, when FXa is in a 1:1 ratio with FVa but becomes increasingly inverse, 0.30 +/- 0.05 and 0.19 +/- 0.04, when FXa/FVa is 1:4, with an increasing FXa and substrate concentration. The rate-determining step changes with the conditions, but the chemical step is not limiting under any circumstance. This corroborates the proposed predominance of the meizothrombin pathway when FXa is well-saturated with the prothrombin complex. In contrast, the FXa-catalyzed hydrolysis of N-alpha-Z-D-Arg-Gly-Arg-pNA.2HCl (S-2765) and H-D-Ile-L-Pro-L-Arg-pNA.HCl (S-2288) is most consistent with two-proton bridges forming at the transition state between Ser195 OgammaH and His57 N(epsilon)2 and His57 Ndelta1 and Asp102 COObeta- at the active site, with transition-state fractionation factors of phi1 = phi2 = 0.57 +/- 0.07 and phiS = 0.78 +/- 0.16 for solvent rearrangement for S-2765 and phi1 = phi2 = 0.674 +/- 0.001 for S-2288 under enzyme saturation with the substrate at pH 8.40 and 25.0 +/- 0.1 degrees C. The rate-determining step(s) in these reactions is most likely the cleavage of the C-N bond and departure of the leaving group.

Full and partial deuterium solvent isotope effect studies of alpha-thrombin-catalyzed reactions of natural substrates.

Proton inventory studies of the thrombin-catalyzed fibrinogen activation to fibrinopeptide A are most consistent with a two-proton bridge forming at the transition state probably between Ser195 OgammaH and His57 Nepsilon2 and His57 Ndelta1 and Asp102 COObeta- at the active site, with fractionation factors 0.66 +/- 0.03 under enzyme saturation with substrate and 0.64 +/- 0.03 at fibrinogen concentration at 0.2 Km, at pH 8.0, pD 8.6, and 25.0 +/- 0.1 degrees C. Strongly inverse solvent isotope effects (SIEs) result from inverse lag times and maximal slopes of blood clotting plots, which are also anion and cation dependent. The blood clot is much coarser in D2O, as indicated in clotting curves with 3-9 times shorter lag time and steeper slopes with respect to H2O. The finer the particles, the weaker the H-bonds interlocking the fibrin mesh and/or in water structure around fibrin. Proton inventories of inverse lag times and maximal slopes of blood clotting curves in buffers containing Na+ and Cl- ions give the best fit to an exponential dependence on deuterium content in the buffer and give fractionation factors 5.6 +/- 0.5 and 7.8 +/- 0.6 at pH 8.0 and 25.0 +/- 0.1 degrees C. The thrombin-catalyzed activation of protein C (PC) to APC is associated with inverse kinetic SIEs (KSIEs) of 0.75 +/- 0.09 and 1.02 +/- 0.06 in 0.3 M NaCl and 0.3 M choline chloride, respectively, at substrate concentrations = 0.2 Km. In comparison, thrombin-catalyzed hydrolysis of chromogenic substrates gives greater KSIEs (Enyedy, E. I.; Kovach. I. M J. Am. Chem. Soc. 2004, 126, 6017-6024) and more complex proton inventories than the ones reported here for the first time for natural substrates. The present study illuminates differences in the character of the rate-determining transition state for the initial phase of the two physiological reactions catalyzed by thrombin.

Proton inventory studies of alpha-thrombin-catalyzed reactions of substrates with selected P and P' sites.

Deuterium kinetic solvent isotope effects for the human alpha-thrombin-catalyzed hydrolysis of (1) substrates with selected P(1)-P(3) sites, Z-Pro-Arg-7-amido-4-methylcoumarin (7-AMC), N-t-Boc-Val-Pro-Arg-7-AMC, Bz-Phe-Val-Arg-4-nitroanilide (pNA), and H-D-Phe-L-Pip-Arg-pNA, are (DOD)k(cat) = (2.8-3.3) +/- 0.1 and (DOD)(k(cat)/K(m)) = (0.8-2.1) +/- 0.1 and (2) internally fluorescence-quenched substrates (a) (AB)Val-Phe-Pro-Arg-Ser-Phe-Arg-Leu-Lys(DNP)-Asp-OH, an optimal sequence, and (b) (AB)Val-Ser-Pro-Arg-Ser-Phe-Gln-Lys(DNP)-Asp-OH, recognition sequence for factor VIII, are (DOD)k(cat) = 2.2 +/- 0.2 and (DOD)(k(cat)/K(m)) = (0.8-0.9) +/- 0.1, at the pL (L = H, D) maximum, 8.4-9.0, and (25.0-26.0) +/- 0.1 degrees C. The most plausible models fitting the partial isotope effect (proton inventory) data have been selected on the basis of lowest values of the reduced chi squared and consistency of fractionation factors at all substrate concentrations, assuming rate-determining acylation. The data for Z-Pro-Arg-7-AMC are consistent with a single-proton bridge at the transition state phi(TS) = 0.39 +/- 0.05 and components for solvent reorganization phi(S) = 0.8 +/- 0.1 and phi(S) = 1.22 for k(cat) and k(cat)/K(m), respectively. The data for tripeptide amides fit bowl-shaped curves; an example is N-t-Boc-Val-Pro-Arg-7-AMC: phi(TS)(1) = phi(TS)(2) = 0.57 +/- 0.01 and phi(S) = 1 for k(cat) and 1.6 +/- 0.1 for k(cat)/K(m). Proton inventories for the nonapeptide (2b) are linear. The data for k(cat) for H-D-Phe-L-Pip-Arg-pNA and the decapeptide (2a) are most consistent with two identical fractionation factors for catalytic proton bridging, phi(TS)(1) = phi(TS)(2) = 0.68 +/- 0.02 and a large inverse component (phi(S) = 3.1 +/- 0.5) for the latter, indicative of substantial solvent reorganization upon leaving group departure. Proton inventory curves for k(cat)/K(m) for nearly all substrates are dome-shaped with an inverse isotope effect component (phi(S) = 1.2-2.4) originating from solvent reorganization during association of thrombin with substrate. These large contributions from medium effects are in full accord with the conformational adjustments required for the fulfillment of the dual, hemostatic and thrombolytic, functions of thrombin.