Evaluation of hydrogen-bonding and enantiomeric P2-S2 hydrophobic contacts in dynamic aspects of molecular recognition by papain

Cysteine proteases: mode of action and role in epidermal differentiation

Desquamation or cell shedding in mammalian skin is known to involve serine proteases, aspartic proteases and glycosidases. In addition, evidence continues to accumulate that papain-like cysteine proteases and an inhibitor cystatin M/E largely confined to the cutaneous epithelia also play key roles in the process. This involves the complete proteolysis of cell adhesive structures of the stratum corneum, the corneodesmosomes and notably of the desmogleins. Continual cell replacement in the epidermis is the result of the balance between the loss of the outer squames and mitosis of the cells in the basal cell layer. This article provides a brief account of the salient features of the characteristics and catalytic mechanism of cysteine proteases, followed by a discussion of the relevant epidermal biology. The proteases include the asparaginyl endopeptidase legumain, which exerts a strict specificity for the hydrolysis of asparaginyl bonds, cathepsin-V and cathepsin-L. The control of these enzymes by cystatin M/E regulates the processing of transglutaminases and is crucial in the biochemical pathway responsible for regulating the cross-linking and desquamation of the stratum corneum. In addition, caspase-14 has now been shown to play a major part in epidermal maturation. Uncontrolled proteolytic activity leads to abnormal hair follicle formation and deleterious effects on the skin barrier function.

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.

Temperature-dependences of the kinetics of reactions of papain and actinidin with a series of reactivity probes differing in key molecular recognition features

The temperature-dependences of the second-order rate constants (k) of the reactions of the catalytic site thiol groups of two cysteine peptidases papain (EC 3.4.22.2) and actinidin (EC 3.4.22.14) with a series of seven 2-pyridyl disulphide reactivity probes (R-S-S-2-Py, in which R provides variation in recognition features) were determined at pH 6.7 at temperatures in the range 4-30 degrees C by stopped-flow methodology and were used to calculate values of DeltaS++, DeltaH++ and DeltaG++. The marked changes in DeltaS++ from negative to positive in the papain reactions consequent on provision of increase in the opportunities for key non-covalent recognition interactions may implicate microsite desolvation in binding site-catalytic site signalling to provide a catalytically relevant transition state. The substantially different behaviour of actinidin including apparent masking of changes in DeltaH++ by an endothermic conformational change suggests a difference in mechanism involving kinetically significant conformational change.

Energetics and Cooperativity of the Hydrogen Bonding and Anchor Interactions that Bind Peptides to MHC Class II Protein

The complexity of the interaction between major histocompatibility complex class II (MHC II) proteins and peptide ligands has been revealed through structural studies and crystallographic characterization. Peptides bind through side-chain "anchor" interactions with MHC II pockets and an extensive array of genetically conserved hydrogen bonds to the peptide backbone. Here we quantitatively investigate the kinetic hierarchy of these interactions. We present results detailing the impact of single side-chain mutations of peptide anchor residues on dissociation rates, utilizing two I-A(d)-restricted peptides, one of which has a known crystal structure, and 24 natural and non-natural amino acid mutant variants of these peptides. We find that the N-terminal P1, P4 and P6 anchor-pocket interactions can make significant contributions to binding stability. We also investigate the interactions of these peptides with four I-A(d) MHC II proteins, each mutated to disrupt conserved hydrogen bonds to the peptide backbone. These complexes exhibit kinetic behavior suggesting that binding energy is disproportionately invested near the peptide N terminus for backbone hydrogen bonds. We then evaluate the effects of simultaneously modifying both anchor and hydrogen bonding interactions. A quantitative analysis of 71 double mutant cycles reveals that there is little apparent cooperativity between anchor residue interactions and hydrogen bonds, even when they are directly adjacent (<5A).

Probing the Cruzain S2 Recognition Subsite: A Kinetic and Binding Energy Calculation Study

