cDNA cloning and phylogenetic analysis of pancreatic serine proteases from Japanese flounder,Paralichthys olivaceus

To better understand the digestive physiology and phylogeny of the pancreatic serine proteases of teleosts, we cloned trypsin, chymotrypsin and elastase from flounder (Paralichthys olivaceus). Fifty phage plaques randomly chosen from a flounder pancreatic cDNA library were found to contain three species of trypsin, two species of chymotrypsin and four species of elastase. cDNAs of two species of carboxypeptidase A, one carboxypeptidase B and lipase were also obtained. In total, 23 out of 24 digestive enzyme cDNAs were those of proteolytic enzymes. Such a high ratio of proteolytic enzyme cDNA in the pancreas may reflect the carnivorous feeding habits of flounder. A phylogenetic comparison of the peptide sequences of flounder enzymes with those of other teleosts and mammals suggested that duplication of trypsin, chymotrypsin and elastase occurred before the divergence of the ray finned fish. It is also hypothesized that functional descendants of both duplicated genes of elastase exist in the teleosts and mammals, whereas only one of the genes of trypsin and chymotrypsin gave rise to the functional descendants in the teleosts but not in the mammals.     

Evidence of homology between the β-chain of human haptoglobin and the chymotrypsin family of serine proteases

Amino acid sequences from the β-chain of human haptoglobin are compared with those sequences known for the serine proteases of the chymotrypsin family. In a comparison of some 171 residues of the haptoglobin β-chain (approximately 60% of the protein molecule), approximately 30% of these are identical to residues occurring in sequences of either bovine trypsin, bovine chymotrypsin A, bovine chymotrypsin B, porcine elastase, or bovine thrombin B-chain, and an additional 10% are chemically similar. A combined comparison of the haptoglobin β-chain with the above five serine proteases gave an identity of 56% and a chemical similarity of 11%. Similarity of primary structure is also striking around two of the five half-cystinyl residues so far characterized in long lengths of sequence. These data provide substantial evidence that the β-chain of haptoglobin is homologous to the chymotrypsin family of serine proteases. Proposals are also presented to explain the occurrence of internal homology in the N-terminal region of the β-chain.     

Amino acid sequence of human D of the alternative complement pathway

The primary structure of human D, the serine protease activating the C3 convertase of the alternative complement pathway, has been deduced by sequencing peptides derived from various chemical (CNBr and o-iodosobenzoic acid) and enzymatic (trypsin, lysine protease, Staphylococcus aureus V8 protease, and chymotrypsin) cleavages. Carboxypeptidase A was also used to confirm the COOH-terminal sequence. The peptides were purified by high-pressure liquid chromatography. The proposed sequence of human D contains 222 amino acids and has a calculated molecular weight of 23 748. It exhibits a high degree of homology with other serine proteases, especially around the NH2-terminus as well as the three residues corresponding to the active-site His-57, Asp-102, and Ser-195 (chymotrypsinogen numbering). This sequence homology is highest (40%) with plasmin, intermediate (35%) with pancreatic serine proteases, such as elastase, trypsin, chymotrypsin, and kallikrein, and least (30%) with the serum enzymes thrombin and factor X. D, however, exhibits only minimal amino acid homology with the other sequenced complement serine proteases, Clr (25%) and Bb (20%). The substitution of a basic lysine for a neutral amino acid three residues NH2-terminal to the active-site serine as well as a small serine residue for a bulky aromatic amino acid at position 215 (chymotrypsinogen numbering) in the binding pocket may be important in determining the exquisite substrate specificity of D. The presence of His-40 which interacts with Asp-194 (chymotrypsinogen numbering) to stabilize other serine protease zymogens [Freer, S. T., Kraut, J., Robertus, J. D., Wright, H. T., & Xuong, N. H. (1970) Biochemistry 9, 1997] argues in favor of such a D precursor molecule.     

Identification and expression analysis of the Steinernema carpocapsae elastase-like serine protease gene during the parasitic stage

