The genetic code as a clue to understanding of molecular evolution

The genetic code is comprised of a system concerning the distribution of doublets of the first two codon bases among amino acids. According to this system a definite order in the relative distribution of the first and the second codon bases coincides with a definite order among the common amino acids and their distribution for the number of hydrogen atoms per molecule (an unexpected parameter). The pattern of the relative distribution of the first and the second codon bases suggests it originated from a crystalline-like structure in which the set of bases AUGC served as an elementary structural unit and the base doublets played the role of structural analogs to the amino acids. These hypothetical crystalline-like aggregates are composed of the free molecules of amino acids and bases, and although different in their composition, should have an even number of hydrogen atoms per standard structural module.     

Experimental studies of molecular interactions between nitrogen bases of nucleic acids.

New experimental results concerning molecular interactions between the nitrogen bases of nucleic acids in the crystalline phase and in vacuo are reported. The temperature dependence of the evaporation rate is measured for solid species. The sensitivity of conventional methods of sublimation heat measurements was improved essentially using a quartz resonator serving as a precise sensor of evaporation rate. Sublimation heats were found for both canonical bases and a number of their derivatives. The in vacuo formation of base associates interacting through hydrogen bonds was observed with a field mass spectrometer. The dimer formation enthalpies, which are indicative of a stronger attraction in complementary pairs compared with noncomplementary ones, were derived from the temperature dependence of ionic currents. Hydrogen-bound complexes of more intricate associates (base trimers and aqueous molecules associates) were studied. The energy gain in the formation of trimers of identical molecules was shown to be larger (per base molecule) than that for dimers.     

Comparisons of the polypeptide chains of globins

Summary Vertebrate myoglobins and hemoglobins each consist of units containing a heme group and an associated polypeptide chain. The polypeptide chains form an homologous series and can be compared with each other to measure amino acid differences and minimum base differences per codon. These differences are the result of mutations which have been incorporated during divergent evolution from a common ancestral gene. Each such base replacement is termed an evolutionary event. Each amino acid replacement is the result of one or more evolutionary events. However there can be only one amino acid difference between two sites. The minimum mutational difference between two sites is perceptible only as one, two or three base differences per codon, but there may be more evolutionary events than base differences, because of revertants and multiplicity of base replacements at the same site. When all three bases in a codon are changed, the result is recognizable in only about five per cent of cases. Therefore such recognizable three-base changes indicate a large number of evolutionary events.     

Evolution of the Genetic Triplet Code via Two Types of Doublet Codons

Explaining the apparent non-random codon distribution and the nature and number of amino acids in the ‘standard’ genetic code remains a challenge, despite the various hypotheses so far proposed. In this paper we propose a simple new hypothesis for code evolution involving a progression from singlet to doublet to triplet codons with a reading mechanism that moves three bases each step. We suggest that triplet codons gradually evolved from two types of ambiguous doublet codons, those in which the first two bases of each three-base window were read (‘prefix’ codons) and those in which the last two bases of each window were read (‘suffix’ codons). This hypothesis explains multiple features of the genetic code such as the origin of the pattern of four-fold degenerate and two-fold degenerate triplet codons, the origin of its error minimising properties, and why there are only 20 amino acids. Reviewing Editor: Dr. Laura Landweber An erratum to this article can be found at .     

不具有3-碱基周期性的编码序列初探

对120个较短编码序列(＜1 200 bp)的Fourier频谱进行分析表明,3-碱基周期性在短编码序列中并不是绝对存在的.统计分析提示,编码序列有无3-碱基周期性与序列的碱基组成和分布、所编码蛋白质氨基酸的选用和顺序以及同义密码子的使用都有一定的关系.一般地,非周期-3序列中A+U含量高于G+C含量,周期-3序列的情况则相反;非周期-3序列中碱基在密码子三个位点上的分布比周期-3序列中的分布均匀;非周期-3序列密码子和氨基酸的使用偏向没有周期-3序列的大.在利用Fourier分析方法预测DNA序列中的基因和外显子时,应充分考虑到这些现象.     

On the Structural Regularity in Nucleobases and Amino Acids and Relationship to the Origin and Evolution of the Genetic Code

