Code dependent conservation of the physico-chemical properties in amino acid substitutions

The frequency of amino acid replacements in families of typical proteins has been elegantly analyzed by Argyle (1980) showing that the most frequent replacements involve a conservation of the amino acid chemical properties. The cyclic arrangement of the twenty amino acids resulting from the most frequent replacements has been described as an amino acid chemical ring. In this work, a novel amino acid replacement frequency ring is proposed, for which a conservation of over 90% of the most general physico-chemical properties can be deduced. The amino acid chemical similarity ring is also analyzed in terms of the genetic code base probability changes, showing that the discrepancy that exists between the standard deviation value of the amino acid replacement frequency matrix and its respective ideal value is almost equal to that deduced from the corresponding base codon replacement probability matrices. These differences are finally evaluated and discussed in terms of the restrictions imposed by the structure of the genetic code and the physico-chemical dissimilarities between some codons of amino acids which are chemically similar.     

The nature of amino acid substitution in antigenic drift of hemagglutinin N3 and neuraminidase N2 from the influenza virus

The nature of amino acid replacements in 16 drift variants of hemagglutinin H3 subtype and 5 drift variants of neuraminidase N2 subtype of the influenza A virus were studied. The dependences of relative replacement frequencies and relative quantities of frequent replacements upon differences of properties of substituted residues are plotted. In contrast to most of the known proteins, amino acid replacements in hemagglutinin and neuraminidase depend weakly on the physico-chemical parameters of amino acids. For the antigenic determinants studied the replacement frequencies were compared to those calculated according to two models: one for conservative replacements and the other for accidental mutation of the genetic code. The differences in the nature of amino acid replacements are found in four antigenic determinants of hemagglutinin. The replacements in experimentally selected proteins are shown to go beyond limitations of natural variants. The explanations of the reasons of low epidemicity of some strains and ineffective attempt to imitate the natural antigenic drift of viruses by using experimental selection are proposed. The causes of time-limited circulation of H3N2 influenza virus subtype are discussed.     

Conservation of the secondary structure of protein during evolution and the role of the genetic code

In this communication we demonstrate, in a group of modern proteins, following an algorithm described by Argyle (1980), that the ordination of the amino acids in terms of the most frequent substitutions agrees with the conservation of the-helix,-sheet, and-turn formation tendencies during evolution. The same correspondence has been demonstrated for the conservation of the physico-chemical properties in the amino acid substitutions. Both parameters are similar in showing higher correlation with the most frequent amino acid substitutions than with the feasibility of changes at the level of the respective codons.Some kind of restrictions for the expression of the genomic changes, due to the conservation of the secondary structure of proteins and/or the physicochemical properties of the substituted amino acids, could account for the differences found between the distribution of the amino acid substitutions and the most probable codon changes.This work has been partially supported by Departamento de Investigación y Bibliotecas, Universidad de Chile y Fondo Nacional de Investigación Científica y Tecnológica.     

The optimality of the standard genetic code assessed by an eight-objective evolutionary algorithm

Background

The standard genetic code (SGC) is a unique set of rules which assign amino acids to codons. Similar amino acids tend to have similar codons indicating that the code evolved to minimize the costs of amino acid replacements in proteins, caused by mutations or translational errors. However, if such optimization in fact occurred, many different properties of amino acids must have been taken into account during the code evolution. Therefore, this problem can be reformulated as a multi-objective optimization task, in which the selection constraints are represented by measures based on various amino acid properties.

Results

To study the optimality of the SGC we applied a multi-objective evolutionary algorithm and we used the representatives of eight clusters, which grouped over 500 indices describing various physicochemical properties of amino acids. Thanks to that we avoided an arbitrary choice of amino acid features as optimization criteria. As a consequence, we were able to conduct a more general study on the properties of the SGC than the ones presented so far in other papers on this topic. We considered two models of the genetic code, one preserving the characteristic codon blocks structure of the SGC and the other without this restriction. The results revealed that the SGC could be significantly improved in terms of error minimization, hereby it is not fully optimized. Its structure differs significantly from the structure of the codes optimized to minimize the costs of amino acid replacements. On the other hand, using newly defined quality measures that placed the SGC in the global space of theoretical genetic codes, we showed that the SGC is definitely closer to the codes that minimize the costs of amino acids replacements than those maximizing them.

