Optimality of the genetic code with respect to protein stability and amino-acid frequencies

Background

The genetic code is known to be efficient in limiting the effect of mistranslation errors. A misread codon often codes for the same amino acid or one with similar biochemical properties, so the structure and function of the coded protein remain relatively unaltered. Previous studies have attempted to address this question quantitatively, by estimating the fraction of randomly generated codes that do better than the genetic code in respect of overall robustness. We extended these results by investigating the role of amino-acid frequencies in the optimality of the genetic code.

Results

We found that taking the amino-acid frequency into account decreases the fraction of random codes that beat the natural code. This effect is particularly pronounced when more refined measures of the amino-acid substitution cost are used than hydrophobicity. To show this, we devised a new cost function by evaluating in silico the change in folding free energy caused by all possible point mutations in a set of protein structures. With this function, which measures protein stability while being unrelated to the code's structure, we estimated that around two random codes in a billion (10⁹) are fitter than the natural code. When alternative codes are restricted to those that interchange biosynthetically related amino acids, the genetic code appears even more optimal.

Conclusions

These results lead us to discuss the role of amino-acid frequencies and other parameters in the genetic code's evolution, in an attempt to propose a tentative picture of primitive life.     

An Analysis of the Metabolic Theory of the Origin of the Genetic Code

A computer program was used to test Wong's coevolution theory of the genetic code. The codon correlations between the codons of biosynthetically related amino acids in the universal genetic code and in randomly generated genetic codes were compared. It was determined that many codon correlations are also present within random genetic codes and that among the random codes there are always several which have many more correlations than that found in the universal code. Although the number of correlations depends on the choice of biosynthetically related amino acids, the probability of choosing a random genetic code with the same or greater number of codon correlations as the universal genetic code was found to vary from 0.1% to 34% (with respect to a fairly complete listing of related amino acids). Thus, Wong's theory that the genetic code arose by coevolution with the biosynthetic pathways of amino acids, based on codon correlations between biosynthetically related amino acids, is statistical in nature. Received: 8 August 1996 / Accepted: 26 December 1996     

The Genetic Code: What Is It Good For? An Analysis of the Effects of Selection Pressures on Genetic Codes

How did the ``universal' genetic code arise? Several hypotheses have been put forward, and the code has been analyzed extensively by authors looking for clues to selection pressures that might have acted during its evolution. But this approach has been ineffective. Although an impressive number of properties has been attributed to the universal code, it has been impossible to determine whether selection on any of these properties was important in the code's evolution or whether the observed properties arose as a consequence of selection on some other characteristic. Therefore we turned the question around and asked, what would a genetic code look like if it had evolved in response to various different selection pressures? To address this question, we constructed a genetic algorithm. We found first that selecting on a particular measure yields codes that are similar to each other. Second, we found that the universal code is far from minimized with respect to the effects of mutations (or translation errors) on the amino acid compositions of proteins. Finally, we found that the codes that most closely resembled real codes were those generated by selecting on aspects of the code's structure, not those generated by selecting to minimize the effects of amino acid substitutions on proteins. This suggests that the universal genetic code has been selected for a particular structure—a structure that confers an important flexibility on the evolution of genes and proteins—and that the particular assignments of amino acids to codons are secondary. Received: 29 December 1998 / Accepted: 8 July 1999     

On the Evolution of Redundancy in Genetic Codes

We simulate a deterministic population genetic model for the coevolution of genetic codes and protein-coding genes. We use very simple assumptions about translation, mutation, and protein fitness to calculate mutation-selection equilibria of codon frequencies and fitness in a large asexual population with a given genetic code. We then compute the fitnesses of altered genetic codes that compete to invade the population by translating its genes with higher fitness. Codes and genes coevolve in a succession of stages, alternating between genetic equilibration and code invasion, from an initial wholly ambiguous coding state to a diversified frozen coding state. Our simulations almost always resulted in partially redundant frozen genetic codes. Also, the range of simulated physicochemical properties among encoded amino acids in frozen codes was always less than maximal. These results did not require the assumption of historical constraints on the number and type of amino acids available to codes nor on the complexity of proteins, stereochemical constraints on the translational apparatus, nor mechanistic constraints on genetic code change. Both the extent and timing of amino-acid diversification in genetic codes were strongly affected by the message mutation rate and strength of missense selection. Our results suggest that various omnipresent phenomena that distribute codons over sites with different selective requirements—such as the persistence of nonsynonymous mutations at equilibrium, the positive selection of the same codon in different types of sites, and translational ambiguity—predispose the evolution of redundancy and of reduced amino acid diversity in genetic codes. Received: 21 December 2000 / Accepted: 12 March 2001     

A Novel Theory on the Origin of the Genetic Code: A GNC-SNS Hypothesis

