The joint evolution of defence and inducibility against natural enemies

The aminoacyl-tRNA synthetases (aaRSs) ensure the fidelity of the translation of the genetic code, covalently attaching appropriate amino acids to the corresponding nucleic acid adaptor molecules-tRNA. The fundamental role of aminoacylation reaction catalysed by aaRSs implies that representatives of the family are thought to be among the earliest proteins to appear. Based on sequence analysis and catalytic domain structure, aaRSs have been partitioned into two classes of 10 enzymes each. However, based on the structural and sequence data only, it will not be easily understood that the present partitioning is not governed by chance. Our findings suggest that organization of amino acid biosynthetic pathways and clustering of aaRSs into different classes are intimately related to one another. A plausible explanation for such a relationship is dictated by early link between aaRSs and amino acids biosynthetic proteins. The aaRSs catalytic cores are highly relevant to the ancient metabolic reactions, namely, amino acids and cofactors biosynthesis. In particular we show that class II aaRSs mostly associated with the primordial amino acids, while class I aaRSs are usually related to amino acids evolved lately. Reasoning from this we propose a possible chronology of genetic code evolution.     

Relationship of the CCA sequence of tRNA with the early evolutional aspect of aminoacyl-tRNA synthetases.

The CCA sequence is common to the 3'-ends of all tRNAs. We investigated the requirement of the CCA sequence in aminoacylation with the cognate aminoacyl-tRNA synthetases (aaRSs) and several interesting conclusions could be drawn. In tRNAs belonging to the class I aaRSs, decreased aminoacylation activities resulted from the substitution of A76 with a pyrimidine, whereas in tRNAs belonging to the class II aaRSs, decreased aminoacylation activities resulted from the substitution with guanine. The results suggest that aminoacylation of proto-tRNA might have started through the direct hydrophobic (or stacking) interaction between the large, hydrophobic amino acid residue (now utilizing a class I aaRS) of aminoacyl-AMP and the 3'-terminal adenine. The shorter distance between the adenine and the 2'-OH position than the 3'-OH position, and the bulkiness and hydrophobicity of amino acids may be important reasons why class I aaRSs select the 2'-OH position in aminoacylation. Molecular mechanics-based conformation modeling also indicated that the resulting positioning of the adenine and the amino acid residue of 2'-aminoacyl-adenosine for large amino acid is in the vicinity. In contrast, in the case of small amino acids (with class II aaRSs) which would not be able to use the hydrophobic interaction, a protein enzyme might have participated in the aminoacylation reaction from an early stage. The active-site folds of aaRSs belonging to each class reflect the history of evolution: typical nucleotide-binding fold (Rossman fold) in the case of class I aaRSs, and primitive fold which is found also among the family of nonribosomal peptide synthetases in the case of class II aaRSs.     

Peripheral neuropathy via mutant tRNA synthetases: Inhibition of protein translation provides a possible explanation

Recent evidence indicates that inhibition of protein translation may be a common pathogenic mechanism for peripheral neuropathy associated with mutant tRNA synthetases (aaRSs). aaRSs are enzymes that ligate amino acids to their cognate tRNA, thus catalyzing the first step of translation. Dominant mutations in five distinct aaRSs cause Charcot‐Marie‐Tooth (CMT) peripheral neuropathy, characterized by length‐dependent degeneration of peripheral motor and sensory axons. Surprisingly, loss of aminoacylation activity is not required for mutant aaRSs to cause CMT. Rather, at least for some mutations, a toxic‐gain‐of‐function mechanism underlies CMT‐aaRS. Interestingly, several mutations in two distinct aaRSs were recently shown to inhibit global protein translation in Drosophila models of CMT‐aaRS, by a mechanism independent of aminoacylation, suggesting inhibition of translation as a common pathogenic mechanism. Future research aimed at elucidating the molecular mechanisms underlying the translation defect induced by CMT‐mutant aaRSs should provide novel insight into the molecular pathogenesis of these incurable diseases.     

