Interaction of adjacent amino acid residues in the structured regions of globular proteins

Distribution of the pairs of amino acids i, i + 1 in alpha-helical, beta-sheet and random coil regions from 46 globular proteins comprising 8115 amino acid residues was analyzed. Statistical analysis of the data excludes null hypothesis about random pairing of the amino acid residues i, i + 1 in beta-sheet and random coil configurations. The distribution of the amino acid pairs, i, i + 1 in alpha-helical configurations does not differ from the random pairing.     

Analysis of distributions of amino acids and amino acid pairs in human tumour necrosis factor precursor and its eight mutations according to a random mechanism

In this study, we use the random principle to analyse the distributions of amino acids and amino acid pairs in human tumour necrosis factor precursor (TNF-!) and its eight mutations, to compare the measured distribution probability with the theoretical distribution probability and to rank the measured distribution probability against the theoretical distribution probability. In this way, we can suggest that distributions with a high random rank should not be deliberately evolved and conserved and those with a low random rank should be deliberately evolved and conserved in human TNF-!. An increased distribution probability in a mutation means probabilistically that the mutation is more likely to occur spontaneously, whereas a decreased distribution probability in a mutation means probabilistically that the mutation is less likely to occur spontaneously and perhaps is more related to a certain cause. The results, for example, show that the distributions of 30% of the amino acids are identical with their probabilistic simplest distributions, and the distributions of some of the remaining amino acids are very close to their probabilistic simplest distributions. With respect to probabilities of distributions of amino acids in mutations, the results show that mutations lead to an increase in eight probabilities, which are thus more likely to occur. Eight probabilities decrease and are thus less likely to occur. With respect to the random ranks against the theoretical probabilities of distributions of amino acids, the results show that mutations lead to an increase in seven and a decrease in seven probabilities, with two probabilities unchanged.     

Changes in neutral amino acid transport activity in myeloid leukemia cells differentiated by lipopolysaccharide

M1 cells derived from mouse myeloid leukemia have been reported to differentiate to macrophage-like cells upon treatment with substances such as lipopolysaccharide. Previously we found that in mouse peritoneal macrophages most of the neutral amino acids were taken up through a unique Na+-independent system. In this paper we have investigated the neutral amino acid transport in M1 cells and in those treated with lipopolysaccharide. In M1 cells serine, alanine and proline were taken up mainly by Na+-dependent transport systems, and leucine was largely transported by a Na+-independent system. By treating the cells with lipopolysaccharide, the activities of the Na+-dependent systems markedly decreased, whereas the activity of the Na+-independent system was little affected. The amino acid concentrations in the cells and the culture medium were measured. As a whole, the intracellular to extracellular distribution ratios for neutral amino acids that are preferred substrates for Na+-dependent systems were decreased on lipopolysaccharide treatment, whereas those for amino acids that are mainly transported by a Na+-independent system were slightly increased. From these results we conclude that M1 cells treated with lipopolysaccharide tend to differentiate to macrophage-like cells with respect to the neutral amino acid transport.     

Patterns of temperature adaptation in proteins from the bacteria Deinococcus radiodurans and Thermus thermophilus

Asymmetrical patterns of amino acid substitution in proteins of organisms living at moderate and high temperatures (mesophiles and thermophiles, respectively) are generally taken to indicate selection favoring different amino acids at different temperatures due to their biochemical properties. If that were the case, comparisons of different pairs of mesophilic and thermophilic taxa would exhibit similar patterns of substitutional asymmetry. A previous comparison of mesophilic versus thermophilic Methanococcus with mesophilic versus thermophilic Bacillus revealed several pairs of amino acids for which one amino acid was favored in thermophilic Bacillus and the other was favored in thermophilic Methanococcus. Most of this could be explained by the higher G+C content of the DNA of thermophilic Bacillus, a phenomenon not seen in the Methanococcus comparison. Here, I compared the mesophilic bacterium Deinococcus radiodurans and its thermophilic relative Thermus thermophilus, which are similar in G+C content. Of the 190 pairs of amino acids, 83 exhibited significant substitutional asymmetry, consistent with the pervasive effects of selection. Most of these significantly asymmetrical pairs of amino acids were asymmetrical in the direction predicted from the Methanococcus data, consistent with thermal adaptation resulting from universal biochemical properties of the amino acids. However, 12 pairs of amino acids exhibited asymmetry significantly different from and in the opposite direction of that found in the Methanococcus comparison, and 21 pairs of amino acids exhibited asymmetry that was significantly different from that found in the Bacillus comparison and could not be explained by the greater G+C content in thermophilic Bacillus. This suggests that selection due to universal biochemical properties of the amino acids and differences in G+C content are not the only causes of substitutional asymmetry between mesophiles and thermophiles. Instead, selection on taxon-specific properties of amino acids, such as their metabolic cost, may play a role in causing asymmetrical patterns of substitution.     