We have previously proposed an SNS hypothesis on the origin of the genetic code (Ikehara and Yoshida 1998). The hypothesis predicts that the universal genetic code originated from the SNS code composed of 16 codons and 10 amino acids (S and N mean G or C and either of four bases, respectively). But, it must have been very difficult to create the SNS code at one stroke in the beginning. Therefore, we searched for a simpler code than the SNS code, which could still encode water-soluble globular proteins with appropriate three-dimensional structures at a high probability using four conditions for globular protein formation (hydropathy, α-helix, β-sheet, and β-turn formations). Four amino acids (Gly [G], Ala [A], Asp [D], and Val [V]) encoded by the GNC code satisfied the four structural conditions well, but other codes in rows and columns in the universal genetic code table do not, except for the GNG code, a slightly modified form of the GNC code. Three three-amino acid systems ([D], Leu and Tyr; [D], Tyr and Met; Glu, Pro and Ile) also satisfied the above four conditions. But, some amino acids in the three systems are far more complex than those encoded by the GNC code. In addition, the amino acids in the three-amino acid systems are scattered in the universal genetic code table. Thus, we concluded that the universal genetic code originated not from a three-amino acid system but from a four-amino acid system, the GNC code encoding [GADV]-proteins, as the most primitive genetic code. Received: 11 June 2001 / Accepted: 11 October 2001     

Revisiting the Physico-Chemical Hypothesis of Code Origin: An Analysis Based on Code-Sequence Coevolution in a Finite Population

The origin of the genetic code marked a major transition from a plausible RNA world to the world of DNA and proteins and is an important milestone in our understanding of the origin of life. We examine the efficacy of the physico-chemical hypothesis of code origin by carrying out simulations of code-sequence coevolution in finite populations in stages, leading first to the emergence of ten amino acid code(s) and subsequently to 14 amino acid code(s). We explore two different scenarios of primordial code evolution. In one scenario, competition occurs between populations of equilibrated code-sequence sets while in another scenario; new codes compete with existing codes as they are gradually introduced into the population with a finite probability. In either case, we find that natural selection between competing codes distinguished by differences in the degree of physico-chemical optimization is unable to explain the structure of the standard genetic code. The code whose structure is most consistent with the standard genetic code is often not among the codes that have a high fixation probability. However, we find that the composition of the code population affects the code fixation probability. A physico-chemically optimized code gets fixed with a significantly higher probability if it competes against a set of randomly generated codes. Our results suggest that physico-chemical optimization may not be the sole driving force in ensuring the emergence of the standard genetic code.     

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.     

Selenocysteine biosynthesis and mechanism of incorporation into growing proteins

The universal genetic code codes for the 20 canonical amino acids, while selenocysteine (Sec) is encoded by UGA, one of the three well-known stop codons. Selenocysteine is of particular interest of molecular biology, principally differing in the mechanism of incorporation into growing polypeptide chains from the other 20 amino acids. The process involves certain cis- and trans-active factors, such as the Sec insertion sequence (SECIS). The SECIS is in the 3′-untranslated mRNA region in eukaryotes and within the open reading frame located immediately downstream of the Sec UGA codon in bacteria, the difference leading to differences in the mechanism of Sec incorporation between the two domains of life. The trans-active factors include Sec-tRNA^[Ser]Sec, which is synthesized by a unique system; the Sec-specific elongation factor EFsec; and a SECIS-binding protein (SBP2). Thus, many additional molecules are to be synthesized in the cell to allow Sec incorporation during translation. The fact makes Sec-containing proteins rather “expensive” and emphasizes their crucial role in metabolism.     

Physicochemical Optimization in the Genetic Code Origin as the Number of Codified Amino Acids Increases

We have assumed that the coevolution theory of genetic code origin (Wong JT, Proc Natl Acad Sci USA 72:1909–1912, 1975) is essentially correct. This theory makes it possible to identify at least 10 evolutionary stages through which genetic code organization might have passed prior to reaching its current form. The calculation of the minimization level of all these evolutionary stages leads to the following conclusions. (1) The minimization percentages increased linearly with the number of amino acids codified in the codes of the various evolutionary stages when only the sense changes are considered in the analysis. This seems to favor the physicochemical theory of genetic code origin even if, as discussed in the paper, this observation is also compatible with the coevolution theory. (2) For the first seven evolutionary stages of the genetic code, this trend is less clear and indeed is inverted when we consider the global optimisation of the codes due to both sense changes and synonymous changes. This inverse correlation between minimization percentages and the number of amino acids codified in the codes of the intermediate stages seems to favor neither the physicochemical nor the stereochemical theories of genetic code origin, as it is in the early and intermediate stages of code development that these theories would expect minimization to have played a crucial role, and this does not seem to be the case. However, these results are in agreement with the coevolution theory, which attributes a role to the physicochemical properties of amino acids that, while important, is nevertheless subordinate to the mechanism which concedes codons from the precursor amino acids to the product amino acids as the primary factor determining the evolutionary structuring of the genetic code. The results are therefore discussed in the context of the various theories proposed to explain genetic code origin. Received: 25 October 1998 / Accepted: 19 February 1999     

