Differential effect of amino acid residues on the stability of double helices formed from polyribonucleotides and its possible relation to the evolution of the genetic code

The interaction of amino acid residues with polyribonucleotides was characterized by measurements of melting temperatures (tm) for poly(A).poly(U) and poly(I).poly(C) as functions of the concentrations of various amino acid amides. The amides of hydrophilic amino acids lead to a continuous increase of tm with increasing concentration, whereas amides of hydrophobic amino acids induce a decrease of tm at low concentrations (approximately 1 mM) followed by an increase at higher concentrations. Analysis of the data by a simple site model provides the affinity of each ligand for the double helix relative to that for the single strands. This parameter decreases in the order Ala greater than Gly greater than Ser greater than Asn greater than Pro greater than Met, Val greater than Ile, Leu for poly(A).poly(U) and Ala, Gly, Ser greater than Asn greater than Pro greater than Val greater than Ile, Met, Leu for poly(I).poly(C). The special effects of hydrophobic amino acids may be related to the similarity of the codons for these amino acids. A simple model for assignment of codons to amino acids is proposed.     

The regularities of the changes of amino acid physico-chemical properties within the genetic code

Summary All the codons of the genetic code can be arranged into the closed one-step mutation ring, containing three periods of the same sequence of mutations (2,3,3,3,1,3,3,3,1,3,3,3,1,3,3,3,2,3,3,3). The codons of Gly play a role of the connecting element between the end of the third, and the beginning of the first period of the genetic code. The reactivity of amino acids, expressed by the reaction rates of aminolysis reaction of N-hydroxysuccinimide esters of protected amino acids with p-anisidine, changes periodically with the respect to the mutation periods of the genetic code. Chou-Fasman P as well as P conformational parameters of amino acids, and also the compositional frequencies of amino acids in proteins, demonstrate the pseudosymmetry pattern with respect to the center of one-step mutation ring, which is situated between Thr ACY and ACR codons.     

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.     

Local stability and evolution of the genetic code

The standard genetic code is known to be robust to translation errors and point mutations. We studied how small modifications of the standard code affect its robustness. The robustness was assessed in terms of a proper stability function, the negative variations of which correspond to a more robust code. The fraction of more robust codes obtained under small modifications appeared to be unexpectedly high, about 0.1-0.4 depending on the choice of stability function and code modifications, yet significantly lower than the corresponding fraction in the random codes (about a half). In this sense the standard code ought to be considered distinctly non-random in accordance with previous observations. The distribution of the negative variations of stability function revealed very abrupt drop beyond one standard deviation, much sharper than for Gaussian distribution or for the random codes with the same number of codons in the sets coding for amino acids or stop-codons. This behavior holds for both the standard code as a whole and its binary NRN-NYN, NWN-NSN, and NMN-NKN blocks. Previously, it has been proved that such binary block structure is necessary for the robustness of a code and is inherent to the standard genetic code. The modifications of the standard code corresponding to more robust coding may be related to the different variants of the code. These effects may also contribute to the rates of replacements of amino acids. The observed features demonstrate the joint impact of random factors and natural selection during evolution of the genetic code.     

Structure of the genetic code suggested by the hydropathy correlation between anticodons and amino acid residues

The correlation between hydropathies of anticodons and amino acids, detected by other authors utilizing scales of amino acid molecules in solution, was improved with the utilization of scales of amino acid residues in proteins. Three partitions were discerned in the correlation plot with the principal dinucleotides of anticodons (pDiN, excluding the wobble position). (a) The set of outliers of the correlation: Gly-CC, Pro-GG, Ser-GA and Ser-CU. The amino acids are consistently small, hydro-apathetic, stabilizers of protein N-ends, preferred in aperiodic protein conformations and belong to synthetases class II. The pDiN sequences are representative of the homogeneous sector (triplets NRR and NYY), distinguished from the mixed sector (triplets NRY and NYR), that depict a 70% correspondence to the synthetases class II and I, respectively. The triplet pairs proposed to be responsible for the coherence in the set of outliers are of the palindromic kind, where the lateral bases are the same, CCC: GGG and AGA: UCU. This suggests that UCU previously belonged to Ser, adding to other indications that the attribution of Arg to YCU was due to an expansion of the Arg-tRNA synthetase specificity. The other attributions produced two correlation sets. (b) One corresponds to the remaining pDiN of the homogeneous sector, containing both synthetase classes; its regression line overlapped the one formed by the remaining attributions to class II. (c) The other contains the pDiN of the mixed sector and produced steeper slopes, especially with the class I attributions. It is suggested that the correlation was established when the amino acid composition of the protein synthetases became progressively enriched and that the set of outliers were the earliest to have been fixed.     

