An unequivocal example of cysteine proteinase activity affected by multiple electrostatic interactions

Understanding Cysteine Reactivity in Protein Environments with Electric Fields

The role cysteine residues play in proteins is mediated by their protonation state, whereby the thiolate form of the side chain is highly reactive while the thiol form is more inert. However, the pKa of cysteine residues is hard to predict as it can differ widely from its reference value in solution, an effect that is accentuated by local effects in the heterogeneous protein environment. Here, we present a new approach to the prediction of cysteine reactivity based on electric field calculations at the thiol/thiolate group. We validated our approach by predicting the protonation state of cysteine residues in different protein environments (in the active site, at the protein surface, and buried within the protein interior), including Cys-25 in papaya protease omega, which was proven problematic for the more traditional constant pH molecular dynamics (MD) technique. We predict pKa shifts consistent with experimental observations, and the decomposition of the electric fields into contributions from molecular fragments provides a direct handle to rationalize local pH and pKa effects in proteins without introducing parameters other than those of the force field used for MD simulations.

An Unusual Member of the Papain Superfamily: Mapping the Catalytic Cleft of the Marasmius oreades agglutinin (MOA) with a Caspase Inhibitor

Papain-like cysteine proteases (PLCPs) constitute the largest group of thiol-based protein degrading enzymes and are characterized by a highly conserved fold. They are found in bacteria, viruses, plants and animals and involved in a number of physiological and pathological processes, parasitic infections and host defense, making them interesting targets for drug design. The Marasmius oreades agglutinin (MOA) is a blood group B-specific fungal chimerolectin with calcium-dependent proteolytic activity. The proteolytic domain of MOA presents a unique structural arrangement, yet mimicking the main structural elements in known PLCPs. Here we present the X-ray crystal structure of MOA in complex with Z-VAD-fmk, an irreversible caspase inhibitor known to cross-react with PLCPs. The structural data allow modeling of the substrate binding geometry and mapping of the fundamental enzyme-substrate interactions. The new information consolidates MOA as a new, yet strongly atypical member of the papain superfamily. The reported complex is the first published structure of a PLCP in complex with the well characterized caspase inhibitor Z-VAD-fmk.

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.

Structural insights into the substrate specificity and activity of ervatamins, the papain-like cysteine proteases from a tropical plant, Ervatamia coronaria

Multiple proteases of the same family are quite often present in the same species in biological systems. These multiple proteases, despite having high homology in their primary and tertiary structures, show deviations in properties such as stability, activity, and specificity. It is of interest, therefore, to compare the structures of these multiple proteases in a single species to identify the structural changes, if any, that may be responsible for such deviations. Ervatamin-A, ervatamin-B and ervatamin-C are three such papain-like cysteine proteases found in the latex of the tropical plant Ervatamia coronaria, and are known not only for their high stability over a wide range of temperature and pH, but also for variations in activity and specificity among themselves and among other members of the family. Here we report the crystal structures of ervatamin-A and ervatamin-C, complexed with an irreversible inhibitor 1-[l-N-(trans-epoxysuccinyl)leucyl]amino-4-guanidinobutane (E-64), together with enzyme kinetics and molecular dynamic simulation studies. A comparison of these results with the earlier structures helps in a correlation of the structural features with the corresponding functional properties. The specificity constants (k(cat)/K(m)) for the ervatamins indicate that all of these enzymes have specificity for a branched hydrophobic residue at the P2 position of the peptide substrates, with different degrees of efficiency. A single amino acid change, as compared to ervatamin-C, in the S2 pocket of ervatamin-A (Ala67-->Tyr) results in a 57-fold increase in its k(cat)/K(m) value for a substrate having a Val at the P2 position. Our studies indicate a higher enzymatic activity of ervatamin-A, which has been subsequently explained at the molecular level from the three-dimensional structure of the enzyme and in the context of its helix polarizibility and active site plasticity.

The catalytic mode of cysteine proteinases of papain (C1) family

The Proton Inventory (PI) method has been applied in the hydrolysis of synthetic substrates by papain, chymopapain and stem bromelain, comparing also their corresponding pH-(k(cat)/K(m)) profiles, and it was found: (a) k(cat)/K(m)=k(1), and thus K(S)=k(2)/k(1) is a dynamic equilibrium constant, (b) bowed-downward PI for k(cat)/K(m) exhibiting large inverse SIE, and (c) linear PI exhibiting large normal SIE for K(S), k(2) and k(3). A novel finding of this work is that the association of substrates onto all three studied cysteine proteinases proceeds via a stepwise pathway, in contrast to purely concerted pathways found previously for both acylation and deacylation. A hydrogen bond, which seems more likely to be developed across a pK(a)-value close to 4.00, connecting [see text] (papain/chymopapain or bromelain numbering), constitutes another novelty of this work.

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

Modulation of the electrostatic charge at the active site of foot-and-mouth-disease-virus leader proteinase, an unusual papain-like enzyme

