Modeling and structural analysis of evolutionarily diverse S8 family serine proteases

Serine proteases are an abundant class of enzymes that are involved in a wide range of physiological processes and are classified into clans sharing structural homology. The active site of the subtilisin-like clan contains a catalytic triad in the order Asp, His, Ser (S8 family) or a catalytic tetrad in the order Glu, Asp and Ser (S53 family). The core structure and active site geometry of these proteases is of interest for many applications. The aim of this study was to investigate the structural properties of different S8 family serine proteases from a diverse range of taxa using molecular modeling techniques. In conjunction with 12 experimentally determined three-dimensional structures of S8 family members, our predicted structures from an archaeon, protozoan and a plant were used for analysis of the catalytic core. Amino acid sequences were obtained from the MEROPS database and submitted to the LOOPP server for threading based structure prediction. The predicted structures were refined and validated using PROCHECK, SCRWL and MODELYN. Investigation of secondary structures and electrostatic surface potential was performed using MOLMOL. Encompassing a wide range of taxa, our structural analysis provides an evolutionary perspective on S8 family serine proteases. Focusing on the common core containing the catalytic site of the enzyme, the analysis presented here is beneficial for future molecular modeling strategies and structure-based rational drug design.     

Cleavage specificity of the serine protease of Aeromonas sobria, a member of the kexin family of subtilases

Subtilisin-like proteases have been grouped into six families based on a sequence of the catalytic domain. One of the six is the kexin family, of which furin is a representative protease. All members of the kexin family, except one, are from eukaryotes. The one prokaryotic protease is a serine protease of Aeromonas sorbria (ASP). Here, we examined the substrate specificity of ASP based on the cleavage of short peptides. The results showed that ASP preferentially cleaves the peptide bond following two basic residues, one of which is Lys, but not the bond following a single basic residue. This indicates that the tertiary structure around the catalytic domain of ASP resembles, but is not identical to that of furin. Prekallikrein was cleaved into four fragments by ASP, indicating that the protein must be cleaved at specific sequences.     

Activation of prothrombin by two subtilisin-like serine proteases from Acremonium sp

Two novel subtilisin-like serine proteases (AS-E1 and -E2) that activate prothrombin have been identified in a culture of the fungus Acremonium sp. The enzymes were purified through repeated hydrophobic interaction chromatography. The N-terminal sequences of AS-E1 (34.4 kDa) and AS-E2 (32 kDa) showed high similarity to the internal sequences of two distinct subtilisin-like hypothetical proteins from Chaetomium globosum. Both enzymes proteolytically activated prothrombin to meizothrombin(desF1)-like molecules, while the activation cleavage seemed to occur at a site (Tyr(316)-Ile(317)) that is four residues proximal to the canonical Xa cleavage site (Arg(320)-Ile(321)). Both enzymes inhibited plasma clotting, possibly due to extensive degradation of fibrinogen and production of meizothrombin(desF1)-like molecule.     

Recombinant hirustasin: production in yeast,crystallization, and interaction with serine proteases.

A synthetic gene coding for the 55-amino acid protein hirustasin, a novel tissue kallikrein inhibitor from the leech Hirudo medicinalis, was generated by polymerase chain reaction using overlapping oligonucleotides, fused to the yeast alpha-factor leader sequence and expressed in Saccharomyces cerevisiae. Recombinant hirustasin was secreted mainly as incompletely processed fusion protein, but could be processed in vitro using a soluble variant of the yeast yscF protease. The processed hirustasin was purified to better than 97% purity. N-terminal sequence analysis and electrospray ionization mass spectrometry confirmed a correctly processed N-terminus and the expected amino acid sequence and molecular mass. The biological activity of recombinant hirustasin was identical to that of the authentic leech protein. Crystallized hirustasin alone and in complex with tissue kallikrein diffracted beyond 1.4 A and 2.4 A, respectively. In order to define the reactive site of the inhibitor, the interaction of hirustasin with kallikrein, chymotrypsin, and trypsin was investigated by monitoring complex formation in solution as well as proteolytic cleavage of the inhibitor. During incubation with high, nearly equimolar concentration of tissue kallikrein, hirustasin was cleaved mainly at the peptide bond between Arg 30 and Ile 31, the putative reactive site, to yield a modified inhibitor. In the corresponding complex with chymotrypsin, mainly uncleaved hirustasin was found and cleaved hirustasin species accumulated only slowly. Incubation with trypsin led to several proteolytic cleavages in hirustasin with the primary scissile peptide bond located between Arg 30 and Ile 31. Hirustasin appears to fall into the class of protease inhibitors displaying temporary inhibition.     