Natural history and experimental evolution of the genetic code

The standard genetic code is a set of rules that relates the 20 canonical amino acids in proteins to groups of three bases in the mRNA. It evolved from a more primitive form and the attempts to reconstruct its natural history are based on its present-day features. Genetic code engineering as a new research field was developed independently in a few laboratories during the last 15 years. The main intention is to re-program protein synthesis by expanding the coding capacities of the genetic code via re-assignment of specific codons to un-natural amino acids. This article focuses on the question as to which extent hypothetical scenarios that led to codon re-assignments during the evolution of the genetic code are relevant for its further evolution in the laboratory. Current attempts to engineer the genetic code are reviewed with reference to theoretical works on its natural history. Integration of the theoretical considerations into experimental concepts will bring us closer to designer cells with target-engineered genetic codes that should open not only tremendous possibilities for the biotechnology of the twenty-first century but will also provide a basis for the design of novel life forms.     

On the origin and evolution of the genetic code

Two ideas have essentially been used to explain the origin of the genetic code: Crick's frozen accident and Woese's amino acid-codon specific chemical interaction. Whatever the origin and codon-amino acid correlation, it is difficult to imagine the sudden appearance of the genetic code in its present form of 64 codons coding for 20 amino acids without appealing to some evolutionary process. On the contrary, it is more reasonable to assume that it evolved from a much simpler initial state in which a few triplets were coding for each of a small number of amino acids. Analysis of genetic code through information theory and the metabolism of pyrimidine biosynthesis provide evidence that suggests that the genetic code could have begun in an RNA world with the two letters A and U grouped in eight triplets coding for seven amino acids and one stop signal. This code could have progressively evolved by making gradual use of letters G and C to end with 64 triplets coding for 20 amino acids and three stop signals. According to proposed evidence, DNA could have appeared after the four-letter structure was already achieved. In the newborn DNA world, T substituted U to get higher physicochemical and genetic stability.     

Use of Chou-Fasman amino acid conformational parameters to analyze the organization of the genetic code and to construct protein genealogies

Summary Chou-Fasman parameters, measuring preferences of each amino acid for different conformational regions in proteins, were used to obtain an amino acid difference index of conformational parameter distance (CPD) values. CPD values were found to be significantly lower for amino acid exchanges representing in the genetic code transitions of purines, GA than for exchanges representing either transitions of pyrimidines, CU, or transversions of purines and pyrimidines. Inasmuch as the distribution of CPD values in these non GA exchanges resembles that obtained for amino acid pairs with double or triple base differences in their underlying codons, we conclude that the genetic code was not particularly designed to minimize effects of mutation on protein conformation. That natural selection minimizes these changes, however, was shown by tabulating results obtained by the maximum parsimony method for eight protein genealogies with a total occurrence of 4574 base substitutions. At the beginning position of the codons GA transitions were in very great excess over other base substitutions, and, conversely, CU transitions were deficient. At the middle position of the codons only fast evolving proteins showed an excess of GA transitions, as though selection mainly preserved conformation in these proteins while weeding out mutations affecting chemical properties of functional sites in slow evolving proteins. In both fast and slow evolving proteins the net direction of transitions and transversions was found to be from G beginning codons to non-G beginning codons resulting in more commonly occurring amino acids, especially alanine with its generalized conformational properties, being replaced at suitable sites by amino acids with more specialized conformational and chemical properties. Historical circumstances pertaining to the origin of the genetic code and the nature of primordial proteins could account for such directional changes leading to increases in the functional density of proteins.In order to further explore the course of protein evolution, a modified parsimony algorithm was developed for constructing protein genealogies on the basis of minimum CPD length. The algorithm's ability to judge with finer discrimination that in protein evolution certain pathways of amino acid substitution should occur more readily than others was considered a potential advantage over strict maximum parsimony. In developing this CPD algorithm, the path of minimum CPD length through intermediate amino acids allowed by the genetic code for each pair of amino acids was determined. It was found that amino acid exchanges representing two base changes have a considerably lower average CPD value per base substitution than the amino acid exchanges representing single base changes. Amino acid exchanges representing three base changes have yet a further marked reduction in CPD per base change. This shows how extreme constraining effects of stabilizing selection can be circumvented, for by way of intermediate amino acids almost any amino acid can ultimately be substituted for another without damage to an evolving protein's conformation during the process.     

