Searching for local similarities between HIV-1 and human proteins. Application to vaccines

A method for identification of fragments with a high local similarity to human proteins within potentially immunopathogenic regions of HIV-1 proteins was developed. The method is based on the use of an original matrix of antigenic similarity of amino acids. The regions, whose fragments are frequent in human proteins, and regions, exhibiting a high similarity to the proteins responsible for important physiological functions, were identified in HIV-1 proteins. A possibility of participation of such regions in immunopathogenesis of HIV-infection due either to induction of cross-reacting effectors of immune system or through molecular mimicry of physiologically important human proteins, leading to an alteration of homeostasis of the organism, is discussed. Most of regions, identified in HIV-1 proteins, contain either T-cell (CD8+ CTL or CD4+ Th) or B-cell epitopes, or both of them simultaneously. The criteria for the design of safe polyepitopic antiviral vaccines which enable the exclusion of epitopes, having the (immuno)pathogenic potential, are discussed. According to these criteria, polyepitopic immunogens should be free of the virus protein regions, whose fragments are frequent in human proteins, as well as of regions exhibiting a pronounced local similarity to proteins that provide for important physiological functions.     

Amino acid similarity matrix for detecting the antigenically similar sequences in non-related proteins

A special matrix of amino acid antigenic similarity for computer detection of the potential antigenic proximity of unrelated proteins is proposed. The matrix was built using the data concerning affinities of amino acid residue interactions between subunits in oligomeric proteins. The diagonal elements of the matrix characterize the recognition of amino acid residues and the non-diagonal ones represent the relative similarity measure of antibody--amino acid residue interactions specificity. The application of the new matrix for comparing proteins allows the hydrophilic potentially immunologically active regions of sequences to be picked out as similar fragments. When the influenza virus hemagglutinin was compared with 116 human proteins, eight fragments were picked out, that could not be determined by means of the routinely used MDM78 matrix. The antigenic similarity matrix for defining the forbidden structures is proposed to be used for preparing the peptidic antiviral vaccines.     

Searching for Local Similarities between HIV-1 and Human Proteins. Application to Vaccines

A method for identification of fragments with high local similarity to human proteins within potentially immunopathogenic regions of HIV-1 proteins was developed. The method is based on the use of an original matrix of antigenic similarity of amino acids. The regions whose fragments are frequent in human proteins, and regions exhibiting high similarity to the proteins responsible for important physiological functions, were identified in HIV-1 proteins. A possibility of participation of such regions in immunopathogenesis of HIV infection either through induction of cross-reacting effectors of the immune system or through molecular mimicry of physiologically important human proteins, leading to alteration of homeostasis of the organism, is discussed. Most of the regions identified in HIV-1 proteins contain either T-cell (CD8⁺ CTL or CD4⁺ Th) or B-cell epitopes, or both of them simultaneously. The criteria for the design of safe polyepitopic antiviral vaccines that would allow exclusion of epitopes with (immuno)pathogenic potential are discussed. According to these criteria, polyepitopic immunogens should be free of the virus protein regions whose fragments are frequent in human proteins, as well as of regions exhibiting pronounced local similarity to proteins that provide for important physiological functions.     

Modification of the Bacillus sphaericus 51- and 42-kilodalton mosquitocidal proteins: effects of internal deletions, duplications, and formation of hybrid proteins

The 51- and 42-kDa proteins which constitute the binary mosquitocidal toxin of Bacillus sphaericus 2362 have a low overall sequence similarity but share several regions of near identity (L. Baumann, A. H. Broadwell, and P. Baumann, J. Bacteriol. 170:2045-2050, 1988). By using site-directed mutagenesis, deletions of 6 to 16 amino acids in three of these regions of the 51- and 42-kDa proteins were made, and the modified proteins were expressed in Bacillus subtilis. Deletions in both of these proteins resulted in a loss of toxicity for mosquito larvae. Hybrid proteins containing exchanged fragments of the 51- and 42-kDa proteins were inactive when tested in a variety of combinations, thereby indicating that potentially analogous fragments of these two proteins were not functionally equivalent. An internal duplication of 73 amino acids in the 51-kDa protein and 72 amino acids in the 42-kDa protein resulted in a major reduction in toxicity. These results indicate that the conserved regions of the 51- and 42-kDa proteins are necessary for toxicity to larvae and that the 51- and 42-kDa proteins, despite their sequence similarity, are unique, differing from each other by at least one essential attribute.     