A cDNA encoding elastase was isolated from Steinernema carpocapsae by suppression subtractive hybridization and rapid amplification of 5′ cDNA ends. The predicted protein contained a 19-aa signal peptide, a 44-aa N-terminal propeptide, and a 264-aa mature protein with a predicted molecular mass of 28,949 Da and a theoretical pI of 8.88. BLAST analysis showed 27-35% amino acid sequence identity to serine proteases from insects, mammals, fish and other organisms. The Sc-ela gene contains three exons and two introns with at least two copies in the S. carpocapsae genome. Expression analysis indicated that the Sc-ela gene was upregulated during the initial parasitic stage. Sequence comparison and evolutionary marker analysis revealed that Sc-ELA was a member of the elastase serine protease family with potential degradative, developmental and fibrinolytic activities. Homology modeling showed that Sc-ELA adopts a two β-barrel fold typical of trypsin-like serine proteases, and phylogenetic analysis indicates that Sc-ELA branched off early during elastase evolution.     

The inactive subunit of ruminant procarboxypeptidase A-S6 complexes. Structural basis of inactivity and physiological role

Subunit III has so far been found only in the pancreas of ruminants in a non-covalent association (procarboxypeptidase A-S6) with two different proteins: the procarboxypeptidase A itself (subunit I) and a C-type chymotrypsinogen (subunit II). In contrast with these latter two proteins, which are zymogens of pancreatic proteases, subunit III seems to be devoid of any activity towards specific substrates of pancreatic proteases. However, it possesses a weakly functional active site which allows it to hydrolyze a non-specific ester, p-nitrophenyl acetate, and to react with several active-site titrants. The binding of proflavin to subunit III shows that this protein owns a non-polar binding site with a very high Kd compared to that of chymotrypsin. The comparison of the amino acid sequences of subunit III and some serine proteases showed that subunit III is closely related to an elastase. Models of the tertiary structure of subunit III suggest a conformational modification that affects the substrate binding and could explain the lack of specific enzymatic activity. The presence of subunit III in the ternary complex is not related to an enzymatic function. This protein does not participate in the activation process of subunit I but prevents the denaturation of this subunit at low pH. This may represent its biological role in the acidic environment of the duodenum in ruminants.     

Homology modelling of the catalytic domain of early mammalian protein C: evolution of structural features

By means of a novel cDNA-based strategy employing the maximum parsimony principle, we have previously deduced probable amino acid sequences for the catalytic domains of the early mammalian ancestors of each of the five extant vitamin K-dependent serine proteases of coagulation, and for their common ancestor from a still earlier stage of vertebrate evolution. In the present study, we employed one of these sequences to construct a molecular model of the catalytic domain of early mammalian protein C and to explore its functional architecture. Following the domain’s progression from the common ancestor of the vitamin K-dependent serine proteases toward extant human protein C, this novel application of homology modelling to a reconstructed amino acid sequence has allowed us to trace the evolution of structural features in a vital coagulation protein. Received: 23 May 1997 / Accepted: 23 July 1997     

A serine alkaline protease from the fungusConidiobolus coronatus with a distinctly different structure than the serine protease subtilisin Carlsberg

In view of the functional similarities between subtilisin Carlsberg and the alkaline protease fromConidiobolus coronatus, the biochemical and structural properties of the two enzymes were compared. In spite of their similar biochemical properties, e.g., pH optima, heat stability, molecular mass, pI, esterase activity, and inhibition by diisopropyl fluorophosphate and phenylmethlysulfonylfluoride, the proteases were structurally dissimilar as revealed by (1) their amino acid compositions, (2) their inhibition by subtilisin inhibitor, (3) their immunological response to specific anti-Conidiobolus protease antibody, and (4) their tryptic peptide maps. Our results demonstrate that although they are functionally analogous, theConidiobolus protease is structurally distinct from subtilisin Carlsberg. TheConidiobolus protease was also different from other bacterial and animal proteases (e.g. pronase, protease K, trypsin, and chymotrypsin) as evidenced by their lack of response to anti-Conidiobolus protease antibody in double diffusion and in neutralization assays. TheConidiobolus serine protease fails to obey the general rule that proteins with similar functions have similar primary sequences and, thus, are evolutionarily related. Our results strengthen the concept of convergent evolution for serine proteases and provide basis for research in evolutionary relationships among fungal, bacterial, and animal proteases.     

