Two types of amino acid substitutions in protein evolution

Summary The frequency of amino acid substitutions, relative to the frequency expected by chance, decreases linearly with the increase in physico-chemical differences between amino acid pairs involved in a substitution. This correlation does not apply to abnormal human hemoglobins. Since abnormal hemoglobins mostly reflect the process of mutation rather than selection, the correlation manifest during protein evolution between substitution frequency and physico-chemical difference in amino acids can be attributed to natural selection. Outside of abnormal proteins, the correlation also does not apply to certain regions of proteins characterized by rapid rates of substitution. In these cases again, except for the largest physico-chemical differences between amino acid pairs, the substitution frequencies seem to be independent of the physico-chemical parameters. The limination of the substituents involving the largest physicochemical differences can once more be attributed to natural selection. For smaller physico-chemical differences, natural selection, if it is operating in the polypeptide regions, must be based on parameters other than those examined.     

Accumulation pattern of amino acid substitutions in protein evolution

Summary A simple method for the evolutionary analysis of amino acid sequence data is presented and used to examine whether the number of variable sites (NVS) of a protein is constant during its evolution. The NVSs for hemoglobin and for mitochondrial cytochrome c are each found to be almost constant, and the ratio between the NVSs is close to the ratio between the unit evolutionary periods. This indicates that the substitution rate per variable site is almost uniform for these proteins, as the neutral theory claims. An advantage of the present analysis is that it can be done without knowledge of paleontological divergence times and can be extended to bacterial proteins such as bacterial c-type cytochromes. It is suggested that the NVS of cytochrome c has been almost constant even over the long period (ca. 3.0 billion years) of bacterial evolution but that at least two different substitution rates are necessary to describe the accumulated changes in the sequence. This two clock interpretation is consistent with fossil evidence for the appearance times of photosynthetic bacteria and eukaryotes.     

Effect of conserved intersubunit amino acid substitutions on Hfq protein structure and stability

Hfq is a thermostable RNA-binding bacterial protein that forms a uniquely shaped homohexamer. Based on sequence and structural similarity, Hfq belongs to the like-Sm (LSm) protein family. In spite of a rather high degree of homology between archaeal and eukaryotic LSm proteins, their quaternary structure is different, usually consisting of five to eight monomers. In this work, the importance of conserved intersubunit hydrogen bonds for the Hfq spatial organization was tested. The structures and stabilities for the Gln8Ala, Asn28Ala, Asp40Ala, and Tyr55Ala Hfq mutants were determined. All these proteins have the same hexamer organization, but their stability is different. Elimination of a single intersubunit hydrogen bond due to Gln8Ala, Asp40Ala, and Tyr55Ala substitutions results in decreased stability of the Hfq hexamer. Tyr55Ala Hfq as well as the earlier studied His57Ala Hfq has reduced protein thermostability, which seems to correspond to an opening of the protein hydrophobic core.     

Prediction of protein disorder on amino acid substitutions

Intrinsically disordered regions of proteins are known to have many functional roles in cell signaling and regulatory pathways. The altered expression of these proteins due to mutations is associated with various diseases. Currently, most of the available methods focus on predicting the disordered proteins or the disordered regions in a protein. On the other hand, methods developed for predicting protein disorder on mutation showed a poor performance with a maximum accuracy of 70%. Hence, in this work, we have developed a novel method to classify the disorder-related amino acid substitutions using amino acid properties, substitution matrices, and the effect of neighboring residues that showed an accuracy of 90.0% with a sensitivity and specificity of 94.9 and 80.6%, respectively, in 10-fold cross-validation. The method was evaluated with a test set of 20% data using 10 iterations, which showed an average accuracy of 88.9%. Furthermore, we systematically analyzed the features responsible for the better performance of our method and observed that neighboring residues play an important role in defining the disorder of a given residue in a protein sequence. We have developed a prediction server to identify disorder-related mutations, and it is available at http://www.iitm.ac.in/bioinfo/DIM_Pred/.     

Volume and polarity changes accompanied by amino acid substitutions in protein evolution.

We evaluated the volume and polarity changes accompanied by amino acid substitutions along branches of the phylogenetic trees of cytochrome c, myoglobin and hemoglobin alpha and beta chains. In most cases the volume changes accompanied by the substitutions were found to be much larger than the volume of cavities existing in the interior of X-ray-analysed proteins. This implies that the interior of the proteins is very flexible and the necessary space for a larger amino acid residue substitution can be provided by adjusting nearby structures. Also, the volume and polarity changes are not particularly dependent on whether the substituted site is located in the exterior or interior of the proteins. This result supports the concept of the covarions by Fitch and Markowitz, when combined with the known fact that the exterior sites are more variable than the interior ones during protein evolution.     

Identification of amino acid latent periodicity within 94 protein families.

