Theoretical possibilities and limitations of protease primary specificity determination by statistical analysis of MALDI mass-spectra of proteolysis products

Possibilities and limitations of method examination of proteolytic enzymes' primary specificity by statistical analysis of MALDI (matrix-assisted laser desorption/ionization) mass spectra of products obtained by protein substrates proteolysis without direct determination of their amino acid sequences were investigated theoretically. The optimum ranges given by the errors of the peptides masses measuring for the fabrication of statistical set of the events and the form of statistical data presentation were chosen. It was shown that the proposed method can be applied only for proteases with a relatively narrow primary specificity (two or three amino acids). The influence of protein substrate molecular weight and amino acid composition on the efficiency of specific to a particular protease amino acids display under statistical treatment of the set of proteolysis products masses was studied on the model of trypsin, chymotrypsin, glutamylendopeptidase, pepsin (pH 1.3).     

Analysis of Protein Three-Dimension Structure Using Amino Acids Depths

The issue of amino acid depth in proteins gives important insights to our understanding of protein’s three-dimensional structure. There has already been much research done in mathematical and statistical sciences regarding the general definitions, properties and algorithms describing the particle depth of spatially extended systems. We constructed a method of calculating the amino acids depths and applied it to a set of 527 protein structures. We propose the introduction of amino acid depth tendency factors for three-dimensional structures of proteins. The depth tendency factors relate not only to the hydrophobicity indices but also to the electrostatic charge. We found a relationship between the protein size and the number of residues using the distance between the deepest residue and surface residues. We made a prediction regarding the number of residues on the surface of a protein, the deepest amino acid, and the average depth, all of which are fitted well to a linear functional relationship with the length of the protein. Finally, we have predicted the depths of multiple peptides in protein’s three-dimension structure. Electronic supplementary material The online version of this article () contains supplementary material, which is available to authorized users.     

The preparation of matrix-bound proteases and their use in the hydrolysis of proteins

Four proteases, crude acid protease from Aspergillus, pronase, amino-peptidase M, and prolidase, have been covalently attached to activated agarose and to amino propyl glass beads. The matrix-bound enzymes have been tested as catalysts for the complete hydrolysis of protein substrates, with the primary goal to isolate unstable amino acid derivatives present in the substrate protein. Under conditions used in the present work, the total amino acid release from the protease-catalyzed hydrolysis of four substrate proteins (pancreatic ribonuclease, egg white lysozyme, yeast enolase, and bovine insulin) was 95–103% of that observed in standard acid hydrolysis. Recovery of individual amino acids showed greater deviation from the theoretical values, but cystine was the only amino acid recovered in low yields (42–77%) from all four proteins. Derivatized amino acids, such as methionine sulfoxide, O-(butylcarbamoyl)-serine, and N-glycosyl asparagine have been obtained from chemically modified proteins or from unmodified glycoprotein in good yield, and normal amino acid constituents of proteins which cannot be quantified after acid hydrolysis (tryptophan, asparagine, and glutamine) have also been determined either directly after proteolysis or after proteolysis in conjunction with acid hydrolysis.     

An o-phthalaldehyde spectrophotometric assay for proteinases

A rapid and convenient spectrophotometric assay has been devised to measure proteolysis. The assay is based on the reaction of o-phthalaldehyde (OPA) and 2-mercaptoethanol with amino groups released during proteolysis of a protein substrate. The reaction is specific for primary amines in amino acids, peptides, and proteins, approaches completion within 1 to 2 min at 25 degrees C (half-times of approx 10-15 s), and requires no preliminary heating or separation of the hydrolyzed products from the undegraded protein substrate prior to performing the assay. The OPA assay was relatively as successful as a 2,4,6-trinitrobenzenesulfonic acid (TNBS) procedure in predicting the extent of hydrolysis of a protein substrate. The utility of the OPA method was demonstrated by measuring the degree of proteolytic degradation caused by trypsin, subtilisin, Pronase, and chymotrypsin of various soluble protein substrates. Ethanethiol (instead of 2-mercaptoethanol) or 50% of dimethyl sulfoxide can be included in the assay solution to stabilize certain OPA-amine products. The present method approaches the sensitivity of ninhydrin and TNBS procedures, is more convenient and rapid, and could substitute for these reagents in most assay systems.     

Length recognition at the N-terminal tail for the initiation of FtsH-mediated proteolysis