α -Aminoalkylphosphonate DI(Chlorophenyl) Esters as Inhibitors of Serine Proteases

Abstract

α -Aminoalkylphosphonate di(chlorophenyl) esters and (α -aminoalkyl)phenylphosphinate phenylesters have been tested as irreversible inhibitors of human neutrophil elastase, porcine pancreatic elastase and chymotrypsin, serine proteases important in biochemical processes. Peptidyl derivatives of diphenyl (α -aminoalkyl) phosphonates have previously been shown to be potent and specific inhibitors of serine proteases at low concentrations. Addition of a halogen to the phenoxy group of the inhibitors should make the leaving group more electrophilic, and thus more reactive. Peptide phosphonate inhibitors with chlorine in the meta- or para- positions of the phenoxy ester moiety were synthesized and shown to be potent inhibitors of elastase. Tripeptide phosphonates are more potent inhibitors than dipeptide phosphonates, however, addition of the halogen did not increase the inhibitory potency of these phosphonates with elastase compared to the non-halogenated phosphonates. In the case of chymotrypsin, the halogenated phenoxy esters were more reactive, possibly due to an alternate binding mode. The novel (α -aminoalkyl)phenylphosphinate phenylesters were poor inhibitors of serine proteases.     

Biochemical characterization of digestive proteases and carbohydrases of the carob moth,Ectomyelois ceratoniae (Zeller) (Lepidoptera: Pyralidae)

Digestive proteinases and carbohydrases of Ectomyelois ceratoniae (Zeller) larvae were investigated using appropriate substrates and inhibitors. Midgut pH in larvae was determined to be slightly alkaline. Midgut extracts showed optimum activity for proteolysis of hemoglobin at pH 9–12. Midgut proteinases also hydrolyzed the synthetic substrates of trypsin, chymotrypsin, and elastase at pH 8–11. Maximum digestive α-amylase activity was also observed at pH 8–11. However, optimum activity for α- and β-glucosidase occurred at pH 5–8. Alpha- and β-galactosidases optimum activities occurred at pH 5 and pH 6, respectively. Inhibitors of serine proteases were effective on midgut serine proteases (trypsin and chymotrypsin proteases). Zymogram analyses revealed at least five bands of total proteolytic activity in the larval midgut. Protease-specific zymogram analyses revealed at least four, two, and one isozymes for trypsin-, chymotrypsin-, and elastase-like activities respectively. Two α-amylase isozymes were found in the midgut of fifth instar larvae and in the whole bodies of 1st through 5th instar larvae. Zymogram studies also revealed the presence of one and two bands of activity for β- and α-glucosidase, respectively. Recycling of α-amylase and proteases in the larval midgut was not complete. At least one isozyme of trypsin, chymotrypsin, elastase, and α-amylase were not recycled and were observed in the larval hindgut.     

Crystal structure of bovine procarboxypeptidase A-S6 subunit III, a highly structured truncated zymogen E.

Subunit III, a defective serine endopeptidase lacking the typical N-terminal hydrophobic dipeptide is secreted by the pancreas of ruminant species as part of the bovine ternary complex procarboxypeptidase A-S6. Two monoclinic crystal forms were obtained and subsequently used to solve its X-ray structure. The highest resolution model of subunit III was refined at 1.7 A resolution to a crystallographic R-factor of 18.4%, with r.m.s. bond deviations from ideality of 0.012 A. About 80% of the model presents the characteristic architecture of trypsin-like proteases. The remaining zones, however, have well-defined, unique conformations. The regions from residues 70 to 80 and from 140 to 155 present maximum distances of 16 and 18 A relative to serine proteases and zymogens. Comparisons with the structures of porcine elastase 1 and chymotrypsinogen A indicate that the specific binding pocket of subunit III adopts a zymogen-like conformation and thus provide a basis for its inactivity. In general, the structural analysis of subunit III strongly suggests that it corresponds to a truncated version of a new class of highly structured elastase-like zymogen molecules. Based on the structures of subunit III and elastase 1, it is concluded that large concerted movements are necessary for the activation of zymogen E.     

