Structural basis of the unusual stability and substrate specificity of ervatamin C, a plant cysteine protease from Ervatamia coronaria

Virtual screening and molecular dynamics simulation study of plant protease inhibitors against SARS-CoV-2 envelope protein

Due to the outbreak of a new strain of pandemic coronavirus, there is a huge loss of economy and health. In 2021, some vaccines are recommended as emergency licensed vaccines to protect against the virus, and efforts are continuously ongoing to evaluate the vaccine safety measures for licensed vaccines. Recently, there was an increase in the cases of a new variant of coronavirus (omicron). Envelope protein plays an important role in virus packaging and assembly. If viral assembly is blocked, there is less chance of spreading the infection to another cell.In the present study, the plant protease inhibitors (PPIs) were screened against the envelope protein of SARS CoV 2. The structures were downloaded from the protein data bank. The plant protease inhibitors cystatin-I, Eravatmin, squash, Kunitz, Bowman-Birk, Alpha-amylase inhibitors, and potato serine protease inhibitors were screened and out of them Kunitz, alpha-amylase, and squash protease inhibitors have shown maximum binding energy. The molecular dynamics simulation was performed for docked complexes showing the lowest binding energy by NMA (normal mode analysis) to visualize the motion and stability of complexes. These plant-based protease inhibitors are a good target to fight against the new emerging strain of coronavirus because plant extracted compounds are natural and there is fewer side effect than synthetic compounds.

Directed-Mutagenesis of Flavobacterium meningosepticum Prolyl-Oligopeptidase and a Glutamine-Specific Endopeptidase From Barley

Wheat gluten proteins are the known cause of celiac disease. The repetitive tracts of proline and glutamine residues in these proteins make them exceptionally resilient to digestion in the gastrointestinal tract. These indigested peptides trigger immune reactions in susceptible individuals, which could be either an allergic reaction or celiac disease. Gluten exclusion diet is the only approved remedy for such disorders. Recently, a combination of a glutamine specific endoprotease from barley (EP-B2), and a prolyl endopeptidase from Flavobacterium meningosepticum (Fm-PEP), when expressed in the wheat endosperm, were shown to reasonably detoxify immunogenic gluten peptides under simulated gastrointestinal conditions. However useful, these "glutenases" are limited in application due to their denaturation at high temperatures, which most of the food processes require. Variants of these enzymes from thermophilic organisms exist, but cannot be applied directly due to their optimum activity at temperatures higher than 37°C. Though, these enzymes can serve as a reference to guide the evolution of peptidases of mesophilic origin toward thermostability. Therefore, a sequence guided site-saturation mutagenesis approach was used here to introduce mutations in the genes encoding Fm-PEP and EP-B2. A thermostable variant of Fm-PEP capable of surviving temperatures up to 90°C and EP-B2 variant with a thermostability of up 60°C were identified using this approach. However, the level of thermostability achieved is not sufficient; the present study has provided evidence that the thermostability of glutenases can be improved. And this pilot study has paved the way for more detailed structural studies in the future to obtain variants of Fm-PEP and EP-B2 that can survive temperatures ~100°C to allow their packing in grains and use of such grains in the food industry.

Survey of the rubber tree genome reveals a high number of cysteine protease-encoding genes homologous to Arabidopsis SAG12

Arabidopsis thaliana SAG12, a senescence-specific gene encoding a cysteine protease, is widely used as a molecular marker for the study of leaf senescence. To date, its potential orthologues have been isolated from several plant species such as Brassica napus and Nicotiana tabacum. However, little information is available in rubber tree (Hevea brasiliensis), a rubber-producing plant of the Euphorbiaceae family. This study presents the identification of SAG12-like genes from the rubber tree genome. Results showed that an unexpected high number of 17 rubber orthologues with a single intron were found, contrasting the single copy with two introns in Arabidopsis. The gene expansion was also observed in another two Euphorbiaceae plants, castor bean (Ricinus communis) and physic nut (Jatropha curcas), both of which contain 8 orthologues. In accordance with no occurrence of recent whole-genome duplication (WGD) events, most duplicates in castor and physic nut were resulted from tandem duplications. In contrast, the duplicated HbSAG12H genes were derived from tandem duplications as well as the recent WGD. Expression analysis showed that most HbSAG12H genes were lowly expressed in examined tissues except for root and male flower. Furthermore, HbSAG12H1 exhibits a strictly senescence-associated expression pattern in rubber tree leaves, and thus can be used as a marker gene for the study of senescence mechanism in Hevea.

