The catalytic mechanism of serine proteases: single proton versus double proton transfer

The two possible mechanisms of proton transfer on the catalytic process of serine proteases (single or double proton transfer) have been analysed. Intermediate neglect of differential overlap calculations have been performed in the absence and in the presence of the substrate molecule and one water molecule localized in the active site. It is shown that, in the absence of the substrate and water, double proton transfer seems to be the most feasible mechanism. However, when these molecules are introduced in the calculation, the role played by them is to facilitate the formation of the zwitterionic structure (single proton transfer) and to destabilize the intermediate structure which leads to double proton transfer. All calculations were made in vacuo.     

Quantum chemical study of the "catalytic triad" of serine proteases

Using the semi-empirical MNDO/H method several systems simulating the reaction of tetrahedral intermediate formation in the active site of serine proteases have been studied. The role played by elements of the "catalytic triad" in increasing the reactivity of serine hydroxyl has been discussed. The formation of a strong hydrogen bond between His and Asp was shown to be important in lowering the activation energy in the reaction of Ser with substrate. The change in position of the proton located between Ser and His and between His and Asp was analysed. The influence of substrate distortion on the energy of intermediate formation has been considered.     

Stereochemical aspects of peptide hydrolysis catalyzed by serine proteases of the chymotrypsin type

The steric course of peptide hydrolysis catalyzed by serine proteases has been studied on the basis of the available, extensive structural data and taking into account the stereoelectronic theory of Deslongchamps (Heterocycles, 7, 1271 (1977)). These studies allowed elucidation of the structure of intermediates, in particular of the tetrahedral intermediate, and of the main structural events taking place during catalysis. They reveal a difficulty inherent in the generally accepted mechanism of peptide hydrolysis: protonation of the leaving nitrogen in the configuration arising from nucleophilic attack of Ser-195 on the carbonyl carbon cannot take place internally from His-57. Two alternative mechanisms are discussed which are compatible with all implications of the stereoelectronic theory. The main features of the more probable mechanism are: (i) a conformational change allowing the imidazole ring of His-57 to occupy two distinct positions; in one position a proton is abstracted from O^γ of Ser-195, and in the other this proton is donated to the leaving nitrogen; (ii) a configurational change (inversion) of the pyramidal leaving nitrogen reorienting the lone-pair orbital developed during nucleophilic attack; in one orientation CO bond breaking, and in the other CN bond breaking, is allowed. This inversion process confers on the nitrogen the property of a switch controlling the breakdown of the tetrahedral intermediate.     

The 1.8 A structure of the complex between chymostatin and Streptomyces griseus protease A. A model for serine protease catalytic tetrahedral intermediates

The naturally occurring serine protease inhibitor, chymostatin, forms a hemiacetal adduct with the catalytic Ser195 residue of Streptomyces griseus protease A. Restrained parameter least-squares refinement of this complex to 1.8 A resolution has produced an R index of 0 X 123 for the 11,755 observed reflections. The refined distance of the carbonyl carbon atom of the aldehyde to O gamma of Ser195 is 1 X 62 A. Both the R and S configurations of the hemiacetal occur in equal populations, with the end result resembling the expected configuration for a covalent tetrahedral product intermediate of a true substrate. This study strengthens the concept that serine proteases stabilize a covalent, tetrahedrally co-ordinated species and elaborates those features of the enzyme responsible for this effect. We propose that a major driving force for the hydrolysis of peptide bonds by serine proteases is the non-planar distortion of the scissile bond by the enzyme, which thereby lowers the activation energy barrier to hydrolysis by eliminating the resonance stabilization energy of the peptide bond.     

Crystal structure at 1.63 A resolution of the native form of porcine beta-trypsin: revealing an acetate ion binding site and functional water network

The active center of a serine protease is the catalytic triad composed of His-57, Ser-195 and Asp-102. The existing crystal structure data on serine proteases have not fully answered a number of fundamental questions relating to the catalytic activity of serine proteases. The new high resolution native porcine beta-trypsin (BPT) structure is aimed at extending the knowledge on the conformation of the active site and the ordered water structure within and around the active site. The crystal structure of BPT has been determined at 1.63 A resolution. An acetate ion bound at the active site of a trypsin molecule by both classical hydrogen bonds and C-HellipsisO hydrogen bonds has been identified for the first time. A large network of water molecules extending from the recognition amino acid Asp-184 to the entry of the active site has been observed in the BPT structure. A detailed comparison with inhibitor complexes and autolysates indicates that the sulfate ion and the acetate ion bind at the same site of the trypsin molecule. The Ser-195 Cbeta-Ogamma-His-57 Nepsilon angle in the catalytic triad of BPT is intermediate between the corresponding values of the complex and native structure due to acetate ion binding. The network of waters from the recognition amino acid to the active site entry is probably the first ever complete picture of functional waters around the active site. Structural comparisons show that the functional waters involved in the binding of small molecule inhibitors and protease inhibitors are distinctly different.     