Here, we have applied information decomposition, cyclic profile alignment, and noise decomposition techniques to search for latent repeats within protein families of various functions. We have identified 94 protein families with a family-specific periodicity. In each case, the periodic element was found in greater than 70% of family members. Latent periodicity profiles with specific length and signature were obtained in each case. The possible relationship between the periodic elements thus identified and the evolutionary development of the protein families are discussed with specific reference to the possibility that there is a correlation between the periodic elements and protein function.     

The evolution of the structures of amino acid families]

Natural amino acids possessing common antiamino acids are divided into groups and families according to the genetic code algorithm a-n-n-a (amino acid-codon-anticodon-antiamino acid). In an attempt to study structural evolution of amino acid families, artificial genetic code models were constructed. It is suggested that after inclusion of asparaginase and glutamine into the coding system, one of the two natural amino acid families is split into two parts ("half-families").     

Reciprocal amino acid substitutions in the evolution of homologous peptides

Amino acid sequences of peptides are often inferred from their amino acid compositions by comparison with homologous peptides of known sequence. The probabilities are considered that by such an approach errors are made due to the occurrence of balanced double changes, i.e. reciprocal substitutions, between two homologous peptides of identical compositions. Formulae are derived for the calculation of these probabilities, depending on peptide length and evolutionary distance. However, such calculations requiring too much computer time, the probabilities for reciprocal substitutions are estimated by simulation of evolutionary changes in peptides. It can be concluded from the resulting data that for many purposes the possible errors in amino acid sequences partially inferred from amino acid compositions are acceptably small.     

Approaches to predicting effects of single amino acid substitutions on the function of a protein

The relative activities of 313 mutants of the gene V protein of bacteriophage f1, assayed in vivo, have been used to evaluate two approaches to predicting the effects of single amino acid substitutions on the function of a protein. First, we tested methods that only depend on the properties of the wild-type and substituting amino acids. None of the properties or measures of the functional equivalence of amino acids we tested, including the frequency of exchange of amino acids among homologous proteins as well as changes in side-chain size, hydrophobicity, and charge, were found to be more than weakly correlated with the activities of mutants. The principal reason for this poor correlation was found to be that the effect of a particular substitution varies considerably from site to site. We then tested an approach using the activities of several mutants with substitutions at a site to predict the activity of another mutant, and we find that this is a relatively good indicator of whether the other mutant at that site will be functional. A predictive scheme was developed that combines the weak information from the models depending on the properties of the wild-type and substituting amino acids with the stronger information from the tolerance of a site to substitution. Although this scheme requires no knowledge of the structure of a mutant protein, it is useful in predicting the activities of mutants.     

Conformational effects of selected cancer-related amino acid substitutions in the p53 protein.

The tumor suppressor gene p53 has been identified as the most frequent target of genetic alterations in human cancers. Cancer-related mutations in the human p53 protein tend to cluster in four of the five highly conserved domains of the protein, and, in particular, in the central region of domain IV from residues 241 to 253. Using conformational energy analysis based on ECEPP (Empirical Conformational Energies for Polypeptides Program), we have determined the preferred three dimensional structures for this tridecapeptide sequence for the human wild-type p53 protein and four cancer-related mutant p53 proteins (Ala 245, Ile 246, Trp 248, Ser 249). The results show that the mutant peptides adopt conformations that are distinctly different from that of the wild-type peptide. These results are consistent with experimental conformational studies demonstrating altered detectability of antigenic epitopes in murine wild-type and mutant p53 proteins. These results suggest that the oncogenic effects of human mutant p53 proteins may be mediated by distinct local conformational changes in the protein.     

A new method to model membrane protein structure based on silent amino acid substitutions.

The importance of accurately modeling membrane proteins cannot be overstated, in lieu of the difficulties in solving their structures experimentally. Often, however, modeling procedures (e.g., global searching molecular dynamics) generate several possible candidates rather then pointing to a single model. Herein we present a new approach to select among candidate models based on the general hypothesis that silent amino acid substitutions, present in variants identified from evolutionary conservation data or mutagenesis analysis, do not affect the stability of a native structure but may destabilize the non-native structures also found. The proof of this hypothesis has been tested on the alpha-helical transmembrane domains of two homodimers, human glycophorin A and human CD3-zeta, a component of the T-cell receptor. For both proteins, only one structure was identified using all the variants. For glycophorin A, this structure is virtually identical to the structure determined experimentally by NMR. We present a model for the transmembrane domain of CD3-zeta that is consistent with predictions based on mutagenesis, homology modeling, and the presence of a disulfide bond. Our experiments suggest that this method allows the prediction of transmembrane domain structure based only on widely available evolutionary conservation data.     

Nature of amino acid substitutions in homologous proteins during evolution

