Search

Proteases play a fundamental role in the control of intra- and extracellular processes by binding and cleaving specific amino acid sequences. Identifying these targets is extremely challenging. Current computational attempts to predict cleavage sites are limited, representing these amino acid sequences as patterns or frequency matrices. Here we present PoPS, a publicly accessible bioinformatics tool (http://pops.csse.monash.edu.au/) which provides a novel method for building computational models of protease specificity that, while still being based on these amino acid sequences, can be built from any experimental data or expert knowledge available to the user. PoPS specificity models can be used to predict and rank likely cleavages within a single substrate, and within entire proteomes. Other factors, such as the secondary or tertiary structure of the substrate, can be used to screen unlikely sites. Furthermore, the tool also provides facilities to infer, compare and test models, and to store them in a publicly accessible database.

The classical complement pathway, which plays a vital role in preventing infection, is initiated by the action of the serine proteases C1r and C1s. We have examined the hydrolysis of substrates representing cleavage sequences in the physiological substrates for C1s, C2 and C4. These studies showed that the P(1)'-P(4)' substrate residues of C2 and C4 conferred greater affinity of substrate for enzyme and also induced a sigmoidal dependence of enzyme velocity on substrate concentration. This indicates that the substrate gave rise to homotropic positive cooperative behavior in the enzyme. When C1s was in complex with C1q and C1r, as would occur under physiological conditions, the same behavior was observed, indicating that this mechanism is relevant in the complement pathway in vivo. We further investigated the requirements of C1s for prime side amino acids by examining a substrate library in which each of the P(1)'-P(4)' positions had been substituted by different classes of amino acids. This revealed that the P(1)' position was a major determinant of the selectivity of the enzyme, while certain substitutions at the P(1)'-P(4)' positions abolished the allosteric behavior, indicating that contact residues at these positions in the C1s enzyme must mediate the cooperativity. The studies reported here highlight the importance of prime subsites in C1s for interaction with its cognate substrates in the complement pathway and therefore yield greater understanding of the mechanism of interaction between this vital protease and its physiological substrates.

MENT is a developmentally regulated heterochromatin-associated protein that condenses chromatin in terminally differentiated avian blood cells. Its homology to the serpin protein family suggests that the conserved serpin reactive center loop (RCL) and the unique M-loop are important for its function. To examine the role of these domains, we studied the interaction of wild-type and mutant MENT with naked DNA and biochemically defined nucleosome arrays reconstituted from 12-mer repeats containing nucleosome positioning sequences. Wild-type MENT folded the naked DNA duplexes into closely juxtaposed parallel structures ("tramlines"). Deletion of the M-loop, but not inactivation of the RCL, prevented tramline formation and the cooperative interaction of MENT with DNA. Reconstitution of wild-type MENT with nucleosome arrays caused their tight folding and self-association. M-loop deletion inhibited nucleosome array folding, whereas the inactive RCL mutant was competent to fold the nucleosome arrays, but had a significantly impaired ability to cause their self-association. Bifunctional chemical cross-linking of MENT revealed oligomerization of wild-type MENT in the presence of chromatin and DNA. This oligomerization was severely reduced in the RCL mutant. We propose that the mechanism of MENT-induced heterochromatin formation involves two independent events: bringing together nucleosome linkers within a chromatin fiber and formation of protein bridges between chromatin fibers. Ordered binding of MENT to linker DNA via its unique M-loop domain promotes the folding of chromatin, whereas bridging of chromatin fibers is facilitated by MENT oligomerization mediated by the RCL.

Porphyromonas gingivalis is a pathogen associated with periodontal disease, and arginine-specific proteases (gingipains-R) from the bacterium are important virulence factors. The specificity of two forms of gingipain-R, HRgpA and RgpB, for substrate positions C-terminal to the cleavage site was analyzed, and notable differences were observed between the enzymes. Molecular modeling of the HRgpA catalytic domain, based on the structure of RgpB, revealed that there are four amino acid substitutions around the active site of HRgpA relative to RgpB that may explain their different specificity. Previously, differences in the ability of these two gingipain-R forms to cleave a number of proteins were attributed to additional adhesins on HRgpA mediating increased interaction with the substrates. Here, purified RgpA(cat), the catalytic domain of HRgpA, which like RgpB also lacks adhesin subunits, was used to show that the differences between HRgpA and RgpB are probably due to the amino acid substitutions at the active site. The kinetics of cleavage of fibrinogen, a typical protein substrate for the gingipain-R enzymes, which is bound by HRgpA but not RgpA(cat) or RgpB, were evaluated, and it was shown that there was no difference in the cleavage of the fibrinogen Aalpha-chain between the different enzyme forms. HRgpA degraded the fibrinogen Bbeta-chain more efficiently, generating distinct cleavage products. This indicates that while the adhesin domain(s) play(s) a minor role in the cleavage of protein substrates, the major effect is still provided by the amino acid substitutions at the active site of rgpA gene products versus those of the rgpB gene.

The sequence of a genome contains the plans of the possible life of an organism, but implementation of genetic information depends on the functions of the proteins and nucleic acids that it encodes. Many individual proteins of known sequence and structure present challenges to the understanding of their function. In particular, a number of genes responsible for diseases have been identified but their specific functions are unknown. Whole-genome sequencing projects are a major source of proteins of unknown function. Annotation of a genome involves assignment of functions to gene products, in most cases on the basis of amino-acid sequence alone. 3D structure can aid the assignment of function, motivating the challenge of structural genomics projects to make structural information available for novel uncharacterized proteins. Structure-based identification of homologues often succeeds where sequence-alone-based methods fail, because in many cases evolution retains the folding pattern long after sequence similarity becomes undetectable. Nevertheless, prediction of protein function from sequence and structure is a difficult problem, because homologous proteins often have different functions. Many methods of function prediction rely on identifying similarity in sequence and/or structure between a protein of unknown function and one or more well-understood proteins. Alternative methods include inferring conservation patterns in members of a functionally uncharacterized family for which many sequences and structures are known. However, these inferences are tenuous. Such methods provide reasonable guesses at function, but are far from foolproof. It is therefore fortunate that the development of whole-organism approaches and comparative genomics permits other approaches to function prediction when the data are available. These include the use of protein-protein interaction patterns, and correlations between occurrences of related proteins in different organisms, as indicators of functional properties. Even if it is possible to ascribe a particular function to a gene product, the protein may have multiple functions. A fundamental problem is that function is in many cases an ill-defined concept. In this article we review the state of the art in function prediction and describe some of the underlying difficulties and successes.

Factor Xa is a central protease in the coagulation cascade and the target for many anticoagulant compounds currently under development. The preferences of the enzyme for substrates incorporating residues N-terminal to the cleavage site (P1, P2, etc.) have been elucidated, but little is known of its preferences for residues C-terminal to the cleavage site (P1', P2', etc.). The preferences of bovine factor Xa for substrate residues in the P1', P2' and P3' positions were mapped using fluorescence-quenched substrates. Bovine factor Xa, often used as a model for factor Xa, was most selective for the P2' position, less selective at the P1' position and almost completely non-selective at the P3' position. It appears that while the prime side subsites of factor Xa impose some selectivity towards substrates, the influence of these sites on factor Xa cleavage specificity is relatively low in comparison to related enzymes such as thrombin.

Members of the serpin (serine proteinase inhibitor) superfamily have been identified in higher multicellular eukaryotes (plants and animals) and viruses but not in bacteria, archaea, or fungi. Thus, the ancestral serpin and the origin of the serpin inhibitory mechanism remain obscure. In this study we characterize 12 serpin-like sequences in the genomes of prokaryotic organisms, extending this protein family to all major branches of life. Notably, these organisms live in dramatically different environments and some are evolutionarily distantly related. A sequence-based analysis suggests that all 12 serpins are inhibitory. Despite considerable sequence divergence between the proteins, in four of the 12 sequences the region of the serpin that determines proteinase specificity is highly conserved, indicating that these inhibitors are likely to share a common target. Inhibitory serpins are typically prone to polymerization upon heating; thus, the existence of serpins in the moderate thermophilic bacterium Thermobifida fusca, the thermophilic bacterium Thermoanaerobacter tengcongensis, and the hyperthermophilic archaeon Pyrobaculum aerophilum is of particular interest. Using molecular modeling, we predict the means by which heat stability in the latter protein may be achieved without compromising inhibitory activity.

The control of coagulation enzymes by antithrombin is vital for maintenance of normal hemostasis. Antithrombin requires the co-factor, heparin, to efficiently inhibit target proteinases. A specific pentasaccharide sequence (H5) in high affinity heparin induces a conformational change in antithrombin that is particularly important for factor Xa (fXa) inhibition. Thus, synthetic H5 accelerates the interaction between antithrombin and fXa 100-fold as compared with only 2-fold versus thrombin. We built molecular models and identified residues unique to the active site of fXa that we predicted were important for interacting with the reactive center loop of H5-activated antithrombin. To test our predictions, we generated the mutants E37A, E37Q, E39A, E39Q, Q61A, S173A, and F174A in human fXa and examined the rate of association of these mutants with antithrombin in the presence and absence of H5. fXa(Q61A) interacts with antithrombin alone with a nearly normal k(ass); however, we observe only a 4-fold increase in k(ass) in the presence of H5. The x-ray crystal structure of fXa reveals that Gln(61) forms part of the S1' and S3' pocket, suggesting that the P' region of the reactive center loop of antithrombin is crucial for mediating the acceleration in the rate of inhibition of fXa by H5-activated antithrombin.

MENT (Myeloid and Erythroid Nuclear Termination stage-specific protein) is a developmentally regulated chromosomal serpin that condenses chromatin in terminally differentiated avian blood cells. We show that MENT is an effective inhibitor of the papain-like cysteine proteinases cathepsins L and V. In addition, ectopic expression of MENT in mammalian cells is apparently sufficient to inhibit a nuclear papain-like cysteine proteinase and prevent degradation of the retinoblastoma protein, a major regulator of cell proliferation. MENT also accumulates in the nucleus, causes a strong block in proliferation, and promotes condensation of chromatin. Variants of MENT with mutations or deletions within the M-loop, which contains a nuclear localization signal and an AT-hook motif, reveal that this region mediates nuclear transport and morphological changes associated with chromatin condensation. Non-inhibitory mutants of MENT were constructed to determine whether its inhibitory activity has a role in blocking proliferation. These mutations changed the mode of association with chromatin and relieved the block in proliferation, without preventing transport to the nucleus. We conclude that the repressive effect of MENT on chromatin is mediated by its direct interaction with a nuclear protein that has a papain-like cysteine proteinase active site.

Serpins form loop-sheet polymers through the formation of a partially folded intermediate. Through mutagenesis and biophysical analysis, we have probed the conformational stability of the F-helix, demonstrating that it is almost completely unfolded in the intermediate state. The replacement of Tyr160 on the F-helix of alpha1-antitrypsin to alanine results in the loss of a conserved hydrogen bond that dramatically reduces the stability of the protein to both heat and solvent denaturation, indicating the importance of Tyr160 in the stability of the molecule. The mutation of Tyr160 to a tryptophan residue, within a fluorescently silent variant of alpha1-antitrypsin, results in a fully active, stable serpin. Fluorescence analysis of the equilibrium unfolding behavior of this variant indicates that the F-helix is highly disrupted in the intermediate conformation. Iodide quenching experiments demonstrate that the tryptophan residue is exposed to a similar extent in both the intermediate and unfolded states. Cumulatively, these data indicate that the F-helix plays an important role in controlling the early conformational changes involved in alpha1-antitrypsin unfolding. The implications of these data on both alpha1-antitrypsin function and misfolding are discussed.

The N-terminal one-third of the NS3 protein of Dengue virus type 2 (DEN-2) complexes with co-factor NS2B to form an active serine proteinase which cleaves the viral polyprotein. To identify sites within NS3 that may interact with NS2B, seven regions within the NS3 proteinase outside the conserved flavivirus enzyme motifs were mutated by alanine replacement. Five sites contained clusters of charged residues and were hydrophilic. Two sites were hydrophobic and highly conserved among flaviviruses. The effects of five mutations on NS2B/3 processing were examined using a COS cell expression system. Four retained significant proteinase activity. Three of these mutations and two more were introduced into genomic-length cDNA and tested for their effects on virus replication. The five mutant viruses showed reduced plaque size and two of the five showed significantly reduced titres. All seven mutations were mapped on the X-ray crystal structure of the DEN-2 NS3 proteinase: three were located at the N terminus and two at the C terminus of the NS2B-binding cleft. Two mutations were at the C terminus of the proteinase domain and one was solvent-exposed. The study demonstrated that charged-to-alanine mutagenesis in the viral proteinase can be used to produce growth-restricted flaviviruses that may be useful in the production of attenuated vaccine strains.

Myxoma virus (MV), a member of the family Poxviridae, is the causative agent of myxomatosis, a fatal disease of the European rabbit. The MV genome is a linear, double-stranded DNA molecule that encodes several factors important for evasion of the host immune system. Sequencing the right-end region of the MV genome identified an 801 bp open reading frame (ORF) encoding a polypeptide that belongs to the serpin superfamily. To date, two MV-encoded serpins have been characterized: SERP-1 binds to several targets and is an anti-inflammatory molecule, whereas Serp2 is essential for virus virulence and has both anti-inflammatory and anti-apoptotic effects. Thus, Serp3 is the third MV-encoded serpin. DNA sequence analysis of Serp3 indicated a similarity to poxvirus late promoters, which was confirmed by mRNA expression analysis. Serp3 has an atypical serpin motif and has significant sequence deletions as compared to most cellular and viral serpins. However, molecular modelling studies suggested that Serp3 can retain the overall serpin fold. Insertional inactivation of the serp3 ORF led to a significant attenuation of virulence in vivo (as measured by the increase in survival of infected rabbits) and limited dissemination of the virus to secondary sites of infection. In rabbits infected with a Serp3 deletion mutant (MV-Serp3(-)), the main histopathological feature is the absence of secondary myxomas. Both wild-type MV and MV-Serp3(-) replicate at comparable levels in vivo. Serp3 may represent a significant virulence factor of MV and probably acts in synergy with other viral proteins.