Mutational separation of two pathways for editing by a class I tRNA synthetase

Aminoacyl tRNA synthetases (aaRSs) catalyze the first step in protein biosynthesis, establishing a connection between codons and amino acids. To maintain accuracy, aaRSs have evolved a second active site that eliminates noncognate amino acids. Isoleucyl tRNA synthetase edits valine by two tRNA(Ile)-dependent pathways: hydrolysis of valyl adenylate (Val-AMP, pretransfer editing) and hydrolysis of mischarged Val-tRNA(Ile) (posttransfer editing). Not understood is how a single editing site processes two distinct substrates--an adenylate and an aminoacyl tRNA ester. We report here distinct mutations within the center for editing that alter adenylate but not aminoacyl ester hydrolysis, and vice versa. These results are consistent with a molecular model that shows that the single editing active site contains two valyl binding pockets, one specific for each substrate.     

Exclusive cytosolic localization and broad tRNA<Superscript>Ser</Superscript> specificity of <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis> seryl-tRNA synthetase

Aminoacyl-tRNA synthetases (aaRSs) decipher the genetic code, covalently linking amino acids to cognate tRNAs, thus preparing substrates for the process of translation. Although aaRSs funtion primarily in translation and are localized in cytosol, mitochondria and chloroplasts there are many reports on their additional functions and subcellular destinations beyond translation. However, data on plant aaRSs are scarce. Initial analysis of amino acid sequence of Arabidopsis thaliana seryl-tRNA synthetase (SerRS) suggested that protein contains putative nuclear localization signals. GFP-localization experiments in transiently transformed epidermal onion cells and Arabidopsis protoplasts gave ambiguous results because in some cells SerRS appeared to be dually localized to both cytosol and nucleus. However, data obtained on transgenic lines expressing SerRS-TAP and GFP-SerRS revealed exclusive cytosolic location of SerRS. Subcellular distribution of SerRS did not change during stress. Cytosolic Arabidopsis SerRS was expressed and purified. The enzyme efficiently aminoacylated eukaryotic and bacterial tRNAs^Ser, that are structurally very different. Given the fact that the same behavior was previously shown for monocot maize SerRS, it seems that plant SerRSs exhibit unusually broad tRNA^Ser specificity, unlike SerRSs from other organisms. Possible functional implications of this unique characteristic of plant SerRSs are discussed.     

Intra-protein Compensatory Mutations Analysis Highlights the tRNA Recognition Regions in Aminoacyl-tRNA Synthetases

Abstract

The aminoacyl-tRNA synthetases (aaRSs) covalently attach amino acids to their corresponding nucleic acid adapter molecules, tRNAs. The interactions in the tRNA-aaRSs complexes are mostly non-specific, and largely electrostatic. Tracing a way of aaRS-tRNA mutual adaptation throughout evolution offers a clearer view of understanding how aaRS-tRNA systems preserve patterns of tRNA recognition and binding. In this study, we used the compensatory mutations analysis to explore adaptation of aaRSs in respond to random mutations that can occur in the tRNA-recognition area. We showed that the frequency of compensatory mutations among residues that belong to the recognition region is 1.75-fold higher than that of the exposed residues. The highest frequencies of compensatory mutations are observed for pairs of charged residues, wherein one residue is located within the tRNA-recognition area, while the second is placed outside of the area, and contributes to the formation of the aaRS electrostatic landscape. Given charged residues are compensated by buried charge residues in more than 60% of the analyzed mutations. The cytoplasmatic and mitochondrial aaRSs preserve similar patterns of compensatory mutations in the tRNA recognition areas. Moreover, we found that mitochondrial aaRSs demonstrate a significant increase in the frequency of compensatory mutations in the area. Our findings shed light on the physical nature of compensatory mutations in aaRSs, thereby keeping unchanged tRNA-recognition patterns.     