Modification of the Bacillus sphaericus 51- and 42-kilodalton mosquitocidal proteins: effects of internal deletions, duplications, and formation of hybrid proteins.

The 51- and 42-kDa proteins which constitute the binary mosquitocidal toxin of Bacillus sphaericus 2362 have a low overall sequence similarity but share several regions of near identity (L. Baumann, A. H. Broadwell, and P. Baumann, J. Bacteriol. 170:2045-2050, 1988). By using site-directed mutagenesis, deletions of 6 to 16 amino acids in three of these regions of the 51- and 42-kDa proteins were made, and the modified proteins were expressed in Bacillus subtilis. Deletions in both of these proteins resulted in a loss of toxicity for mosquito larvae. Hybrid proteins containing exchanged fragments of the 51- and 42-kDa proteins were inactive when tested in a variety of combinations, thereby indicating that potentially analogous fragments of these two proteins were not functionally equivalent. An internal duplication of 73 amino acids in the 51-kDa protein and 72 amino acids in the 42-kDa protein resulted in a major reduction in toxicity. These results indicate that the conserved regions of the 51- and 42-kDa proteins are necessary for toxicity to larvae and that the 51- and 42-kDa proteins, despite their sequence similarity, are unique, differing from each other by at least one essential attribute.     

A program to find regions of similarity between homologous protein sequences using dotmatrix analysis

MATRIX is a program designed primarily to enable the user to visualize all regions of similarity between two proteins at a glance. The program helps the user to see where they are similar—at what relative positions in the amino acid sequences of the two proteins in question does the similarity exist; how they are similar—what functional characteristics the two similar sequences have in common; and to what extent they are similar—is the similarity significant, if so how significant relative to other similar sequences in the protein. This is achieved by constructing a diagram in which quantitative parameters of amino acids are used to compare every amino acid residue of the first protein with every amino acid residue of the second.Another function of the program is, given two sets of atomic coordinates—either of different proteins or for the same protein (for self-comparison)—to demonstrate which residues of the two proteins, when the two proteins are superimposed upon each other, appear in the same space (or are close to each other).     

Protein structure, neighbor effect, and a new index of amino acid dissimilarities.

Amino acids interact with each other, especially with neighboring amino acids, to generate protein structures. We studied the pattern of association and repulsion of amino acids based on 24,748 protein-coding genes from human, 11,321 from mouse, and 15,028 from Escherichia coli, and documented the pattern of neighbor preference of amino acids. All amino acids have different preferences for neighbors. We have also analyzed 7,342 proteins with known secondary structure and estimated the propensity of the 20 amino acids occurring in three of the major secondary structures, i.e., helices, sheets, and turns. Much of the neighbor preference can be explained by the propensity of the amino acids in forming different secondary structures, but there are also a number of intriguing association and repulsion patterns. The similarity in neighbor preference among amino acids is significantly correlated with the number of amino acid substitutions in both mitochondrial and nuclear genes, with amino acids having similar sets of neighbors replacing each other more frequently than those having very different sets of neighbors. This similarity in neighbor preference is incorporated into a new index of amino acid dissimilarities that can predict nonsynonymous codon substitutions better than the two existing indices of amino acid dissimilarities, i.e., Grantham's and Miyata's distances.     

Similarity between average distance maps of structurally homologous proteins