Replacement substitutions of mitochondrial cytochrome $\text{c}$ and α- and β-chains of haemoglobin have been studied by considering the structural similarity among amino acid residues at the secondary and tertiary structural levels. Secondary structural similarity explains ～70% while tertiary structural similarity explains ～50% of observed replacements for most of the cases. These structural similarities could not account for all the replacement substitutions. The study was extended to consider the composition of codons, and the chemical nature and polarity of the replacing and replaced residues. These also could not individually account for all the affected replacements. In general, no property of amino acid residues is conserved for substitutions occurring at any single position during evolution of proteins.     

Residual structure in large fragments of staphylococcal nuclease: effects of amino acid substitutions

In an attempt to develop a model of the denatured state of staphylococcal nuclease that can be analyzed experimentally under physiological conditions, a series of four large fragments of this small protein which extend from residues 1 to 103, 1 to 112, 1 to 128, and 1 to 136 have been generated through the overexpression of nuclease genes containing stop codons at defined positions. Large amounts of protein fragments were accumulated in induced cells and were purified by carrying out all fractionation steps in the presence of 6 M urea. The far-ultraviolet circular dichroism spectra of all four fragments suggested the presence of small to moderate amounts of residual structure. When the CD spectra were monitored as a function of concentrations of the tight-binding ligands Ca2+ and thymidine 3',5'-bisphosphate and the known affinity constants for wild-type nuclease (1-149) were used, apparent equilibrium constants of 160 and 2000 for the reversible denaturation reaction for fragments 1-136 and 1-128, respectively, were estimated. Four single and two double mutations, all of which exhibit unusual behavior in the full-length protein on solvent denaturation [Shortle, D., & Meeker, A. K. (1986) Proteins: Struct., Funct., Genet. 1, 81-89] and thermal denaturation [Shortle, D., Meeker, A. K., & Freire, E. (1988) Biochemistry 27, 4761-4768], were recombined into the 1-136 and 1-128 fragment expression vectors, and purified mutant fragments were characterized.(ABSTRACT TRUNCATED AT 250 WORDS)     

Effect of amino acid substitutions on conformational stability of a protein

This paper reviews studies on thermostable proteins from thermophilic bacteria and on mutant proteins of human hemoglobin, tryptophan synthase α-subunit of E. coli, T₄ phage lysozyme, and phage λ repressor with respect to the role of the consisting amino acid residues in stabilization of conformation. The stability of a protein is easily affected by single amino acid substitutions, by which the protein undergoes change(s) of one or more of the following: a hydrogen bond, a salt bridge, a hydrophobic interaction, the volume of the residue, a disulfide bond, or the relative position of two aromatic rings.     

Pervasive adaptive protein evolution apparent in diversity patterns around amino acid substitutions in Drosophila simulans

In Drosophila, multiple lines of evidence converge in suggesting that beneficial substitutions to the genome may be common. All suffer from confounding factors, however, such that the interpretation of the evidence-in particular, conclusions about the rate and strength of beneficial substitutions-remains tentative. Here, we use genome-wide polymorphism data in D. simulans and sequenced genomes of its close relatives to construct a readily interpretable characterization of the effects of positive selection: the shape of average neutral diversity around amino acid substitutions. As expected under recurrent selective sweeps, we find a trough in diversity levels around amino acid but not around synonymous substitutions, a distinctive pattern that is not expected under alternative models. This characterization is richer than previous approaches, which relied on limited summaries of the data (e.g., the slope of a scatter plot), and relates to underlying selection parameters in a straightforward way, allowing us to make more reliable inferences about the prevalence and strength of adaptation. Specifically, we develop a coalescent-based model for the shape of the entire curve and use it to infer adaptive parameters by maximum likelihood. Our inference suggests that ～13% of amino acid substitutions cause selective sweeps. Interestingly, it reveals two classes of beneficial fixations: a minority (approximately 3%) that appears to have had large selective effects and accounts for most of the reduction in diversity, and the remaining 10%, which seem to have had very weak selective effects. These estimates therefore help to reconcile the apparent conflict among previously published estimates of the strength of selection. More generally, our findings provide unequivocal evidence for strongly beneficial substitutions in Drosophila and illustrate how the rapidly accumulating genome-wide data can be leveraged to address enduring questions about the genetic basis of adaptation.     

Nonrandom amino acid substitution and estimation of the number of nucleotide substitutions in evolution

Summary A method of estimating the number of nucleotide substitutions from amino acid sequence data is developed by using Dayhoff's mutation probability matrix. This method takes into account the effect of nonrandom amino acid substitutions and gives an estimate which is similar to the value obtained by Fitch's counting method, but larger than the estimate obtained under the assumption of random substitutions (Jukes and Cantor's formula). Computer simulations based on Dayhoff's mutation probability matrix have suggested that Jukes and Holmquist's method of estimating the number of nucleotide substitutions gives an overestimate when amino acid substitution is not random and the variance of the estimate is generally very large. It is also shown that when the number of nucleotide substitutions is small, this method tends to give an overestimate even when amino acid substitution is purely at random.