Comparative model-building of the mammalian serine proteases

Proteins have been classified into families based upon sequence homology. An accurate, systematic comparative model-building procedure for a homologous family of proteins would be very valuable scientifically. This paper presents such a procedure and applies it to the mammalian serine proteases, which are ubiquitous and involved in many important biological functions. Eleven proteins of this family are considered here, including a variety of blood serum, intestinal and pancreatic proteins as well as a closely related bacterial enzyme.The modeling method capitalizes upon the availability of three experimentally determined structures for mammalian serine proteases. These structures show that the molecule is divided into structurally conserved regions, which contain the strong sequence homology, and structurally variable regions, which include all the additions and deletions. We show that by applying this structural distinction to new sequences, erroneous alignments of the sequences are greatly minimized.For each aligned new sequence, the structurally conserved regions can be constructed from any of the known structures. In examining the variable regions, we have found that a variable region that has the same length and residue character in two different known structures usually has the same conformation in both. Thus, when the eight structurally unknown proteins are modeled, most of the variable regions can be constructed directly from the known structures. A minority of the variable regions require more sophisticated analysis to evaluate the relative merits of a small number of possible conformations. Only a very few are so different that modeling by homology is entirely ruled out. We demonstrate, therefore, that by this modeling procedure, the maximum of each of these mammalian serine proteases is constructed directly from the experimentally determined structures and the necessity to build from intuition or from energy considerations is greatly reduced.     

The HtrA family of serine proteases

HtrA, also known as DegP and probably identical to the Do protease, is a heat shock-induced serine protease that is active in the periplasm of Escherichia coli . Homologues of HtrA have been described in a wide range of bacteria and in eukaryotes. Its chief role is to degrade misfolded proteins in the periplasm. Substrate recognition probably involves the recently described PDZ domains in the C-terminal half of HtrA and, we suspect, has much in common with the substrate recognition system of the tail-specific protease, Prc (which also possesses a PDZ domain). The expression of htrA is regulated by a complex set of signal transduction pathways, which includes an alternative sigma factor, RpoE, an anti-sigma factor, RseA, a two-component regulatory system, CpxRA, and two phosphoprotein phosphatases, PrpA and PrpB. Mutations in the htrA genes of Salmonella , Brucella and Yersinia cause decreased survival in mice and/or macrophages, and htrA mutants can act as vaccines, as cloning hosts and as carriers of heterologous antigens.     

Granzymes: a family of lymphocyte granule serine proteases

Granzymes, a family of serine proteases, are expressed exclusively by cytotoxic T lymphocytes and natural killer (NK) cells, components of the immune system that protect higher organisms against viral infection and cellular transformation. Following receptor-mediated conjugate formation between a granzyme-containing cell and an infected or transformed target cell, granzymes enter the target cell via endocytosis and induce apoptosis. Granzyme B is the most powerful pro-apoptotic member of the granzyme family. Like caspases, cysteine proteases that play an important role in apoptosis, it can cleave proteins after acidic residues, especially aspartic acid. Other granzymes may serve additional functions, and some may not induce apoptosis. Granzymes have been well characterized only in human and rodents, and can be grouped into three subfamilies according to substrate specificity: members of the granzyme family that have enzymatic activity similar to the serine protease chymotrypsin are encoded by a gene cluster termed the 'chymase locus'; granzymes with trypsin-like specificities are encoded by the 'tryptase locus'; and a third subfamily cleaves after unbranched hydrophobic residues, especially methionine, and is encoded by the 'Met-ase locus'. All granzymes are synthesized as zymogens and, after clipping of the leader peptide, maximal enzymatic activity is achieved by removal of an amino-terminal dipeptide. They can all be blocked by serine protease inhibitors, and a new group of inhibitors has recently been identified - serpins, some of which are specific for granzymes. Future studies of serpins may bring insights into how cells that synthesize granzymes are protected from inadvertent cell suicide.     

ClpP: a distinctive family of cylindrical energy-dependent serine proteases

Processes maintaining protein homeostasis in the cell are governed by the activities of molecular chaperones that mainly assist in the folding of polypeptide chains and by a large class of proteases that regulate protein levels through degradation. ClpP proteases define a distinctive family of cylindrical, energy-dependent serine proteases that are highly conserved throughout bacteria and eukaryota. They typically interact with ATP-dependent AAA+ chaperones that bind and unfold target substrates and then translocate them into ClpP for degradation. Structural and functional studies have provided a detailed view of the mechanism of function of this class of proteases.     

Clp P represents a unique family of serine proteases