The crystal structure of guamerin in complex with chymotrypsin and the development of an elastase-specific inhibitor

Guamerin, a canonical serine protease inhibitor from Hirudo nipponia, was identified as an elastase-specific inhibitor and has potential application in various diseases caused by elevated elastase concentration. However, the application of guamerin is limited because it also shows inhibitory activity against other proteases. To improve the selectivity of guamerin as an elastase inhibitor, it is essential to understand the binding mode of the inhibitor to elastase and to other proteases. For this purpose, we determined the crystal structure of guamerin in complex with chymotrypsin at 2.5 Å resolution. The binding mode of guamerin on elastase was explored from the model structure of guamerin/elastase. Guamerin binds to the hydrophobic pocket of the protease in a substrate-like manner using its binding loop. In order to improve the binding selectivity of guamerin to elastase, several residues in the binding loop were mutated and the inhibitory activities of the mutants against elastase and chymotrypsin were monitored. The substitution of the Met36 residue for Ala in the P1 site increased the inhibitory activity against elastase up to 14-fold, while the same mutant showed 7-fold decreased activity against chymotrypsin compared to the wild-type guamerin. Furthermore, the M36A guamerin mutant more effectively protected endothelial cells against cell damage caused by elastase than the wild-type guamerin.     

Alignment statistic for identifying related protein sequences

Summary Closely related proteins show an obvious kinship by having numerous matching amino acids in their aligned sequences. Kinship between anciently separated proteins requires a statistical evaluation to rule out fortuitous similarities. A simple statistic is developed which assumes equal probability for all codon pairs, and a table of critical values for amino acid sequence alignments of length 200 or less is presented. Applying this statistic toV andC regions of immunoglobulin chains, aligned on the basis of shared features of three-dimensional structure, provides evidence that theV andC sequences descended from a common ancestor. Similarly the distant evolutionary relationship of dehydrogenases, flavdoxin, and subtilisin, suggested by structural alignments, is verified. On the other hand, the statistic does not verify a common evolutionary origin for the heme binding pocket in globins and cytochromeb ₅. Empirical evidence from the distribution of MMD values of amino acid pairs in comparisons of misaligned polypeptide chains and from Monte Carlo trials of sequences aligned with arbitrary gaps supports the validity of the statistic.     

A practical method for calculating evolutionary trees from sequence data

In a previous paper (Klotz et a1., 1979) we described a method for determining evolutionary trees from sequence data when rates of evolution of the sequences might differ greatly. It was shown theoretically that the method always gave the correct topology and root when the exact number of mutation differences between sequences and from their common ancestor was known. However, the method is impractical to use in most situations because it requires some knowledge of the ancestor. In this present paper we describe another method, related to the previous one, in which a present-day sequence can serve temporarily as an ancestor for purposes of determining the evolutionary tree regardless of the rates of evolution of the sequences involved. This new method can be carried out with high precision without the aid of a computer, and it does not increase in difficulty rapidly as the number of sequences involved in the study increases, unlike other methods.     

Distinct contributions of residue 192 to the specificity of coagulation and fibrinolytic serine proteases