A similarity between average distance maps (Kikuchiet al., 1988a)—that is, predicted contact maps of two tertiary structurally homologous proteins—is examined. Comparisons of shapes of average distance maps (we refer to this as ADM) are made by superpositions of ADMs for two homologous proteins. Also, we compare shapes of actual contact maps for the pair of proteins. We search a optimal superposition mode of each pair of maps showing that two proteins are most similar. It is concluded that two ADMs are also similar when actual tertiary structures between two proteins show similarity. A criterion for similarity of maps is also proposed. The possibility of application of this method to detect weak homology between protein structures is discussed.     

Sequencing of bovine enamelin ("tuftelin") a novel acidic enamel protein

Enamelins are a major group of 28-70-kDa acidic proteins rich in aspartic acid, glutamic acid, serine, and glycine found in developing and mature extracellular enamel; a unique and highly mineralized ectodermal tissue covering vertebrate teeth. They have been associated with the mineralization and structural organization of this tissue. In an attempt to elucidate the primary structure of enamelin, a 2674-base pair cDNA isolated from a bovine ameloblast-enriched, lambda Zap 2 expression library, was sequenced. The identity and localization of the deduced protein was confirmed by amino acid composition, enzyme-linked immunosorbent assay, Western blotting, indirect immunohistochemistry, and high resolution protein-A gold immunocytochemistry. The immunological techniques were employed using antibodies directed against synthetic peptides corresponding to the protein sequence deduced from the cloned cDNA sequence. The results reveal the deduced protein to be a novel acidic enamel protein. It contains 389 amino acids and has a calculated molecular weight of 43,814. Its amino acid composition is similar to that of "tuft" proteins (enamel matrix protein fragments remaining in the mature tissue). It contains one potential N-glycosylation site and 5 cysteine residues. Southern hybridization of the cloned cDNA with genomic bovine DNA indicated the existence of a single gene with one or more introns.     

Computer prediction of peptide maps: assignment of polypeptides to human and mouse mitochondrial DNA genes by analysis of two-dimensional-proteolytic digest gels.

We have prepared a computer program that predicts complete and partial peptide maps from amino acid sequences. The program fragments amino acid sequences at designated cleavage sites and calculates the molecular weight and relative labeling of each peptide. These data are graphed as log molecular weight of the original protein (X-axis) vs. log molecular weight of the component peptides (Y-axis). The program is interactive, permitting adjustment of a number of graphic parameters and alteration of the position of proteins in the first dimension to accommodate aberrations in protein mobility. The program has been used to predict the V8 protease peptide maps of the 13 open reading frames (ORFs) identified in the human and the mouse mitochondrial DNA (mtDNA) sequences. The results were compared to the V8 protease peptide maps obtained for mouse and human mitochondrially synthesized proteins by two-dimensional proteolytic digest gels. A high correlation was observed between the predicted and observed peptide maps. These results suggest the assignment of several proteins to mtDNA genes.     

Alignment of protein sequences using secondary structure: a modified dynamic programming method

A method for comparison of protein sequences based on their primary and secondary structure is described. Protein sequences are annotated with predicted secondary structures (using a modified Chou and Fasman method). Two lettered code sequences are generated (Xx, where X is the amino acid and x is its annotated secondary structure). Sequences are compared with a dynamic programming method (STRALIGN) that includes a similarity matrix for both the amino acids and secondary structures. The similarity value for each paired two-lettered code is a linear combination of similarity values for the paired amino acids and their annotated secondary structures. The method has been applied to eight globin proteins (28 pairs) for which the X-ray structure is known. For protein pairs with high primary sequence similarity (greater than 45%), STRALIGN alignment is identical to that obtained by a dynamic programming method using only primary sequence information. However, alignment of protein pairs with lower primary sequence similarity improves significantly with the addition of secondary structure annotation. Alignment of the pair with the least primary sequence similarity of 16% was improved from 0 to 37% 'correct' alignment using this method. In addition, STRALIGN was successfully applied to seven pairs of distantly related cytochrome c proteins, and three pairs of distantly related picornavirus proteins.     