To explore how chemical structures of both nucleobases and amino acids may have played a role in shaping the genetic code, numbers of sp² hybrid nitrogen atoms in nucleobases were taken as a determinative measure for empirical stereo-electronic property to analyze the genetic code. Results revealed that amino acid hydropathy correlates strongly with the sp² nitrogen atom numbers in nucleobases rather than with the overall electronic property such as redox potentials of the bases, reflecting that stereo-electronic property of bases may play a role. In the rearranged code, five simple but stereo-structurally distinctive amino acids (Gly, Pro, Val, Thr and Ala) and their codon quartets form a crossed intersection “core”. Secondly, a re-categorization of the amino acids according to their β-carbon stereochemistry, verified by charge density (at β-carbon) calculation, results in five groups of stereo-structurally distinctive amino acids, the group leaders of which are Gly, Pro, Val, Thr and Ala, remarkably overlapping the above “core”. These two lines of independent observations provide empirical arguments for a contention that a seemingly “frozen” “core” could have formed at a certain evolutionary stage. The possible existence of this codon “core” is in conformity with a previous evolutionary model whereby stereochemical interactions may have shaped the code. Moreover, the genetic code listed in UCGA succession together with this codon “core” has recently facilitated an identification of the unprecedented icosikaioctagon symmetry and bi-pyramidal nature of the genetic code.     

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.     

A novel complex mutation in the LDL receptor gene probably caused by the simultaneous occurrence of deletion and insertion in the same region

A novel complex mutation with the presence of both deletion and insertion in very close proximity in the same region was detected in exon 8 of the LDL receptor gene from two apparently unrelated Japanese families with familial hypercholesterolemia (FH). In this mutant LDL receptor gene, the nine bases from nucleotide (nt) 1115 to nt 1123 (AGGGTGGCT) were replaced by six different bases (CACTGA), and consequently the four amino acids from codon 351 to 354, Glu-Gly-Gly-Tyr, were replaced by three amino acids, Ala-Leu-Asn, in the conserved amino acid region of the growth factor repeat B of the LDL receptor. The nature of the amino acid substitution and data on the families suggest that this mutation is very likely to affect the LDL receptor function and cause FH. The generation of this complex mutation can be explained by the simultaneous occurrence of deletion and insertion through the formation of a hairpin-loop structure mediated by inverted repeat sequences. Thus this mutation supports the hypothesis that inverted repeat sequences influence the stability of a given gene and promote human gene mutations.     

Codon bias in Escherichia coli: the influence of codon context on mutation and selection.

The codon bias in Escherichia coli for all two-fold degenerate amino acids was studied as dependent on the context from the six bases in the nearest surrounding codons. By comparing the results in genes at different expression levels, effects that are due to differences in mutation rates can be distinguished from those that are due to selection. Selective effects on the codon bias is found mostly from the first neighbouring base in the 3'direction, while neighbouring bases further away influence mostly the mutational bias. In some cases it is also possible to identify specific molecular processes, repair or avoidance of frame shift, that lead to the context dependence of the bias.     

Codon substitution models based on residue similarity and their applications

Codon models are now widely used to draw evolutionary inferences from alignments of homologous sequence data. Incorporating physicochemical properties of amino acids into codon models, two novel codon substitution models describing the evolution of protein-coding DNA sequences are presented based on the similarity scores of amino acids. To describe substitutions between codons a continue-time Markov process is used. Transition/transversion rate bias and nonsynonymous codon usage bias are allowed in the models. In our implementation, the parameters are estimated by maximum-likelihood (ML) method as in previous studies. Furthermore, instantaneous mutations involving more than one nucleotide position of a codon are considered in the second model. Then the two suggested models are applied to five real data sets. The analytic results indicate that the new codon models considering physicochemical properties of amino acids can provide a better fit to the data comparing with existing codon models, and then produce more reliable estimates of certain biologically important measures than existing methods.     

简并性还是多态性?——遗传密码子设定的进化特点的剖析(英文)

遗传密码子的设定表现出令人困惑的多态性特点 :不同氨基酸拥有的密码子的数目 ,除 5个外 ,从 1个到 6个都有 .这种特点显示出密码子无论在翻译行为还是进化轨迹上 ,都存在诸多的异质性 .因此 ,简并性一词的收敛含义 ,并不能表征这种多态性的进化内涵 .没有同义密码子的AUG(Met)和UGG (Trp)并无简并现象 .其余的密码子则可分为两大类 :一类是 ,4个同义密码子为 1组 ,具有相同的第 1、2位碱基 ,并遵循“3中读 2”的读出规则 .同组的 4个同义密码子 ,不过是来自同一个双字母原始密码子 (XYN)的孑遗物 ,从这个意义上讲 ,也不宜视为简并现象 ;另一类则主要是 ,2个同义密码子为一组 ,并遵循“3中读 3”读出规则 .它们是由编码 2个氨基酸的双义原始密码子 ,第 3位的未定碱基N进一步设定形成 .至于有 6个同义密码子的 ,特别令人困感不解的组别 ,实际上是 4 + 2个 ,这启示它们可能源于上述两大类 .遗传密码子多态性的起源 ,可能始于最初阶段 ,氨基酸同某类寡核苷酸的起始二联体的相互作用 ,而完成于所有的双义原始密码子的第 3位碱基的分化 .这种进化轨迹被传统的简并性一词所模糊 ,并导致鉴定各有关理论可信性的坚实依据和令不同观点取得共识的基础被掩盖起来 .这可能就是在遗传密码子起源领域里 ,长期存在着众     