Purification and molecular characterization of subtilisin-like alkaline protease BPP-A from Bacillus pumilus strain MS-1

AIMS: The present study was conducted by screening zein-degrading bacteria in an attempt to obtain zein-degrading protease. METHODS AND RESULTS: Soil bacteria were screened by formation of a clear zone on zein plates. Characterization of a zein-degrading bacterium indicated a taxonomic affiliation to Bacillus pumilus, and was named MS-1 strain. The strain produced two different types of extracellular proteases, BPP-A and BPP-B. In this study, we purified and characterized BPP-A because it exhibited a higher ability to hydrolyze zein than BPP-B. When casein was used as the substrate, the optimal pH for BPP-A was 11.0. In BPP-A, zein was better substrate than casein at pH 13.0, whereas casein was better one than zein at pH 11.0. The bppA gene encoded a 383-amino acid pre-pro form of BPP-A, and mature BPP-A contained 275 amino acid residues. It was concluded that BPP-A belonged to the subtilisin family. CONCLUSION: A zein-degrading bacterium assigned to B. pumilus produced two different types of extracellular proteases, BPP-A and BPP-B. BPP-A exhibited an ability to hydrolyze zein in an extreme alkaline condition. SIGNIFICANCE AND IMPACT OF THE STUDY: This is a first report on screening for zein-degrading micro-organisms. The subtilisin-like protease BPP-A is possible to utilize as an industrial enzyme for the production of zein hydrolysates.     

The 0.93A crystal structure of sphericase: a calcium-loaded serine protease from Bacillus sphaericus

We have previously isolated sphericase (Sph), an extracellular mesophilic serine protease produced by Bacillus sphaericus. The Sph amino acid sequence is highly homologous to two cold-adapted subtilisins from Antarctic bacilli S39 and S41 (76% and 74% identity, respectively). Sph is calcium-dependent, 310 amino acid residues long and has optimal activity at pH 10.0. S41 and S39 have not as yet been structurally analysed.In the present work, we determined the crystal structure of Sph by the Eu/multiwavelength anomalous diffraction method. The structure was extended to 0.93A resolution and refined to a crystallographic R-factor of 9.7%. The final model included all 310 amino acid residues, one disulfide bond, 679 water molecules and five calcium ions. Although Sph is a mesophilic subtilisin, its amino acid sequence is similar to that of the psychrophilic subtilisins, which suggests that the crystal structure of these subtilisins is very similar.The presence of five calcium ions bound to a subtilisin molecule, as found here for Sph, has not been reported for the subtilisin superfamily. None of these calcium-binding sites correlates with the well-known high-affinity calcium-binding site (site I or site A), and only one site has been described previously. This calcium-binding pattern suggests that a reduction in the flexibility of the surface loops of Sph by calcium binding may be responsible for its adaptation to mesophilic organisms.     

Specificity and reactive loop length requirements for crmA inhibition of serine proteases