Archetypal members of the chymotrypsin family of serine proteases, such as trypsin, chymotrypsin, and elastase, exhibit relatively broad substrate specificity. However, the successful development of efficient proteolytic cascades, such as the blood coagulation and fibrinolytic systems, required the evolution of proteases that displayed restricted specificity. Tissue-type plasminogen activator (t-PA), for example, possesses exquisitely stringent substrate specificity, and the molecular basis of this important biochemical property of t-PA remains obscure. Previous investigations of related serine proteases, which participate in the blood coagulation cascade, have focused attention on the residue that occupies position 192 (chymotrypsin numbering system), which plays a pivotal role in determining both the inhibitor and substrate specificity of these enzymes. Consequently, we created and characterized the kinetic properties of new variants of t-PA that contained point mutations at position 192. These studies demonstrated that, unlike in coagulation serine proteases, Gln-192 does not contribute significantly to the substrate or inhibitor specificity of t-PA in physiologically relevant reactions. Replacement of Gln-192 with a glutamic acid residue did, however, decrease the catalytic efficiency of mature, two-chain t-PA toward plasminogen in the absence of a fibrin co-factor.     

Primary structure of two distinct rat pancreatic preproelastases determined by sequence analysis of the complete cloned messenger ribonucleic acid sequences

The mRNA sequences for two rat pancreatic elastolytic enzymes have been cloned by recombinant DNA technology and their nucleotide sequences determined. Rat elastase I mRNA is 1113 nucleotides in length, plus a poly(A) tail, and encodes a preproelastase of 266 amino acids. The amino acid sequence of the predicted active form of rat elastase I is 84% homologous to porcine elastase 1. Key amino acid residues involved in determining substrate specificity of porcine elastase 1 are retained in the rat enzyme. The activation peptide of the zymogen does not appear related to that of other mammalian pancreatic serine proteases. The mRNA for elastase I is localized in the rough endoplasmic reticulum of acinar cells, as expected for the site of synthesis of an exocrine secretory enzyme. Rat elastase II mRNA is 910 nucleotides in length, plus a poly(A) tail, and encodes a preproenzyme of 271 amino acids. The amino acid sequence is more closely related to porcine elastase 1 (58% sequence identity) than to the other pancreatic serine proteases (33-39% sequence identity). Predictions of substrate preference based upon key amino acid residues that define the substrate binding cleft are consistent with the broad specificity observed for mammalian pancreatic elastase 2. The activation peptide is similar to that of the chymotrypsinogens and retains an N-terminal cysteine available to form a disulfide link to an internal conserved cysteine residue.     

Inhibitory spectrum of mouse contrapsin and alpha-1-antitrypsin against mouse serine proteases

Contrapsin and alpha-1-antitrypsin have been recently characterized as major protease inhibitors in mouse plasma (Takahara, H. & Sinohara, H. (1982) J. Biol. Chem. 257, 2438-2446). We have studied the effects of the two inhibitors upon various serine proteases prepared from mouse tissues. Trypsin, plasmin and trypsin-like proteases of the submaxillary gland were inhibited by contrapsin but not by alpha-1-antitrypsin. On the other hand, chymotrypsin, elastase, and thrombin were inactivated by alpha-1-antitrypsin but not by contrapsin. Thus, their inhibitory spectra did not overlap each other in spite of their broad specificities. The inhibition of trypsin, chymotrypsin, and elastase was rapid and stoichiometric, whereas the inhibition of the other proteases was relatively slow. Contrapsin accounted for almost the total capacities of mouse plasma to inhibit both trypsin and submaxillary gland trypsin-like proteases, whereas alpha-1-antitrypsin was responsible for nearly all the capacities of plasma to inhibit both chymotrypsin and elastase.     

Functional and Structural Characterization of Vibrio cholerae Extracellular Serine Protease B,VesB