Quantum chemistry analysis of the mechanism of action of proteolytic enzymes. II. Nucleophilic attack

Dynamic systems, modelling the elementary act of nucleophilic attack on the carbonyl carbon atom of amide (N-methyl-acetamide) and ester (methylacetate) substrates by some compounds, simulating the nucleophilic group of various types proteases active site were calculated and analysed by the CNDO/2 method, namely: 1) methoxyanion and the methanol (serine proteases), 2) water molecule in the presence of formate anion (carboxylic proteases) and 3) mercaptide anion (CH3S-) (thiolic proteases). The formation of productive enzyme-substrate complex was shown not only to orient reactive groups of the enzyme and substrate, but also to activate it by induction of a certain degree of cleavable pyramidalization, as a result of the partial resonance stabilization energy loss.     

Structural analyses on intermediates in serine protease catalysis

Although the subject of many studies, detailed structural information on aspects of the catalytic cycle of serine proteases is lacking. Crystallographic analyses were performed in which an acyl-enzyme complex, formed from elastase and a peptide, was reacted with a series of nucleophilic dipeptides. Multiple analyses led to electron density maps consistent with the formation of a tetrahedral species. In certain cases, apparent peptide bond formation at the active site was observed, and the electron density maps suggested production of a cis-amide rather than a trans-amide. Evidence for a cis-amide configuration was also observed in the noncovalent complex between elastase and an alpha1-antitrypsin-derived tetrapeptide. Although there are caveats on the relevance of the crystallographic data to solution catalysis, the results enable detailed proposals for the pathway of the acylation step to be made. At least in some cases, it is proposed that the alcohol of Ser-195 may preferentially attack the carbonyl of the cis-amide form of the substrate, in a stereoelectronically favored manner, to give a tetrahedral oxyanion intermediate, which undergoes N-inversion and/or C-N bond rotation to enable protonation of the leaving group nitrogen. The mechanistic proposals may have consequences for protease inhibition, in particular for the design of high energy intermediate analogues.     

QM/MM study of the active site of free papain and of the NMA-papain complex.

Hybrid quantum mechanical/molecular mechanical (QM/MM) calculations using restricted and unrestricted Hartree-Fock and B3LYP ab initio (QM) and Amber force field (MM), respectively, have been applied to study the catalytic site of papain in both free and substrate bonded forms. Ab initio geometry optimizations have been performed for the active site of papain and the N-methyl-acetamide (NMA)-papain complex within the molecular mechanical treatment of the protein environment. A covalent tetrahedral intermediate structure could be obtained only when the amide N atom of the substrate molecule was protonated through a proton transfer from the His-159 in the catalytic site. Our results support the previous assumption that a proton transfer from His-159 to the amide N atom of the substrate occurs prior to or concerted with the nucleophilic attack of the Cys-25 sulfur atom to the carbonyl group of the substrate. The electron correlation effect will reduce the proton transfer barrier. Therefore, this proton transfer can be easily observed in the B3LYP/6-31G* calculations. The HF/6-31G* method overestimates the reaction barrier against this proton transfer. The sulfur atom of Cys-25 and the imidazole ring of His-159 are found to be coplanar in the free form of the enzyme. However, the rotation of the imidazole ring of His-159 was observed during the formation of the tetrahedral intermediate. Without the papain environment, the coplanar thiolate-imidazolium ion pair RS-...ImH+ is much less stable than the neutral form of RSH....Im. Within the protein environment, however, the thiolate-imidazolium ion pair becomes more stable than its neutral form by 4.1 and 0.4 kcal/mol in HF/6-31G* and B3LYP/6-31G* calculations, respectively. The barrier of proton transfer from S-H group of Cys-25 to the imidazole ring of His-159 was reduced from 22.0 kcal/mol to 15.2 kcal/mol by the protein environment in HF/6-31G* calculations. This barrier is found to be much smaller (2.5 kcal/mol) in B3LYP/6-31G* calculations.     