The amino acid sequence of Clp P, the proteolytic subunit of the ATP-dependent Clp protease of Escherichia coli, closely resembles a protein encoded by chloroplast DNA, which is well conserved between chloroplasts of different plant species. The homology extends over almost the full length of the sequences of both proteins and consists of approximately 46% identical and approximately 70% similar amino acids. Antibodies against E. coli Clp P cross-reacted with proteins with Mr of 20,000-30,000 in bacteria, lower eukaryotes, plants, and animal cells. Since the regulatory subunit of Clp protease, Clp A, also has a homolog in plants, as well as in other bacteria and in lower eukaryotes, it is likely that ATP-dependent proteolysis in chloroplasts is catalyzed in part by a Clp-like protease and that both components of Clp-like proteases are widespread in living cells. We have identified Ser-111 as the active site serine in E. coli Clp P modified by diisopropyl fluorophosphate. Mutational alteration of Ser-111 or His-136 eliminates proteolytic activity of Clp P. Both residues are found in highly conserved regions of the protein. The sequences around the active site residues suggest that Clp P represents a unique class of serine protease. Amino-terminal processing of cloned Clp P mutated at either Ser-111 or His-136 occurs efficiently when wild-type clpP is present in the chromosome but is blocked in clpP- hosts. Processing of Clp P appears, therefore, to involve an intermolecular autocatalytic cleavage reaction. Since processing of Clp P occurs in clpA- cells, the autoprocessing activity of Clp P is independent of Clp A.     

A despecialization step underlying evolution of a family of serine proteases

In the trypsin superfamily of serine proteases, non-trypsin-like primary specificities have arisen in only two monophyletic descendent subbranches. We have recreated an ancestor to one of these subbranches (granzyme) using phylogenetic inference, gene synthesis, and protein expression. This ancestor has two unusual properties. First, it has broad primary specificity encompassing the entire repertoire of novel primary specificities found in its descendents. Second, unlike extant members that have narrow primary specificities, the ancestor exhibits tolerance to mutational changes in primary specificity-conferring residues-that is, structural plasticity. Molecular modeling and mutagenesis studies indicate that these unusual properties are due to a particularly wide substrate binding pocket. These two crucial properties of the ancestor not only distinguish it from its extant descendents but also from the trypsin-like proteases that preceded it. This indicates that a despecialization step, characterized by broad specificity and structural plasticity, underlies evolution of new primary specificities in this protease superfamily.     

Collagenolytic activity of crustacean midgut serine proteases: comparison with the bacterial and mammalian enzymes.

1. We have investigated the collagenolytic activity of the following serine proteases: proteinase K, subtilisin Novo, Staphylococcal endoproteinase Glu-C, Streptomyces pronases, the trypsins and chymotrypsins from shrimp midgut and bovine pancreas. 2. By assays on both the insoluble 3H-collagen fibrils and the soluble type I collagen, it was demonstrated that the shrimp midgut serine proteases, and less efficiently, the pronases from Streptomyces griseus, could hydrolyze collagen while the other serine proteases tested could not. 3. Our data indicate that the trypsins and chymotrypsins of shrimp (Penaeus monodon) directly and indirectly digest native collagen, and that the indirect pathway probably involves activation of procollagenase in the native collagen by these serine proteases.     

Comparative study of the active site caging of serine proteases: thrombin and factor Xa

Bovine thrombin and human factor Xa were acylated at their active site selectively with inhibitors derived from the parent compound 4-guanidinophenyl (E)-4-diethylamino-2-hydroxy-alpha-methylcinnamate hydrochloride, 1b. Peptidyl side chains were attached to the phenol ring via amide connection, which served as a recognition motif in inhibiting different serine proteases. Upon irradiation with 366 nm light, the trans-cinnamate attached to the active-site serine isomerizes to the cis isomer which then rapidly lactonizes to release the free enzyme. The peptidyl side chain sequences specific for each serine protease were revealed via constructing and screening a library of homologous compounds. This methodology may be applied to other proteases. One application based on enzyme-specific, photoactivatable inhibitors is to isolate a designated active protease from a mixture of several proteases. Thus, a cinnamate inhibitor with a biotin moiety, 1d, was synthesized. A solution of enzyme-specific, biotinylated inhibitor was added into a mixture of proteases containing a target enzyme. The target enzyme was acylated at the active site and subsequently bore a biotin tail. An avidin column was used to separate the biotinylated enzyme from the unmodified ones, by a strong binding between biotin and avidin. After a brief irradiation on the avidin column, the retained enzymes were released from the biotin tag and eluted off the column. To demonstrate the idea, thrombin and factor Xa have been separated from each other by this strategy.     