Stability of the genetic code and optimal parameters of amino acids

The standard genetic code is known to be much more efficient in minimizing adverse effects of misreading errors and one-point mutations in comparison with a random code having the same structure, i.e. the same number of codons coding for each particular amino acid. We study the inverse problem, how the code structure affects the optimal physico-chemical parameters of amino acids ensuring the highest stability of the genetic code. It is shown that the choice of two or more amino acids with given properties determines unambiguously all the others. In this sense the code structure determines strictly the optimal parameters of amino acids or the corresponding scales may be derived directly from the genetic code. In the code with the structure of the standard genetic code the resulting values for hydrophobicity obtained in the scheme “leave one out” and in the scheme with fixed maximum and minimum parameters correlate significantly with the natural scale. The comparison of the optimal and natural parameters allows assessing relative impact of physico-chemical and error-minimization factors during evolution of the genetic code. As the resulting optimal scale depends on the choice of amino acids with given parameters, the technique can also be applied to testing various scenarios of the code evolution with increasing number of codified amino acids. Our results indicate the co-evolution of the genetic code and physico-chemical properties of recruited amino acids.     

Challenges of the genetic code for exploring sequence space in directed protein evolution

Directed protein evolution is the most versatile method for studying protein structure–function relationships, and for tailoring a protein's properties to the needs of industrial applications. In this review, we performed a statistical analysis on the genetic code to study the extent and consequence of the organization of the genetic code on amino acid substitution patterns generated in directed evolution experiments. In detail, we analyzed amino acid substitution patterns caused by (a) a single nucleotide (nt) exchange at each position of all 64 codons, and (b) two subsequent nt exchanges (first and second nt, first and third nt, second and third nt). Additionally, transitions and transversions mutations were compared at the level of amino acid substitution patterns. The latter analysis showed that single nucleotide substitution in a codon generates only 39.5% of the natural diversity on the protein level with 5.2–7 amino acid substitutions per codon. Transversions generate more complex amino acid substitution patterns (increased number and chemically more diverse amino acid substitutions) than transitions. Simultaneous nt exchanges at both first and second nt of a codon generates very diverse amino acid substitution patterns, achieving 83.2% of the natural diversity. The statistical analysis described in this review sets the objectives for novel random mutagenesis methods that address the consequences of the organization of the genetic code. Random mutagenesis methods that favor transversions or introduce consecutive nt exchanges can contribute in this regard.     

Challenges of the genetic code for exploring sequence space in directed protein evolution

Directed protein evolution is the most versatile method for studying protein structure-function relationships, and for tailoring a protein's properties to the needs of industrial applications. In this review, we performed a statistical analysis on the genetic code to study the extent and consequence of the organization of the genetic code on amino acid substitution patterns generated in directed evolution experiments. In detail, we analyzed amino acid substitution patterns caused by (a) a single nucleotide (nt) exchange at each position of all 64 codons, and (b) two subsequent nt exchanges (first and second nt, first and third nt, second and third nt). Additionally, transitions and transversions mutations were compared at the level of amino acid substitution patterns. The latter analysis showed that single nucleotide substitution in a codon generates only 39.5% of the natural diversity on the protein level with 5.2-7 amino acid substitutions per codon. Transversions generate more complex amino acid substitution patterns (increased number and chemically more diverse amino acid substitutions) than transitions. Simultaneous nt exchanges at both first and second nt of a codon generates very diverse amino acid substitution patterns, achieving 83.2% of the natural diversity. The statistical analysis described in this review sets the objectives for novel random mutagenesis methods that address the consequences of the organization of the genetic code. Random mutagenesis methods that favor transversions or introduce consecutive nt exchanges can contribute in this regard.     

Block structure and stability of the genetic code