Role of catalytic residues in the formation of a tetrahedral adduct in the acylation reaction of bovine β-trypsin: A molecular orbital study

In the acylation reaction of serine proteases the effect of amino acid residues on the geometrical change of the catalytic site from Michaelis to tetrahedral state was studied by using ab initio molecular orbital calculations. Amino acid residues in the catalytic site and the peptide substrate were calculated as a quantum mechanical region, and all the other amino acid residues and the calcium ion were included in the calculation as the electrostatic effects. The effects of Asp102, Asp194, N-terminus and the oxyanion binding site are large. The oxyanion binding site directly stabilizes the tetrahedral substrate. Asp102 stabilizes the enzyme intermediate, interacting with the protonated His57 residue. In order to elucidate the roles of Asp102 and the oxyanion binding site, energy decomposition analyses were done for the intermolecular interactions. The contribution of Asp102 and the oxyanion binding site to the decrease of energy in the geometrical change is due to the electrostatic effect. The energies of the proton shuttle from Ser195 O^γ to the leaving group of the substrate were calculated for amide and ester substrate models.     

Catalysis by dienelactone hydrolase: A variation on the protease mechanism

Dienelactone hydrolase (DLH), an enzyme from the β-ketoadipate pathway, catalyzes the hydrolysis of dienelactone to maleylacetate. Our inhibitor binding studies suggest that its substrate, dienelactone, is held in the active site by hydrophobic interactions around the lactone ring and by the ion pairs between its carboxylate and Arg-81 and Arg-206. Like the cysteine/serine proteases, DLH has a catalytic triad (Cys-123, His-202, Asp-171) and its mechanism probably involves the formation of covalently bound acyl intermediate via a tetrahedral intermediate. Unlike the proteases, DLH seems to protonate the incipient leaving group only after the collapse of the first tetrahedral intermediate, rendering DLH incapable of hydrolyzing amide analogues of its ester substrate. In addition, the triad His probably does not protonate the leaving group (enolate) or deprotonate the water for deacylation; rather, the enolate anion abstracts a proton from water and, in doing so, supplies the hydroxyl for deacylation. © 1993 Wiley-Liss, Inc.     

Mechanism-based inactivation of human leukocyte elastase via an enzyme-induced sulfonamide fragmentation process

We describe herein the design and in vitro biochemical evaluation of a novel class of mechanism-based inhibitors of human leukocyte elastase (HLE) that inactivate the enzyme via an unprecedented enzyme-induced sulfonamide fragmentation cascade. The inhibitors incorporate in their structure an appropriately functionalized saccharin scaffold. Furthermore, the inactivation of the enzyme by these inhibitors was found to be time-dependent and to involve the active site. Biochemical, HPLC, and mass spectrometric studies show that the interaction of these inhibitors with HLE results in the formation of a stable acyl complex and is accompanied by the release of (L) phenylalanine methyl ester. The data are consistent with initial formation of a Michaelis-Menten complex and subsequent formation of a tetrahedral intermediate with the active site serine (Ser(195)). Collapse of the tetrahedral intermediate with tandem fragmentation results in the formation of a highly reactive conjugated sulfonyl imine which can either react with water to form a stable acyl enzyme and/or undergo a Michael addition reaction with an active site nucleophilic residue (His(57)). It is also demonstrated herein that this class of compounds can be used in the design of inhibitors of serine proteases having either a neutral or basic primary substrate specificity. Thus, the results suggest that these inhibitors constitute a potential general class of mechanism-based inhibitors of (chymo)trypsin-like serine proteases.     

Proton inventory of the water-catalyzed hydrolysis of p-nitrotrifluoroacetanilide

A proton inventory study was made of the water-catalyzed hydrolysis of p-nitrotrifluoro-acetanilide at pH 4.0, 70°C. Multiple protons in the transition state were demonstrated; although two or three protons were shown to be involved, they were not distinguishable. Neither imidazolyl cation nor acetic acid catalyzed the water-catalyzed hydrolysis. The water-catalyzed hydrolysis proceeds through acid catalysis by a water molecule of the breakdown of the tetrahedral intermediate between the anilide and another water molecule. Acid catalysis by the first water molecule is probably assisted by proton transfer from the second water molecule.     

Energetic and structural consequences of perturbing Gly-193 in the oxyanion hole of serine proteases

