Roles of Tyr13 and Phe219 in the unique substrate specificity of pepsin B

Pepsin B is known to be distributed throughout mammalia, including carnivores. In this study, the proteolytic specificity of canine pepsin B was clarified with 2 protein substrates and 37 synthetic octapeptides and compared with that of human pepsin A. Pepsin B efficiently hydrolyzed gelatin but very poorly hydrolized hemoglobin. It was active against only a group of octapeptides with Gly at P2, such as KPAGF/LRL and KPEGF/LRL (arrows indicate cleavage sites). In contrast, pepsin A hydrolyzed hemoglobin but not gelatin and showed high activity against various types of octapeptides, such as KPAEF/FRL and KPAEF/LRL. The specificity of pepsin B is unique among pepsins, and thus, the enzyme provides a suitable model for analyzing the structure and function of pepsins and related aspartic proteinases. Because Tyr13 and Phe219 in/around the S2 subsites (Glu/Ala13 and Ser219 are common in most pepsins) appeared to be involved in the specificity of pepsin B, site-directed mutagenesis was undertaken to replace large aromatic residues with small residues and vice versa. The Tyr13Ala/Phe219Ser double mutant of pepsin B was found to demonstrate broad activity against hemoglobin and various octapeptides, whereas the reverse mutant of pepsin A had significantly decreased activity. According to molecular modeling of pepsin B, Tyr13 OH narrows the substrate-binding space and a peptide with Gly at P2 might be preferentially accommodated because of its high flexibility. The hydroxyl can also make a hydrogen bond with nitrogen of a P3 residue and fix the substrate main chain to the active site, thus restricting the flexibility of the main chain and strengthening preferential accommodation of Gly at P2. The phenyl moiety of Phe219 is bulky and narrows the S2 substrate space, which also leads to a preference for Gly at P2, while lowering the catalytic activity against other peptide types without making a hydrogen-bonding network in the active site.     

Primary structure of a precursor to the aspartic proteinase from Rhizomucor miehei shows that the enzyme is synthesized as a zymogen

In order to characterize the zymogen of the milk-clotting enzyme from Rhizomucor miehei, we constructed a cDNA library on pBR327 in Escherichia coli. Aspartic proteinase-specific recombinants were isolated by colony hybridization to a specific oligonucleotide mixture, and the cDNA sequence corresponding to a precursor form of the enzyme was determined. The deduced amino acid sequence shows that this secreted fungal proteinase is synthesized as a precursor. The first 22 amino acid residues in this precursor constitute a typical signal peptide. The amino acid sequence of the following 47-amino-acid-long prosegment shows homology to the prosegments from both the extracellular and intracellular vertebrate aspartic proteinases, and to the prosegments from the yeast and Mucor pusillus aspartic proteinases as well. These observations suggest that all aspartic proteinases are synthesized with a prosegment and that this prosegment is essential for the correct folding of all the mature enzymes. The active Rhizomucor miehei enzyme consists of 361 amino acid residues with a total molecular weight of 38,701. Clusters of identities around the active site cleft support the assumption that these proteinases have a common folding of their peptide chains. The disulphide bridges were localized in the fungal enzyme, and 2 N-glycosylation sites were identified.     

Chemical modification studies of Rhizomucor miehei protease: evidence for the role of basic amino acids in enzyme catalysis.

The effect of chemical modification on milk clotting and proteolytic activities of aspartyl protease obtained from Rhizomucor miehei NRRL 3500 was examined in the absence and the presence of its specific inhibitor pepstatin A. The effect on the ratio of milk clotting activity (MC) to proteolytic activity (PA), an index of the quality of milk clotting proteases was also determined. Modification of the enzyme with trinitrobenzenesulfonic acid, diethylpyrocarbonate and phenylglyoxal produced an increase in the ratio of MC/PA, while modification with 2- hydroxy-5-nitrobenzyl bromide did not affect the ratio. Modification with N-acetylimidazole resulted in a marginal increase in MC/PA ratio. Protection using pepstatin A during modification with phenylglyoxal, N-acetylimidazole and 2-hydroxy-5-nitrobenzyl bromide, protected both MC and PA. In the case of modification by diethylpyrocarbonate, pepstatin A protected only MC. Pepstatin A did not protect both the activities on the modification of the enzyme by trinitrobenzene sulfonic acid. These observations indicate the presence of arginine, tyrosine and tryptophan at the catalytic site of the enzyme, for eliciting MC and PA of the enzyme. In general, modification of the positively charged residues increases the MC/PA ratio of the enzyme. In addition the modified lysine residues responsible for the inactivation of the enzyme were not involved in the active site of the enzyme. Thus the lysine residues might have a secondary role in enzyme catalysis. Further, histidine at the catalytic site was found to be exclusively involved in milk clotting activity. The enzyme with modified histidine residues were more susceptible to autocatalysis, indicating that histidine residues protect the enzyme against autolysis.     

