Structuring of the genetic code took place at acidic pH

I have observed that in multiple regression the number of codons specifying amino acids in the genetic code is positively correlated with the isoelectric point of amino acids and their molecular weight. Therefore basic amino acids are, on average, codified in the genetic code by a larger number of codons, which seems to imply that the genetic code originated in an acidic 'intracellular' environment. Moreover, I compare the proteins from Picrophilus torridus and Thermoplasma volcanium, which have different intracellular pH and I define the ranks of acidophily for the amino acids. A simple index of acidophily (AI), which can be easily obtained from acidophily ranks, can be associated to any protein and, therefore, can also be associated to the genetic code if the number of synonymous codons attributed to the amino acids in the code is assumed to be the frequency with which the amino acids appeared in ancestral proteins. Finally, the sampling of the variable AI among organisms having an intracellular pH less than or equal to 6.6 and those having a non-acidic intracellular pH leads to the conclusion that the value of the genetic code's AI is not typical of proteins of the latter organisms. As the genetic code's AI value is also statistically not different from that of proteins of the organisms having an acidic intracellular pH, this supports the hypothesis that the structuring of the genetic code took place in acidic pH conditions.     

The universal ancestor and the ancestors of Archaea and Bacteria were anaerobes whereas the ancestor of the Eukarya domain was an aerobe

The use of an oxyphobic index (OI) based on the propensity of amino acids to enter more frequently the proteins of anaerobes makes it possible to make inferences on the environment in which the last universal common ancestor (LUCA) lived. The reconstruction of the ancestral sequences of proteins using a method based on maximum likelihood and their attribution by means of the OI to the set of aerobe or anaerobe sequences has led to the following conclusions: the LUCA was an anaerobic 'organism', as were the ancestors of Archaea and Bacteria, whereas the ancestor of Eukarya was an aerobe. These observations seem to falsify the hypothesis that the LUCA was an aerobe and help to identify better the environment in which the first organisms lived.     

A methanogen hosted the origin of the genetic code

A comparison is made between orthologous proteins from a methanogen (Methanopyrus kandleri) and from a non-methanogen (Pyrococcus abyssi) in order to determine the amino acid substitution pattern. This analysis makes it possible to establish which amino acids are significantly and asymmetrically utilised by these two organisms. A methanophily index (MI) based on this asymmetry makes it possible for any protein to be associated with a numerical value which, when calculated for the same orthologous protein from methanogenic and non-methanogenic organisms, turns out to have the power to discriminate between these two groups of organisms, even if only for about 20% of the analysed proteins. The MI can also be associated to the genetic code under the assumption that the frequency of synonymous codons specifying the amino acids in the genetic code also reflects the frequency with which amino acids appeared in ancestral proteins. Finally a t-test shows that the MI value associated to the genetic code is not different from the mean value of the MI deriving from methanogen proteins, but it differs from the mean MI of non-methanogen proteins. This might indicate that the genetic code evolved in a methanogenic ‘organism’.     

The late stage of genetic code structuring took place at a high temperature

The correlation between the optimal growth temperature of organisms and a thermophily index based on the propensity of amino acids to enter more frequently into (hyper)thermophile proteins is used to conduct an analysis aiming to establish whether genetic code structuring took place at a low or a high temperature. If the number of codons attributed to the various amino acids in the genetic code constitutes an estimate of the mean amino acid composition of proteins produced when the genetic code was definitively structured, then the thermophily index can also be associated to the genetic code. This value and the sampling of the variable thermophily index of different alignments of protein sequences from mesophile, thermophile and hyperthermophile species make it possible to establish, with an extremely high statistical confidence, that the late stage of genetic code structuring took place in a hyperthermophile (or thermophile) 'organism'. Moreover the 95% confidence interval of the temperature at which the genetic code was fixed turned out to be 91+/-24 degrees C. These observations seem to support the hypothesis that the origin of life might have taken place at a high temperature.     

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.     

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.     

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.     

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.     

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.     

Applications of Genetic Code Expansion in Studying Protein Post-translational Modification