Albizia lebbeck Seed Coat Proteins Bind to Chitin and Act as a Defense against Cowpea Weevil Callosobruchus maculatus

The seed coat is an external tissue that participates in defense against insects. In some nonhost seeds, including Albizia lebbeck, the insect Callosobruchus maculatus dies during seed coat penetration. We investigated the toxicity of A. lebbeck seed coat proteins to C. maculatus. A chitin-binding protein fraction was isolated from seed coat, and mass spectrometry showed similarity to a C1 cysteine protease. By ELM program an N-glycosylation interaction motif was identified in this protein, and by molecular docking the potential to interact with N-acetylglucosamine (NAG) was shown. The chitin-binding protein fraction was toxic to C. maculatus and was present in larval midgut and feces but not able to hydrolyze larval gut proteins. It did not interfere, though, with the intestinal cell permeability. These results indicate that the toxicity mechanism of this seed coat fraction may be related to its binding to chitin, present in the larvae gut, disturbing nutrient absorption.

Enhancement of proteolytic activity of a thermostable papain-like protease by structure-based rational design

Ervatamins (A, B and C) are papain-like cysteine proteases from the plant Ervatamia coronaria. Among Ervatamins, Ervatamin-C is a thermostable protease, but it shows lower catalytic efficiency. In contrast, Ervatamin-A which has a high amino acid sequence identity (∼90%) and structural homology (Cα rmsd 0.4 Å) with Ervatamin-C, has much higher catalytic efficiency (∼57 times). From the structural comparison of Ervatamin-A and -C, two residues Thr32 and Tyr67 in the catalytic cleft of Ervatamin-A have been identified whose contributions for higher activity of Ervatamin-A are established in our earlier studies. In this study, these two residues have been introduced in Ervatamin-C by site directed mutagenesis to enhance the catalytic efficiency of the thermostable protease. Two single mutants (S32T and A67Y) and one double mutant (S32T/A67Y) of Ervatamin-C have been generated and characterized. All the three mutants show ∼ 8 times higher catalytic efficiency (k cat/K m) than the wild-type. The thermostability of all the three mutant enzymes remained unchanged. The double mutant does not achieve the catalytic efficiency of the template enzyme Ervatamin-A. By modeling the structure of the double mutant and probing the role of active site residues by docking a substrate, the mechanistic insights of higher activity of the mutant protease have been addressed. The in-silico study demonstrates that the residues beyond the catalytic cleft also influence the substrate binding and positioning of the substrate at the catalytic centre, thus controlling the catalytic efficiency of an enzyme.

The origin and functional transition of P34

P34, a storage protein and major soybean allergen, has undergone a functional transition from a cysteine peptidase to a syringolide receptor. An exploration of the evolutionary mechanism of this functional transition is made. To identify homologous genes of P34, syntenic network was constructed using syntenic relationships from the Plant Genome Duplication Database. The collected homologous genes, along with SPE31, a highly homologous protein to P34 from the seeds of Pachyrhizus erosus, were used to construct a phylogenetic tree. The results show that multiple gene duplications, exon shuffling and following granulin domain loss and some critical point mutations are associated with the functional transition. Although some tests suggested the existence of positive selection, the possibility that random fixation under relaxation of purifying selection results in the functional transition is also supported. In addition, the genes Glyma08g12340 and Medtr8g086470 may belong to a new group within the papain family.

The structure of a thermostable mutant of pro-papain reveals its activation mechanism