Folding of peptide chains during formation of the spatial structure of protein molecules is determined by the code of amino acid codon roots]

Basing on the analysis of a large number of protein sequences (Cserzo M., Simon I., 1989), the structure of the amino acid nearest neighbour pair whose occurrence has a maximal positive deviation from the mean statistical value, is shown to correspond in most cases to the code of the amino acid codon roots. It reveals particularly amino acid pairs in n and n+5 positions of polypeptide chains. Amino acids belonging to A/U family contribute mostly to the folding of peptide chains.     

Sequence analysis of an artificial family of RNA-binding peptides

Diverse peptide sequences recognizing the lambda boxB RNA hairpin were previously isolated from a library encoding the 22-residue lambda N peptide with random amino acids at positions 13-22 using mRNA display. We have statistically analyzed amino acid distributions in 65 unique sequences from rounds 11 and 12 of this selection and evaluated the resulting structural and functional predictions by alanine-scanning mutagenesis and circular dichroism spectrometry. This artificial sequence family has a consensus structure that continues the bent alpha helix of lambda N up to position 17 when bound to lambda boxB. A charge pair (E(14)R(15)) and hydrophobic patch (A(21)L(22) or V(21)L(22)) have important functional roles in this context. Notably, amino acid covariance reveals six specific pairs of random region positions with >95% significant linkage and strong overall helical (i+1, i+3, and i+4) couplings. The covariance analysis suggests that (1) the sequence context of every residue in each insert has been optimized, (2) selected sequences are local optima on a rugged fitness landscape, and (3) it is possible to detect more subtle structural features with artificial protein sequence families than natural homologs. Our results provide a framework for investigating the structures of in vitro selected proteins by functional minimization, reselection, and covariance analysis.     

Amino acid stimulation of ATP cleavage by two Ehrlich cell membrane preparations in the presence of ouabain.

Two membrane fractions prepared from the Ehrlich ascites-tumor cell show non-identical stimulatory responses to certain amino acids in their Mg+2 -dependent activity to cleave ATP, despite the presence of ouabain and the absence of Na+ or K+. The first of these, previously described, shows little (Na+ + K+)-ATPase activity, and is characteristicallly stimulated by the presence of certain diamino acids with low pK2, and at pH values suggesting that the cationic forms of these amino acids are effective. The evidence indicates that these effects are not obtained through occupation of the kinetically discernible receptor site serving characteristically for the uphill transport of these amino acids into the Ehrlich cell. The second membrane preparation was purified with the goal of concentrating the (Na+ +K+)-ATPase activity. It also is stimulated by the model diamino acid, 4-amino-1-methylpiperidine-4-carboxylic acid, and several ordinary amino acids. The diamino acids were most effective at pH values where the neutral zwitterionic forms might be responsible. Among the optically active amino acids tested, the effects of ornithine and leucine were substantially stronger for the L than for the D isomers. The list of stimulatory amino acids again corresponds poorly to any single transport system, although the possibility was not excluded that stimulation might occur for both preparations by occupation of a membrane site which ordinarily is kinetically silent in the transport sequence. The high sensitivity to deoxycholate and to dicyclohexylcarbodiimide of the hydrolytic activity produced by the presence of L-ornithine and 4-amino-1-methyl-piperidine-4-carboxylic acid suggests that the stimulatory effect is not merely a general intensification of the background Mg+ -dependent hydrolytic activity.     

A survey of amino acid side-chain interactions in 21 proteins

Based on the atomic co-ordinate data for 21 representative proteins, the frequencies of long-range interactions between side-chain groups and 15 different types of side-chain atoms have been determined. The observed frequencies are compared to the results expected for random association in order to define a scale of relative affinities. Thirty-five residue-atom pairs exhibit frequencies of interaction that differ by at least 50% from the expected values. The amino acids tend to fall into three classes: non-polar, neutral and polar amino acids. The data are regrouped in a different way to determine the average affinity of each amino acid side-chain group for all other types of side-chain groups. Fourteen side-chain pairs have at least 50% fewer interactions than expected, while 21 side-chain pairs have at least 50% more interactions than expected. Unusual patterns of association are discussed and compared with current ideas about the organization of protein structure.     