The leader proteinase (L(pro)) of foot-and-mouth-disease virus is an unusual papain-like cysteine proteinase. Synthesized without an N-terminal pro precursor region, it frees itself from the growing polypeptide chain by cleavage at its own C-terminus. It also possesses a unique electrostatic environment around the active site, essentially due to Asp(163), which orients the catalytic histidine residue, and Asp(164); the equivalent residues in papain are Asn(175) and Ser(176). The importance of these residues for L(pro) activity was examined by site-directed mutagenesis. Replacement of Asp(163) with asparagine reduced activity by five-fold towards a hexapeptide substrate and slightly delayed self-processing when expressed in rabbit reticulocyte lysates. However, no effect on the cleavage of the only known cellular substrate of L(pro), eukaryotic initiation factor 4GI (eIF4GI), was observed. In contrast, replacement of Asp(164) by either alanine, asparagine or lysine abrogated activity towards the hexapeptide. Furthermore, in all cases, the onset of both self-processing and eIF4GI cleavage were significantly delayed; the reaction rates were also diminished compared with those of the wild-type enzyme. The alanine-substituted enzyme was least affected, followed by those substituted with asparagine and lysine. The double mutant protein in which both aspartate residues were replaced by asparagine was most severely affected; it failed to complete either self-processing or eIF4GI cleavage within 3 h, compared with the 8 min required by the wild-type enzyme. Hence, we propose that the electrostatic charge of Asp(164), and to a lesser extent that of Asp(163), is extremely important for L(pro) to attain full activity upon synthesis.

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 two sides of enzyme-substrate specificity: lessons from the aspartic proteinases

Like most proteolytic enzymes, the aspartic proteinases bind substrates and most inhibitors within an extended active site cleft. Bound ligands typically adopt a beta-strand conformation. Interactions with groups on both sides of the cleft determine the primary as well as secondary specificity of the enzymes. We have pursued the discovery of the sometimes subtle distinctions between members of the aspartic proteinase family by two routes. In the first case, we have constructed sets of oligopeptide substrates with systematic variation in each position to assess interactions at one position at a time. In the second type of experiment, we have altered residues of the enzymes in order to test theories of selectivity. The combination of the two approaches has provided a better understanding of the forces involved in determining specificity of enzyme action.

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.

The prosequence of procaricain forms an alpha-helical domain that prevents access to the substrate-binding cleft

BACKGROUND

Cysteine proteases are involved in a variety of cellular processes including cartilage degradation in arthritis, the progression of Alzheimer's disease and cancer invasion: these enzymes are therefore of immense biological importance. Caricain is the most basic of the cysteine proteases found in the latex of Carica papaya. It is a member of the papain superfamily and is homologous to other plant and animal cysteine proteases. Caricain is naturally expressed as an inactive zymogen called procaricain. The inactive form of the protease contains an inhibitory proregion which consists of an additional 106 N-terminal amino acids; the proregion is removed upon activation.

RESULTS

The crystal structure of procaricain has been refined to 3.2 A resolution; the final model consists of three non-crystallographically related molecules. The proregion of caricain forms a separate globular domain which binds to the C-terminal domain of mature caricain. The proregion also contains an extended polypeptide chain which runs through the substrate-binding cleft, in the opposite direction to that of the substrate, and connects to the N terminus of the mature region. The mature region does not undergo any conformational change on activation.

CONCLUSIONS

We conclude that the rate-limiting step in the in vitro activation of procaricain is the dissociation of the prodomain, which is then followed by proteolytic cleavage of the extended polypeptide chain of the proregion. The prodomain provides a stable scaffold which may facilitate the folding of the C-terminal lobe of procaricain.

Crystal structure of a caricain D158E mutant in complex with E-64

The structure of the D158E mutant of caricain (previously known as papaya protease omega) in complex with E-64 has been determined at 2.0 A resolution (overall R factor 19.3%). The structure reveals that the substituted glutamate makes the same pattern of hydrogen bonds as the aspartate in native caricain. This was not anticipated since in the native structure there is insufficient room to accommodate the glutamate side chain. The glutamate is accommodated in the mutant by a local expansion of the structure demonstrating that small structural changes are responsible for the change in activity.

Inhibition of cysteine proteinases by Carica papaya cystatin produced in Escherichia coli

A papaya cystatin (Cst)-encoding cDNA clone was isolated from a papaya leaf cDNA library and the active protein produced in Escherichia coli. The amino-acid sequence reveals a protein of 11,262 Da with over 40% identity to other published plant Cst. Unique features of the papaya Cst include a single Cys residue, variation in the papain-binding region, and the first reported inhibition of papaya proteinase IV by a Cst.

Theory of electrostatic interactions in macromolecules

In the past year, substantial progress has been made in the modeling of electrostatic interactions in biomolecules. This review highlights advances in the following areas: first, the efficient computation of long-range electrostatic interactions in detailed molecular simulations; second, the application of the Poisson-Boltzmann electrostatic model in conformational analysis; third, the application of the Poisson-Boltzmann model in quantum chemistry calculations; fourth, the development of atomic parameters; and finally, the modeling of ionization equilibria in proteins.