Penicillopepsin-JT2, a recombinant enzyme from Penicillium janthinellum and the contribution of a hydrogen bond in subsite S3 to k(cat)

The nucleotide sequence of the gene (pepA) of a zymogen of an aspartic proteinase from Penicillium janthinellum with a 71% identity in the deduced amino acid sequence to penicillopepsin (which we propose to call penicillopepsin-JT1) has been determined. The gene consists of 60 codons for a putative leader sequence of 20 amino acid residues, a sequence of about 150 nucleotides that probably codes for an activation peptide and a sequence with two introns that codes for the active aspartic proteinase. This gene, inserted into the expression vector pGPT-pyrG1, was expressed in an aspartic proteinase-free strain of Aspergillus niger var. awamori in high yield as a glycosylated form of the active enzyme that we call penicillopepsin-JT2. After removal of the carbohydrate component with endoglycosidase H, its relative molecular mass is between 33,700 and 34,000. Its kinetic properties, especially the rate-enhancing effects of the presence of alanine residues in positions P3 and P2' of substrates, are similar to those of penicillopepsin-JT1, endothiapepsin, rhizopuspepsin, and pig pepsin. Earlier findings suggested that this rate-enhancing effect was due to a hydrogen bond between the -NH- of P3 and the hydrogen bond accepting oxygen of the side chain of the fourth amino acid residue C-terminal to Asp215. Thr219 of penicillopepsin-JT2 was mutated to Ser, Val, Gly, and Ala. Thr219Ser showed an increase in k(cat) when a P3 residue was present in the substrate, which was similar to that of the wild-type, whereas the mutants Thr219Val, Thr219Gly, and Thr219Ala showed no significant increase when a P3 residue was added. The results show that the putative hydrogen bond alone is responsible for the increase. We propose that by locking the -NH- of P3 to the enzyme, the scissile peptide bond between P1 and P1' becomes distorted toward a tetrahedral conformation and becomes more susceptible to nucleophilic attack by the catalytic apparatus without the need of a conformational change in the enzyme.     

Synthesis of ethyl isovalerate using Rhizomucor miehei lipase: optimization

Immobilized lipase from Rhizomucor miehei (Lipozyme IM-20) was used to catalyze the esterification reaction between isovaleric acid and ethanol to synthesize ethyl isovalerate in n-hexane. Response surface methodology based on five-level four-variable central composite rotatable design was employed to optimize four important reaction variables such as enzyme/substrate E/S ratio, substrate concentration, incubation time, and temperature affecting the synthesis of ethyl isovalerate. The optimum conditions predicted for achieving maximum ester yield (500 mM) are as follows: E/S ratio, 48.41 g/mol; substrate concentration, 1 M; reaction time, 60 h; temperature, 60 degrees C. The predicted value matched well with experimentally obtained value of 487 mM.     

Roles of catalytic residues in alpha-amylases as evidenced by the structures of the product-complexed mutants of a maltotetraose-forming amylase.

The crystal structures of the four product-complexed single mutants of the catalytic residues of Pseudomonas stutzeri maltotetraose-forming alpha-amylase, E219G, D193N, D193G and D294N, have been determined. Possible roles of the catalytic residues Glu219, Asp193 and Asp294 have been discussed by comparing the structures among the previously determined complexed mutant E219Q and the present mutant enzymes. The results suggested that Asp193 predominantly works as the base catalyst (nucleophile), whose side chain atom lies in close proximity to the C1-atom of Glc4, being involved in the intermediate formation in the hydrolysis reaction. While Asp294 works for tightly binding the substrate to give a twisted and a deformed conformation of the glucose ring at position -1 (Glc4). The hydrogen bond between the side chain atom of Glu219 and the O1-atom of Glc4, that implies the possibility of interaction via hydrogen, consistently present throughout these analyses, supports the generally accepted role of this residue as the acid catalyst (proton donor).     