Papain is the archetype of a broad class of cysteine proteases (clan C1A) that contain a pro-peptide in the zymogen form which is required for correct folding and spatio-temporal regulation of proteolytic activity in the initial stages after expression. This study reports the X-ray structure of the zymogen of a thermostable mutant of papain at 2.6 Å resolution. The overall structure, in particular that of the mature part of the protease, is similar to those of other members of the family. The structure provides an explanation for the molecular basis of the maintenance of latency of the proteolytic activity of the zymogen by its pro-segment at neutral pH. The structural analysis, together with biochemical and biophysical studies, demonstrated that the pro-segment of the zymogen undergoes a rearrangement in the form of a structural loosening at acidic pH which triggers the proteolytic activation cascade. This study further explains the bimolecular stepwise autocatalytic activation mechanism by limited proteolysis of the zymogen of papain at the molecular level. The possible factors responsible for the higher thermal stability of the papain mutant have also been analyzed.

Structural insight of dopamine β-hydroxylase, a drug target for complex traits, and functional significance of exonic single nucleotide polymorphisms

Background

Human dopamine β-hydroxylase (DBH) is an important therapeutic target for complex traits. Several single nucleotide polymorphisms (SNPs) have also been identified in DBH with potential adverse physiological effect. However, difficulty in obtaining diffractable crystals and lack of a suitable template for modeling the protein has ensured that neither crystallographic three-dimensional structure nor computational model for the enzyme is available to aid rational drug design, prediction of functional significance of SNPs or analytical protein engineering.

Principal Findings

Adequate biochemical information regarding human DBH, structural coordinates for peptidylglycine alpha-hydroxylating monooxygenase and computational data from a partial model of rat DBH were used along with logical manual intervention in a novel way to build an in silico model of human DBH. The model provides structural insight into the active site, metal coordination, subunit interface, substrate recognition and inhibitor binding. It reveals that DOMON domain potentially promotes tetramerization, while substrate dopamine and a potential therapeutic inhibitor nepicastat are stabilized in the active site through multiple hydrogen bonding. Functional significance of several exonic SNPs could be described from a structural analysis of the model. The model confirms that SNP resulting in Ala318Ser or Leu317Pro mutation may not influence enzyme activity, while Gly482Arg might actually do so being in the proximity of the active site. Arg549Cys may cause abnormal oligomerization through non-native disulfide bond formation. Other SNPs like Glu181, Glu250, Lys239 and Asp290 could potentially inhibit tetramerization thus affecting function.

Conclusions

The first three-dimensional model of full-length human DBH protein was obtained in a novel manner with a set of experimental data as guideline for consistency of in silico prediction. Preliminary physicochemical tests validated the model. The model confirms, rationalizes and provides structural basis for several biochemical data and claims testable hypotheses regarding function. It provides a reasonable template for drug design as well.

Selective and reversible thiol-pegylation, an effective approach for purification and characterization of five fully active ficin (iso)forms from Ficus carica latex

The latex of Ficus carica constitutes an important source of many proteolytic components known under the general term of ficin (EC 3.4.22.3) which belongs to the cysteine proteases of the papain family. So far, no data on the purification and characterization of individual forms of these proteases are available. An effective strategy was used to fractionate and purify to homogeneity five ficin forms, designated A, B, C, D1 and D2 according to their sequence of elution from a cation-exchange chromatographic support. Following rapid fractionation on a SP-Sepharose Fast Flow column, the different ficin forms were chemically modified by a specific and reversible monomethoxypolyethylene glycol (mPEG) reagent. In comparison with their un-derivatized counterparts, the mPEG-protein derivatives behaved differently on the ion-exchanger, allowing us for the first time to obtain five highly purified ficin molecular species titrating 1mol of thiol group per mole of enzyme. The purified ficins were characterized by de novo peptide sequencing and peptide mass fingerprinting analyzes, using mass spectrometry. Circular dichroism measurements indicated that all five ficins were highly structured, both in term of secondary and tertiary structure. Furthermore, analysis of far-UV CD spectra allowed calculation of their secondary structural content. Both these data and the molecular masses determined by MS reinforce the view that the enzymes belong to the family of papain-like proteases. The five ficin forms also displayed different specific amidase activities against small synthetic substrates like dl-BAPNA and Boc-Ala-Ala-Gly-pNA, suggesting some differences in their active site organization. Enzymatic activity of the five ficin forms was completely inhibited by specific cysteine and cysteine/serine proteases inhibitors but was unaffected by specific serine, aspartic and metallo proteases inhibitors.