The viral serpin, crmA, is distinguished by its small size and ability to inhibit both serine and cysteine proteases utilizing a reactive loop shorter than most other serpins. Here, we characterize the mechanism of crmA inhibition of serine proteases and probe the reactive loop length requirements for inhibition with two crmA reactive loop variants. P1 Arg crmA inhibited the trypsin-like proteases, thrombin, and factor Xa, with moderate efficiencies (approximately 10(2)-10(4) M(-1)sec(-1)), near equimolar inhibition stoichiometries, and formation of SDS-stable complexes which were resistant to dissociation (k(diss) approximately 10(-7) sec(-1)), consistent with a serpin-type inhibition mechanism. Trypsin was not inhibited, but efficiently cleaved the variant crmA as a substrate (k(cat)/K(M) of approximately 10(6) M(-1) sec(-1)). N-terminal sequencing confirmed that the P1 Arg-P1'Cys bond was the site of cleavage. Altering the placement of the Arg in a double mutant P1 Gly-P1'Arg crmA resulted in minimal ability to inhibit any of the trypsin family proteases. This variant was cleaved by the proteases approximately 10-fold less efficiently than P1 Arg crmA. Surprisingly, pancreatic elastase was rapidly inhibited by wild-type and P1 Arg crmAs (10(5)-10(6) M(-1)sec(-1)), although with elevated inhibition stoichiometries and higher rates of complex dissociation. N-terminal sequencing showed that elastase attacked the P1'Cys-P2'Ala bond, indicating that crmA can inhibit proteases using a reactive loop length similar to that used by other serpins, but with variations in this inhibition arising from different effective P2 residues. These results indicate that crmA inhibits serine proteases by the established serpin conformational trapping mechanism, but is unusual in inhibiting through either of two adjacent reactive sites.     

NETosis occurs independently of neutrophil serine proteases

Neutrophils are primary host innate immune cells defending against pathogens. One proposed mechanism by which neutrophils prevent the spread of pathogens is NETosis, the extrusion of cellular DNA resulting in neutrophil extracellular traps (NETs). The protease neutrophil elastase (NE) has been implicated in the formation of NETs through proteolysis of nuclear proteins leading to chromatin decondensation. In addition to NE, neutrophils contain three other serine proteases that could compensate if the activity of NE was neutralized. However, whether they do play such a role is unknown. Thus, we deployed recently described specific inhibitors against all four of the neutrophil serine proteases (NSPs). Using specific antibodies to the NSPs along with our labeled inhibitors, we show that catalytic activity of these enzymes is not required for the formation of NETs. Moreover, the NSPs that decorate NETs are in an inactive conformation and thus cannot participate in further catalytic events. These results indicate that NSPs play no role in either NETosis or arming NETs with proteolytic activity.     

The two main endoproteases present in dark-induced senescent wheat leaves are distinct subtilisin-like proteases

We have previously reported the occurrence of two serine endoproteases (referred to as P1 and P2) in dark-induced senescent wheat (Triticum aestivum L.) leaves. P1 enzyme was already purified and identified as a subtilisin-like serine endoprotease (Roberts et al. in Physiol Plant 118:483–490, 2003). In this paper, we demonstrate by Western blot analysis of extracts obtained from dark-induced senescent leaves that an antiserum raised against P1 was able to recognise a second protein band of 78 kDa which corresponded to P2 activity. This result suggested that both enzymes must be structurally related. Therefore, we purified and characterised P2 activity. According to its biochemical and physical properties (inhibition by chymostatin and PMSF, broad pH range of activity, thermostability and ability to hydrolyse Suc-AAPF-pNA) P2 was classified as a serine protease with chymotrypsin-like activity. In addition, P2 was identified by mass spectrometry as a subtilisin-like protease distinct from P1. Western blot analysis demonstrated that P1 appeared in extracts from non-detached dark-induced senescent leaves but was undetectable in leaves senescing after nitrogen (N) deprivation. In contrast, P2 was already present in non-senescent leaves and showed increased levels in leaves senescing after N starvation or incubation in darkness. P1 signal was detected at late stages of ethephon or methyl jasmonate-induced senescence but was undetectable in senescent leaves from plants treated with abscisic acid. None of the three hormones have any effect on P2 protein levels. These results indicate that despite their biochemical and structural similarities, both enzymes are probably involved in different physiological roles.Electronic Supplementary Material Supplementary material is available to authorised users in the online version of this article at .     