It is known that different codons may be unified into larger groups related to the hierarchical structure, approximate hidden symmetries, and evolutionary origin of the universal genetic code. Using a simplified evolutionary motivated two-letter version of genetic code, the general principles of the most stable coding are discussed. By the complete enumeration in such a reduced code it is strictly proved that the maximum stability with respect to point mutations and shifts in the reading frame needs the fixation of the middle letters within codons in groups with different physico-chemical properties, thus, explaining a key feature of the universal genetic code. The translational stability of the genetic code is studied by the mapping of code onto de Bruijn graph providing both the compact visual representation of mutual relationships between different codons as well as between codons and protein coding DNA sequence and a powerful tool for the investigation of stability of protein coding. Then, the results are extended to four-letter codes. As is shown, the universal genetic code obeys mainly the principles of optimal coding. These results demonstrate the hierarchical character of optimization of universal genetic code with strictly optimal coding being evolved at the earliest stages of molecular evolution. Finally, the universal genetic code is compared with the other natural variants of genetic codes.     

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

Summary The 20 naturally occurring amino acids are characterized by 20 variables: pK_{NH 2}, pK_COOH, pI, molecular weight, substituent van der Waals volume, seven¹H and¹³C nuclear magnetic resonance shift variables, and eight hydrophobicity-hydrophilicity scales. The 20-dimensional data set is reduced to a few new dimensions by principal components analysis. The three first principal components reveal relationships between the properties of the amino acids and the genetic code. Thus the amino acids coded for by adenosine (A), uracil (U), or cytosine (C) in their second codon position (corresponding to U, A, or G in the second anticodon position) are grouped in these components. No grouping was detected for the amino acids coded for by guanine (G) in the second codon position (corresponding to C in the second anticodon position). The results show that a relationship exists between the physical-chemical properties of the amino acids and which of the A (U), U (A), or C (G) nucleotide is used in the second codon (anticodon) position. The amino acids coded for by G (C) in the second codon (anticodon) position do not participate in this relationship.     

Globular proteins,GU wobbling,and the evolution of the genetic code

Summary It has previously been shown that the formation of GU base pairs in RNA copying processes leads to an accumulation of G and U in both strands of the replicating RNA, which results in a non-random distribution of base triplets. In the present paper, this distribution is calculated, and, using the ²-test, a correlation between the distribution of triplets and the amino acid composition of the evolutionarily conservative interior regions of selected globular proteins is established.It is suggested that GU wobbling in early replication of RNA could have led to the observed amino acid composition of present-day protein interiors. If this hypothesis is correct, the GU wobbling must have been very extensive in the imprecisely replicating RNA, even reaching values close to the critical for stability of its double-helical structure. Implications of the hypothesis both for the evolution of the genetic code and of proteins are discussed.     

Use of fuzzy clustering technique and matrices to classify amino acids and its impact to Chou's pseudo amino acid composition

In this paper we present a study of classification of the 20 amino acids via a fuzzy clustering technique. In order to calculate distances among the various elements we employ two different distance functions: the Minkowski distance function and the NTV metric. In the clustering procedure we take into account several physical properties of the amino acids. We examine the effect of the number and nature of properties taken into account to the clustering procedure as a function of the degree of similarity and the distance function used. It turns out that one should use the properties that determine in the more important way the behavior of the amino acids and that the use of the appropriate metric can help in defining the separation into groups.     

Structure,function and evolution of seryl-tRNA synthetases: Implications for the evolution of aminoacyl-tRNA synthetases and the genetic code

Two aspects of the evolution of aminoacyl-tRNA synthetases are discussed. Firstly, using recent crystal structure information on seryl-tRNA synthetase and its substrate complexes, the coevolution of the mode of recognition between seryl-tRNA synthetase and tRNA^ser in different organisms is reviewed. Secondly, using sequence alignments and phylogenetic trees, the early evolution of class 2 Amnoacyl-tRNA synthetases is traced. Arguments are presented to suggest that synthetases are not the oldest of protein enzymes, but survived as RNA enzymes during the early period of the evolution of protein catalysts. In this view, the relatedness of the current synthetases, as evidenced by the division into two classes with their associated subclasses, reflects the replacement of RNA synthetases by protein synthetases. This process would have been triggered by the acquisition of tRNA 3 end charging activity by early proteins capable of activating small molecules (e.g., amino acids) with ATP. If these arguments are correct, the genetic code was essentially frozen before the protein synthetases that we know today came into existence. Correspondence to: S. CusackBased on a presentation made at a workshop-Aminoacyl-tRNA Synthetases and the Evolution of the Genetic Code-held at Berkeley, CA, July 17–20, 1994     