氨酰tRNA合成酶作用机制及其应用

氨酰tRNA合成酶（aminoacyl tRNA synthetases, aaRSs）通过催化氨基酸与相应tRNA的氨酰化以保证遗传信息翻译的准确性,在生物体内具有重要作用。近年来,随着对aaRS催化机制理解的不断加深,aaRS的应用逐渐成为研究热点。在细菌中,aaRS活性被抑制后会导致其生命活动发生紊乱,根据aaRS在人体与病原菌内不同的催化特点设计针对病原体的特异性aaRS抑制剂,将有助于开发以aaRS为靶标的新型抗生素。另外,通过突变aaRS可以在蛋白质序列中定点掺入非天然氨基酸,扩展蛋白质工程。本文简述了aaRS的分类、结构与功能的特点,并在此基础上综述了aaRS在研发新型抑制剂,设计改造特殊蛋白质等方面的应用。     

Molecular Evolution of Aminoacyl tRNA Synthetase Proteins in the Early History of Life

Aminoacyl-tRNA synthetases (aaRS) consist of several families of functionally conserved proteins essential for translation and protein synthesis. Like nearly all components of the translation machinery, most aaRS families are universally distributed across cellular life, being inherited from the time of the Last Universal Common Ancestor (LUCA). However, unlike the rest of the translation machinery, aaRS have undergone numerous ancient horizontal gene transfers, with several independent events detected between domains, and some possibly involving lineages diverging before the time of LUCA. These transfers reveal the complexity of molecular evolution at this early time, and the chimeric nature of genomes within cells that gave rise to the major domains. Additionally, given the role of these protein families in defining the amino acids used for protein synthesis, sequence reconstruction of their pre-LUCA ancestors can reveal the evolutionary processes at work in the origin of the genetic code. In particular, sequence reconstructions of the paralog ancestors of isoleucyl- and valyl- RS provide strong empirical evidence that at least for this divergence, the genetic code did not co-evolve with the aaRSs; rather, both amino acids were already part of the genetic code before their cognate aaRSs diverged from their common ancestor. The implications of this observation for the early evolution of RNA-directed protein biosynthesis are discussed.     

The case for an error minimizing set of coding amino acids

The fidelity of the translation machinery largely depends on the accuracy by which the tRNAs within the living cells are charged. Aminoacyl-tRNA synthetases (aaRSs) attach amino acids to their cognate tRNAs ensuring the fidelity of translation in coding sequences. Based on the sequence analysis and catalytic domain structure, these enzymes are classified into two major groups of 10 enzymes each. In this study, we have generally tackled the role of aaRSs in decreasing the effects of mistranslations and consequently the evolution of the translation machinery. To this end, a fitness function was introduced in order to measure the accuracy by which each tRNA is charged with its cognate amino acid. Our results suggest that the aaRSs are very well optimized in "load minimization" based on their classes and their mechanisms in distinguishing the correct amino acids. Besides, our results support the idea that from an evolutionary point, a selectional pressure on the translational fidelity seems to be responsible in the occurrence of the 20 coding amino acids.     

Redesigning the stereospecificity of tyrosyl‐tRNA synthetase