Three serine proteases from the latex of Euphorbia cyparissias

Three serine-centred proteolytic enzymes, euphorbains y-1, ?2 and ?3, were isolated from the latex of Euphorbia cyparissias. These proteases have different specific activities to azocoll and CBZ glycine p-nitrophenyl ester. The pIs and M_rs of y-1, ?2 and ?3 are 5.2 and 67 000, 5.2 and 33 000, and 6.3 and 67 000, respectively. The enzymes, which are glycoproteins, are immunologically distinct from euphorbain 1, but clearly related to that enzyme in amino acid composition.     

Structure and energetics of protein-protein interactions: the role of conformational heterogeneity in OMTKY3 binding to serine proteases

Proteins with flexible binding surfaces can interact with numerous binding partners. However, this promiscuity is more difficult to understand in "rigid-body" proteins, whose binding results in little, or no, change in the position of backbone atoms. The binding of Kazal inhibitors to serine proteases is considered a classic case of rigid-body binding, although they bind to a wide range of proteases. We have studied the thermodynamics of binding of the Kazal serine protease inhibitor, turkey ovomucoid third domain (OMTKY3), to the serine protease subtilisin Carlsberg using isothermal titration calorimetry and have determined the crystal structure of the complex at very high resolution (1.1A). Comparison of the binding energetics and structure to other OMTKY3 interactions demonstrates that small changes in the position of side-chains can make significant contributions to the binding thermodynamics, including the enthalpy of binding. These effects emphasize that small, "rigid-body" proteins are still dynamic structures, and these dynamics make contributions to both the enthalpy and entropy of binding interactions.     

Induction of cysteine and serine proteases during xylogenesis in Zinnia elegans

The terminal process of xylogenesis, autolysis, is essential for the formation of a tubular system for conduction of water and solutes throughout the whole plant. Several hydrolase types are implicated in autolysis responsible for the breakdown of cytoplasm. Here, we characterize p48h-17 cDNA from in vitro tracheary elements (TEs) of Zinnia elegans which encodes a preproprotein similar to papain. The putative mature protein, a cysteine protease, has a molecular mass of 22,699 Da with a pI of 5.7. DNA gel blot analysis indicated that p48h-17 is likely encoded by one or two genes. The p48h-17 mRNA accumulated markedly in in vitro differentiating TEs, whereas it appeared not to be induced in response to senescence and wounding in the leaves or H₂O₂ challenge in the cultured mesophyll cells. In stems, the expression of the p48h-17 gene was preferentially associated with differentiating xylem. Activity gel assays demonstrated that a cysteine and a serine protease, which had apparent molecular masses of 20 kDa and 60 kDa, respectively, were markedly induced during in vitro TE differentiation. The cysteine protease activity was also preferentially present in the xylem of Zinnia stems. Transient expression of the p48h-17 cDNA in tobacco protoplasts resulted in the production of a 20 kDa cysteine protease. Taken together, the results indicate that the p48h-17 gene appears to be preferentially associated with xylogenesis, and both the cysteine and serine proteases might be involved in autolysis during xylogenesis.     

A subtilisin-like serine protease specifically expressed in reproductive organs in rice

To identify genes specifically expressed in flowering pistils and that are related to reproductive phenomena, simplified differential display was performed with cDNA obtained from pistils and ovaries at several stages. One clone preferentially expressed in pistils at flowering and 1 day before flowering was identified as a rice subtilisin-like serine protease (RSP1). Sequence comparisons revealed that RSP1 has several characteristics in common with preproproteins. RT-PCR, northern blot, and in situ hybridization analysis revealed that RSP1 mRNA accumulates in pistils and in the filaments of stamens, whereas mRNA was undetectable in tissue from leaves, roots, panicles, and embryos. The mRNA levels in pistils increased slightly at flowering and decreased afterwards. Possible roles of the subtilisin-like serine protease in plant reproduction are discussed. Received: 27 February 2000 / Revision accepted: 24 October 2000     