Conclusions

The standard genetic code represents most likely only partially optimized systems, which emerged under the influence of many different factors. Our findings can be useful to researchers involved in modifying the genetic code of the living organisms and designing artificial ones.

     

Molecular Variation and Possible Lines Of Evolution of Peptide and Protein Hormones

The widespread distribution of certain steroids and amino acidderivatives with hormonal properties is considered evidencein support of the dictum that "it is not the hormones that change,but rather the uses to which they are put." However, analysesof the distributions, biological activities, immunological cross-reactivities,and sequences of amino acids of five representative peptideand protein hormones or groups of hormones—lactogenichormone, growth hormone, the corticotropin-MSH-ß lipotropinfamily, insulin, and the neurohypophysial hormones—supporta concept of change and of molecular evolution of these polypeptidicmolecules. When analyzed in terms of the genetic code, the aminoacid interchanges which have been revealed by determinationof sequences of amino acids can, most often, be explained bysingle base mutations in the appropriate codons. In two instanceswhere two base mutations within a single codon are required,intermediate replacements of amino acid have been suggested;one of these would lead to a 2-ALA-ß MSH, and theother to a 4-PRO, 8-ILE oxytocin.     

Stereochemical origins of the genetic code

The origin of the genetic code may be attributed to a postulated prebiological stereochemistry in which amino acid dimers, the trans -R,R'-diketopiperazines, interacted with prototype codon and anticodon nucleotide sequences. An intricately coupled stereochemistry is formulated which displays a binary logic for amino acid-codon recognition. It is shown that the diketopiperazine ring system can be inserted between any terminal pair of base paired nucleotides in a codon-anticodon structure with exact registration of complementary hydrogen bonding functional groups. This yields a codon-dimer-anticodon structure in which each amino acid residue is projected towards and interacts with a particular sequence of vicinal nucleotides on either codon or anticodon. The projection direction and the sequence of nucleotides encountered is a strongly coupled function of the choice of codon terminal nucleotide and the handedness of the amino acid. The reciprocal chemical nature of the complementary base pairs drives the selection of dimers containing quite dissimilar and chirally opposed amino acids. Application of the stereochemical model to the in vivo system leads to a general correlation for amino acid-codon assignments. The genetic code is restated in terms of the dimers selected. The profound symmetry of the code is elucidated and this proves useful for correlative and predictive purposes.     

Logic of the genetic code: Conservation of long-range interactions among amino acids as a prime factor

Any statement on the optimality of the existing code ought to imply that this code is ideal for conserving a certain hierarchy of properties while implying that other codes may have been better suited for conservation of other hierarchies of properties. We have evaluated the capability of mutations in the genetic code to convert one amino acid into another in relation to the consequent changes in physical properties of those amino acids. A rather surprising result emerging from this analysis is that the genetic code conserves long-range interactions among amino acids and not their short-range stereochemical attributes. This observation, based directly on the genetic code itself and the physical properties of the 20 amino acids, lends credibility to the idea that the genetic code has not originated by a frozen accident (the null hypothesis rejected by these studies) nor are stereochemical attributes particularly useful in our understanding of what makes the genetic code ‘tick’. While the argument that replacement of, say, an aspartate by a glutamate is less damaging than replacement by arginine makes sense, in order to subject such statements to rigorous statistical tests it is essential to define what constitutes a random sample for the genetic code. The present investigation describes one possible specification. In addition to obvious statistical considerations of testing hypotheses, this procedure points to the more exciting notion that alternative codes may have existed.     

