Characterization of the electrostatic perturbation of a catalytic site (Cys)-S-/(His)-Im+H ion-pair in one type of serine proteinase architecture by kinetic and computational studies on chemically mutated subtilisin variants

Identification of interactions involved in the generation of nucleophilic reactivity and of catalytic competence in the catalytic site Cys/His ion pair of papain

Understanding the roles of noncovalent interactions within the enzyme molecule and between enzyme and substrate or inhibitor is an essential goal of the investigation of active center chemistry and catalytic mechanism. Studies on members of the papain family of cysteine proteinases, particularly papain (EC 3.4.22.2) itself, continue to contribute to this goal. The historic role of the catalytic site Cys/His ion pair now needs to be understood within the context of multiple dynamic phenomena. Movement of Trp177 may be necessary to expose His159 to solvent with consequent decrease in its degree of electrostatic solvation of (Cys25)-S(-). Here we report an investigation of this possibility using computer modeling of quasi-transition states and pH-dependent kinetics using 3,3'-dipyridazinyl disulfide, its n-propyl and phenyl derivatives, and 4,4'-dipyrimidyl disulfide as reactivity probes that differ in the location of potential hydrogen-bonding acceptor atoms. Those interactions that influence ion pair geometry and thereby catalytic competence, including by transmission of the modulatory effect of a remote ionization with pK(a) 4, were identified. A key result is the correlation between the kinetic influence of the modulatory trigger of pK(a) 4 and disruption of the hydrogen bond donated by the indole N-H of Trp177, the hydrophobic shield of the initial "intimate" ion pair. This hydrogen bond is accepted by the amide O of Gln19-a component of the oxyanion hole that binds the tetrahedral species formed from the substrate during the catalytic act. The disruption would be expected to contribute to the mobility of Trp177 and possibly to the effectiveness of the binding of the developing oxyanion.

Combination of a Dinuclear Zn2+ Complex and a Medium Effect Exerts a 1012-Fold Rate Enhancement of Cleavage of an RNA and DNA Model System

The catalytic ability of a dinuclear Zn2+ complex of 1,3-bis-N1-(1,5,9-triazacyclododecyl)propane (3) in promoting the cleavage of an RNA model, 2-hydroxypropyl-p-nitrophenyl phosphate (HPNPP, 1), and a DNA model, methyl p-nitrophenyl phosphate (MNPP, 4), was studied in methanol solution in the presence of added CH3O- at 25 degrees C. The di-Zn2+ complex (Zn2 :3), in the presence of 1 equiv of added methoxide, exhibits a second-order rate constant of (2.75 +/- 0.10) x 10(5) M(-1) s(-1) for the reaction with 1 at s(s)pH 9.5, this being 10(8)-fold larger than the k2 value for the CH3O- promoted reaction (kOCH3 = (2.56 +/- 0.16) x 10(-3) M(-1) s(-1)). The complex is also active toward the DNA model 4, exhibiting Michaelis-Menten kinetics with a KM and kmax of 0.37 +/- 0.07 mM and (4.1 +/- 0.3) x 10(-2) s(-1), respectively. Relative to the background reactions at s(s)pH 9.5, Zn2 :3 accelerates cleavage of each phosphate diester by a remarkable factor of 1012-fold. A kinetic scheme common to both substrates is discussed. The study shows that a simple model system comprising a dinuclear Zn2+ complex and a medium effect of the alcohol solvent achieves a catalytic reactivity that approaches enzymatic rates and is well beyond anything seen to date in water for the cleavage of these phosphate diesters.

Electrostatic properties in the catalytic site of papain: A possible regulatory mechanism for the reactivity of the ion pair