Cloning and sequence analysis of the glyceraldehyde-3-phosphate dehydrogenase gene from the zygomycetes fungus Rhizomucor miehei

Rhizomucor miehei is important from a biotechnological aspect in consequence of its content of aspartic proteinase, which has high milk-clotting activity. A genomic library of R. miehei NRRL 5901 has been constructed in a phage (Lambda Fix II) vector. The glyceraldehyde-3-phosphate dehydrogenase (gpd) gene was isolated from this library by hybridization of the recombinant phage clones with a gpd-specific gene probe generated by polymerase chain reaction. The complete nucleotide sequence encodes a putative polypeptide chain of 336 amino acids interrupted by 5 introns. The predicted amino acid sequence of this gene shows a high degree of sequence similarity to the glyceraldehyde-3-phosphate dehydrogenase proteins from yeast and filamentous fungi. The promoter region, containing a consensus TATA box, and 246-bp downstream from the putative stop codon were also determined. The possibility of using the gpd promoter in the construction of new transformation vectors is discussed.     

Specificity of milk-clotting enzymes towards bovine kappa-casein

Limited proteolysis of bovine kappa-casein has been investigated with porcine pepsin A and C, and with the 2 microbial proteinases Mucor miehei proteinase and Endothia parasitica proteinase. The liberated C-terminal glycomacropeptide of kappa-casein was isolated after precipitation in 3% trichloroacetic acid followed by high-performance gel-permeation chromatography on a TSK G3000 SW column. From amino acid analyses and N-terminal sequencing of the liberated peptide it is concluded that porcine pepsin A, C and Mucor miehei proteinase cleave the same bond as chymosin: Phe-105-Met-106 whereas Endothia parasitica proteinase cleaves the bond Ser-104-Phe-105.     

Stucture of the complex between <Emphasis Type="Italic">Mucor pusillus</Emphasis> pepsin and the key domain of κ-casein for site-directed mutagenesis: a combined molecular modeling and docking approach

In the structural-based mutagenesis of Mucor pusillus pepsin (MPP), understanding how κ-casein interacts with MPP is a great challenge for us to explore. Chymosin-sensitive peptide is the key domain of κ-casein that regulates milk clotting through the specific proteolytic cleavage of its peptide bond (Phe¹⁰⁵-Met¹⁰⁶) by MPP to produce insoluble para-κ-casein. Here, we built the model of this large peptide using molecular modeling technique. Docking study showed that MPP can accommodate the designed model with a favorable binding energy and the docked complex has proven to locally resemble the inhibitor-chymosin complex. The catalytic mechanism for the peptide model binding with MPP was explored in terms of substrate-enzyme interaction and property of contacting surface. Some critical amino acid residues in the substrate binding cleft have been identified as an important guide for further site-directed mutagenesis. Glu13 and Leu11 in the S3 region of MPP, predicted as the special mutation sites, were confirmed to retain clotting activity and decrease the proteolytic activity. These novel mutants may provide a promising application for improving cheese flavor.     

Substrate specificity of the alpha-L-arabinofuranosidase from Rhizomucor pusillus HHT-1

The alpha-L-arabinofuranosidase (AF) from the fungus Rhizomucor pusillus HHT-1 released arabinose at appreciable rates from (1-->5)-alpha-L-arabinofuranooligosaccharides, sugar beet arabinan and debranched arabinan. This enzyme preferentially hydrolyzed the terminal arabinofuranosyl residue [alpha-(1-->5)-linked] of the arabinan backbone rather than the arabinosyl side chain [alpha-(1-->3)-linked residues]. The enzyme-hydrolyzed arabinan reacted at and debranched the arabinan almost at the same rate, and the degree of conversion for both cases was 65%. Methylation analysis of arabinan showed that the arabinosyl-linkage proportions were 2:2:2:1, respectively, for (1-->5)-Araf, T-Araf, (1-->3, 5)-Araf and (1-->3)-Araf, while the ratios for the AF-digested arabinan shifted to 3:1:2:1. Enzyme digestion resulted in an increase in the proportion of (1-->5)-linked arabinose and a decrease in the proportion of terminal arabinose indicated this AF cleaved the terminal arabinosyl residue of the arabinan back bone [alpha-(1-->5)-linked residues]. Peak assignments in the 13C NMR spectra also confirmed this linkage composition of four kinds of arabinose residues. Both 1H and 13C NMR spectra are dominated by signals of the alpha-anomeric configuration of the arabinofuranosyl moieties. No signals were recorded for arabinopyranosyl moieties in the NMR spectra. Methylation and NMR analysis of native and AF-digested arabinan revealed that this alpha-L-arabinofuranosidase can only hydrolyse alpha-L-arabinofuranosyl residues of arabinan.     