Cysteine proteases are relevant to several aspects of the parasite life cycle and the parasite-host relationship. Moreover, they appear as promising targets for antiparasite chemotherapy. Here, a quantitative investigation on the catalytic properties of cruzain, the papain-like cysteine protease from epimastigotes of Trypanosoma cruzi, is reported. The results indicate that kinetics for the cruzain catalyzed hydrolysis of N-alpha-benzyloxycarbonyl-l-arginyl-l-alanine-(7-amino-4-methylcoumarin), N-alpha-benzyloxycarbonyl-l-phenylalanyl-l-alanine-(7-amino-4-methylcoumarin), and N-alpha-benzyloxycarbonyl-l-tyrosyl-l-alanine-(7-amino-4-methylcoumarin) can be consistently fitted to the minimum three-step mechanism of cysteine proteases involving the acyl.enzyme intermediate E.P; the deacylation step is rate-limiting in enzyme catalysis. Remarkably, these substrates show identical catalytic parameters. This reflects the ability of the cruzain Glu205 residue, located at the bottom of the S(2) subsite, to neutralize the substrate/inhibitor polar P(2) residues (e.g., Arg or Tyr) and to be solvent-exposed when substrate/inhibitor nonpolar P(2) residues (e.g., Phe) fit the S(2) subsite. More complex catalytic mechanisms are also discussed. Binding free-energy calculation provides a quantitative framework for the interpretation of these results; in particular, direct evidence for the compensatory effect between Coulomb interaction(s) and solvation effect(s) is reported. These results appear of general significance for a deeper understanding of (macro)molecular recognition and for the rational design of novel inhibitors of parasitic cysteine proteases.

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

Binding interactions between peptides and proteins of the class II major histocompatibility complex

The activation of helper T cells by peptides bound to proteins of the class II Major Histocompatibility Complex (MHC II) is pivotal to the initiation of an immune response. The primary functional requirement imposed on MHC II proteins is the ability to efficiently bind thousands of different peptides. Structurally, this is reflected in a unique architecture of binding interactions. The peptide is bound in an extended conformation within a groove on the membrane distal surface of the protein that is lined with several pockets that can accommodate peptide side-chains. Conserved MHC II protein residues also form hydrogen bonds along the length of the peptide main-chain. Here we review recent advances in the study of peptide-MHC II protein reactions that have led to an enhanced understanding of binding energetics. These results demonstrate that peptide-MHC II protein complexes achieve high affinity binding from the array of hydrogen bonds that are energetically segregated from the pocket interactions, which can then add to an intrinsic hydrogen bond-mediated affinity. Thus, MHC II proteins are unlike antibodies, which utilize cooperativity among binding interactions to achieve high affinity and specificity. The significance of these observations is discussed within the context of possible mechanisms for the HLA-DM protein that regulates peptide presentation in vivo and the design of non-peptide molecules that can bind MHC II proteins and act as vaccines or immune modulators.

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.

Subsite specificity of trypanosomal cathepsin L-like cysteine proteases. Probing the S2 pocket with phenylalanine-derived amino acids

The S2 subsite of mammalian cysteine proteinases of the papain family is essential for specificity. Among natural amino acids, all these enzymes prefer bulky hydrophobic residues such as phenylalanine at P2. This holds true for their trypanosomal counterparts: cruzain from Trypanosoma cruzi and congopain from T. congolense. A detailed analysis of the S2 specificity of parasitic proteases was performed to gain information that might be of interest for the design of more selective pseudopeptidyl inhibitors. Nonproteogenic phenylalanyl analogs (Xaa) have been introduced into position P2 of fluorogenic substrates dansyl-Xaa-Arg-Ala-Pro-Trp, and their kinetic constants (Km, kcat/Km) have been determined with congopain and cruzain, and related host cathepsins B and L. Trypanosomal cysteine proteases are poorly stereoselective towards D/L-Phe, the inversion of chirality modifying the efficiency of the reaction but not the Km. Congopain binds cyclohexylalanine better than aromatic Phe derivatives. Another characteristic feature of congopain compared to cruzain and cathepsins B and L was that it could accomodate a phenylglycyl residue (kcat/Km = 1300 mM-1.s-1), while lengthening of the side chain by a methylene group only slightly impaired the specificity constant towards trypanosomal cysteine proteases. Mono- and di-halogenation or nitration of Phe did not affect Km for cathepsin L-like enzymes, but the presence of constrained Phe derivatives prevented a correct fitting into the S2 subsite. A model of congopain has been built to study the fit of Phe analogs within the S2 pocket. Phe analogs adopted a positioning within the S2 pocket similar to that of the Tyr of the cruzain/Z-Tyr-Ala-fluoromethylketone complex. However, cyclohexylalanine has an energetically favorable chair-like conformation and can penetrate deeper into the subsite. Fitting of modeled Phe analogs were in good agreement with kinetic parameters. Furthermore, a linear relationship could be established with logP, supporting the suggestion that fitting into the S2 pocket of trypanosomal cysteine proteases depends on the hydrophobicity of Phe analogs.