We present an analysis of the electrostatic properties in the catalytic site of papain (EC 3.4.22.2), an archetype enzyme of the C1 cysteine proteinase family, and we investigate their possible role in the formation, stabilization and regulation of the Cys25((-))...His159((+)) catalytic ion pair. The electrostatic properties were computed using a reassociation method based in multicentered multipolar expansions obtained from ab initio quantum calculations of overlapping protein fragments. Solvent effects were introduced by coupling the use of multicentered multipolar expansions to two continuum boundary element methods to solve the Poisson and the linearized Poisson-Boltzmann equations. The electrostatic profile found in the proton transfer region of papain showed that this enzyme has a well-defined electrostatic environment to favor the formation and stabilization of the catalytic ion pair. The papain catalytic site electrostatic profile can be considered as an electrostatic fingerprint of the papain family with the following characteristics: (i) the presence of a net electric field highly aligned in the (Cys25)-SG-->(His159)-ND1 direction; (ii) the electrostatic profile has a saddle-point character; (iii) it is basically a local environmental effect. Furthermore, our analysis describes a possible regulatory mechanism (the E(SG-->ND1) attenuation effect) controlling the ion pair reactivity and permits to infer the Asp57 acidic residue as the most probable candidate to act as the electrostatic modulator.

Variation in the pH-dependent pre-steady-state and steady-state kinetic characteristics of cysteine-proteinase mechanism: evidence for electrostatic modulation of catalytic-site function by the neighbouring carboxylate anion

The acylation and deacylation stages of the hydrolysis of N -acetyl-Phe-Gly methyl thionoester catalysed by papain and actinidin were investigated by stopped-flow spectral analysis. Differences in the forms of pH-dependence of the steady-state and pre-steady-state kinetic parameters support the hypothesis that, whereas for papain, in accord with the traditional view, the rate-determining step is the base-catalysed reaction of the acyl-enzyme intermediate with water, for actinidin it is a post-acylation conformational change required to permit release of the alcohol product and its replacement in the catalytic site by the key water molecule. Possible assignments of the kinetically influential p K (a) values, guided by the results of modelling, including electrostatic-potential calculations, and of the mechanistic roles of the ionizing groups, are discussed. It is concluded that Asp(161) is the source of a key electrostatic modulator (p K (a) 5.0+/-0.1) in actinidin, analogous to Asp(158) in papain, whose influence is not detected kinetically; it is always in the 'on' state because of its low p K (a) value (2.8+/-0.06).

Direct proton magnetic resonance determination of the pKa of the active center histidine in thiolsubtilisin

The serine proteases constitute a group of endopeptidases whose members owe their catalytic activity to the presence of a catalytic triad of amino acids consisting of a serine, a histidine and an aspartate. The pK(a) values for this histidine have been determined for several cases in which there is a negative charge installed at the serine to mimic the oxyanionic intermediate and related transition state for the catalytic pathway. Instances from this laboratory include (1) replacement of the serine by a cysteine in subtilisin to create a thiolate; (2) formation of monoisopropylphosphoryl-Ser 195 monoanionic phosphodiesters (in trypsin and chymotrypsin, Ser 221 in subtilisins); and (3) tetrahedral boronates formed with peptide boronic acids. The nuclear magnetic resonance (NMR) signals pertinent to this histidine, or signals indirectly reflecting the state of ionization of this histidine, have been used effectively to monitor changes in the active center ionization state. In every case studied, there is elevation of the pK(a) at the histidine when the negative charge is installed at the serine position. Herein is reported the first NMR measurement of the active center His 63 pK(a) in thiolsubtilisin Carlsberg; it is elevated by 3 units compared with the parent enzyme. Using a numerical solution (finite difference) of the Poisson-Boltzmann equation, a protein dielectric constant of 4 provides a good estimate of the experimentally observed pK(a) elevations. Very significantly, a very low protein dielectric constant (epsilon(p) = 3-5) is required in all of the comparisons, and for all three enzymes used (chymotrypsin, trypsin, and subtilisin). Finally, we discuss why the electrostatic perturbation sensed at His of the active center is more amplified by a negative charge on the Ser side than the same charge on the Asp side. A plausible explanation is that the positive charge on the imidazolium ring of the His is localized, with the N(delta 1) carrying a smaller fraction, the N(epsilon 2) carrying the bulk of the positive charge.