Aldocyanoin microspheres: Partial amino acid analysis of the microparticulates formed from simple reactants under various conditions

Summary The work of Kenyon and Nissenbaum on aldocyanoin microspheres was repeated and extended. It was determined that the microspheres contained amino acids and that specific amino acids could be incorporated into the microspheres by adding the requisite aldehyde or ketone precursor to the model mixture. Microsphere formation was found to be dependent on the availability of oxygen. Under anaerobic conditions of synthesis, no micro-spheres formed in the time allotted and the amino acid composition of the macromolecular material was simple. Microparticulate material synthesized by C. Folsome using a quenched spark technique was analyzed and found to contain amino acids that had a qualitative composition similar to both a Miller-Urey discharge and the Kenyon-Nissenbaum microspheres.     

Quantifying the effect of burial of amino acid residues on protein stability

The average contribution of individual residue to folding stability and its dependence on buried accessible surface area (ASA) are obtained by two different approaches. One is based on experimental mutation data, and the other uses a new knowledge-based atom-atom potential of mean force. We show that the contribution of a residue has a significant correlation with buried ASA and the regression slopes of 20 amino acid residues (called the buriability) are all positive (pro-burial). The buriability parameter provides a quantitative measure of the driving force for the burial of a residue. The large buriability gap observed between hydrophobic and hydrophilic residues is responsible for the burial of hydrophobic residues in soluble proteins. Possible factors that contribute to the buriability gap are discussed.     

Conformational fluctuations versus constraints in amino acid side chains: the evolution of information content from free amino acids to proteins

Like all other complex biological systems, proteins exhibit properties not found in free amino acids (i.e., emergent properties). Here, we explore top-down constraints experienced by the residue side chains in proteins compared to amino acids in increasingly complex molecular environments: free amino acids, end-capped amino acids, and the central residue in an alpha-helical nonapeptide. The crystalline structure of the contractile protein profilin Ib and the enzyme trypsin were chosen as objects of study, and submitted to 10 ns molecular dynamics (MD) simulations. The results revealed increased conformational constraints on the side chains when going from the simpler to the more complex compounds. A Shannon entropy (SE) analysis of the conformational behavior of the side chains showed in most cases a progressive and marked decrease in the SE of the chi1 and chi2 dihedral angles. This is equivalent to stating that conformational constraints on the side chain of residues increase their information content and, hence, recognition specificity compared to free amino acids. In other words, the vastly increased information content of a protein relative to its free monomers is embedded not only in the tertiary structure of the backbone, but also in the conformational behavior of the side chains. The postulated implication is that both backbone and side chains, by virtue of being conformationally constrained, contribute to the protein's recognition specificity toward other macromolecules and ligands.     

Determination of amino acid tissue concentrations by microdialysis: method evaluation and relation to plasma values

Summary. Microdialysis is an in vivo technique to monitor tissue concentrations of low molecular weight substances by means of a continuously perfused artificial capillary with a semipermeable membrane placed into the region of interest. The suitability of microdialysis to determine tissue concentrations of amino acids was evaluated in vitro by placing the catheter into Ringer buffer or into a plasma protein (50 g/l) solution containing 32 different amino acids (150 μmol/l each). All amino acids tested crossed freely the microdialysis membrane with recoveries close to 100%. Microdialysis fluid was sampled from subcutaneous tissue of five newborns and amino acid content analysed. Total and non protein bound amino acids were determined in the patients plasma by acid precipitation or ultrafiltration, respectively. Mean subcutaneous tissue concentrations were lower as compared to plasma for taurine, serine, alanine, aspartate, glutamate and ornithine and higher for valine, isoleucine, leucine, methionine, phenylalanine, tyrosine and arginine, indicating net uptake or release of amino acids from subcutaneous tissue. Thus, microdialysis offers a convenient and minimal invasive way to study tissue amino acid composition and appears to be a promising analytical tool for the study of amino acid metabolism in vivo. Received August 7, 2000 Accepted January 7, 2001