Representation of amino acids as five-bit or three-bit patterns for filtering protein databases

MOTIVATION: We propose representing amino acids by bit-patterns so they may be used in a filter algorithm for similarity searches over protein databases, to rapidly eliminate non-homologous regions of database sequences. The filter algorithm would be based on dynamic programming optimization. It would have the advantage over previous filter algorithms that its substitution scoring function distinguishes between conservative and non-conservative amino acid substitutions. RESULTS: Simulated annealing was used to search for the best five-bit or three-bit patterns to represent amino acids, where similar amino acids were given similar bit-patterns. The similarity between amino acids was estimated from the BLOSUM45 matrix. Representing amino acids by these five-bit and three-bit patterns, the Escherichia coli PhoE precursor and the bacteriophage PA2 LC precursor were aligned. The alignments were nearly the same as that obtained when BLOSUM45 was used to score substitutions. AVAILABILITY: The C code of the optimization algorithm for searching for the optimal bit-pattern representation of amino acids is available from the authors upon request.     

Use of Chou-Fasman amino acid conformational parameters to analyze the organization of the genetic code and to construct protein genealogies

Summary Chou-Fasman parameters, measuring preferences of each amino acid for different conformational regions in proteins, were used to obtain an amino acid difference index of conformational parameter distance (CPD) values. CPD values were found to be significantly lower for amino acid exchanges representing in the genetic code transitions of purines, GA than for exchanges representing either transitions of pyrimidines, CU, or transversions of purines and pyrimidines. Inasmuch as the distribution of CPD values in these non GA exchanges resembles that obtained for amino acid pairs with double or triple base differences in their underlying codons, we conclude that the genetic code was not particularly designed to minimize effects of mutation on protein conformation. That natural selection minimizes these changes, however, was shown by tabulating results obtained by the maximum parsimony method for eight protein genealogies with a total occurrence of 4574 base substitutions. At the beginning position of the codons GA transitions were in very great excess over other base substitutions, and, conversely, CU transitions were deficient. At the middle position of the codons only fast evolving proteins showed an excess of GA transitions, as though selection mainly preserved conformation in these proteins while weeding out mutations affecting chemical properties of functional sites in slow evolving proteins. In both fast and slow evolving proteins the net direction of transitions and transversions was found to be from G beginning codons to non-G beginning codons resulting in more commonly occurring amino acids, especially alanine with its generalized conformational properties, being replaced at suitable sites by amino acids with more specialized conformational and chemical properties. Historical circumstances pertaining to the origin of the genetic code and the nature of primordial proteins could account for such directional changes leading to increases in the functional density of proteins.In order to further explore the course of protein evolution, a modified parsimony algorithm was developed for constructing protein genealogies on the basis of minimum CPD length. The algorithm's ability to judge with finer discrimination that in protein evolution certain pathways of amino acid substitution should occur more readily than others was considered a potential advantage over strict maximum parsimony. In developing this CPD algorithm, the path of minimum CPD length through intermediate amino acids allowed by the genetic code for each pair of amino acids was determined. It was found that amino acid exchanges representing two base changes have a considerably lower average CPD value per base substitution than the amino acid exchanges representing single base changes. Amino acid exchanges representing three base changes have yet a further marked reduction in CPD per base change. This shows how extreme constraining effects of stabilizing selection can be circumvented, for by way of intermediate amino acids almost any amino acid can ultimately be substituted for another without damage to an evolving protein's conformation during the process.     

Independent genes coding for three acidic proteins of the large ribosomal subunit from Saccharomyces cerevisiae