d ‐Amino acids are largely excluded from protein synthesis, yet they are of great interest in biotechnology. Unnatural amino acids have been introduced into proteins using engineered aminoacyl‐tRNA synthetases (aaRSs), and this strategy might be applicable to d ‐amino acids. Several aaRSs can aminoacylate their tRNA with a d ‐amino acid; of these, tyrosyl‐tRNA synthetase (TyrRS) has the weakest stereospecificity. We use computational protein design to suggest active site mutations in Escherichia coli TyrRS that could increase its d ‐Tyr binding further, relative to l ‐Tyr. The mutations selected all modify one or more sidechain charges in the Tyr binding pocket. We test their effect by probing the aminoacyl‐adenylation reaction through pyrophosphate exchange experiments. We also perform extensive alchemical free energy simulations to obtain l ‐Tyr/d ‐Tyr binding free energy differences. Agreement with experiment is good, validating the structural models and detailed thermodynamic predictions the simulations provide. The TyrRS stereospecificity proves hard to engineer through charge‐altering mutations in the first and second coordination shells of the Tyr ammonium group. Of six mutants tested, two are active towards d ‐Tyr; one of these has an inverted stereospecificity, with a large preference for d ‐Tyr. However, its activity is low. Evidently, the TyrRS stereospecificity is robust towards charge rearrangements near the ligand. Future design may have to consider more distant and/or electrically neutral target mutations, and possibly design for binding of the transition state, whose structure however can only be modeled. Proteins 2016; 84:240–253. © 2015 Wiley Periodicals, Inc.     

Mechanism of the activation step of the aminoacylation reaction: a significant difference between class I and class II synthetases