FtsH-mediated proteolysis against membrane proteins is processive, and presumably involves dislocation of the substrate into the cytosol where the enzymatic domains of FtsH reside. To study how such a mode of proteolysis is initiated, we manipulated N-terminal cytosolic tails of three membrane proteins. YccA, a natural substrate of FtsH was found to require the N-terminal tail of 20 amino acid residues or longer to be degraded by FtsH in vivo. Three unrelated sequences of this segment conferred the FtsH sensitivity to YccA. An artificially constructed TM9-PhoA protein, derived from SecY, as well as the SecE protein, were sensitized to FtsH by addition of extra amino acid sequences to their N-terminal cytosolic tails. Thus, FtsH recognizes a cytosolic region of sufficient length (~20 amino acids) to initiate the processive proteolysis against membrane proteins. Such a region is typically at the N-terminus and can be diverse in amino acid sequences.     

Peptide phosphorylation by calcium-dependent protein kinase from maize seedlings.

Ca2+-dependent protein kinase (CDPK-1) was purified from maize seedlings, and its substrate specificity studied using a set of synthetic peptides derived from the phosphorylatable sequence RVLSRLHS15VRER of maize sucrose synthase 2. The decapeptide LARLHSVRER was found to be efficiently phosphorylated as a minimal substrate. The same set of peptides were found to be phosphorylated by mammalian protein kinase Cbeta (PKC), but showed low reactivity with protein kinase A (PKA). Proceeding from the sequence LARLHSVRER, a series of cellulose-membrane-attached peptides of systematically modified structure was synthesised. These peptides had hydrophobic (Ala, Leu) and ionic (Arg, Glu) amino acids substituted in each position. The phosphorylation of these substrates by CDPK-1 was measured and the substrate specificity of the maize protein kinase characterised by the consensus sequence motif A/L-5X-4R-3X-2X-1SX+1R+2Z+3R+4, where X denotes a position with no strict amino acid requirements and Z a position strictly not tolerating arginine compared with the other three varied amino acids. This motif had a characteristic sequence element RZR at positions +2 to +4 and closely resembled the primary structure of the sucrose synthase phosphorylation site. The sequence surrounding the phosphorylatable serine in this consensus motif was similar to the analogous sequence K/RXXS/TXK/R proposed for mammalian PKC, but different from the consensus motif RRXS/TX for PKA.     

Domain organization and biochemical features of Sulfolobus solfataricus DNA polymerase

DNA polymerase from Sulfolobus solfataricus, strain MT4 (Sso DNA pol), was one of the first archaeal DNA polymerases to be isolated and characterized. Its encoding gene was cloned and sequenced, indicating that Sso DNA pol belongs to family B of DNA polymerases. By limited proteolysis experiments carried out on the recombinant homogeneous protein, we were able to demonstrate that the enzyme has a modular organization of its associated catalytic functions (DNA polymerase and 3′-5′ exonuclease). Indeed, the synthetic function was ascribed to the enzyme C-terminal portion, whereas the N-terminal half was found to be responsible for the exonucleolytic activity. In addition, partial proteolysis studies were utilized to map conformational changes on DNA binding by comparing the cleavage map in the absence or presence of nucleic acid ligands. This analysis allowed us to identify two segments of the Sso DNA pol amino acid chain affected by structural modifications following nucleic acid binding: region 1 and region 2, in the middle and at the C-terminal end of the protein chain, respectively. Site-directed mutagenesis studies will be performed to better investigate the role of these two protein segments in DNA substrate interaction. Received: January 22, 1998 / Accepted: February 16, 1998     

Proteins as special subsets of polypeptides

A large protein sequence database with over 31,000 sequences and 10 million residues has been analysed. The pair probabilities have been converted to entropies using Boltzmann’s law of statistical thermodynamics. A scoring weight corresponding to “mixing entropy” of the amino acid pairs has been developed from which the entropies of the protein sequences have been calculated. The entropy values of natural sequences are lower than their random counterparts of same length and similar amino acid composition. Based on the results it has been proposed that natural sequences are a special set of polypeptides with additional qualification of biological functionality that can be quantified using the entropy concept as worked out in this paper.     

Cleavage site selection within a folded substrate by the ATP-dependent lon protease