Sol-gel immobilization of serine proteases for application in organic solvents

The serine proteases alpha-chymotrypsin, trypsin, and subtilisin Carlsberg were immobilized in a sol-gel matrix and the effects on the enzyme activity in organic media are evaluated. The percentage of immobilized enzyme is 90% in the case of alpha-chymotrypsin and the resulting specific enzyme activity in the transesterification of N-acetyl-L-phenylalanine ethyl ester with 1-propanol in cyclohexane is 43 times higher than that of a nonimmobilized lyophilized alpha-chymotrypsin. The activities of trypsin and subtilisin Carlsberg are enhanced with 437 and 31 times, respectively. The effect of immobilization on the enzyme activity is highest in hydrophobic solvents.     

Comparative modeling of mammalian aspartate transcarbamylase

Mammalian aspartate transcarbamylase (ATCase) is part of a 243 kDa multidomain polypeptide, called CAD, that catalyzes the first three steps in de novo pyrimidine biosynthesis. The structural organization of the mammalian enzyme is very different from E. coli ATCase, a dodecameric, monofunctional molecule comprised of six copies of separate catalytic and regulatory chains. Nevertheless, sequence similarities and other properties suggested that the mammalian ATCase domain and the E. coli ATCase catalytic chain have the same tertiary fold. A model of mammalian ATCase was built using the X-ray coordinates of the E. coli catalytic chain as a tertiary template. Five small insertions and deletions could be readily accommodated in the model structure. Following energy minimization the RMS difference in the alpha carbon positions of the mammalian and bacterial proteins was 0.93 A. A comparison of the hydrophobic energies, surface accessibility index, and the distribution of hydrophilic and hydrophobic residues of the CAD ATCase structure with correctly and incorrectly folded proteins and with several X-ray structures supported the validity of the model. The mammalian ATCase domain associates to form a compact globular trimer, a prerequisite for catalysis since the active site is comprised of residues from adjacent subunits. Interactions between the clearly defined aspartate and carbamyl phosphate subdomains of the monomer were largely preserved while there was appreciable remodeling of the trimeric interfaces. Several clusters of basic residues are located on the upper surface of the domain which account in part for the elevated isoelectric point (pI = 9.4) and may represent contact regions with other more acidic domains within the chimeric polypeptide. A long interdomain linker connects the monomer at its upper surface to the remainder of the polypeptide. The configuration of active site residues is virtually identical in the mammalian and bacterial enzymes. While the CAD ATCase domain can undergo the local conformational changes that accompany catalysis in the E. coli enzyme, the high activity, closed conformation is probably more stable in the mammalian enzyme.     

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.     

A novel protease from Entamoeba histolytica homologous to members of the family S28 of serine proteases

Serine proteases are one of the biologically most important and widely distributed enzyme families. A protease capable of degrading the substrate Suc-AAF-AMC was isolated from axenically grown trophozoites of Entamoeba histolytica. The enzyme was purified by ion-exchange chromatography and electroelution, and appeared on 2D-PAGE as a spot of 60 kDa and pI of 4.65. Data obtained from zymogram suggest the active protease is present either as homodimer (130 kDa) or homotetramer (250 kDa). The optimal temperature of the enzyme was 37 degrees C, and it exhibited activity over a broad pH range. The protease was strongly inhibited by TPCK and chelating agents. The enzymatic activity was restored upon addition of calcium. BLAST analysis with the sequence of internal peptides of the protein revealed two open reading frames within the genome of E. histolytica, homologous to members of the family S28, clan SC of serine proteases.     

Comparative study of the properties of serine proteases of lower and higher vertebrates

Results of the comparative study of trypsin- and chymotrypsin-like serine proteases from pyloric caeca of salmon fishes and trypsin and chymotrypsin of bulls are presented in the paper. The hydrolytic activity of salmon proteases with respect to methyl ethers of N-benzoyl-L-leucine is 2.4 times higher than that of bull chymotrypsin, but with respect to methyl esters of N-benzoyl-L-tyrosine and N-benzoyl-L-arginine the activity of salmon proteases is 6.5 and 80 times lower than that of bull chymotrypsin and trypsin. Salmon proteases in contrast to bull trypsin and chymotrypsin hydrolyze but slightly N-glutaryl-L-phenylalanine para-nitroanilide. It shown that fish proteases are not absolutely specific to synthetic substrates, which is a result of their less pronounced (than in case of bull trypsin and chymotrypsin) differences in structures of binding centres. The study of the salmon protease interaction with some immobilized ligands has confirmed the higher affinity of enzymes to reagents with two space-separated aromatic rings in their composition. It is supposed that salmon proteases interact with such reagents through two sites: hydrophobic "pockets" and probably additional binding site of the active centre. The salmon protease preparation demonstrates higher resistance to inactivating action of formaldehyde within the range of concentrations 2-16% than bull chymotrypsin does.     