Protein evolution within and between species

Protein evolution can be seen as the successive replacement of amino acids by other amino acids. In general, it is a very slow process which is triggered by point mutations in the nucleotide sequence. These mutations can transform into single nucleotide polymorphisms (SNPs) within populations and diverging proteins between species. It is well known that in many cases amino acids can be replaced by others without impeding the functioning of the protein, even if these are of quite different physico-chemical character. In some cases, however, almost any replacement would result in a functionally deficient protein. Based upon comprehensive published SNP data and applying correlation analysis we quantified the two antagonist factors controlling the process of amino acid replacement and thus protein evolution: First, the degenerate structure of the genetic code which facilitates the exchange of certain amino acids and, second, the physico-chemical forces which limit the range of possible exchanges to maintain a functional protein. We found that the observed frequencies of amino acid exchanges within species are best explained by the genetic code and that the conservation of physico-chemical properties plays a subordinate role, but has nevertheless to be considered as a key factor. Between moderately diverged species genetic code and physico-chemical properties exert comparable influence on amino acid exchanges. We furthermore studied amino acid exchanges in more detail for six species (four mammals, one bird, and one insect) and found that the profiles are highly correlated across all examined species despite their large evolutionary divergence of up to 800 million years. The species specific exchange profiles are also correlated to the exchange profile observed between different species. The currently available huge body of SNP data allows to characterize the role of two major shaping forces of protein evolution more quantitatively than before.     

Evolution of the Genetic Code by Incorporation of Amino Acids that Improved or Changed Protein Function

Fifty years have passed since the genetic code was deciphered, but how the genetic code came into being has not been satisfactorily addressed. It is now widely accepted that the earliest genetic code did not encode all 20 amino acids found in the universal genetic code as some amino acids have complex biosynthetic pathways and likely were not available from the environment. Therefore, the genetic code evolved as pathways for synthesis of new amino acids became available. One hypothesis proposes that early in the evolution of the genetic code four amino acids—valine, alanine, aspartic acid, and glycine—were coded by GNC codons (N = any base) with the remaining codons being nonsense codons. The other sixteen amino acids were subsequently added to the genetic code by changing nonsense codons into sense codons for these amino acids. Improvement in protein function is presumed to be the driving force behind the evolution of the code, but how improved function was achieved by adding amino acids has not been examined. Based on an analysis of amino acid function in proteins, an evolutionary mechanism for expansion of the genetic code is described in which individual coded amino acids were replaced by new amino acids that used nonsense codons differing by one base change from the sense codons previously used. The improved or altered protein function afforded by the changes in amino acid function provided the selective advantage underlying the expansion of the genetic code. Analysis of amino acid properties and functions explains why amino acids are found in their respective positions in the genetic code.     

On the optimality of the genetic code, with the consideration of coevolution theory by comparison of prominent cost measure matrices

Statistical and biochemical studies have revealed non-random patterns in codon assignments. The canonical genetic code is known to be highly efficient in minimizing the effects of mistranslation errors and point mutations, since it is known that when an amino acid is converted to another due to error, the biochemical properties of the resulted amino acid are usually very similar to those of the original one. In this study, using altered forms of the fitness functions used in the prior studies, we have optimized the parameters involved in the calculation of the error minimizing property of the genetic code so that the genetic code outscores the random codes as much as possible. This work also compares two prominent matrices, the Mutation Matrix and Point Accepted Mutations 74-100 (PAM(74-100)). It has been resulted that the hypothetical properties of the coevolution theory of the genetic code are already considered in PAM(74-100), giving more evidence on the existence of bias towards the genetic code in this matrix. Furthermore, our results indicate that PAM(74-100) is biased towards the single base mistranslation occurrences in second codon position as well as the frequency of amino acids. Thus PAM(74-100) is not a suitable substitution matrix for the studies conducted on the evolution of the genetic code.     