Variation in aspects of cysteine proteinase catalytic mechanism deduced by spectroscopic observation of dithioester intermediates, kinetic analysis and molecular dynamics simulations

The possibility of a slow post-acylation conformational change during catalysis by cysteine proteinases was investigated by using a new chromogenic substrate, N-acetyl-Phe-Gly methyl thionoester, four natural variants (papain, caricain, actinidin and ficin), and stopped-flow spectral analysis to monitor the pre-steady state formation of the dithioacylenzyme intermediates and their steady state hydrolysis. The predicted reversibility of acylation was demonstrated kinetically for actinidin and ficin, but not for papain or caricain. This difference between actinidin and papain was investigated by modelling using QUANTA and CHARMM. The weaker binding of hydrophobic substrates, including the new thionoester, by actinidin than by papain may not be due to the well-known difference in their S2-subsites, whereby that of actinidin in the free enzyme is shorter due to the presence of Met211. Molecular dynamics simulation suggests that during substrate binding the sidechain of Met211 moves to allow full access of a Phe sidechain to the S2-subsite. The highly anionic surface of actinidin may contribute to the specificity difference between papain and actinidin. During subsequent molecular dynamics simulations the P1 product, methanol, diffuses rapidly (over<8 ps) out of papain and caricain but 'lingers' around the active centre of actinidin. Uniquely in actinidin, an Asp142-Lys145 salt bridge allows formation of a cavity which appears to constrain diffusion of the methanol away from the catalytic site. The cavity then undergoes large scale movements (over 4.8 A) in a highly correlated manner, thus controlling the motions of the methanol molecule. The changes in this cavity that release the methanol might be those deduced kinetically.

Biochemical and X-ray crystallographic studies on shikimate kinase: the important structural role of the P-loop lysine

Shikimate kinase, despite low sequence identity, has been shown to be structurally a member of the nucleoside monophosphate (NMP) kinase family, which includes adenylate kinase. In this paper we have explored the roles of residues in the P-loop of shikimate kinase, which forms the binding site for nucleotides and is one of the most conserved structural features in proteins. In common with many members of the P-loop family, shikimate kinase contains a cysteine residue 2 amino acids upstream of the essential lysine residue; the side chains of these residues are shown to form an ion pair. The C13S mutant of shikimate kinase was found to be enzymatically active, whereas the K15M mutant was inactive. However, the latter mutant had both increased thermostability and affinity for ATP when compared to the wild-type enzyme. The structure of the K15M mutant protein has been determined at 1.8 A, and shows that the organization of the P-loop and flanking regions is heavily disturbed. This indicates that, besides its role in catalysis, the P-loop lysine also has an important structural role. The structure of the K15M mutant also reveals that the formation of an additional arginine/aspartate ion pair is the most likely reason for its increased thermostability. From studies of ligand binding it appears that, like adenylate kinase, shikimate kinase binds substrates randomly and in a synergistic fashion, indicating that the two enzymes have similar catalytic mechanisms.

Kinetic and titration methods for determination of active site contents of enzyme and catalytic antibody preparations