Production of plant proteases in vivo and in vitro--a review

In the latest two decades, the interest received by plant proteases has increased significantly. Plant enzymes such as proteases are widely used in medicine and the food industry. Some proteases, like papain, bromelain and ficin are used in various processes such as brewing, meat softening, milk-clotting, cancer treatment, digestion and viral disorders. These enzymes can be obtained from their natural source or through in vitro cultures, in order to ensure a continuous source of plant enzymes. The focus of this review will be the production of plant proteases both in vivo and in vitro, with particular emphasis on the different types of commercially important plant proteases that have been isolated and characterized from naturally grown plants. In vitro approaches for the production of these proteases is also explored, focusing on the techniques that do not involve genetic transformation of the plants and the attempts that have been made in order to enhance the yield of the desired proteases.

C-Terminal extension of a plant cysteine protease modulates proteolytic activity through a partial inhibitory mechanism

The amino acid sequence of ervatamin-C, a thermostable cysteine protease from a tropical plant, revealed an additional 24-amino-acid extension at its C-terminus (CT). The role of this extension peptide in zymogen activation, catalytic activity, folding and stability of the protease is reported. For this study, we expressed two recombinant forms of the protease in Escherichia coli, one retaining the CT-extension and the other with it truncated. The enzyme with the extension shows autocatalytic zymogen activation at a higher pH of 8.0, whereas deletion of the extension results in a more active form of the enzyme. This CT-extension was not found to be cleaved during autocatalysis or by limited proteolysis by different external proteases. Molecular modeling and simulation studies revealed that the CT-extension blocks some of the substrate-binding unprimed subsites including the specificity-determining subsite (S2) of the enzyme and thereby partially occludes accessibility of the substrates to the active site, which also corroborates the experimental observations. The CT-extension in the model structure shows tight packing with the catalytic domain of the enzyme, mediated by strong hydrophobic and H-bond interactions, thus restricting accessibility of its cleavage sites to the protease itself or to the external proteases. Kinetic stability analyses (T(50) and t(1/2) ) and refolding experiments show similar thermal stability and refolding efficiency for both forms. These data suggest that the CT-extension has an inhibitory role in the proteolytic activity of ervatamin-C but does not have a major role either in stabilizing the enzyme or in its folding mechanism.

pH-Dependent urea-induced unfolding of stem bromelain: unusual stability against urea at neutral pH

Equilibrium unfolding of stem bromelain (SB) with urea as a denaturant has been monitored as a function of pH using circular dichroism and fluorescence emission spectroscopy. Urea-induced denaturation studies at pH 4.5 showed that SB unfolds through a two-state mechanism and yields DeltaG (free energy difference between the fully folded and unfolded forms) of approximately 5.0 kcal/mol and C(m) (midpoint of the unfolding transition) of approximately 6.5 M at 25 degrees C. Very high concentration of urea (9.5 M) provides unusual stability to the protein with no more structural loss and transition to a completely unfolded state.

Developing novel anthelmintics from plant cysteine proteinases