The yeast ribosome contains three acidic proteins, L44, L44', and L45, closely related from a structural point of view, that seem to play a functional role similar to that of proteins L7 and L12 in the bacterial ribosome. By screening a cDNA bank in lambda gt11 with specific polyclonal and monoclonal antibodies, recombinant phages expressing each one of the acidic proteins have been cloned. A unique copy of each gene is detected using the phage cDNA inserts as probes in nitrocellulose blots of yeast DNA digested with different restriction enzymes. The inserts were subcloned in the plasmid pUC19, and their physical maps and nucleotide sequences were determined. By using the cDNA inserts as probes in genomic DNA banks, DNA fragments carrying the acidic protein genes have been cloned, characterized, and sequenced. The results conclusively show that the three yeast acidic proteins are coded by independent genes and are not the result of a post-translational modification of the product of a unique gene, as in bacteria. Like most ribosomal protein genes, the gene for protein L44' has an intron and two upstream stimulatory boxes (UASrpg) fitting closely to the consensus sequence. The genes coding for proteins L44 and L45 lack introns and seem also exceptional in other characteristics of their sequences. Proteins L44 and L45 have amino acid sequences with about 80% similarity. Protein L44' is only 63% similar to the other two polypeptides. The three proteins have highly conserved carboxyl termini comprising the last 30 amino acids, and the first 10 amino acids of L44 and L45 are identical. The results cast doubts about the possibility of a similar role for the different acidic ribosomal proteins.     

Improved sensitivity of biological sequence database searches

We have increased the sensitivity ofDNA and protein sequencedatabase searches by allowing similar but non-identical aminoacids or nucleotides to match. In addition, one can match k-tuplesor words instead of matching individual residues in order tospeed the search. A matching matrix specifies which k-tuplesmatch each other. The matching matrix can be calculated froma similarity matrix of amino acids and a threshold of similarityrequired for matching. This permits amino acid similarity matricesor replacement matrices (PAM matrices) to be used in the firststep of a sequence comparison rather than in a secondary scoringphase. The concept of matching non-identical k-tuples also increasesthe power ofDNA database searches. For example, a matrix thatspecifies that any 3-tuple in a DNA sequence can match any other3-tuple encoding the same amino acid permits a DNA databasesearch using a DNA query sequence for regions that would encodea similar amino acid sequence. Received on October 10, 1989; accepted on May 1, 1990     

Nucleotide and amino acid sequences of oriT-traM-traJ-traY-traA-traL regions and mobilization of virulence plasmids of Salmonella enterica serovars enteritidis,gallinarum-pullorum,and typhimurium

The virulence plasmid of Salmonella enterica serovar Gallinarum-Pullorum (pSPV) but not those of Salmonella enterica serovars Enteritidis (pSEV) and Typhimurium (pSTV) can be readily mobilized by an F or F-like conjugative plasmid. To investigate the reason for the difference, the oriT-traM-traJ-traY-traA-traL regions of the three salmonella virulence plasmids (pSVs) were cloned and their nucleotide and deduced amino acid sequences were examined. The cloned fragments were generally mobilized more readily than the corresponding full-length pSVs, but the recombinant plasmid containing the oriT of pSPV was, as expected, more readily mobilized, with up to 100-fold higher frequency than the recombinant plasmids containing the oriT of the other two pSVs. The nucleotide sequences of the oriT-traM-traJ-traY-traA-traL region of pSEV and pSTV were almost identical (only 4 bp differences), but differed from that of pSPV. Major nucleotide sequence variations were found in traJ, traY, and the Tra protein binding sites sby and sbm. sby of pSPV showed higher similarity than that of pSEV or pSTV to that of the F plasmid. The reverse was true for sbm: similarity was higher with pSEV and pSTV than with pSPV. In the deduced amino acid sequences of the five Tra proteins, major differences were found in TraY: pSEV's TraY was 75 amino acids, pSTV's was 106 amino acids, and pSPV's was 133 amino acids; and there were duplicate consensus betaalphaalpha fragments in the TraY of pSPV and F plasmid, whereas there was only a single betaalphaalpha fragment in that of pSEV and pSTV.     

Solution structure of biopolymers: a new method of constructing a bead model

