Electrostatic recognition in substrate binding to serine proteases

Serine proteases of the Chymotrypsin family are structurally very similar but have very different substrate preferences. This study investigates a set of 9 different proteases of this family comprising proteases that prefer substrates containing positively charged amino acids, negatively charged amino acids, and uncharged amino acids with varying degree of specificity. Here, we show that differences in electrostatic substrate preferences can be predicted reliably by electrostatic molecular interaction fields employing customized GRID probes. Thus, we are able to directly link protease structures to their electrostatic substrate preferences. Additionally, we present a new metric that measures similarities in substrate preferences focusing only on electrostatics. It efficiently compares these electrostatic substrate preferences between different proteases. This new metric can be interpreted as the electrostatic part of our previously developed substrate similarity metric. Consequently, we suggest, that substrate recognition in terms of electrostatics and shape complementarity are rather orthogonal aspects of substrate recognition. This is in line with a 2‐step mechanism of protein‐protein recognition suggested in the literature.     

Discovery of local packing motifs in protein structures

We present a language for describing structural patterns of residues in protein structures and a method for the discovery of such patterns that recur in a set of protein structures. The patterns impose restrictions on the spatial position of each residue, their order along the amino acid chain, and which amino acids are allowed in each position. Unlike other methods for comparing sets of protein structures, our method is not based on the use of pairwise structure comparisons which is often time consuming and can produce inconsistent results. Instead, the method simultaneously takes into account information from all structures in the search for conserved structure patterns which are potential structure motifs. The method is based on describing the spatial neighborhoods of each residue in each structure as a string and applying a sequence pattern discovery method to find patterns common to subsets of these strings. Finally it is checked whether the similarities between the neighborhood strings correspond to spatially similar substructures. We apply the method to analyze sets of very disparate proteins from the four different protein families: serine proteases, cuprodoxins, cysteine proteinases, and ferredoxins. The motifs found by the method correspond well to the site and motif information given in the annotation of these proteins in PDB, Swiss-Prot, and PROSITE. Furthermore, the motifs are confirmed by using the motif data to constrain the structural alignment of the proteins obtained with the program SAP. This gave the best superposition/alignment of the proteins given the motif assignment.     

Molecular cloning of cDNA for matriptase, a matrix-degrading serine protease with trypsin-like activity.

A major protease from human breast cancer cells was previously detected by gelatin zymography and proposed to play a role in breast cancer invasion and metastasis. To structurally characterize the enzyme, we isolated a cDNA encoding the protease. Analysis of the cDNA reveals three sequence motifs: a carboxyl-terminal region with similarity to the trypsin-like serine proteases, four tandem cysteine-rich repeats homologous to the low density lipoprotein receptor, and two copies of tandem repeats originally found in the complement subcomponents C1r and C1s. By comparison with other serine proteases, the active-site triad was identified as His-484, Asp-539, and Ser-633. The protease contains a characteristic Arg-Val-Val-Gly-Gly motif that may serve as a proteolytic activation site. The bottom of the substrate specificity pocket was identified to be Asp-627 by comparison with other trypsin-like serine proteases. In addition, this protease exhibits trypsin-like activity as defined by cleavage of synthetic substrates with Arg or Lys as the P1 site. Thus, the protease is a mosaic protein with broad spectrum cleavage activity and two potential regulatory modules. Given its ability to degrade extracellular matrix and its trypsin-like activity, the name matriptase is proposed for the protease.     

aspS encoding an unusual aspartyl protease from Sclerotinia sclerotiorum is expressed during phytopathogenesis

The gene aspS encoding an aspartyl protease has been cloned from Sclerotinia sclerotiorum by screening a genomic library with a PCR-amplified fragment of the gene. The open reading frame of 1368 bp interrupted by one intron would encode a preproprotein of 435 amino acids. The catalytic aspartyl residues characteristic of aspartyl proteases are conserved; however, the active-site motif (DSG) in the N-terminal lobe is unusual in that Ser replaced Thr used in the active-site motif (DTG) of the C-terminal lobe and in all other fungal aspartyl proteases. RT-PCR revealed that aspS expression in axenic culture is not subjected to catabolite repression and demonstrated that aspS is expressed from the beginning of infection of sunflower cotyledons.     

Variable context Markov chains for HIV protease cleavage site prediction