Mechanistic studies of ATP-dependent proteolysis demonstrate that substrate unfolding is a prerequisite for processive peptide bond hydrolysis. We show that mitochondrial Lon also degrades folded proteins and initiates substrate cleavage non-processively. Two mitochondrial substrates with known or homology-derived three-dimensional structures were used: the mitochondrial processing peptidase alpha-subunit (MPPalpha) and the steroidogenic acute regulatory protein (StAR). Peptides generated during a time course of Lon-mediated proteolysis were identified and mapped within the primary, secondary, and tertiary structure of the substrate. Initiating cleavages occurred preferentially between hydrophobic amino acids located within highly charged environments at the surface of the folded protein. Subsequent cleavages proceeded sequentially along the primary polypeptide sequence. We propose that Lon recognizes specific surface determinants or folds, initiates proteolysis at solvent-accessible sites, and generates unfolded polypeptides that are then processively degraded.     

Predicting protein crystallization propensity from protein sequence

The high-throughput structure determination pipelines developed by structural genomics programs offer a unique opportunity for data mining. One important question is how protein properties derived from a primary sequence correlate with the protein’s propensity to yield X-ray quality crystals (crystallizability) and 3D X-ray structures. A set of protein properties were computed for over 1,300 proteins that expressed well but were insoluble, and for ~720 unique proteins that resulted in X-ray structures. The correlation of the protein’s iso-electric point and grand average hydropathy (GRAVY) with crystallizability was analyzed for full length and domain constructs of protein targets. In a second step, several additional properties that can be calculated from the protein sequence were added and evaluated. Using statistical analyses we have identified a set of the attributes correlating with a protein’s propensity to crystallize and implemented a Support Vector Machine (SVM) classifier based on these. We have created applications to analyze and provide optimal boundary information for query sequences and to visualize the data. These tools are available via the web site .     

Identification of RNA-binding sites in proteins by integrating various sequence information

RNA–protein interactions play a pivotal role in various biological processes, such as mRNA processing, protein synthesis, assembly, and function of ribosome. In this work, we have introduced a computational method for predicting RNA-binding sites in proteins based on support vector machines by using a variety of features from amino acid sequence information including position-specific scoring matrix (PSSM) profiles, physicochemical properties and predicted solvent accessibility. Considering the influence of the surrounding residues of an amino acid and the dependency effect from the neighboring amino acids, a sliding window and a smoothing window are used to encode the PSSM profiles. The outer fivefold cross-validation method is evaluated on the data set of 77 RNA-binding proteins (RBP77). It achieves an overall accuracy of 88.66% with the Matthew’s correlation coefficient (MCC) of 0.69. Furthermore, an independent data set of 39 RNA-binding proteins (RBP39) is employed to further evaluate the performance and achieves an overall accuracy of 82.36% with the MCC of 0.44. The result shows that our method has good generalization abilities in predicting RNA-binding sites for novel proteins. Compared with other previous methods, our method performs well on the same data set. The prediction results suggest that the used features are effective in predicting RNA-binding sites in proteins. The code and all data sets used in this article are freely available at .     

Dipeptide synthesis by L-amino acid ligase from Ralstonia solanacearum

Despite its utility, dipeptides have not been widely used due to the absence of an efficient manufacturing method. Recently, a novel method for effective production of dipeptides using l-amino acid α-ligase (Lal) is presented. Lal, which is only identified in Bacillus subtilis, catalyzes dipeptide synthesis from unprotected amino acids in an ATP-dependent manner. However, not all the dipeptide can be synthesized by Lal from B. subtilis (BsLal) due to its substrate specificity. Here, we attempted to find a novel Lal exhibiting different substrate specificity from BsLal. By in silico screening based on the amino acid sequence of BsLal, RSp1486a an unknown protein from Ralstonia solanacearum was found to show the Lal activity. RSp1486a exhibited different substrate specificity from BsLal, and preferably synthesized hetero-dipeptides where more bulky amino acid was placed at N terminus and less bulky amino acid was placed at C terminus in opposition to those synthesized by BsLal.     

Hepatocyte plasma membrane ECTO-ATPase (pp120/HA4) is a substrate for tyrosine kinase activity of the insulin receptor

pp120/HA4, a membrane protein found in hepatocyte plasma membranes and a substrate for the insulin receptor tyrosine kinase, was purified to homogeneity, subjected to partial proteolysis, and peptides were sequenced by Edman degradation. Six amino acid sequences were obtained, and they matched the deduced amino acid sequences of six regions of a hepatocyte membrane protein called ecto-ATPase.     