We propose a new, automated method of converting crystallographic data into a bead model used for the calculations of hydrodynamic properties of rigid macromolecules. Two types of molecules are considered: nucleic acids and small proteins. A bead model of short DNA fragments has been constructed in which each nucleotide is represented by two identical, partially overlapping spheres: one for the base and one for the sugar and phosphate group. The optimum radius sigma = 5.0 A was chosen on the basis of a comparison of the calculated translational diffusion coefficients (D(T)) and the rotational relaxation times (tau(R)) with the corresponding experimental data for B-DNA fragments of 8, 12, and 20 basepairs. This value was assumed for the calculation D(T) and tau(R) of tRNA(Phe). Better agreement with the experimental data was achieved for slightly larger sigma = 5.7 A. A similar procedure was applied to small proteins. Bead models were constructed such that each amino acid was represented by a single sphere or a pair of identical, partially overlapping spheres, depending on the amino acid's size. Experimental data of D(T) of small proteins were used to establish the optimum value of sigma = 4.5 A for amino acids. The lack of experimental data on tau(R) for proteins restricted the tests to the translational diffusion properties.     

Structural comparison of Neisseria gonorrhoeae outer membrane proteins.

Outer membranes from opaque colonia variants of Neisseria gonorrhoeae P9 contain a major outer membrane protein (protein I) together with one or more of a series of heat-modifiable proteins (proteins II). Proteins I. II, and IIa have been isolated by detergent extraction of outer membranes. Amino acid analysis showed proteins II and IIa to have a very similar composition. Cyanogen bromide cleavage of proteins II and IIa produced a pair of fragments with identical molecular weight and a pari which differed by an amount (0.5K) equivalent to the difference between the intact proteins. Tryptic peptide maps of 125I-labeled proteins II, IIa, and IIb showed many similarities, with only a few peptides unique to any one protein. Peptide maps of protein IIa from cells which had been surface labeled showed that the unique peptides were exposed on the surface. The heat-modifiable proteins thus appear to form a family of proteins with closely related structure probably differing in that part which is exposed on the bacterial surface.     

On the properties and sequence context of structurally ambivalent fragments in proteins

The goal of this work is to characterize structurally ambivalent fragments in proteins. We have searched the Protein Data Bank and identified all structurally ambivalent peptides (SAPs) of length five or greater that exist in two different backbone conformations. The SAPs were classified in five distinct categories based on their structure. We propose a novel index that provides a quantitative measure of conformational variability of a sequence fragment. It measures the context-dependent width of the distribution of (phi,xi) dihedral angles associated with each amino acid type. This index was used to analyze the local structural propensity of both SAPs and the sequence fragments contiguous to them. We also analyzed type-specific amino acid composition, solvent accessibility, and overall structural properties of SAPs and their sequence context. We show that each type of SAP has an unusual, type-specific amino acid composition and, as a result, simultaneous intrinsic preferences for two distinct types of backbone conformation. All types of SAPs have lower sequence complexity than average. Fragments that adopt helical conformation in one protein and sheet conformation in another have the lowest sequence complexity and are sampled from a relatively limited repertoire of possible residue combinations. A statistically significant difference between two distinct conformations of the same SAP is observed not only in the overall structural properties of proteins harboring the SAP but also in the properties of its flanking regions and in the pattern of solvent accessibility. These results have implications for protein design and structure prediction.     

Amino acid sequences of proteins from Leptospira serovar pomona

This report describes a partial amino acid sequences from three putative outer envelope proteins from Leptospira serovar pomona. In order to obtain internal fragments for protein sequencing, enzymatic and chemical digestion was performed. The enzyme clostripain was used to digest the proteins 32 and 45 kDa. In situ digestion of 40 kDa molecular weight protein was accomplished using cyanogen bromide. The 32 kDa protein generated two fragments, one of 21 kDa and another of 10 kDa that yielded five residues. A fragment of 24 kDa that yielded nineteen residues of amino acids was obtained from 45 kDa protein. A fragment with a molecular weight of 20 kDa, yielding a twenty amino acids sequence from the 40 kDa protein.