Conservation of physico-chemical amino acid properties during the evolution of proteins

Based on a similarity ring constructed from a substitution probability matrix, we have analyzed the conservation of some amino acid properties in the evolution of proteins. Refractive index and bulkiness are highly conserved, hydrophobicity and polarity are fairly well conserved while optical rotation appears to be a less relevant property. On the other hand, the analysis of the correspondence between phenotype and genotype shows that the most frequent amino acid substitutions in proteins do not always correspond to the most feasible codon changes. The apparent disagreement between amino acid substitutions in modern proteins and the primordial amino acid-codon assignment is discussed.     

The complementarity of the chemical parameters of amino acids and anticodons in the genetic code

Chemical language of the genetic code is suggested in which elementary information code units are presented by functional groups of amino acids and nucleotides. Using this language, the existence of correspondence and conformity of chemical parameters of amino acids and of central nucleotides of their anticodons was demonstrated. These findings confirm the idea that the genetic code is determined by chemical properties of amino acids and nucleotides and that this determination is the result of direct specific interactions between amino acids and nucleotide triplets at the stage of the origin of the code. The data obtained reveal primary role of anticodon triplets in the origin of the code. Key role of the central nucleotide in triplets for amino acid coding is confirmed.     

Entropy of the genetic information and evolution.

The entropy of the amino acid sequences coded by DNA is considered as a measure of diversity of variety of proteins, and is taken as a measure of evolution. The DNA or m-RNA sequence is considered as a stationary second-order Markov chain composed of four kinds of bases. Because of the biased nature of the genetic code table, increase of entropy of amino acid sequences is possible with biased nucleotide sequence. Thus the biased DNA base composition and the extreme rarity of the base doublet CpG of higher organisms are explained. It is expected that the amino acid composition was highly biased at the days of the origin of the genetic code table, and the more frequent amino acids have tended to get rarer, and the rarer ones more frequent. This tendency is observed in the evolution of hemoglobin, cytochrome C, fibrinopeptide, immunoglobulin and lysozyme, and protein as a whole.     

A hardware interpretation of the evolution of the genetic code

A quantitative rationale for the evolution of the genetic code is developed considering the principle of minimal hardware. This principle defines an optimal code as one that minimizes for a given amount of information encoded, the product of the number of physical devices used by the average complexity of each device. By identifying the number of different amino acids, number of nucleotide positions per codon and number of base types that can occupy each such position with, respectively, the amount of information, number of devices and the complexity, we show that optimal codes occur for 3, 7 and 20 amino acids with codons having a single, two and three base positions per codon, respectively. The advantage of a code of exactly 4 symbols is deduced, as well as a plausible evolutionary pathway from a code of doublets to triplets. The present day code of 20 amino acids encoded by 64 codons is shown to be the most optimal in an absolute sense. Using a tetraplet code further evolution to a code in which there would be 55 amino acids is in principle possible, but such a code would deviate slightly more than the present day code from the minimal hardware configuration. The change from a triplet code to a tetraplet code would occur at about 32 amino acids. Our conclusions are independent of, but consistent with, the observed physico-chemical properties of the amino acids and codon structures. These correlations could have evolved within the constrains imposed by the minimal hardware principle.     

Ratios of radical to conservative amino acid replacement are affected by mutational and compositional factors and may not be indicative of positive Darwinian selection

The ratio of radical to conservative amino acid replacements is frequently used to infer positive Darwinian selection. This method is based on the assumption that radical replacements are more likely than conservative replacements to improve the function of a protein. Therefore, if positive selection plays a major role in the evolution of a protein, one would expect the radical-conservative ratio to exceed the expectation under neutrality. Here, we investigate the possibility that factors unrelated to selection, i.e., transition-transversion ratio, codon usage, genetic code, and amino acid composition, influence the radical-conservative replacement ratio. All factors that have been studied were found to affect the radical-conservative replacement ratio. In particular, amino acid composition and transition-transversion ratio are shown to have the most profound effects. Because none of the studied factors had anything to do with selection (positive or otherwise) and also because all of them (singly or in combination) affected a measure that was supposed to be indicative of positive selection, we conclude that selectional inferences based on radical-conservative replacement ratios should be treated with suspicion.     