Deciphering the knowledge of HIV protease specificity and developing computational tools for detecting its cleavage sites in protein polypeptide chain are very desirable for designing efficient and specific chemical inhibitors to prevent acquired immunodeficiency syndrome. In this study, we developed a generative model based on a generalization of variable order Markov chains (VOMC) for peptide sequences and adapted the model for prediction of their cleavability by certain proteases. The new method, called variable context Markov chains (VCMC), attempts to identify the context equivalence based on the evolutionary similarities between individual amino acids. It was applied for HIV-1 protease cleavage site prediction problem and shown to outperform existing methods in terms of prediction accuracy on a common dataset. In general, the method is a promising tool for prediction of cleavage sites of all proteases and encouraged to be used for any kind of peptide classification problem as well.     

Proteolytic processing of Ty3 proteins is required for transposition.

Ty3 is a retroviruslike element found in Saccharomyces cerevisiae. It encodes GAG3 and GAG3-POL3 polyproteins which are processed into mature proteins found in the Ty3 viruslike particle. In this study, the region encoding a protease that is homologous to retroviral aspartyl proteases was identified and shown to be required for production of mature Ty3 proteins and transposition. The Ty3 protease has the Asp-Ser-Gly consensus sequence found in copia, Ty1, and Rous sarcoma virus proteases, rather than the Asp-Thr-Gly found in most retroviral proteases. The Asp-Ser-Gly consensus is flanked by residues similar to those which flank the active sites of cellular aspartyl proteases. Mutations were made in the Ty3 active-site sequence to examine the role of the protease in Ty3 particle maturation and to test the functional significance of the Ser active-site variant in the consensus sequence. Mutation of the active-site Asp blocked processing of Gag3 and Gag3-Pol3 and allowed identification of a GAG3-POL3 polyprotein. This protein was turned over rapidly in cells expressing the mutant Ty3. Changing the active-site Ser to Thr caused only a modest reduction in the levels of certain Ty3 proteins. Five putative cleavage sites of this protease in Ty3 GAG3 and GAG3-POL3 polyproteins were defined by amino-terminal sequence analysis. The existence of an additional protein(s) of unknown function, encoded downstream of the protease-coding region, was deduced from the positions of these amino termini and the sizes of known Ty3 proteins. Although Ty3 protease cleavage sites do not correspond exactly to known retroviral protease cleavage sites, there are similarities. Residues P3 through P2' in the regions encompassing each of the five sites are uncharged, and no P1 position is occupied by an amino acid with a branched beta carbon.     

Functional Evolution of Proteins

