Integrated displays of aligned amino acid sequences and protein structures

Programs for the Sun-4 workstation permit the combined displayof a table of aligned amino acid sequences of a family of proteins,and a corresponding three-dimensional fold. Interactive facilitiesinclude the ability, to scroll through the sequences, to rotatethe structure and to connect the examination of the sequencesand the structure by selecting a portion of the sequences andautomatically highlighting the corresponding region in the structureand vice versa. These programs are well suited to support applications suchas the investigation of the structural or functional significanceof conserved patterns of amino acids in the sequences of a familyof proteins. Received on July 3, 1990; accepted on January 25, 1991     

Generation of deviation parameters for amino acid singlets, doublets and triplets from three-dimensional structures of proteins and its implications for secondary structure prediction from amino acid sequences

We present a new method, secondary structure prediction by deviation parameter (SSPDP) for predicting the secondary structure of proteins from amino acid sequence. Deviation parameters (DP) for amino acid singlets, doublets and triplets were computed with respect to secondary structural elements of proteins based on the dictionary of secondary structure prediction (DSSP)-generated secondary structure for 408 selected nonhomologous proteins. To the amino acid triplets which are not found in the selected dataset, a DP value of zero is assigned with respect to the secondary structural elements of proteins. The total number of parameters generated is 15,432, in the possible parameters of 25,260. Deviation parameter is complete with respect to amino acid singlets, doublets, and partially complete with respect to amino acid triplets. These generated parameters were used to predict secondary structural elements from amino acid sequence. The secondary structure predicted by our method (SSPDP) was compared with that of single sequence (NNPREDICT) and multiple sequence (PHD) methods. The average value of the percentage of prediction accuracy for αhelix by SSPDP, NNPREDICT and PHD methods was found to be 57%, 44% and 69% respectively for the proteins in the selected dataset. For Β-strand the prediction accuracy is found to be 69%, 21% and 53% respectively by SSPDP, NNPREDICT and PHD methods. This clearly indicates that the secondary structure prediction by our method is as good as PHD method but much better than NNPREDICT method.     

Analysis of crystal structures of aspartic proteinases: On the role of amino acid residues adjacent to the catalytic site of pepsin-like enzymes

To elucidate the role of amino acid residues adjacent to the catalytic site of pepsin-like enzymes, we analyzed and compared the crystal structures of these enzymes, their complexes with inhibitors, and zymogens in the active site area (a total of 82 structures). In addition to the water molecule (W1) located between the active carboxyls and playing a role of the nucleophile during catalytic reaction, another water molecule (W2) at the vicinity of the active groups was found to be completely conserved. This water molecule plays an essential role in formation of a chain of hydrogen-bonded residues between the active site flap and the active carboxyls on ligand binding. These data suggest a new approach to understanding the role of residues around the catalytic site, which can assist the development of the catalytic reaction. The influence of groups adjacent to the active carboxyls is manifested by pepsin activity at pH 1.0. Some features of pepsin-like enzymes and their mutants are discussed in the framework of the approach.     

An interactive graphics program for comparing and aligning nucleic acid and amino acid sequences.

This paper describes a computer program designed to look for similarities between pairs of nucleic or amino acid sequences. The program looks both for segments of perfect identity or for regions where, using a scoring matrix, a minimum value is exceeded. The results of comparisons are presented as a matrix which is displayed on a simple graphics terminal. Use of a graphics terminal allows the user to display the whole of the two sequences in one screenful or to home-in on regions of interest to examine them in more detail. The program is interactive and so the user can easily see the effect of changes to variables and can use inbuilt editing functions to make insertions to produce alignments of the two sequences. These aligned sequences can then be saved on disk files for further processing.     

Prediction of the secondary and tertiary structures of interferon from four homologous amino acid sequences

The secondary and tertiary structures of interferon were predicted from four homologous amino acid sequences. Three methods of secondary structure prediction gave differing results that were interpreted to suggest that there might be four α-helices that are important in the tertiary fold. The validity of this interpretation was assessed by the application of the methods to predict the secondary structures of two proteins known to consist of four α-helices. A possible tertiary model for interferon is then proposed in which the four α-helices pack into a right-handed bundle similar to that observed in several known protein structures. This model was shown to be stereochemically feasible by an α-helix docking algorithm. One of the resultant structures is shown to be compatible with the known disulphide linkages in interferon. Certain residues that are conserved between the different sequences lie near each other in our model and these residues might form a functional site. In the absence of a crystal structure for interferon, a predicted tertiary model will help further structural and functional studies.     