Various post-translational modifications can naturally occur on proteins, regulating the activity, subcellular localization, interaction, or stability of the proteins. However, it can be challenging to decipher the biological implication or physiological roles of site-specific modifications due to their dynamic and sub-stoichiometric nature. Genetic code expansion method, relying on an orthogonal aminoacyl-tRNA synthetase/tRNA pair, enables site-specific incorporation of non-canonical amino acids. Here we focus on the application of genetic code expansion to study site-specific protein post-translational modification in vitro and in vivo. After a brief introduction, we discuss possibilities of incorporating non-canonical amino acids containing post-translational modifications or their mimics into target proteins. This approach is applicable for Ser/Thr/Tyr phosphorylation, Tyr sulfation/nitration/hydroxylation, Lys acetylation/acylation, Lys/His mono-methylation, as well as Arg citrullination. The next section describes the use of a precursor non-canonical amino acid followed by chemical and/or enzymatic reactions to afford the desired modification, such as Cys/Lys acylation, ubiquitin and ubiquitin-like modifications, as well as Lys/Gln methylation. We also discuss means for functional regulation of enzymes involving in post-translational modifications through genetically incorporated non-canonical amino acids. Lastly, the limitations and perspectives of genetic code expansion in studying protein post-translational modification are described.     

Physico-chemical constraints connected with the coding properties of the genetic system

New insights on the origin of the genetic code, based on the analysis of the physico-chemical properties of its molecular constituents (RNA and amino acids), are reported in this paper. We point out a symmetry in the genetic code table and show that it can be explained by the nature of the anticodon-codon interaction. The importance of the strength of this interaction is examined and a correlation is found between the free-energy change (DeltaG(0)) of anticodon-codon association and the volume of the corresponding amino acids. This correlation is investigated in conjunction with the well-known one linking the hydrophobicity of the anticodons with that of the amino acids. We show that they can be considerated separately and that the energy vs. volume correlation may be explained by the process implicating the peptide bond formation between two successive amino acids during translation. This interpretation is supported by a statistical pattern of bases (purines or pyrimidines), observed in present coding genes, and by considerations involving the availability of the different kinds of amino acids. Finally, we try to explain the hydrophobicity correlation when reconstructing the events at the time of the so-called "RNA World". The whole of our investigation shows that the genetic code might be sufficiently robust to exist without the participation of pre-existing proteins, and that this robustness is a consequence of the physico-chemical properties of the four bases of the genetic system.     

Code dependent conservation of the physicochemical 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.This work was partially supported by OEA and Departamento de Desarrollo de la Investigación.     

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.     

Origins of gene, genetic code, protein and life: comprehensive view of life systems from a GNC-SNS primitive genetic code hypothesis

We have investigated the origin of genes, the genetic code, proteins and life using six indices (hydropathy, α-helix, β-sheet and β-turn formabilities, acidic amino acid content and basic amino acid content) necessary for appropriate three-dimensional structure formation of globular proteins. From the analysis of microbial genes, we have concluded that newly-born genes are products of nonstop frames (NSF) on antisense strands of microbial GC-rich genes [GC-NSF(a)] and from SNS repeating sequences [(SNS)_n] similar to the GC-NSF(a) (S and N mean G or C and either of four bases, respectively). We have also proposed that the universal genetic code used by most organisms on the earth presently could be derived from a GNC-SNS primitive genetic code. We have further presented the [GADV]-protein world hypothesis of the origin of life as well as a hypothesis of protein production, suggesting that proteins were originally produced by random peptide formation of amino acids restricted in specific amino acid compositions termed as GNC-, SNS and GC-NSF(a)-0th order structures of proteins. The [GADV]-protein world hypothesis is primarily derived from the GNC-primitive genetic code hypothesis. It is also expected that basic properties of extant genes and proteins could be revealed by considerations based on the scenario with four stages This review is a modified English version of the paper, which was written in Japanese and published inViva Origino 2001 29 66–85.     

On How Many Fundamental Kinds of Cells are Present on Earth: Looking for Phylogenetic Traits that Would Allow the Identification of the Primary Lines of Descent