The mycosins of Mycobacterium tuberculosis H37Rv: a family of subtilisin-like serine proteases

There is little information regarding the role of proteolysis in Mycobacterium tuberculosis and no studies on the potential involvement of proteases in the pathogenesis of tuberculosis. We identified five M. tuberculosis genes (mycP1-5) that encode a family of serine proteases (mycosins-1 to 5), ranging from 36 to 47% identity. Each protein contains a catalytic triad (Asp, His, Ser) within highly conserved sequences, typical of proteases of the subtilisin family. These genes are also present in M. bovis BCG and other virulent mycobacteria, but only one homologue (mycP3) was detected in M. smegmatis. The mycosins have N-terminal signal sequences and C-terminal transmembrane anchors, and the localisation of the mycosins to the membrane/cell wall was verified by Western blot analysis of heterologously expressed proteins in cellular fractions of M. smegmatis. In M. tuberculosis, all the mycosins were expressed constitutively during growth in broth. Mycosins-2 and 3 were also expressed constitutively in M. bovis BCG, but no expression of mycosin-1 was detected. Mycosin-2 was modified by cleavage in all three mycobacterial species. The multiplicity and constitutive expression of these proteins suggests that they have an important role in the biology of M. tuberculosis.     

Drosophila rhomboid-1 defines a family of putative intramembrane serine proteases.

The polytopic membrane protein Rhomboid-1 promotes the cleavage of the membrane-anchored TGFalpha-like growth factor Spitz, allowing it to activate the Drosophila EGF receptor. Until now, the mechanism of this key signaling regulator has been obscure, but our analysis suggests that Rhomboid-1 is a novel intramembrane serine protease that directly cleaves Spitz. In accordance with the putative Rhomboid active site being in the membrane bilayer, Spitz is cleaved within its transmembrane domain, and thus is, to our knowledge, the first example of a growth factor activated by regulated intramembrane proteolysis. Rhomboid-1 is conserved throughout evolution from archaea to humans, and our results show that a human Rhomboid promotes Spitz cleavage by a similar mechanism. This growth factor activation mechanism may therefore be widespread.     

RNA-based phylogenetic methods: application to mammalian mitochondrial RNA sequences

The PHASE software package allows phylogenetic tree construction with a number of evolutionary models designed specifically for use with RNA sequences that have conserved secondary structure. Evolution in the paired regions of RNAs occurs via compensatory substitutions, hence changes on either side of a pair are correlated. Accounting for this correlation is important for phylogenetic inference because it affects the likelihood calculation. In the present study we use the complete set of tRNA and rRNA sequences from 69 complete mammalian mitochondrial genomes. The likelihood calculation uses two evolutionary models simultaneously for different parts of the sequence: a paired-site model for the paired sites and a single-site model for the unpaired sites. We use Bayesian phylogenetic methods and a Markov chain Monte Carlo algorithm is used to obtain the most probable trees and posterior probabilities of clades. The results are well resolved for almost all the important branches on the mammalian tree. They support the arrangement of mammalian orders within the four supra-ordinal clades that have been identified by studies of much larger data sets mainly comprising nuclear genes. Groups such as the hedgehogs and the murid rodents, which have been problematic in previous studies with mitochondrial proteins, appear in their expected position with the other members of their order. Our choice of genes and evolutionary model appears to be more reliable and less subject to biases caused by variation in base composition than previous studies with mitochondrial genomes.     

Subtilases: the superfamily of subtilisin-like serine proteases.

Subtilases are members of the clan (or superfamily) of subtilisin-like serine proteases. Over 200 subtilases are presently known, more than 170 of which with their complete amino acid sequence. In this update of our previous overview (Siezen RJ, de Vos WM, Leunissen JAM, Dijkstra BW, 1991, Protein Eng 4:719-731), details of more than 100 new subtilases discovered in the past five years are summarized, and amino acid sequences of their catalytic domains are compared in a multiple sequence alignment. Based on sequence homology, a subdivision into six families is proposed. Highly conserved residues of the catalytic domain are identified, as are large or unusual deletions and insertions. Predictions have been updated for Ca(2+)-binding sites, disulfide bonds, and substrate specificity, based on both sequence alignment and three-dimensional homology modeling.