A classical enzyme active center motif lacks catalytic competence until modulated electrostatically

The cysteine proteinase superfamily is a source of natural structural variants of value in the investigation of mechanism. It has long been considered axiomatic that catalytic competence of these enzymes mirrors the generation of the ubiquitous catalytic site imidazolium-thiolate ion pair. We here report definitive evidence from kinetic studies supported by electrostatic potential calculations, however, that at least for some of these enzymes the ion pair state which provides the nucleophilic and acid-base chemistry is essentially fully developed at low pH where the enzymes are inactive. Catalytic competence requires an additional protonic dissociation with a common pKa value close to 4 possibly from the Glu50 cluster to control ion pair geometry. The pH dependence of the second-order rate constant (k) for the reactions of the catalytic site thiol groups with 4,4'-dipyrimidyl disulfide is shown to provide the pKa values for the formation and deprotonation of the (Cys)-S-/(His)-Im+H ion pair state. Analogous study of the reactions with 2,2'-dipyridyl disulfide reveals other kinetically influential ionizations, and all of these pKa values are compared with those observed in the pH dependence of kcat/Km for the catalyzed hydrolysis of N-acetylphenylalanylglycine 4-nitroanilide. The discrepancy between the pKa value for ion pair formation and the common pKa value close to 4 related to generation of catalytic activity is particularly marked for ficin (pKa 2.49 +/- 0.02) and caricain (pKa 2.88 +/- 0.02) but exists also for papain (pKa 3.32 +/- 0.01).

Structure of chymopapain M the late-eluted chymopapain deduced by comparative modelling techniques and active-centre characteristics determined by pH-dependent kinetics of catalysis and reactions with time-dependent inhibitors: the Cys-25/His-159 ion-pair is insufficient for catalytic competence in both chymopapain M and papain

Chymopapain M, the monothiol cysteine proteinase component of the chymopapain band eluted after chymopapains A and B in cation-exchange chromatography, was isolated from the dried latex of Carica papaya and characterized by kinetic and chromatographic analysis. This late-eluted chymopapain is probably a component of the cysteine proteinase fraction of papaya latex discovered by Schack [(1967) Compt. Rend. Trav. Lab. Carlsberg 36, 67-83], named papaya peptidase B by Lynn [(1979) Biochim. Biophys. Acta 569, 193-201] and partially characterized by Polgár [(1981) Biochim. Biophys. Acta 658, 262-269] and is the enzyme with unusual specificity characteristics (papaya proteinase IV) that Buttle, Kembhavi, Sharp, Shute, Rich and Barrett [Biochem. J. (1989) 261, 469-476] claimed to be a previously undetected cysteine proteinase eluted from a cation-exchange column near to the early-eluted chymopapains. A study of the time-dependent chromatographic consequences of thiol-dependent proteolysis of the components of papaya latex is reported. Chymopapain M was isolated by (i) affinity chromatography followed by separation from papain using cation-exchange f.p.l.c. on a Mono S HR5/5 column and (ii) cation-exchange chromatography followed by an unusual variant of covalent chromatography by thiol-disulphide interchange. The existence in chymopapain M of a nucleophilic interactive Cys/His catalytic-site system analogous to those in papain (EC 3.4.22.2) and other cysteine proteinases was deduced from the characteristics shape of the pH-second-order rate constant (k) profiles for its reactions with 2,2'-dipyridyl disulphide and ethyl 2-pyridyl disulphide. Analysis of the pH-k data for the reactions of chymopapain M with the 2-pyridyl disulphides and with 4,4'-dipyridyl disulphide permits the assignment of molecular pKa values of 3.4 and 8.7 to the formation and subsequent dehydronation of the Cys-S-/His-Im+H state of the catalytic site and reveals three other kinetically influential ionizations with pKa values 3.4, 4.3 and 5.6. The pH-dependences of kcat./Km for the hydrolysis of N-acetyl-L-Phe-Gly-4-nitroanilide at 25.0 degrees C and I0.1 M catalysed by chymopapain M and papain were determined. For both enzymes, little catalytic activity (5-7% of the maximal) develops consequent on formation of the catalytic site Cys-S-/His-Im+H ion-pair state (across pKa 3.4 for both enzymes). For papain, full expression of Kcat./Km for the uncharged substrate requires only the additional hydronic dissociation with pKa 4.2. By contrast, full expression of kcat./Km for chymopapain M requires additional hydronic dissociation with pKa values of 4.3 and 5.6.(ABSTRACT TRUNCATED AT 400 WORDS)