Role of an electrostatic network of residues in the enzymatic action of the Rhizomucor miehei lipase family

We have used continuum electrostatic methods to investigate the role of electrostatic interactions in the structure, function, and pH-dependent stability of the fungal Rhizomucor miehei lipase (RmL) family. We identify a functionally important electrostatic network which includes residues S144, D203, H257, Y260, H143, Y28, R80, and D91 (residue numbering is from RmL). This network consists of residues belonging to the catalytic triad (S144, D203, H257), residues located in proximity to the active site (Y260), residues stabilizing the geometry of the active site (Y28, H143), and residues located in the lid (D91) or close to the first hinge (R80). The lid and the first hinge are associated with the interfacial activation of lipases, where an alpha-helical lid opens up by rotating around two hinge regions. All network residues are well conserved in a set of 12 lipase homologues, and 6 of the network residues are located in sequence motifs. We observe that the effects of modeled mutations R86L, D91N, and H257F on the pH-dependent electrostatic free energies differ significantly in the closed and open conformations of RmL. Mutation R86L is especially interesting since it stabilizes the closed conformation but destabilizes the open one. Site-site electrostatic interaction energies reveal that interactions between R86 and D61, D113, and E117 stabilize the open conformation.     

Identification of a glutamic acid and an aspartic acid residue essential for catalytic activity of aspergillopepsin II, a non-pepsin type acid proteinase

Aspergillopepsin II from Aspergillus niger var. macrosporus is a non-pepsin type or pepstatin-insensitive acid proteinase. To identify the catalytic residues of the enzyme, all acidic residues that are conserved in the homologous proteinases of family A4 were replaced with Asn, Gln, or Ala using site-directed mutagenesis. The wild-type and mutant pro-enzymes were heterologously expressed in Escherichia coli and refolded in vitro. The wild-type pro-enzyme was shown to be processed into a two-chain active enzyme under acidic conditions. Most of the recombinant mutant pro-enzymes showed significant activity under acidic conditions because of autocatalytic activation except for the D123N, D123A, E219Q, and E219A mutants. The D123A, E219Q, and E219A mutants showed neither enzymatic activity nor autoprocessing activity under acidic conditions. The circular dichroism spectra of the mutant pro- and mature enzymes were essentially the same as those of the wild-type pro- and mature enzyme, respectively, indicating that the mutant pro-enzymes were correctly folded. In addition, two single and one double mutant pro-enzyme, D123E, E219D, and D123E/E219D, did not show enzymatic activity under acidic conditions. Taken together, Glu-219 and Asp-123 are deduced to be the catalytic residues of aspergillopepsin II.     

Characterization of an alg2 mutant of the zygomycete fungus Rhizomucor pusillus