The code within the codons

For the first time it is shown that each of the three codon bases has a general correlation with a different, predictable amino acid property, depending on position within the codon. In addition to the previously recognized link between the mid-base and the hydrophobic-hydrophilic spectrum, we show that, with the exception of G, the first base is generally invariant within a synthetic pathway. G--coded amino acids show a different order, being found only at the head of the synthetic pathways. The redundancy of the nature of the third base has a previously unrecognised relationship with molecular weight. The bases U and A (transversions) are associated with the most sharply defined or opposite states in both the first and second position, C somewhat less so or intermediate, anf G neutral. The apparently systematic nature of these relationships has profound implications for the origin of the genetic code. It appears to be the remains of the first language of the cell, predating the tRNA/ribosome system, persisting with remarkably little change at a deeper level of organisation than the codon language.     

Analysis of distribution of bases in the coding sequences by a diagrammatic technique.

The frequencies of occurrence of four bases in the first, second and third codon positions and in the total coding sequences have been calculated by the codon usage table published in 1990 by Ikemura et al. The distribution of frequencies are further analysed in detail by a graphic technique presented recently by us. Formulas expressing the frequencies of four bases in the first and second codon positions in terms of frequencies of amino acids have been given. It is shown by the graphic analysis that for 90 species, in the first codon position the purine bases are dominant and in most cases G is the most dominant base. In the second codon position A is the most dominant base, while G is the least dominant base. In the third codon position the G + C content varies from 0.1 to 0.9, keeping the A + C content equal to 1/2 and G content equal to that of C, approximately. If the frequencies for bases A, C, G and U in the total coding sequences are denoted by a, c, g and u, respectively, it is found that the unequal formula: a2 + c2 + g2 + u2 less than 1/3, is valid for each of the 90 species including the human and E.coli etc.     

Codon equilibrium I: Testing for homogeneous equilibrium

Summary We present theoretical considerations that suggest that synonymous-codon usage might be expected to be close to an equilibrium distribution given a very homogeneous process of silent substitution. By homogeneous we mean that substitution depends only on the two bases involved, so that 12 base-substitution rates completely describe the silent substitution process. We have developed a method of statistically testing for such homogeneous equilibrium and applied it to reported data on the codon usages of different classes of organisms. Weakly expressed bacterial sequences and both mammalian and nonmammalian eukaryotic sequences deviate significantly from a random pattern of codon usage, in the direction of homogeneous equilibrium. On the other hand, highly expressed bacterial sequences do not exhibit homogeneous equilibrium, which may be correlated with recent experimental results showing that they are optimized to accept the most abundant tRNAs. To examine the effect of amino acid replacements on the homogeneous model of silent substitution, we divided the amino acids with degenerate codes into two classes, those with high mutabilities and those with low, and performed the same analysis on bacterial and eukaryotic data sets. The codon sets of the highly mutable class of amino acids are not further from homogeneous equilibrium than are the codon sets of the class with low mutabilities. We also found for the eukaryotic data that these independent classes of codon sets show very similar equilibrium patterns. The various results suggest a high level of uniformity in the process of silent fixation in the different synonymous-codon sets, especially in eukaryotes.     

Construction of "small-intelligent" focused mutagenesis libraries using well-designed combinatorial degenerate primers