Structural basis of substrate specificity in the serine proteases.

Structure-based mutational analysis of serine protease specificity has produced a large database of information useful in addressing biological function and in establishing a basis for targeted design efforts. Critical issues examined include the function of water molecules in providing strength and specificity of binding, the extent to which binding subsites are interdependent, and the roles of polypeptide chain flexibility and distal structural elements in contributing to specificity profiles. The studies also provide a foundation for exploring why specificity modification can be either straightforward or complex, depending on the particular system.     

Potential for dramatic improvement in sequence alignment against structures of remote homologous proteins by extracting structural information from multiple structure alignment

A novel method has been developed for acquiring the correct alignment of a query sequence against remotely homologous proteins by extracting structural information from profiles of multiple structure alignment. A systematic search algorithm combined with a group of score functions based on sequence information and structural information has been introduced in this procedure. A limited number of top solutions (15,000) with high scores were selected as candidates for further examination. On a test-set comprising 301 proteins from 75 protein families with sequence identity less than 30%, the proportion of proteins with completely correct alignment as first candidate was improved to 39.8% by our method, whereas the typical performance of existing sequence-based alignment methods was only between 16.1% and 22.7%. Furthermore, multiple candidates for possible alignment were provided in our approach, which dramatically increased the possibility of finding correct alignment, such that completely correct alignments were found amongst the top-ranked 1000 candidates in 88.3% of the proteins. With the assistance of a sequence database, completely correct alignment solutions were achieved amongst the top 1000 candidates in 94.3% of the proteins. From such a limited number of candidates, it would become possible to identify more correct alignment using a more time-consuming but more powerful method with more detailed structural information, such as side-chain packing and energy minimization, etc. The results indicate that the novel alignment strategy could be helpful for extending the application of highly reliable methods for fold identification and homology modeling to a huge number of homologous proteins of low sequence similarity. Details of the methods, together with the results and implications for future development are presented.     

A phylogenomic approach to reconstructing the diversification of serine proteases in fungi

Using a phylogenomic approach with 10 fungi of very different virulence and habitat, we determined that there was substantial diversification of subtilase-type proteases early in ascomycete history (with subsequent loss in many lineages) but with no comparable diversification of trypsins. Patterns of intron loss and the degree of divergence between paralogues demonstrated that the proliferation of proteinase K subtilases and subtilisin type subtilases seen in pathogenic ascomycetes (Metarhizium anisopliae, Magnaporthe grisea, Fusarium graminearum) occurred after the basidiomycete/ascomycete split but predated radiation of ascomycete lineages. This suggests that the early ascomycetes had a lifestyle that selected for multiple proteases, whereas the current disparity in gene numbers between ascomycete lineages results from retention of genes in at least some pathogens that have been lost in other lineages (yeasts, Aspergillus nidulans, Neurospora crassa). A similar prevailing trend towards lineage specific gene loss of trypsins in saprophytes and some pathogens suggests that their phylogenetic breadth will have been much wider in early fungi than currently.     

Gene structure and expression of the gene from Beauveria bassiana encoding bassiasin I, an insect cuticle-degrading serine protease

Genomic and cDNA encoding Beauveria bassiana bassiasin I, a potential cuticle-degrading serine protease, were isolated and analysed. Bassiasin I gene is comprised of 1137 bp (379 aa) and 3 introns which are 69, 62 and 68 bp long. The comparison of a deduced amino acid sequence with Metarhizium anisopliae Pr1, B. bassiana Pr1, and proteinase K showed high homology. When the cDNA including the intact signal peptide was expressed in E. coli, a clear proteolytic-degraded zone on LB-skimmed milk plates was observed.     