Systematic variation of amino acid substitutions for stringent assessment of pairwise covariation

During protein evolution, amino acids change due to a combination of functional constraints and genetic drift. Proteins frequently contain pairs of amino acids that appear to change together (covariation). Analysis of covariation from naturally occurring sets of orthologs cannot distinguish between residue pairs retained by functional requirements of the protein and those pairs existing due to changes along a common evolutionary path. Here, we have separated the two types of covariation by independently recombining every naturally occurring amino acid variant within a set of 15 subtilisin orthologs. Our analysis shows that in this family of subtilisin orthologs, almost all possible pairwise combinations of amino acids can coexist. This suggests that amino acid covariation found in the subtilisin orthologs is almost entirely due to common ancestral origin of the changes rather than functional constraints. We conclude that naturally occurring sequence diversity can be used to identify positions that can vary independently without destroying protein function.     

Determination of amino acid pairs in human p53 protein sensitive to mutations/variants by means of a random approach

This is the continuation of our studies using random approaches to analyse the p53 protein family. In this data-based theoretical analysis, we use the random approach to analyse the amino acid pairs in human p53 protein in order to determine which amino acid pairs are more sensitive to 190 human p53 mutations/variants. The rationale of this study is based on our hypothesis and findings that a harmful mutation is more likely to occur at randomly unpredictable amino acid pairs, and a harmless mutation is more likely to occur at randomly predictable amino acid pairs. This is because we argue that the randomly predictable amino acid pairs should not be deliberately evolved, whereas the randomly unpredictable amino acid pairs should be deliberately evolved with a connection to protein function. The results show, for example, that 93.16% of 190 mutations/variants occur at randomly unpredictable amino acid pairs. Thus, the randomly unpredictable amino acid pairs are more sensitive to mutations/variants in human p53 protein. The results also suggest that the human p53 protein has a tendency for the occurrence of mutation/variants.     

Patterns of temperature adaptation in proteins from Methanococcus and Bacillus

It has long been known that amino acid substitutions in proteins of organisms living at moderate and high temperatures (mesophiles and thermophiles, respectively) are not all symmetrical; for example, more aligned sites have lysine in mesophiles and arginine in thermophiles than have the opposite pattern. This is generally taken to indicate that certain amino acids are favored over others by selection at different temperatures. Previous comparisons of protein sequences from mesophiles and thermophiles have used relatively small numbers of sequences from a diverse array of species, meaning that only the most common amino acid substitutions could be examined and any taxon-specific patterns would be obscured. Here, we compare a large number of proteins between mesophiles and thermophiles in the archaeal genus Methanococcus and the bacterial genus Bacillus. Each genus exhibits dramatically asymmetrical substitution patterns for many pairs of amino acids. There are several pairs of amino acids for which one amino acid is favored in thermophilic Bacillus and the other is favored in thermophilic Methanococcus; this appears to result from the higher G + C content of the DNA of thermophilic Bacillus, a complication not seen in Methanococcus.     

Estimation of amino acid pairs sensitive to variants in human phenylalanine hydroxylase protein by means of a random approach

In this data-based theoretical analysis, we use a random approach to estimate amino acid pairs in human phenylalanine 4-hydroxylase (PAH) protein in order to determine which amino acid pairs are more sensitive to 187 variants in human PAH protein. The rationale of this study is based on our hypothesis and previous findings that the harmful variants are more likely to occur at randomly unpredictable amino acid pairs rather than at randomly predictable pairs. This is reasonable to argue as randomly predictable amino acid pairs are less likely to be deliberately evolved, whereas randomly unpredictable amino acid pairs are probably deliberately evolved in connection with protein function. 94.12% of 187 variants occurred at randomly unpredictable amino acid pairs, which accounted for 71.84% of 451 amino acid pairs in human PAH protein. The chance of a variant occurring is five times higher in randomly unpredictable amino acid pairs than in predictable pairs. Thus, randomly unpredictable amino acid pairs are more sensitive to variance in human PAH protein. The results also suggest that the human PAH protein has a natural tendency towards variants.     

全雌系苦瓜ACC合成酶基因cDNA克隆及序列分析