The zygomycete fungus Rhizomucor pusillus secretes an aspartic proteinase (MPP) that contains asparagine ( N )-linked oligosaccharides at two sites. Mutant strain 1116 defective in N -glycosylation secretes MPP with truncated oligo-saccharide chains. Lipid-linked oligosaccharides in mutant 1116 were labeled with [6-(3)H]glucosamine and [2-(3)H]mannose, prepared by cycles of solvent extraction, and analyzed by gel filtration chromatography on a Bio-Gel P-4 column after mild acid-hydrolysis. Mutant 1116 accumulated an intermediate, Man(1)GlcNAc(2)-dolichol pyrophosphate (PP-Dol), whereas wild-type strain F27 synthesized the fully assembled oligosaccharide precursor Glc(3)Man(9)GlcNAc(2)-PP-Dol. Consistent with this, alg2 encoding a mannosyltransferase in the lipid-linked oligosaccharide biosynthetic pathway in mutant 1116 had a 5 bp insertion that generated a stop codon in the middle of the coding sequence. Transformation of mutant 1116 with the intact alg2 gene on a pUC19-derived plasmid generated transformants that contained multicopies of alg2 at the alg2 locus. Glycosylation of the total proteins in the transformants was recovered to the same level as in strain F27, as determined with peroxidase-concanavalin A. These transformants produced MPP mainly with the same N -linked oligosaccharides as that produced by strain F27, but still with truncated oligosaccharides in small amounts. All of these data show that Alg2 is an alpha-1,3 or alpha-1,6 mannosyltransferase that elongates Man(1)GlcNAc(2)-PP-Dol to Man(2)GlcNAc(2)-PP-Dol. The slower growth of mutant 1116 was significantly recovered on introduction of alg2. The viability of the alg2 mutants of the zygomycete R.pusillus makes a contrast with the lethal effect of ALG2 mutations in the yeast Saccharomyces cerevisiae.     

The study of the bacteriophage T5 deoxynucleoside monophosphate kinase active site by site-directed mutagenesis

The amino acid residues essential for the enzymatic activity of bacteriophage T5 deoxyribonucleoside monophosphate kinase were determined using a computer model of the enzyme active site. By site-directed mutagenesis, cloning, and gene expression in E. coli, a series of proteins were obtained with single substitutions of the conserved active site amino acid residues—S13A, D16N, T17N, T17S, R130K, K131E, Q134A, G137A, T138A, W150F, W150A, D170N, R172I, and E176Q. After purification by ion exchange and affine chromatography electrophoretically homogeneous preparations were obtained. The study of the enzymatic activity with natural acceptors of the phosphoryl group (dAMP, dCMP, dGMP, and dTMP) demonstrated that the substitutions of charged amino acid residues of the NMP binding domain (R130, R172, D170, and E176) caused nearly complete loss of enzymatic properties. It was found that the presence of the OH-group at position 17 was also important for the catalytic activity. On the basis of the analysis of specific activity variations we assumed that arginine residues at positions 130 and 172 were involved in the binding to the donor γ-phosphoryl and acceptor α-phosphoryl groups, as well as the aspartic acid residue at position 16 of the ATP-binding site (P-loop), in the binding to some acceptors, first of all dTMP. Disproportional changes in enzymatic activities of partially active mutants, G137A, T138A, T17N, Q134A, S13A, and D16N, toward different substrates may indicate that different amino acid residues participate in the binding to various substrates.     

High-throughput screening of B factor saturation mutated Rhizomucor miehei lipase thermostability based on synthetic reaction

Conventional lipase screening methods are mostly based on hydrolytic activity, which may not always be the best method to assess the enzyme activity, especially for evaluating synthetic activity. Here we developed a high throughput and visual method to screen clones with high synthetic activity and used it to assess lipases thermostability. All mutants' lipase synthetic activity were identified through esterification of caprylic acid and ethanol with methyl red as the pH indicator adding in the substrates on according to the color change halo around the colony on culture plates since synthetic reaction was often accompanied with a rise in pH. After two rounds operation with the pH indicator screening method, we obtained a double mutant Asn120Lys/Lys131Phe from the Rhizomucor miehei lipase saturation mutated library based on amino acid residue B factors. The mutant's initial synthetic activity was a little higher than wild type and its thermostability in synthetic reaction was enhanced, which remained 63.1% residual activity after being heated at 70°C for 5h comparing to 51.0% of wild type. The double mutant with the two residue replacements balanced well between stability and activity. Yeast surface display technology and the pH indicator method, combined with colony screening were shown to facilitate high-throughput screening for lipase synthetic activity.     

Role of the invariant Asn345 and Asn435 residues in a leucine aminopeptidase from Bacillus kaustophilus as evaluated by site-directed mutagenesis