Kinetic characterization of enzymes and analogous catalysts such as catalytic antibodies requires knowledge of the molarity of functional sites. Various stoichiometric titration methods are available for the determination of active-site concentrations of some enzymes and these are exemplified in the second part of this article. Most of these are not general in that they require the existence of certain types of either intermediate or active-site residues that are susceptible to specific covalent modification. Thus they are not readily applicable to many enzymes and they are rarely available currently for titration of catalytic antibody active sites. In the first part of the article we discuss a general kinetic method for the investigation of active-site availability in preparations of macromolecular catalysts. The method involves steady-state kinetics to provide Vmax and Km and single-turnover first-order kinetics using excess of catalyst over substrate to provide the analogous parameters k(obs)lim and K(m)app. The active-site contents of preparations that contain only active catalyst (Ea) and inert material (Ei) may be calculated as [Ea](T) = Vmax)/k(obs)lim. This is true even if nonproductive binding to E(a) occurs. For polyclonal catalytic antibody preparations, which may contain binding but noncatalytic material (Eb) in addition to Ea and Ei, the significance of Vmax/k(obs)lim is more complex but provides an upper limit to E(a). This can be refined by consideration of the relative values of Km and the equilibrium dissociation constant of EbS. Analysis of the Ea, Eb, Ei system requires the separate determination of Ei. For catalytic antibodies this may be achieved by analytical affinity chromatography using an immobilized hapten or hapten analog and an ELISA procedure to ensure the clean separation of Ei from the Ea + Eb mixture.

Engineering the pH-optimum of a triglyceride lipase: from predictions based on electrostatic computations to experimental results

The optimisation of enzymes for particular purposes or conditions remains an important target in virtually all protein engineering endeavours. Here, we present a successful strategy for altering the pH-optimum of the triglyceride lipase cutinase from Fusarium solani pisi. The computed electrostatic pH-dependent potentials in the active site environment are correlated with the experimentally observed enzymatic activities. At pH-optimum a distinct negative potential is present in all the lipases and esterases that we studied so far. This has prompted us to propose the "The Electrostatic Catapult Model" as a model for product release after cleavage of the ester bond. The origin of the negative potential is associated with the titration status of specific residues in the vicinity of the active site cleft. In the case of cutinase, the role of Glu44 was systematically investigated by mutations into Ala and Lys. Also, the neighbouring Thr45 was mutated into Proline, with the aim of shifting the spatial location of Glu44. All the charge mutants displayed altered titration behaviour of active site electrostatic potentials. Typically, the substitution of the residue Glu44 pushes the onset of the active site negative potential towards more alkaline conditions. We, therefore, predicted more alkaline pH optima, and this was indeed the experimentally observed. Finally, it was found that the pH-dependent computed Coulombic energy displayed a strong correlation with the observed melting temperatures of native cutinase.

The lysosomal cysteine proteases

A significant number of exciting papain-like cysteine protease structures have been determined by crystallographic methods over the last several years. This trove of data allows for an analysis of the structural features that empower these molecules as they efficiently carry out their specialized tasks. Although the structure of the paradigm for the family, papain, has been known for twenty years, recent efforts have reaped several structures of specialized mammalian enzymes. This review first covers the commonalities of architecture and purpose of the papain-like cysteine proteases. From that broad platform, each of the lysosomal enzymes for which there is an X-ray structure (or structures) is then examined to gain an understanding of what structural features are used to customize specificity and activity. Structure-based design of inhibitors to control pathological cysteine protease activity will also be addressed.

The active site of phosphorylating glyceraldehyde-3-phosphate dehydrogenase is not designed to increase the nucleophilicity of a serine residue

Changing a catalytic cysteine into a serine, and vice versa, generally leads to a dramatic decrease in enzymatic efficiency. Except a study done on thiol subtilisin, no extensive study was carried out for determining whether the decrease in activity is due to a low nucleophilicity of the introduced amino acid. In the present study, Cys149 of glyceraldehyde-3-phosphate dehydrogenase from Bacillus stearothermophilus was converted into a Ser residue. This leads to a drastic reduction of the kcat value. The rate-limiting step occurs before the hydride transfer step. Selective, but slow, inactivation is observed with specific, structurally different, inhibitors of serine protease. The esterolytic activity of serine mutant towards activated esters is also strongly decreased. The rate-limiting step of the esterase reaction also shifts from deacylation in the wild type to acylation in the mutant. Altogether, these results strongly suggest that the low catalytic efficiency of the Ser mutant is due to a poor nucleophilicity of the hydroxyl serine group within the active site of the enzyme. The fact that (1) the apo --> holo transition does not change esterolytic and inactivating efficiencies, and (2) Ser149 Asn176 double mutant exhibits the same chemical reactivity and esterolytic catalytic efficiency compared to the Ser149 single mutant indicates that the serine residue is not subject to His176 general base catalysis. A linear relationship between the catalytic dehydrogenase rate, the kcat/KM for esterolysis, and the concentration of OH- is observed, thus supporting the alcoholate entity as the attacking reactive species. Collectively this study shows that the active site environment of GAPDH is not adapted to increase the nucleophilicity of a serine residue. This is discussed in relation to what is known about Ser and Cys protease active sites.