Characterization of the complex between mannose-binding lectin trimer and mannose-binding lectin-associated serine proteases

Mannose-binding lectin (MBL) is an oligomeric serum lectin involved in innate immunity. Human MBL is complexed with three types of serine proteases (MASP-1, MASP-2 and MASP-3) and two types of their truncated forms (sMAP and MAp44). When an MBL complex binds to carbohydrates of pathogens, the complement system is activated via the lectin pathway. Human MBL is a mixture of different sized oligomers that range mainly from trimers to hexamers. It has been suggested that different MBL oligomers may have distinct MASP compositions. In the present study, an MBL trimer (MBL-I) exclusive of other oligomers was isolated from human serum by chromatography. Immunoblot analysis of MBL-I revealed that it had been co-purified with MASP-1 and sMAP. This suggests that MASP-1 and sMAP are bound to each other in MBL-I. The MBL-I complex was found to activate C2, but to lack the ability to activate C4 due to the absence of MASP-2.     

On the modeling of snake venom serine proteinase interactions with benzamidine-based thrombin inhibitors

Pit viper venoms contain a number of serine proteinases that exhibit one or more thrombin-like activities on fibrinogen and platelets, this being the case for the kinin-releasing and fibrinogen-clotting KN-BJ from the venom of Bothrops jararaca. A three-dimensional structural model of the KN-BJ2 serine proteinase was built by homology modeling using the snake venom plasminogen activator TSV-PA as a major template and porcine kallikrein as additional structural support. A set of intrinsic buried waters was included in the model and its behavior under dynamic conditions was molecular dynamics simulated, revealing a most interesting similarity pattern to kallikrein. The benzamidine-based thrombin inhibitors alpha-NAPAP, 3-TAPAP, and 4-TAPAP were docked into the refined model, allowing for a more insightful functional characterization of the enzyme and a better understanding of the reported comparatively low affinity of KN-BJ2 toward those inhibitors.     

A structural model for the prostate disease marker, human prostate-specific antigen.

Prostate-specific antigen (PSA) provides an excellent serum marker for prostate cancer, the most frequent form of cancer in American males. PSA is a 237-residue protease based on sequence homology to kallikrein-like enzymes. To predict the 3-dimensional structure of PSA, homology modeling studies were performed based on sequence and structural alignments with tonin, pancreatic kallikrein, chymotrypsin, and trypsin. The structurally conserved regions of the 4 reference X-ray proteins provided the core structure of PSA, whereas the loop structures were modeled on the loops of tonin and kallikrein. The unique "kallikrein loop" insert, between Ser 95b and Pro 95k of kallikrein, was constructed using molecular mechanics, dynamics, and electrostatics calculations. In the resulting PSA structure, the catalytic triad, involving residues His 57, Asp 102, and Ser 195, and hydrophobic and electrostatic interactions typical of serine proteases were extremely well conserved. Similarly, the 5-disulfide bonds of kallikrein were also conserved in PSA. These results, together with the fact that no major steric clashes arose during the modeling process, provide strong evidence for the validity of the PSA model. Calculation of the electrostatic potential contours of kallikrein and PSA was carried out using the finite difference Poisson-Boltzmann method. The calculations revealed matching areas of negative potential near the catalytic triad, but differences in the positive potential surrounding the active site. The PSA glycosylation site, Asn 61, is fully accessible to the solvent and is enclosed in a positive region of the isopotential map. The bottom of the substrate specificity pocket, residue S1, is a serine (Ser 189) as in chymotrypsin, rather than aspartate (Asp 189) as in tonin, kallikrein, and trypsin. This fact, plus other features of the S1 binding-pocket region, suggest that PSA would prefer substrates with hydrophobic residues at the P1 position. The location of a potential zinc ion binding site involving the side chain of histidines 91, 101, and 233 is also suggested. This PSA model should facilitate the understanding and prediction of structural and functional properties of this important cancer marker.