The chymotrypsin subfamily A of serine proteases consists primarily of eukaryotic proteases, including only a few proteases of bacterial origin. VesB, a newly identified serine protease that is secreted by the type II secretion system in Vibrio cholerae, belongs to this subfamily. VesB is likely produced as a zymogen because sequence alignment with trypsinogen identified a putative cleavage site for activation and a catalytic triad, His-Asp-Ser. Using synthetic peptides, VesB efficiently cleaved a trypsin substrate, but not chymotrypsin and elastase substrates. The reversible serine protease inhibitor, benzamidine, inhibited VesB and served as an immobilized ligand for VesB affinity purification, further indicating its relationship with trypsin-like enzymes. Consistent with this family of serine proteases, N-terminal sequencing implied that the propeptide is removed in the secreted form of VesB. Separate mutagenesis of the activation site and catalytic serine rendered VesB inactive, confirming the importance of these features for activity, but not for secretion. Similar to trypsin but, in contrast to thrombin and other coagulation factors, Na⁺ did not stimulate the activity of VesB, despite containing the Tyr²⁵⁰ signature. The crystal structure of catalytically inactive pro-VesB revealed that the protease domain is structurally similar to trypsinogen. The C-terminal domain of VesB was found to adopt an immunoglobulin (Ig)-fold that is structurally homologous to Ig-folds of other extracellular Vibrio proteins. Possible roles of the Ig-fold domain in stability, substrate specificity, cell surface association, and type II secretion of VesB, the first bacterial multidomain trypsin-like protease with known structure, are discussed.     

The evolutionary origin of proinsulin: Amino acid sequence homology with the trypsin-related serine proteases detected and evaluated by new statistical methods

Proinsulins and pancreatic serine proteases were analyzed for possible amino acid sequence similarity, using an adapted version of the nucleotide sequence alignment technique of Sankoff (1972). The technique allowed us to determine simultaneously the statistical significance of both the sequence alignment and the number of gaps necessary to achieve that alignment. In the course of this work, it was realized that a rigorous analysis required non-parametric statistics.For the B-chain (amino-terminal) of insulin a highly significant gap-free sequence alignment with the serine proteases was found. For the A-chain (carboxy-terminal) of insulin a sequence alignment of modest statistical significance with two gaps could be obtained, while the search for a corresponding alignment for the C-peptide remained unsuccessful. Presumably the rapid evolution of the C-peptide has obscured its origin. Reconstruction of ancestral sequences was of no help. In contrast to the amino acid sequences, three-dimensional structures of the two protein families are quite different.Considering current histophysiological understanding of ontogeny and phylogeny of exocrine and endocrine pancreas, the observed sequence similarity of proinsulins and serine proteases was interpreted to mean that the two protein families have diverged from a common genetic ancestor. Moreoever, from the organismic distribution of these proteins it was concluded that at least one serine protease existed first, and that proinsulin was generated after duplication of a serine protease gene and subsequent drastic modification, such as a large deletion. Thus proinsulin, basically an anabolic hormone, is derived from a serine protease, an enzyme involved in digestion. This constitutes a refinement of a similar proposal by Steiner et al. (1973).The emergence of proinsulin seems to have occurred after coelenterates diverged, and possibly before most other major animal phyla diverged from the line leading to vertebrates, i.e. 520 to 700 million years ago. The evolution of proinsulin seems to have paralleled the evolution of endocrine cells. Homology of the secreted products of endocrine and exocrine cells was most readily reconciled with a common embryological and phylogenetic origin of the two cell types, as considered by Pictet & Rutter (1972).     

Techniques for determining protein sequence evolution applied to primates

Each amino acid in a protein is considered to be an individual, mutable characteristic of the species from which the protein is extracted. For a branching tree representing the evolutionary history of the known sequences in different species, our computer programs use majority logic and parsimony of mutations to determine the most likely ancestral amino acid for each position of the protein at each node of the tree. The number of mutations necessary between the ancestral and present species is summed for each branch and the entire tree. The programs then move branches to make many different configurations, from which we select the one with the minimum number of mutations as the most likely evolutionary history. We used this method to elucidate primate phylogeny from sequences of fibrinopeptides, carbonic anhydrase, and the hemoglobin beta, delta and alpha chains. All available sequences indicate that the early Pongidae had diverged into two lines before the divergence of an ancestor for the human line alone. We have constructed some probable ancestral sequences at major points during primate evolution and have developed tentative trees showing the order of divergences and evolutionary distances among primate groups. Further questions on primate evolution could be answered in the future by the detemination of the appropriate sequences.