A multivariate study of the relationship between the genetic code and the physical-chemical properties of amino acids

Optimizing Nucleotide Mixtures to Encode Specific Subsets of Amino Acids for Semi-Random Mutagenesis

In random mutagenesis, synthesis of an NNN triplet (i.e. equiprobable A, C, G, and T at each of the three positions in the codon) could be considered an optimal nucleotide mixture because all 20 amino acids are encoded. NN(G,C) might be considered a slightly more intelligent "dope" because the entire set of amino acids is still encoded using only half as many codons. Using a general algorithm described herein, it is possible to formulate more complex doping schemes which encode specific subsets of the twenty amino acids, excluding others from the mix. Maximizing the equiprobability of amino acid residues contributing to such a subset is suggested as an optimal basis for performing semi-random mutagenesis. This is important for reducing the nucleotide complexity of combinatorial cassettes so that "sequence space" can be searched more efficiently. Computer programs have been developed to provide tables of optimized dopes compatible with automated DNA synthesizers.

Periodical changes of amino acid reactivity within the genetic code

Enthalpies (delta H++) and entropies (delta S++) of activation for the reaction of 18 N'-hydroxysuccinimide esters of N-protected proteinaceous amino acids with p-anisidine were measured and free enthalpies of activation (delta G++) at 25 degrees C were calculated on this basis. A regular correlation between delta G++s and the corresponding amino acid codons was found. To obtain this correlation all the codons had to be arranged in a closed ring in which the consecutive codons were connected by one-step mutational changes. One-step mutations appeared as a regular series: 2,3,3,3,1,3,3,3,1,3,3,3,1,3,3,3,2,3,3,3. (the numbers denote a codon position in which a change took place). There were three such 'one-step mutation periods' in the ring, each containing 20 codons (in each block of 16 codons with A, U and C, in the central position and 4 codons containing G in the central position). The end of the third period (UG) and the beginning of the first period were bridged by the four codons of glycine with G in the second position. The values of delta G++ change similarly in each period, increasing upon approaching Lys, Pro, and Ile. The periodical relation between the chemical reactivities of the coded amino acids (reflected by delta G++s) and the structure of their codons could be of importance for the origin of the genetic code i.e. for selection of proper codons for the definite amino acids.

A quantitative measure of error minimization in the genetic code

We have calculated the average effect of changing a codon by a single base for all possible single-base changes in the genetic code and for changes in the first, second, and third codon positions separately. Such values were calculated for an amino acid's polar requirement, hydropathy, molecular volume, and isoelectric point. For each attribute the average effect of single-base changes was also calculated for a large number of randomly generated codes that retained the same level of redundancy as the natural code. Amino acids whose codons differed by a single base in the first and third codon positions were very similar with respect to polar requirement and hydropathy. The major differences between amino acids were specified by the second codon position. Codons with U in the second position are hydrophobic, whereas most codons with A in the second position are hydrophilic. This accounts for the observation of complementary hydropathy. Single-base changes in the natural code had a smaller average effect on polar requirement than all but 0.02% of random codes. This result is most easily explained by selection to minimize deleterious effects of translation errors during the early evolution of the code.

The extension reached by the minimization of the polarity distances during the evolution of the genetic code

The level reached by the optimization of the polarity distances during the evolution of the genetic code was investigated. The results, although not conclusive, indicate that this optimization level is higher than the data reported in the literature. The results seem compatible with the reaching of an evolutionary minimum, with respect to the optimization of the polarity distances, by the genetic code during its formation.

Some aspects of the organization and evolution of the genetic code

In this paper, I define a measure of the relative position of each amino acid in the genetic code by means of a 21-dimensional vector describing its potential for mutation, in a single step, to each of the other amino acids, or to a chain termination codon. This measure allows us to make a systematic investigation of the type and number of the physicochemical properties of the amino acids that were involved in evolution. The polar character and size of amino acids are identified in this analysis as properties that played a leading role in the evolutionary history of the genetic code. The application of cluster analysis and discriminant analysis reveals the characteristics of the structural organization of the genetic code. Finally, I suggest the existence of a relationship between the molecular weight of the amino acids and the number of synonymous codons.

Interrelatedness of 5S RNA sequences investigated by correspondence analysis

Correspondence analysis (a form of multivariate statistics) applied to 74 5S ribosomal RNA sequences indicates that the sequences are interrelated in a systematic, nonrandom fashion. Aligned sequences are represented as vectors in a 5N-dimensional space, where N is the number of base positions in the 5S RNA molecule. Mutually orthogonal directions (called factor axes) along which intersequence variance is greatest are defined in this hyperspace. Projection of the sequences onto planes defined by these factorial directions reveals clustering of species that is suggestive of phylogenetic relationships. For each factorial direction, correspondence analysis points to regions of "importance," i.e., those base positions at which the systematic changes occur that define that particular direction. In effect, the technique provides a rapid determination of group-specific signatures. In several instances, similarities between sequences are indicated that have only recently been inferred from visual base-to-base comparisons. These results suggest that correspondence analysis may provide a valuable starting point from which to uncover the patterns of change underlying the evolution of a macromolecule, such as 5S RNA.