Adenovirus proteinases: comparison of amino acid sequences and expression of the cloned cDNA in Escherichia coli

Adenoviruses (Ad) synthesize serine-center endoproteinases (AdEPs) responsible for maturation cleavages within the virus particle. Many questions regarding these enzymes remain unanswered because previous studies utilized crude cells or viral lysates as the enzyme source. Here, we report on the comparison of the amino acid (aa) sequences of several AdEPs and on the expression of the cDNA of the Ad2Ep in Escherichia coli. The AdEPs consist of about 200 aa and their size is around 23 kDa. Among the seven sequences known, 60% of aa were strictly conserved. The usual serine proteinase active site sequence, GDSGG, is absent. The recombinant Ad2EP, produced by an inducible vector as a protein-A fusion product is capable of autocatalytic cleavage, and of cleaving its natural viral substrates as well as foreign proteins. Therefore, other viral proteins or mammalian specific post-translational modifications are not required for enzyme activity.     

Predicting amino acid residues responsible for enzyme specificity solely from protein sequences

We describe a general, modular method for developing protocols to identify the amino acid residues that most likely define the division of a protein superfamily into two subsets. As one possibility, we use PROBE to gather superfamily members and perform an ungapped alignment. We then use a modified BLOSUM62 substitution matrix to determine the discriminating power of each column of aligned residues. The overall method is particularly useful for predicting amino acids responsible for substrate or binding specificity when no structures are available. We apply our method to three pairs of protein classes in three different superfamilies, and present our results, some of which have been experimentally verified. This approach may accelerate the elucidation of enzymic substrate specificity, which is critical for both mechanistic insights into biocatalysis and ultimate application.     

Properties of soluble fusions between mammalian aspartic proteinases and bacterial maltose-binding protein

The mammalian aspartic proteinases procathepsin D and pepsinogen form insoluble inclusion bodies when expressed in bacteria. They become soluble but nonnative when synthesized as fusions to the carboxy terminus of E. coli maltose-binding protein (MBP). Since these nonnative states of the two aspartic proteinases showed no tendency to form insoluble aggregates, their biophysical properties were analyzed. The MBP portions were properly folded as shown by binding to amylose, but the aspartic proteinase moieties failed to bind pepstatin and lacked enzymatic activity, indicating that they were not correctly folded. When treated with proteinase K, only the MBP portion of the fusions was resistant to proteolysis. The fusion between MBP and cathepsin D had increased hydrophobic surface exposure compared to the two unfused partners, as determined by bis-ANS binding. Ultracentrifugal sedimentation analysis of MBP–procathepsin D and MBP–pepsinogen revealed species with very large and heterogeneous sedimentation values. Refolding of the fusions from 8 M urea generated proteins no larger than dimers. Refolded MBP–pepsinogen was proteolytically active, while only a few percent of renatured MBP–procathepsin D was obtained. The results suggest that MBP–aspartic proteinase fusions can provide a source of soluble but nonnative folding states of the mammalian polypeptides in the absence of aggregation.     

Aromatic amino acids and their derivatives as ligands for the isolation of aspartic proteinases

Affinity chromatography was used to study an interaction of aspartic proteinases with immobilized aromatic amino acids and their derivatives. The following ligands were used: L-tyrosine, 3-iodo-L-tyrosine, 3,5-diiodo-L-tyrosine, L-phenylalanine, p-iodo-L-phenylalanine and N-acetyl-L-phenylalanine. With the exception of the last one, ligands were coupled directly to divinyl sulfone activated Sepharose 4B. For the preparation of immobilized N-acetyl-L-phenylalanine, divinyl sulfone activated Sepharose 4-B with linked ethylene diamine was used. Porcine pepsin was used for the evaluation of the capacity of the prepared affinity carriers. The capacity of the immobilized amino acid derivatives significantly increased in comparison with the non-derivatized amino acids. The prepared immobilized ligands were further used for the separation of human pepsinogens.     

The evolution of proteins from random amino acid sequences. I. Evidence from the lengthwise distribution of amino acids in modern protein sequences