Natural reassignment of CUU and CUA sense codons to alanine in Ashbya mitochondria

The discovery of diverse codon reassignment events has demonstrated that the canonical genetic code is not universal. Studying coding reassignment at the molecular level is critical for understanding genetic code evolution, and provides clues to genetic code manipulation in synthetic biology. Here we report a novel reassignment event in the mitochondria of Ashbya (Eremothecium) gossypii, a filamentous-growing plant pathogen related to yeast (Saccharomycetaceae). Bioinformatics studies of conserved positions in mitochondrial DNA-encoded proteins suggest that CUU and CUA codons correspond to alanine in A. gossypii, instead of leucine in the standard code or threonine in yeast mitochondria. Reassignment of CUA to Ala was confirmed at the protein level by mass spectrometry. We further demonstrate that a predicted is transcribed and accurately processed in vivo, and is responsible for Ala reassignment. Enzymatic studies reveal that is efficiently recognized by A. gossypii mitochondrial alanyl-tRNA synthetase (AgAlaRS). AlaRS typically recognizes the G3:U70 base pair of tRNA^Ala; a G3A change in Ashbya abolishes its recognition by AgAlaRS. Conversely, an A3G mutation in Saccharomyces cerevisiae confers tRNA recognition by AgAlaRS. Our work highlights the dynamic feature of natural genetic codes in mitochondria, and the relative simplicity by which tRNA identity may be switched.     

Possible evolutionary steps in the genetic code

In mammalian mitochondrial codes, fourfold codons wobble-pair with UNN anticodons so that U wobbles with U, C, A and G. Twofold pyrimidine-terminated codons pair with GNN and twofold purine-terminated codons pair with UNN. These properties enable a prediction to be made for evolution of the universal genetic code. It was postulated (1) that an archetypal code of 16 quartets coded for 15 amino acids. If this code used UNN anticodons, then duplication of tRNA genes, followed by mutations in the anticodons and aminoacylation sites, would give rise to the present universal code.     

Evolutionary processes possibly limiting the kinds of amino acids in protein to twenty: a review.

Only 20 of more than 250 biosynthetic amino acids are common (coded) constituents of contemporary protein. In this paper, several stages of evolution, both prebiotic and biotic, are examined for means by which other (non-proteinous) amino acids may have been selected against. Simulated prebiotic experiments indicate that some non-proteinous amino acids were present prebiotically, that they could be incorporated during the formation of prebiotic protein, and that they would function in such protein. Biotic selection is thus indicated.Non-proteinous amino acids currently are available via biosynthetic pathways for potential incorporation into bioprotein. Codon-anticodon interaction, peptidyl transferases, and elongation and termination factors of protein synthesis do not show the specificity needed to preclude non-proteinous amino acids. Highly specific recognition among amino acids, tRNAs, and activating enzymes is concluded to be why the kinds of amino acids in contemporary protein are limited to twenty.Some of several theories concerning the origin, nature and evolution of the genetic code can readily accommodate non-proteinous amino acids. Some evidence suggests that such amino acids were eventually eliminated from protein because they were less suitable than related proteinous amino acids. However, deterministic or “direct interaction” theories currently lack sufficient experimental support to answer how non-proteinous amino acids were precluded; such theories, being testable, probably have the most potential for providing an answer.     

Sequence signatures of direct complementarity between mRNAs and cognate proteins on multiple levels

A potential connection between physico-chemical properties of mRNAs and cognate proteins, with implications concerning both the origin of the genetic code and mRNA–protein interactions, is unexplored. We compare pyrimidine content of naturally occurring mRNA coding sequences with the propensity of cognate protein sequences to interact with pyrimidines. The latter is captured by polar requirement, a measure of solubility of amino acids in aqueous solutions of pyridines, heterocycles closely related to pyrimidines. We find that the higher the pyrimidine content of an mRNA, the stronger the average propensity of its cognate protein’s amino acids to interact with pyridines. Moreover, window-averaged pyrimidine profiles of individual mRNAs strongly mirror polar-requirement profiles of cognate protein sequences. For example, 4953 human proteins exhibit a correlation between the two with |R| > 0.8. In other words, pyrimidine-rich mRNA regions quantitatively correspond to regions in cognate proteins containing residues soluble in pyrimidine mimetics and vice versa. Finally, by studying randomized genetic code variants we show that the universal genetic code is highly optimized to preserve these correlations. Overall, our findings redefine the stereo-chemical hypothesis concerning code’s origin and provide evidence of direct complementary interactions between mRNAs and cognate proteins before development of ribosomal decoding, but also presently, especially if both are unstructured.     