以全雌系苦瓜‘X-Hei-d-d’花蕾为材料,根据已报道ACC合成酶(1-aminocyclopropane-1-carboxylic acid synthase,ACS)保守氨基酸序列设计简并引物,采用RT-PCR技术及序列拼接,获得了全雌系苦瓜ACS基因cDNA序列,命名为Mc-ACS4(GenBank登录号:FJ459814)。该序列包含一个1 455 bp的完整开放阅读框,编码484个氨基酸,具有7个保守区;系统进化上与普通苦瓜ACS基因首先聚类,同源性达99%,二者仅有2个氨基酸差异,推测可能与全雌系苦瓜性别分化有关。     

Graphic identification of conservative and variable segments in the amino acid sequences of homologous proteins

An algorithm is presented for localizing variable and constant regions in homologous protein sequences. A set of aligned protein sequences is divided into two groups consisting of m and n sequences. Each group contains sequences of most related species. Value of the position dissimilarity of proteins from different groups of m and n sequences is defined as a number of failures to coincide in comparison with all possible mXn pairs of amino acid residues in the position (each from different group) divided by mXn. The position dissimilarity value of m protein sequences within a group is defined as the number of failures to coincide in comparison with all possible mX X(m-1)/2 pairs of amino acid residues divided by mX(m-1)/2. Ten position average of dissimilarity values is plotted vs. the first position number. Area of the figure included between the profile of dissimilarity values and its mean value line characterizes the overall irregularity of amino acid substitutions along the protein sequences. If the area value is greater than the average area for 1000 random profile by more than two standard deviation units, the profile extrema containing the "surplus" of area are cut off. The cut off stretches are likely to be variable and constant regions. In case of "between groups" comparisons it is found that the overall irregularity of amino acid substitutions is very high for all considered families of proteins; phospholipases A2, aspartate aminotransferases, alpha-subunits of Na+,K(+)-ATPase, L- and M-subunits of photosynthetic bacteria photoreaction centre, human rhodopsins.     

Determination of amino acid pairs sensitive to variants in human beta-glucocerebrosidase by means of a random approach

In this data-based theoretical analysis, we use the random approach to analyse the amino acid pairs in human beta-glucocerebrosidase in order to determine which amino acid pairs are more sensitive to 109 variants from missense mutant human glucocerebrosidase. The rationale of this study is based on our hypothesis and findings that the harmful variants are more likely to occur at randomly unpredictable amino acid pairs and the non-harmful variants are more likely to occur at randomly predictable amino acid pairs. This is because we argue that the randomly predictable amino acid pairs should not be deliberately evolved, whereas the randomly unpredictable amino acid pairs should be deliberately evolved with connection of protein function. The results show, for example, that 93.58% of 109 variants occur at randomly unpredictable amino acid pairs, which account for 71.40% of amino acid pairs in glucocerebrosidase, and the chance of occurrence of the variant is about 4.4 times higher in randomly unpredictable amino acid pairs than in predictable pairs. Hence the randomly unpredictable amino acid pairs are more sensitive to variants in human glucocerebrosidase. The results also suggest that human glucocerebrosidase has a natural tendency to variants.     

The structure of the gene for mouse filaggrin and a comparison of the repeating units

Filaggrins are an important class of intermediate filament-associated proteins that are involved in the organization of keratin filaments in the terminal stages of mammalian epidermal differentiation. Filaggrins are initially synthesized as very large polyprotein precursors consisting of many tandemly arranged repeats that are later liberated by proteolytic processes to yield many copies of the functional protein. We have recently characterized a cDNA clone to mouse filaggrin (Rothnagel, J. A., Mehrel. T., Idler, W. W., Roop, D. R., and Steinert, P. M. (1987) J. Biol. Chem. 262, 15643-15648) which encodes a 750-base pair (250-amino acid) repeating element having properties consistent with a filaggrin molecule. Southern blot analysis of total mouse DNA and the mouse gene isolated from a cosmid library (cosmid clone cFM6.1A2) has also revealed a repeat length of about 750 base pairs. The cosmid clone contains most of the mouse filaggrin gene, but it is missing the 5'-noncoding sequences and possibly some coding sequences as well. We report here that cosmid clone cFM6.1A2 contains 20 filaggrin repeats and 15,213 base pairs of coding sequences. Sequence analysis of this clone has revealed at least two different types of repeating element. Type B has a repeat length of 750 base pairs (250 amino acids), whereas type A is 765 base pairs (255 amino acids) long and contains an additional five amino acids inserted next to an acidic sequence that delineates the amino and carboxyl termini of the filaggrin repeats. It is supposed that these additional five amino acids may alter the proteolytic sensitivity of the acidic linker sequence, thereby affecting the processing of the precursor. The random distribution of the two types of repeats in the precursor indicates that the mouse filaggrin gene arose by a complicated series of duplications and/or rearrangements.     