Summary We examine in this paper one of the expected consequences of the hypothesis that modern proteins evolved from random heteropeptide sequences. Specifically, we investigate the lengthwise distributions of amino acids in a set of 1,789 protein sequences with little sequence identity using the run test statistic (r _o) of Mood (1940,Ann. Math. Stat. 11, 367–392). The probability density ofr _o for a collection of random sequences has mean=0 and variance=1 [the N(0,1) distribution] and can be used to measure the tendency of amino acids of a given type to cluster together in a sequence relative to that of a random sequence. We implement the run test using binary representations of protein sequences in which the amino acids of interest are assigned a value of 1 and all others a value of 0. We consider individual amino acids and sets of various combinations of them based upon hydrophobicity (4 sets), charge (3 sets), volume (4 sets), and secondary structure propensity (3 sets). We find that any sequence chosen randomly has a 90% or greater chance of having a lengthwise distribution of amino acids that is indistinguishable from the random expectation regardless of amino acid type. We regard this as strong support for the random-origin hypothesis. However, we do observe significant deviations from the random expectation as might be expected after billions years of evolution. Two important global trends are found: (1) Amino acids with a strong α-helix propensity show a strong tendency to cluster whereas those with β-sheet or reverse-turn propensity do not. (2) Clustered rather than evenly distributed patterns tend to be preferred by the individual amino acids and this is particularly so for methionine. Finally, we consider the problem of reconciling the random nature of protein sequences with structurally meaningful periodic “patterns” that can be detected by sliding-window, autocorrelation, and Fourier analyses. Two examples, rhodopsin and bacteriorhodopsin, show that such patterns are a natural feature of random sequences.     

Correlation between amino acid residues converted by RNA editing and functional residues in protein three-dimensional structures in plant organelles

Background

In plant organelles, specific messenger RNAs (mRNAs) are subjected to conversion editing, a process that often converts the first or second nucleotide of a codon and hence the encoded amino acid. No systematic patterns in converted sites were found on mRNAs, and the converted sites rarely encoded residues located at the active sites of proteins. The role and origin of RNA editing in plant organelles remain to be elucidated.

Results

Here we study the relationship between amino acid residues encoded by edited codons and the structural characteristics of these residues within proteins, e.g., in protein-protein interfaces, elements of secondary structure, or protein structural cores. We find that the residues encoded by edited codons are significantly biased toward involvement in helices and protein structural cores. RNA editing can convert codons for hydrophilic to hydrophobic amino acids. Hence, only the edited form of an mRNA can be translated into a polypeptide with helix-preferring and core-forming residues at the appropriate positions, which is often required for a protein to form a functional three-dimensional (3D) structure.

Conclusion

We have performed a novel analysis of the location of residues affected by RNA editing in proteins in plant organelles. This study documents that RNA editing sites are often found in positions important for 3D structure formation. Without RNA editing, protein folding will not occur properly, thus affecting gene expression. We suggest that RNA editing may have conferring evolutionary advantage by acting as a mechanism to reduce susceptibility to DNA damage by allowing the increase in GC content in DNA while maintaining RNA codons essential to encode residues required for protein folding and activity.     

Identification of latent periodicity in amino acid sequences of protein families

For detection of the latent periodicity of the protein families responsible for various biological functions, methods of information decomposition, cyclic profile alignment, and the method of noise decomposition have been used. The latent periodicity, being specific to a particular family, is recognized in 94 of 110 analyzed protein families. Family specific periodicity was found for more than 70% of amino acid sequences in each of these families. Based on such sequences the characteristic profile of the latent periodicity has been deduced for each family. Possible relationship between the recognized latent periodicity, evolution of proteins, and their structural organization is discussed.     

How different amino acid sequences determine similar protein structures: The structure and evolutionary dynamics of the globins

To determine how different amino acid sequences form similar protein structures, and how proteins adapt to mutations that change the volume of residues buried in their close-packed interiors, we have analysed and compared the atomic structures of nine different globins. The homology of the sequences in the two most distantly related molecules is only 16%.The principal determinants of three-dimensional structure of these proteins are the approximately 59 residues involved in helix to helix and helix to haem packings. Half of these residues are buried within the molecules. The observed variations in the sequence keep the side-chains of buried residues non-polar, but do not maintain their size: the mean variation of the volume among homologous amino acids is 56 ų.Changes in the volumes of buried residues are accompanied by changes in the geometry of the helix packings. The relative positions and orientations of homologous pairs of helices in the globins differ by rigid body shifts of up to 7 Å and 30 °. In order to retain functional activity these shifts are coupled so that the geometry of the residues forming the haem pocket is very similar in all the globins.We discuss the implications of these results for the mechanism of protein evolution.