In the present work we report, for the first time, a novel difference in the molecular mechanism of the activation step of aminoacylation reaction between the class I and class II aminoacyl tRNA synthetases (aaRSs). The observed difference is in the mode of nucleophilic attack by the oxygen atom of the carboxylic group of the substrate amino acid (AA) to the αP atom of adenosine triphosphate (ATP). The syn oxygen atom of the carboxylic group attacks the α-phosphorous atom (αP) of ATP in all class I aaRSs (except TrpRS) investigated, while the anti oxygen atom attacks in the case of class II aaRSs. The class I aaRSs investigated are GluRS, GlnRS, TyrRS, TrpRS, LeuRS, ValRS, IleRS, CysRS, and MetRS and class II aaRSs investigated are HisRS, LysRS, ProRS, AspRS, AsnRS, AlaRS, GlyRS, PheRS, and ThrRS. The variation of the electron density at bond critical points as a function of the conformation of the attacking oxygen atom measured by the dihedral angle ψ (C(α)-C') conclusively proves this. The result shows that the strength of the interaction of syn oxygen and αP is stronger than the interaction with the anti oxygen for class I aaRSs. This indicates that the syn oxygen is the most probable candidate for the nucleophilic attack in class I aaRSs. The result is further supported by the computation of the variation of the nonbonded interaction energies between αP atom and anti oxygen as well as syn oxygen in class I and II aaRSs, respectively. The difference in mechanism is explained based on the analysis of the electrostatic potential of the AA and ATP which shows that the relative arrangement of the ATP with respect to the AA is opposite in class I and class II aaRSs, which is correlated with the organization of the active site in respective aaRSs. A comparative study of the reaction mechanisms of the activation step in a class I aaRS (Glutaminyl tRNA synthetase) and in a class II aaRS (Histidyl tRNA synthetase) is carried out by the transition state analysis. The atoms in molecule analysis of the interaction between active site residues or ions and substrates are carried out in the reactant state and the transition state. The result shows that the observed novel difference in the mechanism is correlated with the organizations of the active sites of the respective aaRSs. The result has implication in understanding the experimentally observed different modes of tRNA binding in the two classes of aaRSs.     

Evolution of the aminoacyl-tRNA synthetases and the origin of the genetic code

The aminoacyl-tRNA synthetases exist as two enzyme families which were apparently generated by divergent evolution from two primordial synthetases. The two classes of enzymes exhibit intriguing familial relationships, in that they are distributed nonrandomly within the codon-amino acid matrix of the genetic code. For example, all XCX codons code for amino acids handled by class II synthetases, and all but one of the XUX codons code for amino acids handled by class I synthetases. One interpretation of these patterns is that the synthetases coevolved with the genetic code. The more likely explanation, however, is that the synthetases evolved in the context of an already-established genetic code—a code which developed earlier in an RNA world. The rules which governed the development of the genetic code, and led to certain patterns in the coding catalog between codons and amino acids, would also have governed the subsequent evolution of the synthetases in the context of a fixed code, leading to patterns in synthetase distribution such as those observed. These rules are (1) conservative evolution of amino acid and adapter binding sites and (2) minimization of the disruptive effects on protein structure caused by codon meaning changes.     

An Alternative Look at Code Evolution: Using Non-canonical Codes to Evaluate Adaptive and Historic Models for the Origin of the Genetic Code

The canonical code has been shown many times to be highly robust against point mutations; that is, mutations that change a single nucleotide tend to result in similar amino acids more often than expected by chance. There are two major types of models for the origin of the code, which explain how this sophisticated structure evolved. Adaptive models state that the primitive code was specifically selected for error minimization, while historic models hypothesize that the robustness of the code is an artifact or by-product of the mechanism of code evolution. In this paper, we evaluated the levels of robustness in existing non-canonical codes as well as codes that differ in only one codon assignment from the standard code. We found that the level of robustness of many of these codes is comparable or better than that of the standard code. Although these results do not preclude an adaptive origin of the genetic code, they suggest that the code was not selected for minimizing the effects of point mutations.     

An asymmetric underlying rule in the assignment of codons: possible clue to a quick early evolution of the genetic code via successive binary choices

Aminoacyl-tRNA synthetases (aaRSs) are responsible for creating the pool of correctly charged aminoacyl-tRNAs that are necessary for the translation of genetic information (mRNA) by the ribosome. Each aaRS belongs to either one of only two classes with two different mechanisms of aminoacylation, making use of either the 2'OH (Class I) or the 3'OH (Class II) of the terminal A76 of the tRNA and approaching the tRNA either from the minor groove (2'OH) or the major groove (3'OH). Here, an asymmetric pattern typical of differentiation is uncovered in the partition of the codon repertoire, as defined by the mechanism of aminoacylation of each corresponding tRNA. This pattern can be reproduced in a unique cascade of successive binary decisions that progressively reduces codon ambiguity. The deduced order of differentiation is manifestly driven by the reduction of translation errors. A simple rule can be defined, decoding each codon sequence in its binary class, thereby providing both the code and the key to decode it. Assuming that the partition into two mechanisms of tRNA aminoacylation is a relic that dates back to the invention of the genetic code in the RNA World, a model for the assignment of amino acids in the codon table can be derived. The model implies that the stop codon was always there, as the codon whose tRNA cannot be charged with any amino acid, and makes the prediction of an ultimate differentiation step, which is found to correspond to the codon assignment of the 22nd amino acid pyrrolysine in archaebacteria.     

The optimality of the standard genetic code assessed by an eight-objective evolutionary algorithm

Background

The standard genetic code (SGC) is a unique set of rules which assign amino acids to codons. Similar amino acids tend to have similar codons indicating that the code evolved to minimize the costs of amino acid replacements in proteins, caused by mutations or translational errors. However, if such optimization in fact occurred, many different properties of amino acids must have been taken into account during the code evolution. Therefore, this problem can be reformulated as a multi-objective optimization task, in which the selection constraints are represented by measures based on various amino acid properties.

Results

To study the optimality of the SGC we applied a multi-objective evolutionary algorithm and we used the representatives of eight clusters, which grouped over 500 indices describing various physicochemical properties of amino acids. Thanks to that we avoided an arbitrary choice of amino acid features as optimization criteria. As a consequence, we were able to conduct a more general study on the properties of the SGC than the ones presented so far in other papers on this topic. We considered two models of the genetic code, one preserving the characteristic codon blocks structure of the SGC and the other without this restriction. The results revealed that the SGC could be significantly improved in terms of error minimization, hereby it is not fully optimized. Its structure differs significantly from the structure of the codes optimized to minimize the costs of amino acid replacements. On the other hand, using newly defined quality measures that placed the SGC in the global space of theoretical genetic codes, we showed that the SGC is definitely closer to the codes that minimize the costs of amino acids replacements than those maximizing them.

Conclusions

The standard genetic code represents most likely only partially optimized systems, which emerged under the influence of many different factors. Our findings can be useful to researchers involved in modifying the genetic code of the living organisms and designing artificial ones.

     

tRNA genes and the genetic code

The genetic code describes translational assignments between codons and amino acids. tRNAs and aminoacyl-tRNA synthetases (aaRSs) are those molecules by means of which these assignments are established. Any aaRS recognizes its tRNAs according to some of their nucleotides called identity elements (IEs). Let a 1Mut-similarity Sim (1Mut) be the average similarity between such tRNA genes whose codons differ by one point mutation. We showed that: (1) a global maximum of Sim (1Mut) is reached at the standard genetic code 27 times for 4 sets of IEs of tRNA genes of eukaryotic species, while it is so only 5 times for similarities Sim (C&R) between all tRNA genes whose codons lie in the same column or row of the code. Therefore, point mutations of anticodons were tested by nature to recruit tRNAs from one isoaccepting group to another, (2) because plain similarities Sim (all) between tRNA genes of species within any of the three domains of life are higher than between tRNA genes of species belonging to different domains, tRNA genes retained information about early evolution of cells, (3) we searched the order of tRNAs in which they were most probably assigned to their codons and amino acids. The beginning Ala, (Val), Pro, Ile, Lys, Arg, Trp, Met, Asp, Cys, (Ser) of our resulting chronology lies under a plateau on a graph of Sim (1Mut,IE)(univ.ancestors) plotted over this chronology for a set S(IE) of all IEs of tRNA genes, whose universal ancestors were separately computed for each codon. This plateau has remained preserved along the whole line of evolution of the code and is consistent with observations of Ribas de Pouplana and Schimmel [2001. Aminoacy1-tRNA synthetases: potential markers of genetic code development. Trends Biochem. Sci. 26, 591-598] that specific pairs of aaRSs-one from each of their two classes-can be docked simultaneously onto the acceptor stem of tRNA and hence an interaction existed between their ancestors using a reduced code, (4) sharpness of a local maximum of Sim (1Mut) at the standard code is almost 100% along our chronologies.     

Attenuation of the editing activity of the Escherichia coli leucyl-tRNA synthetase allows incorporation of novel amino acids into proteins in vivo

The fidelity of translation is dependent on the specificity of the aminoacyl-tRNA synthetases (aaRSs). The aaRSs that activate the hydrophobic amino acids leucine, isoleucine, and valine employ a proofreading mechanism that hydrolyzes noncognate aminoacyl adenylates and misaminoacylated tRNAs. Discrimination between structurally similar amino acids by these AARSs is believed to operate by a double-sieve principle, wherein a separate editing domain governs hydrolysis on the basis of the size and hydrophilicity of the amino acid side chain. Leucyl-tRNA synthetase (LeuRS) relies on its editing function to correct misaminoacylation of tRNA(Leu) by isoleucine and methionine. Thr252 of Escherichia coli LeuRS has been shown previously to be important in defining the size of the editing cavity. Here we report the isolation and characterization of three LeuRS mutants with point mutations at this position (T252Y, T252L, and T252F). The proofreading activity of the synthetase is significantly impaired when an amino acid bulkier than threonine is introduced. The rate of misaminoacylation of tRNA(Leu) by isoleucine and valine increases with the increasing size of the amino acid substituent at position 252, and the noncognate amino acids norvaline and norleucine are inserted efficiently at the leucine sites of recombinant proteins under conditions of constitutive overexpression of the T252Y mutant in E. coli. In addition, the unsaturated amino acids allylglycine, homoallylglycine, homopropargylglycine, and 2-butynylalanine all support protein synthesis in E. coli hosts harboring the mutant synthetase. These results demonstrate that programmed manipulation of the editing cavity can allow in vivo incorporation of novel protein building blocks.