Modification of orthogonal tRNAs: unexpected consequences for sense codon reassignment

Breaking the degeneracy of the genetic code via sense codon reassignment has emerged as a way to incorporate multiple copies of multiple non-canonical amino acids into a protein of interest. Here, we report the modification of a normally orthogonal tRNA by a host enzyme and show that this adventitious modification has a direct impact on the activity of the orthogonal tRNA in translation. We observed nearly equal decoding of both histidine codons, CAU and CAC, by an engineered orthogonal M. jannaschii tRNA with an AUG anticodon: tRNA^Opt. We suspected a modification of the tRNA^Opt_AUG anticodon was responsible for the anomalous lack of codon discrimination and demonstrate that adenosine 34 of tRNA^Opt_AUG is converted to inosine. We identified tRNA^Opt_AUG anticodon loop variants that increase reassignment of the histidine CAU codon, decrease incorporation in response to the histidine CAC codon, and improve cell health and growth profiles. Recognizing tRNA modification as both a potential pitfall and avenue of directed alteration will be important as the field of genetic code engineering continues to infiltrate the genetic codes of diverse organisms.     

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.     

Modular pathways for editing non-cognate amino acids by human cytoplasmic leucyl-tRNA synthetase

To prevent potential errors in protein synthesis, some aminoacyl-transfer RNA (tRNA) synthetases have evolved editing mechanisms to hydrolyze misactivated amino acids (pre-transfer editing) or misacylated tRNAs (post-transfer editing). Class Ia leucyl-tRNA synthetase (LeuRS) may misactivate various natural and non-protein amino acids and then mischarge tRNA^Leu. It is known that the fidelity of prokaryotic LeuRS depends on multiple editing pathways to clear the incorrect intermediates and products in the every step of aminoacylation reaction. Here, we obtained human cytoplasmic LeuRS (hcLeuRS) and tRNA^Leu (hctRNA^Leu) with high activity from Escherichia coli overproducing strains to study the synthetic and editing properties of the enzyme. We revealed that hcLeuRS could adjust its editing strategy against different non-cognate amino acids. HcLeuRS edits norvaline predominantly by post-transfer editing; however, it uses mainly pre-transfer editing to edit α-amino butyrate, although both amino acids can be charged to tRNA^Leu. Post-transfer editing as a final checkpoint of the reaction was very important to prevent mis-incorporation in vitro. These results provide insight into the modular editing pathways created to prevent genetic code ambiguity by evolution.     

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.     

Genetic code development by stop codon takeover

A novel theoretical consideration of the origin and evolution of the genetic code is presented. Code development is viewed from the perspective of simultaneously evolving codons, anticodons and amino acids. Early code structure was determined primarily by thermodynamic stability considerations, requiring simplicity in primordial codes. More advanced coding stages could arise as biological systems became more complex and precise in their replication. To be consistent with these ideas, a model is described in which codons become permanently associated with amino acids only when a codon-anticodon pairing is strong enough to permit rapid translation. Hence all codons are essentially chain-termination or "stop" codons until tRNA adaptors evolve having the ability to bind tightly to them. This view, which draws support from several lines of evidence, differs from the prevalent thinking on code evolution which holds that codons specifying newer amino acids were derived from codons encoding older amino acids.     

Parallel evolution of the genetic code in arthropod mitochondrial genomes

The genetic code provides the translation table necessary to transform the information contained in DNA into the language of proteins. In this table, a correspondence between each codon and each amino acid is established: tRNA is the main adaptor that links the two. Although the genetic code is nearly universal, several variants of this code have been described in a wide range of nuclear and organellar systems, especially in metazoan mitochondria. These variants are generally found by searching for conserved positions that consistently code for a specific alternative amino acid in a new species. We have devised an accurate computational method to automate these comparisons, and have tested it with 626 metazoan mitochondrial genomes. Our results indicate that several arthropods have a new genetic code and translate the codon AGG as lysine instead of serine (as in the invertebrate mitochondrial genetic code) or arginine (as in the standard genetic code). We have investigated the evolution of the genetic code in the arthropods and found several events of parallel evolution in which the AGG codon was reassigned between serine and lysine. Our analyses also revealed correlated evolution between the arthropod genetic codes and the tRNA-Lys/-Ser, which show specific point mutations at the anticodons. These rather simple mutations, together with a low usage of the AGG codon, might explain the recurrence of the AGG reassignments.