Model calculations of the 1 H/2H, l2C/l?C. and ,60/'l80 kinetic and equilibrium isotope effects at the starred positions for the carbonyl addition reaction of HjCC*H*0* with HO*have been made for two closely related sets of force fields, and for both curved and perpendicular trajectories of nucleophilic approach to carbonyl, for a series of possible transition-state structures at temperatures from 273 to 323 K. In all calculations, the geometrical features and force constants were assumed to change from reactant values toward product values in proportion to the Pauling bond order B of the nucleophile -carbonyl bond. The «-deuterium and /3-deuterium effects are nearly linear functions of B and should be good and consistent probes of transition-state structure. The carbonyl-18# effect is normal and increases steadily with B, as expected. The nucleophile-18# effect is normal for B < 0.5 and inverse for later transition states. The carbonyl-13C is about 1.03 for early transition states and falls to an inverse equilibrium effect. No strikingly anomalous temperature dependences were calculated.Carbonyl-addition reactions, exemplified by the process of eq 1, are of interest in chemistry and biochemistry3•4 and the understanding of their dynamics in terms of transition-state structure is the topic of much current research.5 A very appealing and effective approach to learning the structure of such transition states lies in the use of kinetic isotope effects.5 7 We have employed vibrational analysis of model structures8•9 for the transition state of eq 1 to explore relationships of such CHjCHO + HO CH.j CH.CH(1) OH structures to the following isotopic probes: (1) the «-deuterium secondary isotope effect10 (for the hydrogen attached to the electrophilic carbon); this effect is expected to be inverse (ko > k h), arising from restrictions to motion of the « hydrogen

Molecular mechanics and dynamics combined with semiempirical calculations were carried out for purposes of comparison of the active site characteristics of AChE, trypsin, and chymotrypsin as probed by their diastereomeric adducts with 2-(3,3-dimethylbutyl) methylphosphonofluoridate (soman), methylphosphonate monoester anions, and tetravalent carbonyl intermediates of the reactions of the natural substrates in each case. Glu199 is a key residue in the electrostatic catalytic mechanism of AChE, in removal of the leaving group, and possibly by acting as an alternate general base catalyst. "Pushing" of an alkoxy ligand by Glu199 and the numerous small van der Waals interactions promote dealkylation in phosphonate adducts of AChE much more effectively than any other enzyme. A high concentration of negative charge created by the phosphonate ester monoanion and Glu199 adjacent to it fully accounts for the resistance to the attack of even the strongest nucleophile applied for enzyme reactivation. Stabilization of the developing negative charge on the phosphonates in the soman-inhibited PSCS adducts of serine hydrolases is by electrophilic residues in the oxyanion hole (AChE) and the protonated catalytic His. PR diastereomers of soman-inhibited AChE can be accommodated in an orientation in which the oxyanion hole interactions are lost and for which the stabilizing interactions are 17-26 kcal/mol smaller than in the PS diastereomer. The dealkylation reaction is almost equally likely in all diastereomers of soman-inhibited AChE. The stabilizing interaction energies are approximately 4 kcal/mol greater in the PR than in the PS adducts of the soman-inhibited serine proteases. There is 0.60 unit greater partial negative charge on the phosphonyl fragment in the anion of phosphonate monoesters of Ser than at the oxygens of tetravalent carbonyl transients resulting in approximately 12-22 kcal/mol greater stabilization of the former than the latter.

No abstract

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.

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 o...

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.

The structural profile for the interactions between serine proteases and organophosphorus (OP) compounds can be deduced from recent NMR and X-ray crystallographic data. Using the rationale proposed for serine proteases, dynamic data on the inhibition of acetylcholinesterase by OP compounds is also consistent with structural constraints and an impairment of the proton switch mechanism during phosphorylation.

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.