The oxyanion hole of serine proteases is formed by the backbone N atoms of the catalytic Ser-195 and Gly-193 and engages the backbone O atom of the P1 residue of substrate in an important H-bonding interaction. The energetic contribution of this interaction in the ground and transition states is presently unknown. Measurements of the individual rate constants defining the catalytic mechanism of substrate hydrolysis for wild-type thrombin and trypsin and their G193A and G193P mutants reveal that Gly-193 is required for optimal substrate binding and acylation. Crystal structures of the G193A and G193P mutants of thrombin bound to the active site inhibitor H-d-Phe-Pro-Arg-CH2Cl document the extent of perturbation induced by the replacement of Gly-193. The Ala mutant weakens the H-bonding interaction of the N atom of residue 193, whereas the Pro substitution abrogates it altogether with additional small shifts of the protein backbone. From the kinetic and structural data, we estimate that the H-bonding interaction in the oxyanion hole contributes a stabilization of the ground and transition states of > 1.5 kcal/mol but < 3.0 kcal/mol. These results shed light on a basic aspect of the enzyme-substrate interaction in the entire family of trypsin-like serine proteases.     

Is there a weak H‐bond → LBHB transition on tetrahedral complex formation in serine proteases?

The transformation of a weak hydrogen bond in the free enzyme into a low-barrier hydrogen bond (LBHB) in the tetrahedral intermediate has been suggested as an important factor facilitating catalysis in serine proteases. In this work, we examine the structure of the H-bond in the Asp102-His57 diad of serine proteases in the free enzyme and in a covalent tetrahedral complex (TC) with a trifluoromethylketone inhibitor. We apply ab initio quantum mechanical calculations to models consisting of a large molecular fragment of the enzyme active site, and the combined effect of the rest of the protein body and the solvation by surrounding bulk water was simulated by a self-consistent reaction field method in our novel QM/SCRF(VS) approach. Potential profiles of adiabatic proton transfer in the Asp102-His57 diad in these model systems were calculated. We conclude that the hydrogen bond in both the free enzyme and in the enzyme-inhibitor TC is a strong ionic asymmetric one-well hydrogen bond, in contrast to a previous suggestion that it is a weak H-bond in the former and a double-well LBHB in the latter.     

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.     

Crystal structures of two engineered thiol trypsins

We have determined the three-dimensional structures of engineered rat trypsins which mimic the active sites of two classes of cysteine proteases. The catalytic serine was replaced with cysteine (S195C) to test the ability of sulfur to function as a nucleophile in a serine protease environment. This variant mimics the cysteine trypsin class of thiol proteases. An additional mutation of the active site aspartate to an asparagine (D102N) created the catalytic triad of the papain-type cysteine proteases. Rat trypsins S195C and D102N,S195C were solved to 2.5 and 2.0 A, respectively. The refined structures were analyzed to determine the structural basis for the 10(6)-fold loss of activity of trypsin S195C and the 10(8)-fold loss of activity of trypsin D102N,S195C, relative to rat trypsin. The active site thiols were found in a reduced state in contrast to the oxidized thiols found in previous thiol protease structures. These are the first reported structures of serine proteases with the catalytic centers of sulfhydryl proteases. Structure analysis revealed only subtle global changes in enzyme conformation. The substrate binding pocket is unaltered, and active site amino acid 102 forms hydrogen bonds to H57 and S214 as well as to the backbone amides of A56 and H57. In trypsin S195C, D102 is a hydrogen-bond acceptor for H57 which allows the other imidazole nitrogen to function as a base during catalysis. In trypsin D102N,S195C, the asparagine at position 102 is a hydrogen-bond donor to H57 which places a proton on the imidazole nitrogen proximal to the nucleophile. This tautomer of H57 is unable to act as a base in catalysis.(ABSTRACT TRUNCATED AT 250 WORDS)     

Mechanistic characterization of N-formimino-L-glutamate iminohydrolase from Pseudomonas aeruginosa