The phylogenetic analyses as far as the identification of the number of domains of life is concerned have not reached a clear conclusion. In the attempt to improve this circumstance, I introduce the concept that the amino acids codified in the genetic code might be of markers with outstanding phylogenetic power. In particular, I hypothesise the existence of a biosphere populated, for instance, by three groups of organisms having different genetic codes because codifying at least a different amino acid. Evidently, these amino acids would mark the proteins that are present in the three groups of organisms in an unambiguous way. Therefore, in essence, this mark would not be other than the one that we usually try to make in the phylogenetic analyses in which we transform the protein sequences in phylogenetic trees, for the purpose to identify, for example, the domains of life. Indeed, this mark would allow to classify proteins without performing phylogenetic analyses because proteins belonging to a group of organisms would be recognisable as marked in a natural way by at least a different amino acid among the diverse groups of organisms. This conceptualisation answers the question of how many fundamental kinds of cells have evolved from the Last Universal Common Ancestor (LUCA), as the genetic code has unique proprieties that make the codified amino acids excellent phylogenetic markers. The presence of the formyl-methionine only in proteins of bacteria would mark them and would identify these as domain of life. On the other hand, the presence of pyrrolysine in the genetic code of the euryarchaeota would identify them such as another fundamental kind of cell evolved from the LUCA. Overall, the phylogenetic distribution of formyl-methionine and pyrrolysine would identify at least two domains of life—Bacteria and Archaea—but their number might be actually four; that is to say, Bacteria, Euryarchaeota, archeobacteria that are not euryarchaeota and Eukarya. The usually accepted domains of life represented by Bacteria, Archaea and Eukarya are not compatible with the phylogenetic distribution of these two amino acids and therefore this last classification might be mistaken.     

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.     

The Historical Factor: The Biosynthetic Relationships Between Amino Acids and Their Physicochemical Properties in the Origin of the Genetic Code

Two forces are in general, hypothesized to have influenced the origin of the organization of the genetic code: the physicochemical properties of amino acids and their biosynthetic relationships. In view of this, we have considered a model incorporating these two forces. In particular, we have studied the optimization level of the physicochemical properties of amino acids in the set of amino acid permutation codes that respects the biosynthetic relationships between amino acids. Where the properties of amino acids are represented by polarity and molecular volume we obtain indetermination percentages in the organization of the genetic code of approximately 40%. This indicates that the contingent factor played a significant role in structuring the genetic code. Furthermore, this result is in agreement with the genetic code coevolution hypothesis, which attributes a merely ancillary role to the properties of amino acids while it suggests that it was their biosynthetic relationships that organized the code. Furthermore, this result does not favor the stereochemical models proposed to explain the origin of the genetic code. On the other hand, where the properties of amino acids are represented by polarity alone, we obtain an indetermination percentage of at least 21.5%. This might suggest that the polarity distances played an important role and would therefore provide evidence in favor of the physicochemical hypothesis of genetic code origin. Although, overall, the analysis might have given stronger support to the latter hypothesis, this did not actually occur. The results are therefore discussed in the context of the different theories proposed to explain the origin of the genetic code. Received: 10 September 1996 / Accepted: 3 March 1997     

Neutral adaptation of the genetic code to double-strand coding

Summary We lay new foundations to the hypothesis that the genetic code is adapted to evolutionary retention of information in the antisense strands of natural DNA/RNA sequences. In particular, we show that the genetic code exhibits, beyond the neutral replacement patterns of amino acid substitutions, optimal properties by favoring simultaneous evolution of proteins encoded in DNA/RNA sense-antisense strands. This is borne out in the sense-antisense transformations of the codons of every amino acid which target amino acids physicochemically similar to each other. Moreover, silent mutations in the sense strand generate conservative ones in its antisense counterpart and vice versa. Coevolution of proteins coded by complementary strands is shown to be a definite possibility, a result which does not depend on any physical interaction between the coevolving proteins. Likewise, the degree to which the present genetic code is dedicated to evolutionary sense-antisense tolerance is demonstrated by comparison with many randomized codes. Double-strand coding is quantified from an information-theoretical point of view.     

Recent Technologies for Genetic Code Expansion and their Implications on Synthetic Biology Applications

Genetic code expansion (GCE) enables the site-specific incorporation of non-canonical amino acids as novel building blocks for the investigation and manipulation of proteins. The advancement of genetic code expansion has been benefited from the development of synthetic biology, while genetic code expansion also helps to create more synthetic biology tools. In this review, we summarize recent advances in genetic code expansion brought by synthetic biology progresses, including engineering of the translation machinery, genome-wide codon reassignment, and the biosynthesis of non-canonical amino acids. We highlight the emerging application of this technology in construction of new synthetic biology parts, circuits, chassis, and products.