Automated molecular library generation of proteic fragments by virtual proteolysis for molecular modelling studies

This paper presents a computer aided design method useful for simulation of a set of proteolytic cleavages upon target proteins obtained from the Brookhaven Data Bank. The method was developed by using algorithms that are able to interface themselves with other software environments, in order to assist computer analyses in the molecular modelling field, and allowing the generation of molecular libraries containing protein fragments produced by simulated proteolysis. These libraries include structures that differ for several amino acid deletions upon specified regions of the primary sequence. Target residues chosen for the simulation are compatible with enzymatic proteolysis methods used in conventional laboratory procedures. Furthermore, algorithms were able to identify a set of chemical-physical properties of the starting proteins, leading the simulation to find out the most suitable residues for proteolysis. The goal of these strategies is to generate fragments that are leaded to maintain the native-like condition of starting molecules, avoiding loss of conformational characteristics of the original tertiary structure. Proteins chosen for generating proteolytic libraries were represented by naphthalene 1,2 dioxygenase and Rigidoporus lignosus laccase.     

Prediction of protein structural classes by Chou's pseudo amino acid composition: approached using continuous wavelet transform and principal component analysis

A prior knowledge of protein structural classes can provide useful information about its overall structure, so it is very important for quick and accurate determination of protein structural class with computation method in protein science. One of the key for computation method is accurate protein sample representation. Here, based on the concept of Chou’s pseudo-amino acid composition (AAC, Chou, Proteins: structure, function, and genetics, 43:246–255, 2001), a novel method of feature extraction that combined continuous wavelet transform (CWT) with principal component analysis (PCA) was introduced for the prediction of protein structural classes. Firstly, the digital signal was obtained by mapping each amino acid according to various physicochemical properties. Secondly, CWT was utilized to extract new feature vector based on wavelet power spectrum (WPS), which contains more abundant information of sequence order in frequency domain and time domain, and PCA was then used to reorganize the feature vector to decrease information redundancy and computational complexity. Finally, a pseudo-amino acid composition feature vector was further formed to represent primary sequence by coupling AAC vector with a set of new feature vector of WPS in an orthogonal space by PCA. As a showcase, the rigorous jackknife cross-validation test was performed on the working datasets. The results indicated that prediction quality has been improved, and the current approach of protein representation may serve as a useful complementary vehicle in classifying other attributes of proteins, such as enzyme family class, subcellular localization, membrane protein types and protein secondary structure, etc.     

Partly folded states of members of the lysozyme/lactalbumin superfamily: a comparative study by circular dichroism spectroscopy and limited proteolysis

The partly folded states of protein members of the lysozyme (LYS)/alpha-lactalbumin (LA) superfamily have been analyzed by circular dichroism (CD) measurements and limited proteolysis experiments. Hen, horse, dog, and pigeon LYSs and bovine LA were used in the present study. These are related proteins of 123- to 129-amino-acid residues with similar three-dimensional structures but low similarity in amino acid sequences. Moreover, notable differences among them reside in their calcium-binding properties and capability to adopt partly folded states or molten globules in acid solution (A-state) or on depletion of calcium at neutral pH (apo-state). Far- and near-UV CD measurements revealed that although the structures of hen and dog LYS are rather stable in acid at pH 2.0 or at neutral pH in the absence of calcium, conformational transitions to various extents occur with all other LYS/LA proteins herewith investigated. The most significant perturbation of tertiary structure in acid was observed with bovine LA and LYS from horse milk and pigeon egg-white. Pepsin and proteinase K were used as proteolytic probes, because these proteases show broad substrate specificity, and therefore, their sites of proteolysis are dictated not by the specific amino acid sequence of the protein substrate but by its overall structure and dynamics. Although hen LYS at pH 2.0 was fully resistant to proteolysis by pepsin, the other members of the LYS/LA superfamily were cleaved at different rates at few sites of the polypeptide chain and thus producing rather large protein fragments. The apo-form of bovine LA, horse LYS, and pigeon LYS were attacked by proteinase K at pH 8.3, whereas dog and hen LYSs were resistant to proteolysis when reacted under identical experimental conditions. Briefly, it has been found that the proteolysis data correlate well with the extent of conformational transitions inferred from CD spectra and with existing structural informations regarding the proteins herewith investigated, mainly derived from NMR and hydrogen exchange measurements. The sites of initial proteolytic cleavages in the LYS variants occur at the level of the beta-subdomain (approximately chain region 34-57), in analogy to those observed with bovine LA. Proteolysis data are in agreement with the current view that the molten globule of the LYS/LA proteins is characterized by a structured alpha-domain and a largely disrupted beta-subdomain. Our results underscore the utility of the limited proteolysis approach for analyzing structure and dynamics of proteins, even if adopting an ensemble of dynamic states as in the molten globule.     