Site-saturation mutagenesis is a powerful tool for protein optimization due to its efficiency and simplicity. A degenerate codon NNN or NNS (K) is often used to encode the 20 standard amino acids, but this will produce redundant codons and cause uneven distribution of amino acids in the constructed library. Here we present a novel "small-intelligent" strategy to construct mutagenesis libraries that have a minimal gene library size without inherent amino acid biases, stop codons, or rare codons of Escherichia coli by coupling well-designed combinatorial degenerate primers with suitable PCR-based mutagenesis methods. The designed primer mixture contains exactly one codon per amino acid and thus allows the construction of small-intelligent mutagenesis libraries with one gene per protein. In addition, the software tool DC-Analyzer was developed to assist in primer design according to the user-defined randomization scheme for library construction. This small-intelligent strategy was successfully applied to the randomization of halohydrin dehalogenases with one or two randomized sites. With the help of DC-Analyzer, the strategy was proven to be as simple as NNS randomization and could serve as a general tool to efficiently randomize target genes at positions of interest.     

A model for the prebiotic synthesis of peptides which throws light on the origin of the genetic code and the observed chirality of life

A model is proposed which indicates that peptides could have been synthesized on RNA templates. L-amino acids are selectively trapped and orientated in clefts in an apparently stable RNA structure attached to a solid phase. Each cleft is bounded by three bases which correspond to the codon for that amino acid. Rotation during a dehydration phase orientates the amino acids so that peptide bond formation can occur by condensation. Rehydration releases the formed peptides.     

Effect of temperature and ATP supply on the efficiency of programmed nonsense suppression

Chemical diversity of protein molecules can be expanded through in vitro incorporation of unnatural amino acids in response to a nonsense codon. Chemically misacylated tRNAs are used for tethering unnatural amino acids to a nonsense-mutated target codon (nonsense suppression). In the course of experiments to introduce S-(2-nitrobenzyl)cysteine (NBC) into a targeted location of human erythropoietin, we found that NBC incorporates more efficiently at lower temperatures. In addition, at a fixed reaction temperature, more NBC was incorporated with a reduced supply of ATP. Since the rate of peptide elongation was remarkably higher at the elevated temperature or with enhanced supply of ATP, these results indicate that the efficiency of nonsense suppression is inversely correlated to the peptide elongation rate. Therefore, maximal yield of nonsense-suppressed proteins is obtained at a compromised elongation rate. The present result will offer a primary guideline to optimize the reaction conditions for in vitro production of protein molecules containing unnatural amino acids.     

Codon equilibrium I: Testing for homogeneous equilibrium

We present theoretical considerations that suggest that synonymous-codon usage might be expected to be close to an equilibrium distribution given a very homogeneous process of silent substitution. By homogeneous we mean that substitution depends only on the two bases involved, so that 12 base-substitution rates completely describe the silent substitution process. We have developed a method of statistically testing for such homogeneous equilibrium and applied it to reported data on the codon usages of different classes of organisms. Weakly expressed bacterial sequences and both mammalian and nonmammalian eukaryotic sequences deviate significantly from a random pattern of codon usage, in the direction of homogeneous equilibrium. On the other hand, highly expressed bacterial sequences do not exhibit homogeneous equilibrium, which may be correlated with recent experimental results showing that they are optimized to accept the most abundant tRNAs. To examine the effect of amino acid replacements on the homogeneous model of silent substitution, we divided the amino acids with degenerate codes into two classes, those with high mutabilities and those with low, and performed the same analysis on bacterial and eukaryotic data sets. The codon sets of the highly mutable class of amino acids are not further from homogeneous equilibrium than are the codon sets of the class with low mutabilities. We also found for the eukaryotic data that these independent classes of codon sets show very similar equilibrium patterns. The various results suggest a high level of uniformity in the process of silent fixation in the different synonymous-codon sets, especially in eukaryotes.     

The triplet code from first principles

Temporal order ("chronology") of appearance of amino acids and their respective codons on evolutionary scene is reconstructed. A consensus chronology of amino acids is built on the basis of 60 different criteria each offering certain temporal order. After several steps of filtering the chronology vectors are averaged resulting in the consensus order: G, A, D, V, P, S, E, (L, T), R, (I, Q, N), H, K, C, F, Y, M, W. It reveals two important features: the amino acids synthesized in imitation experiments of S. Miller appeared first, while the amino acids associated with codon capture events came last. The reconstruction of codon chronology is based on the above consensus temporal order of amino acids, supplemented by the stability and complementarity rules first suggested by M. Eigen and P. Schuster, and on the earlier established processivity rule. At no point in the reconstruction the consensus amino-acid chronology was in conflict with these three rules. The derived genealogy of all 64 codons suggested several important predictions that are confirmed. The reconstruction of the origin and evolutionary history of the triplet code becomes, thus, a powerful research tool for molecular evolution studies, especially in its early stages.