Clarification of the pH-dependent kinetic behaviour of papain by using reactivity probes and analysis of alkylation and catalysed acylation reactions in terms of multihydronic state models: implications for electrostatics calculations and interpretation of the consequences of site-specific mutations such as Asp-158-Asn and Asp-158-Glu

1. The complex behaviour of papain (EC 3.4.22.2) in acidic media has been investigated by (a) stopped-flow reactivity probe kinetics using 4,4'-dipyrimidyl disulphide (I) and 2,2'-dipyridyl disulphide (II) as thiol-specific time-dependent inhibitors with markedly different susceptibilities to activation by hydronation (protonation) and (b) using the multitasking application program SKETCHER for the rapid evaluation of pH-dependent kinetic data by means of interactive manipulation of calculated curves. 2. The substantially lower basicity of (I) (pKa 0.91) than that of (II) (pKa 2.45) combined with retention of high reactivity permitted the pKa for the formation of the (Cys-25)-S-/(His-159)-Im+H ion-pair state of papain to be determined kinetically as 3.4, a value close to that (3.3) deduced by potentiometric difference titration [Lewis, Johnson and Shafer (1976) Biochemistry 15, 5009-5017] and lower than the value (approx. 4) often reported from pH-dependent kinetic studies. The higher values are now known to arise from inadequate data analysis that does not take account of other overlapping kinetically influential ionizations. 3. Re-evaluation of the extensive sets of pH-kcat/Km data for the hydrolysis of nine substrates by papain reported by Polgár and Halász (1978) (Eur. J. Biochem. 88, 513-521) by making use of SKETCHER, the known pKa value (3.4) from the reaction with compound (I) and two additional kinetically influential pKa values deduced from the reaction with compound (II) now permits the identification of the pH-dependent events in reactions of papain with inhibitors and substrates. 4. A major conclusion is that, whereas in reactions of simple alkylating agents and compound (I) full nucleophilic character of (Cys-25)-S-/(His-159)-Im+H is provided by hydronic dissociation with pKa 3.3-3.4, in catalysis relatively little catalytic competence is produced consequent upon ion-pair formation. Substantial catalytic competence requires further hydronic dissociation with pKa approx. 4, and for cationic substrates further enhancement is produced by hydronic dissociation with pKa approx. 5. 5. The present work, together with the kinetic analysis of reactions of papain in alkaline media reported by Mellor, Thomas, Topham and Brocklehurst [Biochem. J. (1993) 290, 289-296], defines the kinetically influential ionizations of papain as 3.4, 4.0, 5.0, 8.3 and 10.0 of which 3.4 and 8.3 relate to the formation and subsequent dehydronation of the ion-pair state.(ABSTRACT TRUNCATED AT 400 WORDS)