Cloning of a Novel Omega-6 Desaturase from Flax (<Emphasis Type="Italic">Linum usitatissimum</Emphasis> L.) and Its Functional Analysis in <Emphasis Type="Italic">Saccharomyces cerevisiae</Emphasis>

The Δ12 desaturase represents a diverse gene family in plants and is responsible for conversion of oleic acid (18:1) to linoleic acid (18:2). Several members of this family are known from plants like Arabidopsis and Soybean. Using primers from conserved C- and N-terminal regions, we have cloned a novel Δ12 desaturase gene amplified from flax genomic DNA, denoted as LuFAD2-2. This intron-less gene is 1,149-base pair long encoding 382 amino acids—putative membrane-bound Δ12 desaturase protein. Sequence comparisons show that the novel sequence has 85% similarity with previously reported flax Δ12 desaturase at amino acid level and shows typical features of membrane-bound desaturase such as three conserved histidine boxes along with four membrane-spanning regions that are universally present among plant desaturases. The signature amino acid sequence ‘YNNKL’ was also found to be present at the N terminus of the protein, which is necessary and sufficient for ER localization of enzyme. Neighbor-Joining tree generated from the sequence alignment grouped LuFAD2-2 among the other FAD2 sequences from Ricinus, Hevea, Jatropha, and Vernicia. When LuFAD2-2 and LuFAD2 were expressed in Saccharomyces cerevisiae, they could convert the oleic acid to linoleic acid, with an average conversion rate of 5.25 and 8.85%, respectively. However, exogenously supplied linoleic acid was feebly converted to linolenic acid suggesting that LuFAD2-2 encodes a functional FAD2 enzyme and has substrate specificity similar to LuFAD2.     

A cumulative specificity model for proteases from human immunodeficiency virus types 1 and 2, inferred from statistical analysis of an extended substrate data base.

Statistical analysis of an expanded data base of regions in viral polyproteins and in non-viral proteins that are sensitive to hydrolysis by the protease from human immunodeficiency virus (HIV) type 1 has generated a model which characterizes the substrate specificity of this retroviral enzyme. The model leads to an algorithm for predicting protease-susceptible sites from primary structure. Amino acids in each of the sites from P4 to P4' are tabulated for 40 protein substrates, and the frequency of occurrence for each residue is compared to the natural abundance of that amino acid in a selected data set of globular proteins. The results suggest that the highest stringency for particular amino acid residues is at the P2, P1, and P2' positions of the substrate. The broad specificity of the HIV-1 protease appears to be a consequence of its being able to bind productively substrates in which interactions with only a few Pi or Pi' side-chains need be optimized. The analysis, extended to 22 protein segments cleaved by the HIV-2 protease, delineates marked differences in specificity from that of the HIV-1 enzyme.     

A Combinatorial Approach to Detect Coevolved Amino Acid Networks in Protein Families of Variable Divergence

Communication between distant sites often defines the biological role of a protein: amino acid long-range interactions are as important in binding specificity, allosteric regulation and conformational change as residues directly contacting the substrate. The maintaining of functional and structural coupling of long-range interacting residues requires coevolution of these residues. Networks of interaction between coevolved residues can be reconstructed, and from the networks, one can possibly derive insights into functional mechanisms for the protein family. We propose a combinatorial method for mapping conserved networks of amino acid interactions in a protein which is based on the analysis of a set of aligned sequences, the associated distance tree and the combinatorics of its subtrees. The degree of coevolution of all pairs of coevolved residues is identified numerically, and networks are reconstructed with a dedicated clustering algorithm. The method drops the constraints on high sequence divergence limiting the range of applicability of the statistical approaches previously proposed. We apply the method to four protein families where we show an accurate detection of functional networks and the possibility to treat sets of protein sequences of variable divergence.