The functional evolution of proteins advances through gene duplication followed by functional drift, whereas molecular evolution occurs through random mutational events. Over time, protein active-site structures or functional epitopes remain highly conserved, which enables relationships to be inferred between distant orthologs or paralogs. In this study, we present the first functional clustering and evolutionary analysis of the RCSB Protein Data Bank (RCSB PDB) based on similarities between active-site structures. All of the ligand-bound proteins within the RCSB PDB were scored using our Comparison of Protein Active-site Structures (CPASS) software and database ( http://cpass.unl.edu/ ). Principal component analysis was then used to identify 4431 representative structures to construct a phylogenetic tree based on the CPASS comparative scores ( http://itol.embl.de/shared/jcatazaro ). The resulting phylogenetic tree identified a sequential, step-wise evolution of protein active-sites and provides novel insights into the emergence of protein function or changes in substrate specificity based on subtle changes in geometry and amino acid composition.     

Flexible structural comparison allowing hinge-bending, swiveling motions

We present an efficient method for flexible comparison of protein structures, allowing swiveling motions. In all currently available methodologies developed and applied to the comparisons of protein structures, the molecules are considered to be rigid objects. The method described here extends and generalizes current approaches to searches for structural similarity between molecules by viewing proteins as objects consisting of rigid parts connected by rotary joints. During the matching, the rigid subparts are allowed to be rotated with respect to each other around swiveling points in one of the molecules. This technique straightforwardly detects structural motifs having hinge(s) between their domains. Whereas other existing methods detect hinge-bent motifs by initially finding the matching rigid parts and subsequently merging these together, our method automatically detects recurring substructures, allowing full 3 dimensional rotations about their swiveling points. Yet the method is extremely fast, avoiding the time-consuming full conformational space search. Comparison of two protein structures, without a predefinition of the motif, takes only seconds to one minute on a workstation per hinge. Hence, the molecule can be scanned for many potential hinge sites, allowing practically all C(alpha) atoms to be tried as swiveling points. This algorithm provides a highly efficient, fully automated tool. Its complexity is only O(n2), where n is the number of C(alpha) atoms in the compared molecules. As in our previous methodologies, the matching is independent of the order of the amino acids in the polypeptide chain. Here we illustrate the performance of this highly powerful tool on a large number of proteins exhibiting hinge-bending domain movements. Despite the motions, known hinge-bent domains/motifs which have been assembled and classified, are correctly identified. Additional matches are detected as well. This approach has been motivated by a technique for model based recognition of articulated objects originating in computer vision and robotics.     

RNA-binding protein-related sequence in a malaria antigen, clustered-asparagine-rich protein

Members of the RNA-binding protein superfamily contain RNA binding domains of about 90 amino acids with a highly conserved motif 'GFGF'. Using the conserved motif with some variations G-(F/Y)-(G/A)-(F/Y)-(V/I)-X-(F/Y) as a probe, we screened protein sequences carrying identical amino acids in an NBRF-protein database. It has been found that the C-terminal portion of clustered asparagine-rich protein (CARP), a malaria antigen from Plasmodium falciparum, shows an unexpected sequence similarity with the RNA-binding protein superfamily for the C-terminal half of the RNA-binding domain. Dot matrix comparisons and alignment of these sequences as well as a statistical test have revealed highly significant sequence similarities. From these analyses, we conclude that the malaria antigen CARP belongs to a large family of the RNA-binding proteins. An evolutionary implication of the sequence similarity was also discussed.     

The primary structure of staphylococcal protease.

The amino acid sequence of staphylococcal protease has been determined by analysis of tryptic peptides obtained from cyanogen bromide fragments. Selected peptides obtained from digests with staphylococcal protease, thermolysin, and chymotrypsin provided the information necessary to align the tryptic peptides and the cyanogen bromide fragments. The protease is a single polypeptide chain of some 250 amino acids and is devoid of sulfhydryl groups. The COOH-terminal tryptic peptide of of the protease molecule contains some 43 residues, most of which are aspartic acids, asparagines, and prolines. The amino acid sequence of this peptide was not determined. The primary structure near the active serine residue indicates that staphylococcal protease is related to the pancreatic serine proteases. However, it has little or no additional sequence homologies with these enzymes except for the regions near histidine-50 and aspartic acid - 91. These regions have striking similarities with the corresponding regions of protease B and the trypsin-like enzyme of Streptomyces griseus.     

A protein short motif search tool using amino acid sequence and their secondary structure assignment

We present the development of a web server, a protein short motif search tool that allows users to simultaneously search for a protein sequence motif and its secondary structure assignments. The web server is able to query very short motifs searches against PDB structural data from the RCSB Protein Databank, with the users defining the type of secondary structures of the amino acids in the sequence motif. The output utilises 3D visualisation ability that highlights the position of the motif in the structure and on the corresponding sequence. Researchers can easily observe the locations and conformation of multiple motifs among the results. Protein short motif search also has an application programming interface (API) for interfacing with other bioinformatics tools. AVAILABILITY: The database is available for free at http://birg3.fbb.utm.my/proteinsms.     

A test of proposed rules for helix capping: implications for protein design

alpha-helices within proteins are often terminated (capped) by distinctive configurations of the polypeptide chain. Two common arrangements are the Schellman motif and the alternative alpha(L) motif. Rose and coworkers developed stereochemical rules to identify the locations of such motifs in proteins of unknown structure based only on their amino acid sequences. To check the effectiveness of these rules, they made specific predictions regarding the structural and thermodynamic consequences of certain mutations in T4 lysozyme. We have constructed these mutants and show here that they have neither the structure nor the stability that was predicted. The results show the complexity of the protein-folding problem. Comparison of known protein structures may show that a characteristic sequence of amino acids (a sequence motif) corresponds to a conserved structural motif. In any particular protein, however, changes in other parts of the sequence may result in a different conformation. The structure is determined by sequence as a whole, not by parts considered in isolation.     

Novel protein folds and their nonsequential structural analogs

Newly determined protein structures are classified to belong to a new fold, if the structures are sufficiently dissimilar from all other so far known protein structures. To analyze structural similarities of proteins, structure alignment tools are used. We demonstrate that the usage of nonsequential structure alignment tools, which neglect the polypeptide chain connectivity, can yield structure alignments with significant similarities between proteins of known three-dimensional structure and newly determined protein structures that possess a new fold. The recently introduced protein structure alignment tool, GANGSTA, is specialized to perform nonsequential alignments with proper assignment of the secondary structure types by focusing on helices and strands only. In the new version, GANGSTA+, the underlying algorithms were completely redesigned, yielding enhanced quality of structure alignments, offering alignment against a larger database of protein structures, and being more efficient. We applied DaliLite, TM-align, and GANGSTA+ on three protein crystal structures considered to be novel folds. Applying GANGSTA+ to these novel folds, we find proteins in the ASTRAL40 database, which possess significant structural similarities, albeit the alignments are nonsequential and in some cases involve secondary structure elements aligned in reverse orientation. A web server is available at http://agknapp.chemie.fu-berlin.de/gplus for pairwise alignment, visualization, and database comparison.     

A computer model to dynamically simulate protein folding: studies with crambin

The current work describes a simplified representation of protein structure with uses in the simulation of protein folding. The model assumes that a protein can be represented by a freely rotating rigid chain with a single atom approximating the effect of each side chain. Potentials describing the attraction or repulsion between different types of amino acids are determined directly from the distribution of amino acids in the database of known protein structures. The optimization technique of simulated annealing has been used to dynamically sample the conformations available to this simple model, allowing the protein to evolve from an extended, random coil into a compact globular structure. Many characteristics expected of true proteins, such as the sequence-dependent formation of secondary structure, the partitioning of hydrophobic residues, and specific disulfide pairing, are reproduced by the simulation, suggesting the model may accurately simulate the folding process.     

Targeting of the human adrenoleukodystrophy protein to the peroxisomal membrane by an internal region containing a highly conserved motif

In this study we addressed the targeting requirements of peroxisomal ABC transporters, in particular the human adrenoleukodystrophy protein. This membrane protein is defective or missing in X-linked adrenoleukodystrophy, a neurodegenerative disorder predominantly presenting in childhood. Using adrenoleukodystrophy protein deletion constructs and green fluorescent protein fusion constructs we identified the amino acid regions 1-110 and 67-164 to be sufficient for peroxisomal targeting. However, the minimal region shared by these constructs (amino acids 67-110) is not sufficient for peroxisomal targeting by itself. Additionally, the NH2-terminal 66 amino acids enhance targeting efficiency. Green fluorescent protein-labeled fragments of human peroxisomal membrane protein 69 and Saccharomyces cerevisiae Pxa1 corresponding to the amino acid 67-164 adrenoleukodystrophy protein region were also directed to the mammalian peroxisome. The required region contains a 14-amino-acid motif (71-84) conserved between the adrenoleukodystrophy protein and human peroxisomal membrane protein 69 and yeast Pxa1. Omission or truncation of this motif in the adrenoleukodystrophy protein abolished peroxisomal targeting. The single amino acid substitution L78F resulted in a significant reduction of targeting efficiency. The in-frame deletion of three amino acids (del78-80LLR) within the proposed targeting motif in two patients suffering from X-linked adrenoleukodystrophy resulted in the mislocalization of a green fluorescent protein fusion protein to nucleus, cytosol and mitochondria. Our data define the targeting region of human adrenoleukodystrophy protein containing a highly conserved 14-amino-acid motif.     

Identification of active-site residues in protease 3C of hepatitis A virus by site-directed mutagenesis.

Picornavirus 3C proteases (3Cpro) are cysteine proteases related by amino acid sequence to trypsin-like serine proteases. Comparisons of 3Cpro of hepatitis A virus (HAV) to those of other picornaviruses have resulted in prediction of active-site residues: histidine at position 44 (H44), aspartic acid (D98), and cysteine (C172). To test whether these residues are key members of a putative catalytic triad, oligonucleotide-directed mutagenesis was targeted to 3Cpro in the context of natural polypeptide precursor P3. Autocatalytic processing of the polyprotein containing wild-type or variant 3Cpro was tested by in vivo expression of vaccinia virus-HAV chimeras in an animal cell-T7 hybrid system and by in vitro translation of corresponding RNAs. Comparison with proteins present in HAV-infected cells showed that both expression systems mimicked authentic polyprotein processing. Individual substitutions of H44 by tyrosine and of C172 by glycine or serine resulted in complete loss of the virus-specific proteolytic cascade. In contrast, a P3 polyprotein in which D98 was substituted by asparagine underwent only slightly delayed processing, while an additional substitution of valine (V47) by glycine within putative protein 3A caused a more pronounced loss of processing. Therefore, apparently H44 and C172 are active-site constituents whereas D98 is not. The results, furthermore, suggest that substitution of amino acid residues distant from polyprotein cleavage sites may reduce proteolytic activity, presumably by altering substrate conformation.     

Synthesis of octahydropyrano[3,2-b]pyrrole-2-carboxylic acid derivatives from D-mannose

Bicyclic amino acids are useful building blocks in synthesizing biologically active molecules and peptidomimetics. 2-Carboxy-6-hydroxyloctahydroindole (Choi) is a novel bicyclic amino acid found in the marine natural products aeruginosins. Many compounds in the aeruginosin family exhibit inhibition activities toward serine proteases including thrombin and trypsin. The unique Choi structure is the common feature of this family of oligopeptides and this motif is important for their observed biological activities. To better understand the influence of the stereochemistry of the Choi core structure on the inhibition activities, we have previously synthesized ring-oxygenated variants from glucose. The preparation of octahydro-pyrano[3,2-b]pyrrole 2-carboxylic acids from d-mannose is reported here. These novel bicyclic amino acids can be used in the preparation of aeruginosin analogs, as well as conformationally constrained peptidomimetics or other biologically active molecules.     

Ecotin: lessons on survival in a protease-filled world.

Ecotin, an Escherichia coli periplasmic protein of 142 amino acids, has been shown to be a potent inhibitor of a group of homologous serine proteases with widely differing substrate recognition. It is highly effective against a number of enzymes, including both pancreatic and neutrophil-derived elastases, chymotrypsin, trypsin, factor Xa, and kallikrein. Recent structural and functional studies on ecotin and its interactions with different serine proteases have clarified these initial observations and revealed the remarkable features of this protein in inhibiting a strikingly large subset of the chymotrypsin family of serine proteases. The structures of the ecotin:serine protease complexes provide the first examples of protein-protein recognition where the concept of specificity of interactions needs to be reexamined. The binding sites show a fluidity of protein contacts derived from ecotin's innate flexibility in fitting itself to proteases while strongly interfering with their function.     

Extracting Hydrogen-Bond Signature Patterns from Protein Structure Data

Classification of protein sequences and structures into families is a fundamental task in biology, and it is often used as a basis for designing experiments for gaining further knowledge. Some relationships between proteins are detected by the similarities in their sequences, and many more by the similarities in their structures. Despite this, there are a number of examples of functionally similar molecules without any recognisable sequence or structure similarities, and there are also a number of protein molecules that share common structural scaffolds but exhibit different functions. Newer methods of comparing molecules are required in order to detect similarities and dissimilarities in protein molecules. In this article, it is proposed that the precise 3-dimensional disposition of key residues in a protein molecule is what matters for its function, or what conveys the "meaning" for a biological system, but not what means it uses to achieve this. The concept of comparing two molecules through their intramolecular interaction networks is explored, since these networks dictate the disposition of amino acids in a protein structure. First, signature patterns, or fingerprints, of interaction networks in pre-classified protein structural families are computed using an approach to find structural equivalences and consensus hydrogen bonds. Five examples from different structural classes are illustrated. These patterns are then used to search the entire Protein Data Bank, an approach through which new, unexpected similarities have been found. The potential for finding relationships through this approach is highlighted. The use of hydrogen-bond fingerprints as a new metric for measuring similarities in protein structures is also described.     

蛋白质-核酸复合物界面氨基酸与核苷酸偏好性分析

蛋白质-核酸相互作用机制到目前还不是很清楚,尤其是蛋白质与RNA的相互作用。目前,可得到的蛋白质-核酸复合物结构数据不断增多,作者收集了Protein Data Bank数据库中所有的蛋白质-核酸复合物结构数据,对复合物中结合残基和结合核苷酸的偏好性进行了统计分析。发现：1）不同功能的蛋白质-核酸复合物间的结合残基数量存在显著差异;2）在蛋白 质-DNA和蛋白质-RNA复合物界面,碱性氨基酸都是最受欢迎的;3）氨基酸的极性大小及方向在决定它是否与RNA分子进行结合时起到重要的作用,同时发现氨基酸侧链形成的空间位阻会影响氨基酸残基与RNA分子的相互作用;4）随着定义结合残基距离阈值的增大,其氨基酸使用的特异性降低,而受欢迎与不受欢迎的氨基酸种类均没有变化。