Fitness conferred by replaced amino acids declines with time

The fitness landscape of a locus, the array of fitnesses conferred by its alleles, can be affected by allele replacements at other loci, in the presence of epistatic interactions between loci. In a pair of diverging homologous proteins, the initially high probability that an amino acid replacement in one of them will make it more similar to the other declines with time, implying that the fitness landscapes of homologous sites diverge. Here, we use data on within-population non-synonymous polymorphisms and on amino acid replacements between species to study the dynamics, after an amino acid replacement, of the fitness of the ancestral amino acid, and show that selection against its restoration increases with time. This effect can be owing to increase of fitness conferred by the new amino acid occupying the site, and/or to decline of fitness conferred by the replaced amino acid. We show that the fitness conferred by the replaced amino acid rapidly declines, reaching a new lower steady-state level after approximately 20 per cent of amino acids in the protein get replaced. Therefore, amino acid replacements in evolving proteins are routinely involved in negative epistatic interactions with currently absent amino acids, and chisel off the unused parts of the fitness landscape.     

Evolution of amino acid frequencies in proteins over deep time: inferred order of introduction of amino acids into the genetic code

To understand more fully how amino acid composition of proteins has changed over the course of evolution, a method has been developed for estimating the composition of proteins in an ancestral genome. Estimates are based upon the composition of conserved residues in descendant sequences and empirical knowledge of the relative probability of conservation of various amino acids. Simulations are used to model and correct for errors in the estimates. The method was used to infer the amino acid composition of a large protein set in the Last Universal Ancestor (LUA) of all extant species. Relative to the modern protein set, LUA proteins were found to be generally richer in those amino acids that are believed to have been most abundant in the prebiotic environment and poorer in those amino acids that are believed to have been unavailable or scarce. It is proposed that the inferred amino acid composition of proteins in the LUA probably reflects historical events in the establishment of the genetic code.     

A new hypothesis on the evolution of the genetic code

Summary Starting from the assumption that specific steric and energetic interactions between amino acids and their respective anticodons could exist, the evolution of the genetic code is deduced from purely chemical and physical reasons. In this model the amino acids are intercalated between the two first anticodon bases and their carbon bound hydrogen atoms are assumed to penetrate into the electron clouds of the bases. By these means a gain in energy and a fixation of the amino acid is obtained in such a way that the anticodon nucleotides could be determinant for the nature of the amino acids.     

An alternative model of amino acid replacement

MOTIVATION: The observed correlations between pairs of homologous protein sequences are typically explained in terms of a Markovian dynamic of amino acid substitution. This model assumes that every location on the protein sequence has the same background distribution of amino acids, an assumption that is incompatible with the observed heterogeneity of protein amino acid profiles and with the success of profile multiple sequence alignment. RESULTS: We propose an alternative model of amino acid replacement during protein evolution based upon the assumption that the variation of the amino acid background distribution from one residue to the next is sufficient to explain the observed sequence correlations of homologs. The resulting dynamical model of independent replacements drawn from heterogeneous backgrounds is simple and consistent, and provides a unified homology match score for sequence-sequence, sequence-profile and profile-profile alignment.     

Genetic code preferentially conserves long-range interactions among the amino acids

The physical properties of amino acids were investigated in order to evaluate their possible relationship to the assignment of codons for amino acids in the genetic code. A comparison of the interconversion probability between amino acids and the distances between the amino acids for individual physical properties revealed a striking hierarchy among the physical properties. Surprisingly, it is the long-range/solvent interactions and not the short-range/stereochemical properties which are preferentially conserved in the genetic code.