N-Formimino-l-glutamate iminohydrolase (HutF) from Pseudomonas aeruginosa catalyzes the deimination of N-formimino-l-glutamate in the histidine degradation pathway. An amino acid sequence alignment between HutF and members of the amidohydrolase superfamily containing mononuclear metal centers indicated that residues Glu-235, His-269, and Asp-320 are involved in substrate binding and activation of the nucleophilic water molecule. The purified enzyme contained up to one equivalent of zinc. The metal was removed by dialysis against the metal chelator dipicolinate with the complete loss of catalytic activity. Enzymatic activity was restored by incubation of the apoprotein with Zn2+, Cd2+, Ni2+, or Cu2+. The mutation of Glu-235, His-269, or Asp-320 resulted in the diminution of catalytic activity by two to six orders of magnitude. Bell-shaped profiles were observed for kcat and kcat/Km as a function of pH. The pKa of the group that must be unprotonated for catalytic activity was consistent with the ionization of His-269. This residue is proposed to function as a general base in the abstraction of a proton from the metal-bound water molecule. In the proposed catalytic mechanism, the reaction is initiated by the abstraction of a proton from the metal-bound water molecule by the side chain imidazole of His-269 to generate a tetrahedral intermediate of the substrate. The collapse of the tetrahedral intermediate commences with the abstraction of a second proton via the side chain carboxylate of Asp-320. The C-N bond of the substrate is subsequently cleaved with proton transfer from His-269 to form ammonia and the N-formyl product. The postulated role of the invariant Glu-235 is to ion pair with the positively charged formimino group of the substrate.     

Dynamics of protein conformational fluctuation in enzyme catalysis with special attention to proton transfers in serine proteinases

We have provided a quantum mechanical model for proteinase-catalyzed peptide, amide and ester hydrolysis. The model rests on electron and atom transfer theory, but incorporates the dynamics of conformational nuclear modes as a new element. The model is applied to acylation, but can straightaway be extended to deacylation, and is substantiated by recent structural and kinetic data for proteinase enzyme catalysis. The role of the conformational modes is found to be two-fold. First, the crystallographic distances for the proton transfers involved are far too large for direct transfer. His-57 mobility, handled stochastically, to bring the donor and acceptor groups within suitable reach, is therefore a crucial element of the theory. Secondly, the charge alignment in the Asp-102/His-57/tetrahedral intermediate system implies that the curvature of the potential surface along the conformational coordinates in this state is much lower than in the initial enzyme-substrate and final acyl states. A consequence of this is that the activation energy liberated after the first proton transfer is not dissipated, but stored in the conformational system and used in the second proton transfer step.     

Nitrogen-15 NMR spectroscopy of the catalytic-triad histidine of a serine protease in peptide boronic acid inhibitor complexes

15N NMR spectroscopy was used to examine the active-site histidyl residue of alpha-lytic protease in peptide boronic acid inhibitor complexes. Two distinct types of complexes were observed: (1) Boronic acids that are analogues of substrates form complexes in which the active-site imidazole ring is protonated and both imidazole N-H protons are strongly hydrogen bonded. With the better inhibitors of the class this arrangement is stable over the pH range 4.0-10.5. The results are consistent with a putative tetrahedral intermediate like complex involving a negatively charged, tetrahedral boron atom covalently bonded to O gamma of the active-site serine. (2) Boronic acids that are not substrate analogues form complexes in which N epsilon 2 of the active-site histidine is covalently bonded to the boron atom of the inhibitor. The proton bound to N delta 1 of the histidine in these histidine-boronate adducts remains strongly hydrogen bonded, presumably to the active-site aspartate. Benzeneboronic acid, which falls in this category, forms an adduct with histidine. In both types of complexes the N-H protons of His-57 exchange unusually slowly as evidenced by the room temperature visibility of the low-field 1H resonances and the 15N-H spin couplings. These results, coupled with the kinetic data of the preceding paper [Kettner, C. A., Bone, R., Agard, D. A., & Bachovchin, W. W. (1988) Biochemistry (preceding paper in this issue)], indicate that occupancy of the specificity subsites may be required to fully form the transition-state binding site. The significance of these findings for understanding inhibitor binding and the catalytic mechanism of serine proteases is discussed.     

Model of serine proteases charge relay system - pcilo study

The PCILO (Perturbative Configuration Interaction Using Localized Orbitals) method has been used to determine the electronic structure of the active center of serine proteases. The results show that the carboxyl group of the aspartic acid residue is the ultimate proton acceptor of the catalytic triad (Asp, His, Ser)^?. In the absence of a substrate the negative charge of the active centre is delocalized, causing polarization of the Ser O^γ-H bond and an increase of the nucleophilicity of the O^γ atom. The proton of the O^γ-H bond of the Ser residue is, however, only partially transferred to the N^?2 atom of imidazole His. The hydration of the model charge relay system is also investigated.