The enhancing of a cysteine proteinase activity at acidic pH by protein engineering, the role of glutamic 50 in the enzyme mechanism of caricain

Carica papaya produces four cysteine proteinases. Calculations show that the Cys25, His159 essential ion pair is fully ionised at pH 2.99, where activity cannot be detected, but apparently an additional ionisation with a pKa of 4 is essential for activity (an electrostatic switch). Caricain (EC 3.4.22.30) wt and D158E genetic backgrounds were used to study the contribution of E50A to activity. E50 or E135 are candidates for the switch, E50A would be expected to reduce activity. However, activity increased at pH 5.0 in both backgrounds and at the pH optimum in D158E E50A but decreased slightly in the wt background. This challenges the hypothesis of an electrostatic switch.

Conformational change in the activation of lipase: an analysis in terms of low-frequency normal modes

The interfacial activation of Rhizomucor miehei lipase (RmL) involves the motion of an alpha-helical region (residues 82-96) which acts as a "lid" over the active site of the enzyme, undergoing a displacement from a "closed" to an "open" conformation upon binding of substrate. Normal mode analyses performed in both low and high dielectric media reveal that low-frequency vibrational modes contribute significantly to the conformational transition between the closed and open conformations. In these modes, the lid displacement is coupled to local motions of active site loops as well as global breathing motions. Atomic fluctuations of the first hinge of the lid (residues 83-84) are substantially larger in the low dielectric medium than in the high dielectric medium. Our results also suggest that electrostatic interactions of Arg86 play an important role in terms of both the intrinsic stability of the lid and its displacement, through enhancement of hinge mobility in a high dielectric medium. Additional calculations demonstrate that the observed patterns of atomic fluctuations are an intrinsic feature of the protein structure and not dependent on the nature of specific energy minima.

Characterization of the hydrolytic activity of a polyclonal catalytic antibody preparation by pH-dependence and chemical modification studies: evidence for the involvement of Tyr and Arg side chains as hydrogen-bond donors

The hydrolyses of 4-nitrophenyl 4'-(3-aza-2-oxoheptyl)phenyl carbonate and of a new, more soluble, substrate, 4-nitrophenyl 4'-(3-aza-7-hydroxy-2-oxoheptyl)phenyl carbonate, each catalysed by a polyclonal antibody preparation elicited in a sheep by use of an analogous phosphate immunogen, were shown to adhere closely to the Michaelis-Menten equation, in accordance with the growing awareness that polyclonal catalytic antibodies may be much less heterogeneous than had been supposed. The particular value of studies on polyclonal catalytic antibodies is discussed briefly. Both the kcat and kcat/K(m) values were shown to increase with increase in pH across a pKa of approx. 9. Group-selective chemical modification studies established that the side chains of tyrosine and arginine residues are essential for catalytic activity, and provided no evidence for the involvement of side chains of lysine, histidine or cysteine residues. The combination of evidence from the kinetic and chemical modification studies and from studies on the pH-dependence of binding suggests that catalysis involves assistance to the reaction of the substrate with hydroxide ions by hydrogen-bond donation at the reaction centre by tyrosine and arginine side chains. This combination of hydrogen-bond donors appears to be a feature common to a number of other hydrolytic catalytic antibodies. High-pKa acidic side chains may be essential for the effectiveness of catalytic antibodies that utilize hydroxide ions.