Role of the conserved Asn345 and Asn435 residues of Bacillus kaustophilus leucine aminopeptidase (BkLAP) was investigated by performing computer modeling and site-directed mutagenesis. Replacement of BkLAP Asn345 with Gln or Leu resulted in a dramatic reduction in enzymatic activity. A complete loss of the LAP activity was observed in Asn435 variants. Circular dichroism spectra were nearly identical for wild-type and all mutant enzymes, while measurement of intrinsic tryptophan fluorescence revealed the significant alterations of the microenvironment of aromatic amino acid residues in Asn345 and Asn435 replacements. Except for N435R and N435L, wild-type and other mutant enzymes showed a similar sensitivity towards temperature-induced denaturation. Computer modeling of the active-site structures of wild-type and mutant enzymes exhibits a partial or complete loss of the hydrogen bonding in the variants.     

Cloning and sequence analysis of cDNA encoding Rhizopus niveus lipase.

Complementary DNA encoding Rhizopus niveus lipase (RNL) was isolated from the R. niveus IF04759 cDNA library using a synthetic oligonucleotide corresponding to the amino acid sequence of the enzyme. A clone, which had an insert of 1.0 kilobase pairs, was found to contain the coding region of the enzyme. The lipase gene was expressed in Escherichia coli as a lacZ fusion protein. The mature RNL consisted of 297 amino acid residues with a molecular mass of 32 kDa. The RNL sequence showed significant overall homology to Rhizomucor miehei lipase and the putative active site residues were strictly conserved.     

In vivo and in vitro characterization of TEV protease mutants

Tobacco etch virus protease (TEVp) is frequently applied in the cleavage of fusion protein. However, production of TEV protease in Escherichia coli is hampered by low yield and poor solubility, and auto-cleavage of wild type TEVp gives rise to the loss-of-function. Previously it was reported that TEVp S219V displayed more stability, and TEVp variant containing T17S/N68D/I77V and double mutant L56V/S135G resulted in the enhanced production and solubility, respectively. Here, we introduced T17S/N68D/I77V in TEVp S219V to generate TEVpM1 and combined five amino acid mutations (T17S/L56V/N68D/I77V/S135G) in TEVp S219V to create TEVpM2. Among TEVp S219V, and two constructed variants, TEVpM2 displayed highest solubility and catalytic activity in vivo, using EmGFP as the solubility reporter, and the designed fusion protein as in vivo substrate containing an N-terminal hexahistidine tagged GST, a peptide sequence for thrombin and TEV cut and E. coli diaminopropionate ammonia-lyase. The purified TEVp mutants fused with double hexahistidine-tag at N and C terminus showed highest yield, solubility and cleavage efficiency. Mutations of five amino acid residues in TEVpM2 slightly altered protein secondary structure conformed by circular dichroism assay.     

Stabilization of flavobacterium meningosepticum glycerol kinase by introduction of a hydrogen bond

The thermostability of Flavobacterium meningosepticum glycerol kinase was increased by the change from Ser329 to Asp [Protein Eng., 14, 663-667 (2001)]. Based on a three-dimensional structure model of the mutant, we have postulated that a new charged-neutral hydrogen bond was formed between Asp329 and Ser414, and the formation of the hydrogen bond contributed to the stabilization of the tertiary structure and increased thermostability of the mutant enzyme. If the postulation is the case, FGK thermostabilization would be possible similarly by the single amino acid substitution from Ser414 to another amino acid which could form the hydrogen bond with Ser329. We did a single amino acid substitution of the wild-type enzyme from Ser414 to Asn. As we expected, S414N showed comparable thermostability to that of S329D. On the other hand, a difference in kinetic properties for ATP between S414N and S329D was observed.     

Pepsin as a catalyst for peptide synthesis: formation of peptide bonds not typical for pepsin substrate specificity.

Porcine pepsin in water solutions containing 15-28% of dimethylformamide at pH 5 and 20-37 degrees C catalysed the formation of peptide bonds between Z-Ala-Ala-Phe-OH and various amino acid or peptide derivatives. Substrate binding subsite S1' of pepsin demonstrated broad specificity in these reactions but revealed a certain preference for hydrophobic amino acid residues, including non-proteinous homophenylalanine, p-nitrophenylalanine, S-methylcysteine, as well as for those that contained, in addition to the hydrophobic elements, a group capable of donating a hydrogen bond, e.g. o-nitrotyrosine. This observation increases the range of peptides that might be prepared by pepsin-catalysed synthesis.