Exceptional pairs of amino acid neighbors in alpha-helices

Amino acids seem to have specific preferences for various locations in alpha-helices. These specific preferences, called singlet local propensity (SLP), have been determined by calculating the preference of occurrence of each amino acid in different positions of the alpha-helix. We have studied the occurrence of amino acids, single or pairs, in different positions, singlet or doublet, of alpha-helices in a database of 343 non-homologous proteins representing a unique superfamily from the SCOP database with a resolution better than 2.5 A from the Protein Data Bank. The preference of single amino acids for various locations of the helix was shown by the relative entropy of each amino acid with respect to the background. Based on the total relative entropy of all amino acids occurring in a single position, the N(cap) position was found to be the most selective position in the alpha-helix. A rigorous statistical analysis of amino acid pair occurrences showed that there are exceptional pairs for which, the observed frequency of occurrence in various doublet positions of the alpha-helix is significantly different from the expected frequency of occurrence in that position. The doublet local propensity (DLP) was defined as the preference of occurrences of amino acid pairs in different doublet positions of the alpha-helix. For most amino acid pairs, the observed DLP (DLP(O)) was nearly equal to the expected DLP (DLP(E)), which is the product of the related SLPs. However, for exceptional pairs of amino acids identified above, the DLP(O) and DLP(E) values were significantly different. Based on the relative values of DLP(O) and DLP(E), exceptional amino acid pairs were divided into two categories. Those, for which the DLP(O) values are higher than DLP(E), should have a strong tendency to pair together in the specified position. For those pairs which the DLP(O) values are less than DLP(E), there exists a hindrance in neighboring of the two amino acids in that specific position of the alpha-helix. These cases have been identified and listed in various tables in this paper. The amount of mutual information carried by the exceptional pairs of amino acids was significantly higher than the average mutual information carried by other amino acid pairs. The average mutual information conveyed by amino acid pairs in each doublet position was found to be very small but non-zero.     

Regularities in the primary structure of proteins

In this paper the latest protein database consisting of more than a million amino acids is analyzed to characterize the short range regularities in the primary structure. The amino acid distributions along the polypeptide chain and among the proteins have been studied first. Their influence on the amino acid pair statistics was taken into account. We are primarily interested in the distances of the covalent structure, where the amino acid pair frequencies show non-random characters. The amino acid pairs separated by at least 20 residues in the covalent structure exhibit an exact Gaussian distribution. We found that there is a range of non-random pairing in the covalent structure. We conclude that the pair preference characters are different for each of the 20 x 20 amino acid pairs. The range of the non-random pairing varies from pair to pair, and in most cases it does not extend beyond the 9th neighbour. The preferences of a certain pair in a certain position can not be derived from the character of that pair in another position. The preference values of 400 amino acid pairs are listed for up to the pairs in 9th neighbour position. Some fields of potential application of these data have also been discussed.     

Multi-modality of pI distribution in whole proteome

Multi-modality of pI distribution is a common feature in different whole proteomes. Some researchers considered it relate to the proteins with different subcellular locations, indicating the result of natural selection. We explored the pI distribution of predicted proteomes (including animals, plants, bacterium, archaeans) and random proteome [random protein sequences constructed according to the special amino acid composition and molecular weight (MW) distribution of human predicted proteome]. Our results suggest that the multi-modality is the result of discrete pK(R) values for different amino acids. Amino acid composition and MW distribution of a proteome also contributes to the specific pI distribution. Although protein subcellular location was related to pI value, our analyses revealed that comparing with the random proteome, neither the multi-modality phenomenon nor the distribution bias of pI values is caused by subcellular location. It seems that the multi-modality distribution is just a mathematical fun. The blank region near the neutral pI was caused by the absence of amino acids with neutral pK(R), and suggests that the selection of amino acids with ionizable side chain might be restricted by the requirement for a special pH environment during the origin of life. From this point of view, the special distribution was the result of natural selection.