Intestinal helminth infections of livestock and humans are predominantly controlled by treatment with three classes of synthetic drugs, but some livestock nematodes have now developed resistance to all three classes and there are signs that human hookworms are becoming less responsive to the two classes (benzimidazoles and the nicotinic acetylcholine agonists) that are licensed for treatment of humans. New anthelmintics are urgently needed, and whilst development of new synthetic drugs is ongoing, it is slow and there are no signs yet that novel compounds operating through different modes of action, will be available on the market in the current decade. The development of naturally-occurring compounds as medicines for human use and for treatment of animals is fraught with problems. In this paper we review the current status of cysteine proteinases from fruits and protective plant latices as novel anthelmintics, we consider some of the problems inherent in taking laboratory findings and those derived from folk-medicine to the market and we suggest that there is a wealth of new compounds still to be discovered that could be harvested to benefit humans and livestock.

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.

A thermostable cysteine protease precursor from a tropical plant contains an unusual C-terminal propeptide: cDNA cloning, sequence comparison and molecular modeling studies

We report here the cloning and characterization of the entire cDNA of a papain-like cysteine protease from a tropical flowering plant. The 1098-bp ORF of the cDNA codify a protease precursor having a signal peptide of 19 amino acids, a cathepsin-L like N-terminal proregion of 114 amino acids, a mature enzyme part of 208 amino acids and a C-terminal proregion of 24 amino acids. The derived amino acid sequence of the mature part tallies with the thermostable cysteine protease Ervatamin-C--as was aimed at. The C-terminal proregion of the protease has altogether a different sequence pattern not observed in other members of the family and it contains a negatively charged helical zone. The three-dimensional model of the precursor, based on the homology modeling and X-ray structure, shows that the extended peptide stretch region of the N-terminal propeptide, covering the interdomain cleft, contains protruding side chains of positively charged residues. This study also indicates that the negatively charged zone of C-terminal propeptide may interact with the positively charged zone of the N-terminal propeptide in a cooperative manner in the maturation process of this enzyme.

Crystal structure of a papain-fold protein without the catalytic residue: a novel member in the cysteine proteinase family

A 31kDa cysteine protease, SPE31, was isolated from the seeds of a legume plant, Pachyrizhus erosus. The protein was purified, crystallized and the 3D structure solved using molecular replacement. The cDNA was obtained by RT PCR followed by amplification using mRNA isolated from the seeds of the legume plant as a template. Analysis of the cDNA sequence and the 3D structure indicated the protein to belong to the papain family. Detailed analysis of the structure revealed an unusual replacement of the conserved catalytic Cys with Gly. Replacement of another conserved residue Ala/Gly by a Phe sterically blocks the access of the substrate to the active site. A polyethyleneglycol molecule and a natural peptide fragment were bound to the surface of the active site. Asn159 was found to be glycosylated. The SPE31 cDNA sequence shares several features with P34, a protein found in soybeans, that is implicated in plant defense mechanisms as an elicitor receptor binding to syringolide. P34 has also been shown to interact with vegetative storage proteins and NADH-dependent hydroxypyruvate reductase. These roles suggest that SPE31 and P34 form a unique subfamily within the papain family. The crystal structure of SPE31 complexed with a natural peptide ligand reveals a unique active site architecture. In addition, the clear evidence of glycosylated Asn159 provides useful information towards understanding the functional mechanism of SPE31/P34.

Crystallization and preliminary X-ray diffraction studies of the cysteine protease ervatamin A from Ervatamia coronaria

The ervatamins are highly stable cysteine proteases that are present in the latex of the medicinal plant Ervatamia coronaria and belong to the papain family, members of which share similar amino-acid sequences and also a similar fold comprising two domains. Ervatamin A from this family, a highly active protease compared with others from the same source, has been purified to homogeneity by ion-exchange chromatography and crystallized by the vapour-diffusion method. Needle-shaped crystals of ervatamin A diffract to 2.1 A resolution and belong to space group C222(1), with unit-cell parameters a = 31.10, b = 144.17, c = 108.61 A. The solvent content using an ervatamin A molecular weight of 27.6 kDa is 43.9%, with a VM value of 2.19 A3 Da(-1) assuming one protein molecule in the asymmetric unit. A molecular-replacement solution has been found using the structure of ervatamin C as a search model.