Evolution and relatedness in two aminoacyl-tRNA synthetase families

Cross-Omics Analyses Reveal the Effects of Ambient PM2.5 Exposure on Hepatic Metabolism in Female Mice

Ambient particulate matter (PM2.5) is a potential risk factor for metabolic damage to the liver. Epidemiological studies suggest that elevated PM2.5 concentrations cause changes in hepatic metabolism, but there is a lack of laboratory evidence. Here, we aimed to evaluate the effects of PM2.5 exposure on liver metabolism in C57BL/6j female mice (10 months old) and to explore the mechanisms underlying metabolic alterations and differential gene expressions by combining metabolomics and transcriptomics analyses. The metabolomics results showed that PM2.5 exposure notably affected the metabolism of amino acids and organic acids and caused hepatic lipid and bile acid accumulation. The transcriptomic analyses revealed that PM2.5 exposure led to a series of metabolic pathway abnormalities, including steroid biosynthesis, steroid hormone biosynthesis, primary bile acid biosynthesis, etc. Among them, the changes in the bile acid pathway might be one of the causes of liver damage in mice. In conclusion, this study clarified the changes in liver metabolism in mice caused by PM2.5 exposure through combined transcriptomic and metabolomic analyses, revealed that abnormal bile acid metabolism is the key regulatory mechanism leading to metabolic-associated fatty liver disease (MAFLD) in mice, and provided laboratory evidence for further clarifying the effects of PM2.5 on body metabolism.

Computational Insight into Substrate-Induced Conformational Changes in Methionyl-tRNA Synthetase of Mycobacterium Tuberculosis

Tuberculosis caused by Mycobacterium tuberculosis (M.tb) has killed millions worldwide. Antibiotic resistance leads to the ineffectiveness of the current therapies. Aminoacyl tRNA synthetase (aaRS) class of proteins involved in protein synthesis are promising bacterial targets for developing new therapies. Here, we carried out a systematic comparative study on the aaRS sequences from M.tb and human. We listed important M.tb aaRS that could be explored as potential M.tb targets alongside the detailed conformational space analysis of methionyl-tRNA synthetase (MetRS) in apo- and substrate-bound form, which is among the proposed targets. Understanding the conformational dynamics is central to the mechanistic understanding of MetRS, as the substrate binding leads to the conformational changes causing the reaction to proceed. We performed the most complete simulation study of M.tb MetRS for 6 microseconds (2 systems × 3 runs × 1 microsecond) in the apo and substrate-bound states. Interestingly, we observed differential features, showing comparatively large dynamics for the holo simulations, whereas the apo structures became slightly compact with reduced solvent exposed area. In contrast, the ligand size decreased significantly in holo structures possibly to relax ligand conformation. Our findings correlate with experimental studies, thus validating our protocol. Adenosine monophosphate moiety of the substrate exhibited quite higher fluctuations than the methionine. His21 and Lys54 were found to be the important residues forming prominent hydrogen bond and salt-bridge interactions with the ligand. The ligand-protein affinity decreased during simulations as computed by MMGBSA analysis over the last 500 ns trajectories, which indicates the conformational changes upon ligand binding. These differential features could be further explored for designing new M.tb inhibitors.

Antibiotic hyper-resistance in a class I aminoacyl-tRNA synthetase with altered active site signature motif

Antibiotics target key biological processes that include protein synthesis. Bacteria respond by developing resistance, which increases rapidly due to antibiotics overuse. Mupirocin, a clinically used natural antibiotic, inhibits isoleucyl-tRNA synthetase (IleRS), an enzyme that links isoleucine to its tRNAIle for protein synthesis. Two IleRSs, mupirocin-sensitive IleRS1 and resistant IleRS2, coexist in bacteria. The latter may also be found in resistant Staphylococcus aureus clinical isolates. Here, we describe the structural basis of mupirocin resistance and unravel a mechanism of hyper-resistance evolved by some IleRS2 proteins. We surprisingly find that an up to 103-fold increase in resistance originates from alteration of the HIGH motif, a signature motif of the class I aminoacyl-tRNA synthetases to which IleRSs belong. The structural analysis demonstrates how an altered HIGH motif could be adopted in IleRS2 but not IleRS1, providing insight into an elegant mechanism for coevolution of the key catalytic motif and associated antibiotic resistance.

A Leucyl-tRNA Synthetase Inhibitor with Broad-Spectrum Anti-Mycobacterial Activity

Global infections by non-tuberculous mycobacteria (NTM) are steadily rising. New drugs are needed to treat NTM infections, but the NTM drug pipeline remains poorly populated and focused on repurposing or reformulating approved antibiotics. We sought to accelerate de novo NTM drug discovery by testing advanced compounds with established activity against Mycobacterium tuberculosis 3-aminomethyl 4-halogen benzoxaboroles, a novel class of leucyl-tRNA synthetase inhibitors, were recently discovered as active against M. tuberculosis Here, we report that the benzoxaborole EC/11770 is not only a potent anti-tubercular agent but is active against the M. abscessus and M. avium complexes. Focusing on M. abscessus, which causes the most difficult-to-cure NTM disease, we show that EC/11770 retained potency against drug-tolerant biofilms in vitro and was effective in a mouse lung infection model. Resistant mutant selection experiments showed a low frequency of resistance and confirmed leucyl-tRNA synthetase as the target. This work establishes the benzoxaborole EC/11770 as a novel preclinical candidate for the treatment of NTM lung disease and tuberculosis and validates leucyl-tRNA synthetase as an attractive target for the development of broad-spectrum anti-mycobacterials.

Amino Acid Specificity of Ancestral Aminoacyl-tRNA Synthetase Prior to the Last Universal Common Ancestor Commonote commonote

Extant organisms commonly use 20 amino acids in protein synthesis. In the translation system, aminoacyl-tRNA synthetase (ARS) selectively binds an amino acid and transfers it to the cognate tRNA. It is postulated that the amino acid repertoire of ARS expanded during the development of the translation system. In this study we generated composite phylogenetic trees for seven ARSs (SerRS, ProRS, ThrRS, GlyRS-1, HisRS, AspRS, and LysRS) which are thought to have diverged by gene duplication followed by mutation, before the evolution of the last universal common ancestor. The composite phylogenetic tree shows that the AspRS/LysRS branch diverged from the other five ARSs at the deepest node, with the GlyRS/HisRS branch and the other three ARSs (ThrRS, ProRS and SerRS) diverging at the second deepest node. ThrRS diverged next, and finally ProRS and SerRS diverged from each other. Based on the phylogenetic tree, sequences of the ancestral ARSs prior to the evolution of the last universal common ancestor were predicted. The amino acid specificity of each ancestral ARS was then postulated by comparison with amino acid recognition sites of ARSs of extant organisms. Our predictions demonstrate that ancestral ARSs had substantial specificity and that the number of amino acid types amino-acylated by proteinaceous ARSs was limited before the appearance of a fuller range of proteinaceous ARS species. From an assumption that 10 amino acid species are required for folding and function, proteinaceous ARS possibly evolved in a translation system composed of preexisting ribozyme ARSs, before the evolution of the last universal common ancestor.

Structural Computational Analysis of the Natural History of Class I aminoacyl-tRNA Synthetases Suggests their Role in Establishing the Genetic Code

The evolutionary history of Class I aminoacyl-tRNA synthetases (aaRS) through the reconstruction of ancestral sequences is presented. From structural molecular modeling, we sought to understand its relationship with the acceptor arms and the tRNA anticodon loop, how this relationship was established, and the possible implications in determining the genetic code and the translation system. The results of the molecular docking showed that in 7 out 9 aaRS, the acceptor arm and the anticodon loop bond practically in the same region. Domain accretion process in aaRS and repositioning of interactions between tRNAs and aaRS are illustrated. Based on these results, we propose that the operational code and the anticodon code coexisted, competing for the aaRS catalytic region, while consequently contributed to the stabilization of these proteins.

Inhibitors of aminoacyl-tRNA synthetases as antimycobacterial compounds: An up-to-date review

Aminoacyl-tRNA synthetases (aaRSs) are crucial for the correct assembly of amino acids to cognate tRNA to maintain the fidelity of proteosynthesis. AaRSs have become a hot target in antimicrobial research. Three aaRS inhibitors are already in clinical practice; antibacterial mupirocin inhibits the synthetic site of isoleucyl-tRNA synthetase, antifungal tavaborole inhibits the editing site of leucyl-tRNA synthetase, and antiprotozoal halofuginone inhibits proline-tRNA synthetase. According to the World Health Organization, tuberculosis globally remains the leading cause of death from a single infectious agent. The rising incidence of multidrug-resistant tuberculosis is alarming and urges the search for new antimycobacterial compounds, preferably with yet unexploited mechanism of action. In this literature review, we have covered the up-to-date state in the field of inhibitors of mycobacterial aaRSs. The most studied aaRS in mycobacteria is LeuRS with at least four structural types of inhibitors, followed by TyrRS and AspRS. Inhibitors of MetRS, LysRS, and PheRS were addressed in a single significant study each. In many cases, the enzyme inhibition activity translated into micromolar or submicromolar inhibition of growth of mycobacteria. The most promising aaRS inhibitor as an antimycobacterial compound is GSK656 (compound 8), the only aaRS inhibitor in clinical trials (Phase IIa) for systemic use against tuberculosis. GSK656 is orally available and shares the oxaborole tRNA-trapping mechanism of action with antifungal tavaborole.

The structural basis of the genetic code: amino acid recognition by aminoacyl-tRNA synthetases

Storage and directed transfer of information is the key requirement for the development of life. Yet any information stored on our genes is useless without its correct interpretation. The genetic code defines the rule set to decode this information. Aminoacyl-tRNA synthetases are at the heart of this process. We extensively characterize how these enzymes distinguish all natural amino acids based on the computational analysis of crystallographic structure data. The results of this meta-analysis show that the correct read-out of genetic information is a delicate interplay between the composition of the binding site, non-covalent interactions, error correction mechanisms, and steric effects.

Protein Sequences Recapitulate Genetic Code Evolution

Several hypotheses predict ranks of amino acid assignments to genetic code's codons. Analyses here show that average positions of amino acid species in proteins correspond to assignment ranks, in particular as predicted by Juke's neutral mutation hypothesis for codon assignments. In all tested protein groups, including co- and post-translationally folding proteins, 'recent' amino acids are on average closer to gene 5' extremities than 'ancient' ones. Analyses of pairwise residue contact energies matrices suggest that early amino acids stereochemically selected late ones that stablilize residue interactions within protein cores, presumably producing 5'-late-to-3'-early amino acid protein sequence gradients. The gradient might reduce protein misfolding, also after mutations, extending principles of neutral mutations to protein folding. Presumably, in self-perpetuating and self-correcting systems like the genetic code, initial conditions produce similarities between evolution of the process (the genetic code) and 'ontogeny' of resulting structures (here proteins), producing apparent teleonomy between process and product.

Genetic validation of aminoacyl-tRNA synthetases as drug targets in Trypanosoma brucei

Human African trypanosomiasis (HAT) is an important public health threat in sub-Saharan Africa. Current drugs are unsatisfactory, and new drugs are being sought. Few validated enzyme targets are available to support drug discovery efforts, so our goal was to obtain essentiality data on genes with proven utility as drug targets. Aminoacyl-tRNA synthetases (aaRSs) are known drug targets for bacterial and fungal pathogens and are required for protein synthesis. Here we survey the essentiality of eight Trypanosoma brucei aaRSs by RNA interference (RNAi) gene expression knockdown, covering an enzyme from each major aaRS class: valyl-tRNA synthetase (ValRS) (class Ia), tryptophanyl-tRNA synthetase (TrpRS-1) (class Ib), arginyl-tRNA synthetase (ArgRS) (class Ic), glutamyl-tRNA synthetase (GluRS) (class 1c), threonyl-tRNA synthetase (ThrRS) (class IIa), asparaginyl-tRNA synthetase (AsnRS) (class IIb), and phenylalanyl-tRNA synthetase (α and β) (PheRS) (class IIc). Knockdown of mRNA encoding these enzymes in T. brucei mammalian stage parasites showed that all were essential for parasite growth and survival in vitro. The reduced expression resulted in growth, morphological, cell cycle, and DNA content abnormalities. ThrRS was characterized in greater detail, showing that the purified recombinant enzyme displayed ThrRS activity and that the protein localized to both the cytosol and mitochondrion. Borrelidin, a known inhibitor of ThrRS, was an inhibitor of T. brucei ThrRS and showed antitrypanosomal activity. The data show that aaRSs are essential for T. brucei survival and are likely to be excellent targets for drug discovery efforts.

Evidence for late resolution of the aux codon box in evolution

Recognition strategies for tRNA aminoacylation are ancient and highly conserved, having been selected very early in the evolution of the genetic code. In most cases, the trinucleotide anticodons of tRNA are important identity determinants for aminoacylation by cognate aminoacyl-tRNA synthetases. However, a degree of ambiguity exists in the recognition of certain tRNA(Ile) isoacceptors that are initially transcribed with the methionine-specifying CAU anticodon. In most organisms, the C34 wobble position in these tRNA(Ile) precursors is rapidly modified to lysidine to prevent recognition by methionyl-tRNA synthetase (MRS) and production of a chimeric Met-tRNA(Ile) that would compromise translational fidelity. In certain bacteria, however, lysidine modification is not required for MRS rejection, indicating that this recognition strategy is not universally conserved and may be relatively recent. To explore the actual distribution of lysidine-dependent tRNA(Ile) rejection by MRS, we have investigated the ability of bacterial MRSs from different clades to differentiate cognate tRNACAU(Met) from near-cognate tRNACAU(Ile). Discrimination abilities vary greatly and appear unrelated to phylogenetic or structural features of the enzymes or sequence determinants of the tRNA. Our data indicate that tRNA(Ile) identity elements were established late and independently in different bacterial groups. We propose that the observed variation in MRS discrimination ability reflects differences in the evolution of genetic code machineries of emerging bacterial clades.

The crystal structures of the α-subunit of the α(2)β (2) tetrameric Glycyl-tRNA synthetase

Aminoacyl-tRNA synthetases (AARSs) are ligases (EC.6.1.1.-) that catalyze the acylation of amino acids to their cognate tRNAs in the process of translating genetic information from mRNA to protein. Their amino acid and tRNA specificity are crucial for correctly translating the genetic code. Glycine is the smallest amino acid and the glycyl-tRNA synthetase (GlyRS) belongs to Class II AARSs. The enzyme is unusual because it can assume different quaternary structures. In eukaryotes, archaebacteria and some bacteria, it forms an α(2) homodimer. In some bacteria, GlyRS is an α(2)β(2) heterotetramer and shows a distant similarity to α(2) GlyRSs. The human pathogen eubacterium Campylobacter jejuni GlyRS (CjGlyRS) is an α(2)β(2) heterotetramer and is similar to Escherichia coli GlyRS; both are members of Class IIc AARSs. The two-step aminoacylation reaction of tetrameric GlyRSs requires the involvement of both α- and β-subunits. At present, the structure of the GlyRS α(2)β(2) class and the details of the enzymatic mechanism of this enzyme remain unknown. Here we report the crystal structures of the catalytic α-subunit of CjGlyRS and its complexes with ATP, and ATP and glycine. These structures provide detailed information on substrate binding and show evidence for a proposed mechanism for amino acid activation and the formation of the glycyl-adenylate intermediate for Class II AARSs.

Inferring the ancient history of the translation machinery and genetic code via recapitulation of ribosomal subunit assembly orders

Universally conserved positions in ribosomal proteins have significant biases in amino acid usage, likely indicating the expansion of the genetic code at the time leading up to the most recent common ancestor(s) (MRCA). Here, we apply this principle to the evolutionary history of the ribosome before the MRCA. It has been proposed that the experimentally determined order of assembly for ribosomal subunits recapitulates their evolutionary chronology. Given this model, we produce a probabilistic evolutionary ordering of the universally conserved small subunit (SSU) and large subunit (LSU) ribosomal proteins. Optimizing the relative ordering of SSU and LSU evolutionary chronologies with respect to minimizing differences in amino acid usage bias, we find strong compositional evidence for a more ancient origin for early LSU proteins. Furthermore, we find that this ordering produces several trends in specific amino acid usages compatible with models of genetic code evolution.

Nucleotide triplet based molecular phylogeny of class I and class II aminoacyl t-RNA synthetase in three domain of life process: bacteria, archaea, and eukarya

The aminoacyl-tRNA synthetases are one of the major protein components in the translation machinery. These essential proteins are found in all forms of life and are responsible for charging their cognate tRNAs with the correct amino acid. These important enzymes have been the subject of intense scientific inquiry for nearly half a century, but their complete evolutionary history has yet to emerge. Amino acids sequence based phylogeny has some limitation due to very low sequence similarity amongst the different tRNA synthetases and structure based phylogeny has also its limitation. In our study, tRNA nucleotide sequences of E. coli K12 (Bacteria), Saccharomyces cerevisiae (Eukarya), Thermococcus kodakaraensis KOD1, and Archaeoglobus fulgidus DSM 4304 (Archaea) were used for phylogenetic analysis. Our results complement the observation with the earlier studies based on multiple sequence alignment and structural alignment. We observed that relationship between archaeal tRNA synthetases are different that of bacteria and eucarya. Violation of Class rule of LysRS is observed here also. The uniqueness of this method is that it does not employ sequence alignment of complete nucleotide sequence of the corresponding gene.

Aminoacyl-tRNA synthetase classes and groups in prokaryotes

Knowledge on the evolution of aminoacyl-tRNA synthetases is crucial to studies on the origins of life. The relationships between the different aminoacyl-tRNA synthetase specificities in prokaryotic organisms are studied in this work. We reconstructed the ancestor sequences and the phylogenetic relationships utilizing the Maximum Likelihood method. The results suggest that in class I the evolution of the N-terminal segment was strongly influenced by the amino acid hydropathy in both domains of prokaryotes. The results for the C-terminal segments of class I were different in the two domains, indicating that its evolution was strongly influenced by the specific types of tRNA modification in each domain. The class II groups in Archaea were more heterogeneous with respect to the hydropathy of amino acids, indicating the interference of other influences. In bacteria, the configuration was also complex but the overall consensual division in two groups was maintained, group IIa forming a single branch with the five hydroapathetic amino acid specificities and group IIb containing the specificities for the moderately hydrophobic together with the hydrophilic amino acids. It is indicated that the aminoacyl-tRNA synthetase in both domains were subjected to different selective forces in diverse parts of the proteins, resulting in complex phylogenetic patterns.

In silico detection of tRNA sequence features characteristic to aminoacyl-tRNA synthetase class membership

Aminoacyl tRNA synthetases (aaRS) are grouped into Class I and II based on primary and tertiary structure and enzyme properties suggesting two independent phylogenetic lineages. Analogously, tRNA molecules can also form two respective classes, based on the class membership of their corresponding aaRS. Although some aaRS–tRNA interactions are not extremely specific and require editing mechanisms to avoid misaminoacylation, most aaRS–tRNA interactions are rather stereospecific. Thus, class-specific aaRS features could be mirrored by class-specific tRNA features. However, previous investigations failed to detect conserved class-specific nucleotides. Here we introduce a discrete mathematical approach that evaluates not only class-specific ‘strictly present’, but also ‘strictly absent’ nucleotides. The disjoint subsets of these elements compose a unique partition, named extended consensus partition (ECP). By analyzing the ECP for both Class I and II tDNA sets from 50 (13 archaeal, 30 bacterial and 7 eukaryotic) species, we could demonstrate that class-specific tRNA sequence features do exist, although not in terms of strictly conserved nucleotides as it had previously been anticipated. This finding demonstrates that important information was hidden in tRNA sequences inaccessible for traditional statistical methods. The ECP analysis might contribute to the understanding of tRNA evolution and could enrich the sequence analysis tool repertoire.

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.

Co-repressor induced order and biotin repressor dimerization: a case for divergent followed by convergent evolution

BirA catalyzes the adenylation and subsequent covalent attachment of biotin to the biotin carboxyl carrier protein (BCCP). In the absence of apo-BCCP, biotin-5'-AMP acts as a co-repressor that induces BirA dimerization and binding to the bio operator to repress biotin biosynthesis. The crystal structures of apo-BirA, and BirA in complex with biotin have been reported. We here describe the 2.8A resolution crystal structure of BirA in complex with the co-repressor analog biotinol-5'-AMP. It was previously shown that the structure of apo-BirA is monomeric and that binding of biotin weakly induces a dimeric structure in which three disordered surface loops become organized to form the dimer interface. The structure of the co-repressor complex is also a dimer, clearly related to the BirA.biotin structure, but with several significant conformational changes. A hitherto disordered "adenylate binding loop" forms a well-defined structure covering the co-repressor. The co-repressor buttresses the dimer interface, resulting in improved packing and a 12 degrees change in the hinge-bending angle along the dimer interface relative to the BirA.biotin structure. This helps explain why the binding of the co-repressor is necessary to optimize the binding of BirA to the bioO operator. The structure reveals an unexpected use of the nucleotide-binding motif GXGXXG in binding adenylate and controlling the repressor function. Finally, based on structural analysis we propose that the class of adenylating enzymes represented by BirA, lipoate protein ligase and class II tRNA synthetases diverged early and were selected based on their ability to sequester co-factors or amino acid residues, and adenylation activity arose independently through functional convergence.

Evolution of different oligomeric glycyl-tRNA synthetases

There are two oligomeric types of glycyl-tRNA synthetases (GlyRSs) in genome, the alpha2beta2 tetramer and alpha2 dimer. Here, we showed that the anticodon-binding domains (ABDs) of dimeric and tetrameric GlyRSs are non-homologous, although their catalytic central domains (CCDs) are homologous. The dimeric GlyRS_ABD is fused to the C-terminal of CCD in alpha-subunit, but the tetrameric GlyRS_ABD is to the C-terminal in beta-subunit during evolution. Generally, one species only contains one oligomeric type of GlyRS, but the both oligomeric GlyRSs with the multiple homologous domains can be observed in Magnetospirillum magnetotacticum genome, nevertheless, these homologous domains are probably from different genomes.

Amino acid biogenesis, evolution of the genetic code and aminoacyl-tRNA synthetases

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.

On the evolution of structure in aminoacyl-tRNA synthetases

The aminoacyl-tRNA synthetases are one of the major protein components in the translation machinery. These essential proteins are found in all forms of life and are responsible for charging their cognate tRNAs with the correct amino acid. The evolution of the tRNA synthetases is of fundamental importance with respect to the nature of the biological cell and the transition from an RNA world to the modern world dominated by protein-enzymes. We present a structure-based phylogeny of the aminoacyl-tRNA synthetases. By using structural alignments of all of the aminoacyl-tRNA synthetases of known structure in combination with a new measure of structural homology, we have reconstructed the evolutionary history of these proteins. In order to derive unbiased statistics from the structural alignments, we introduce a multidimensional QR factorization which produces a nonredundant set of structures. Since protein structure is more highly conserved than protein sequence, this study has allowed us to glimpse the evolution of protein structure that predates the root of the universal phylogenetic tree. The extensive sequence-based phylogenetic analysis of the tRNA synthetases (Woese et al., Microbiol. Mol. Biol. Rev. 64:202-236, 2000) has further enabled us to reconstruct the complete evolutionary profile of these proteins and to make connections between major evolutionary events and the resulting changes in protein shape. We also discuss the effect of functional specificity on protein shape over the complex evolutionary course of the tRNA synthetases.

Crystal structures that suggest late development of genetic code components for differentiating aromatic side chains

Early forms of the genetic code likely generated "statistical" proteins, with similar side chains occupying the same sequence positions at different ratios. In this scenario, groups of related side chains were treated by aminoacyl-tRNA synthetases as a single molecular species until a discrimination mechanism developed that could separate them. The aromatic amino acids tryptophan, tyrosine, and phenylalanine likely constituted one of these groups. A crystal structure of human tryptophanyl-tRNA synthetase was solved at 2.1 A with a tryptophanyl-adenylate bound at the active site. A cocrystal structure of an active fragment of human tyrosyl-tRNA synthetase with its cognate amino acid analog was also solved at 1.6 A. The two structures enabled active site identifications and provided the information for structure-based sequence alignments of approximately 45 orthologs of each enzyme. Two critical positions shared by all tyrosyl-tRNA synthetases and tryptophanyl-tRNA synthetases for amino acid discrimination were identified. The variations at these two positions and phylogenetic analyses based on the structural information suggest that, in contrast to many other amino acids, discrimination of tyrosine from tryptophan occurred late in the development of the genetic code.

On the evolution of structure in aminoacyl-tRNA synthetases

The aminoacyl-tRNA synthetases are one of the major protein components in the translation machinery. These essential proteins are found in all forms of life and are responsible for charging their cognate tRNAs with the correct amino acid. The evolution of the tRNA synthetases is of fundamental importance with respect to the nature of the biological cell and the transition from an RNA world to the modern world dominated by protein-enzymes. We present a structure-based phylogeny of the aminoacyl-tRNA synthetases. By using structural alignments of all of the aminoacyl-tRNA synthetases of known structure in combination with a new measure of structural homology, we have reconstructed the evolutionary history of these proteins. In order to derive unbiased statistics from the structural alignments, we introduce a multidimensional QR factorization which produces a nonredundant set of structures. Since protein structure is more highly conserved than protein sequence, this study has allowed us to glimpse the evolution of protein structure that predates the root of the universal phylogenetic tree. The extensive sequence-based phylogenetic analysis of the tRNA synthetases (Woese et al., Microbiol. Mol. Biol. Rev. 64:202-236, 2000) has further enabled us to reconstruct the complete evolutionary profile of these proteins and to make connections between major evolutionary events and the resulting changes in protein shape. We also discuss the effect of functional specificity on protein shape over the complex evolutionary course of the tRNA synthetases.

Leucyl-tRNA synthetase consisting of two subunits from hyperthermophilic bacteria Aquifex aeolicus

In a hyperthermophilic bacterium, Aquifex aeolicus, leucyl-tRNA synthetase (LeuRS) consists of two non-identical polypeptide subunits (alpha and beta), different from the canonical LeuRS, which has a single polypeptide chain. By PCR, using genome DNA of A. aeolicus as a template, genes encoding the alpha and beta subunits were amplified and cloned in Escherichia coli. The alpha subunit could not be expressed stably in vivo, whereas the beta subunit was overproduced and purified by a simple procedure. The beta subunit was inactive in catalysis but was able to bind tRNA(Leu). Interestingly, the heterodimer alphabeta-LeuRS could be overproduced in E. coli cells containing both genes and was purified to 95% homogeneity as a hybrid dimer. The kinetics of A. aeolicus LeuRS in pre-steady and steady states and cross-recognition of LeuRS and tRNA(Leu) from A. aeolicus and E. coli were studied. Magnesium concentration, pH value, and temperature aminoacylation optima were determined to be 12 mm, 7.8, and 70 degrees C, respectively. Under optimal conditions, A. aeolicus alphabeta-LeuRS is stable up to 65 degrees C.

Monophyly of class I aminoacyl tRNA synthetase, USPA, ETFP, photolyase, and PP-ATPase nucleotide-binding domains: implications for protein evolution in the RNA

Protein sequence and structure comparisons show that the catalytic domains of Class I aminoacyl-tRNA synthetases, a related family of nucleotidyltransferases involved primarily in coenzyme biosynthesis, nucleotide-binding domains related to the UspA protein (USPA domains), photolyases, electron transport flavoproteins, and PP-loop-containing ATPases together comprise a distinct class of alpha/beta domains designated the HUP domain after HIGH-signature proteins, UspA, and PP-ATPase. Several lines of evidence are presented to support the monophyly of the HUP domains, to the exclusion of other three-layered alpha/beta folds with the generic "Rossmann-like" topology. Cladistic analysis, with patterns of structural and sequence similarity used as discrete characters, identified three major evolutionary lineages within the HUP domain class: the PP-ATPases; the HIGH superfamily, which includes class I aaRS and related nucleotidyltransferases containing the HIGH signature in their nucleotide-binding loop; and a previously unrecognized USPA-like group, which includes USPA domains, electron transport flavoproteins, and photolyases. Examination of the patterns of phyletic distribution of distinct families within these three major lineages suggests that the Last Universal Common Ancestor of all modern life forms encoded 15-18 distinct alpha/beta ATPases and nucleotide-binding proteins of the HUP class. This points to an extensive radiation of HUP domains before the last universal common ancestor (LUCA), during which the multiple class I aminoacyl-tRNA synthetases emerged only at a late stage. Thus, substantial evolutionary diversification of protein domains occurred well before the modern version of the protein-dependent translation machinery was established, i.e., still in the RNA world.

On the relative content of G,C bases in codons of amino acids corresponding to class I and II aminoacyl-tRNA synthetases

We have analyzed the relative G,C content from protein coding regions of 530 organisms and found that the ratio of the G,C content of the codons of the amino acids corresponding to Class II and Class I aminoacyl-tRNA synthetases decreases in a statistically significant way from prokaryotes to animals. This can be interpreted assuming that an initial asymmetry between the G,C content of codons of Class I and II amino acids existed and has decreased in the course of evolution.

Seryl-tRNA synthetase is not responsible for the evolution of CUG codon reassignment in Candida albicans

A number of Candida species translate the standard leucine-CUG codon as serine using a novel ser-tRNA(CAG). This tRNA, which has an unusual anticodon stem-loop structure, has been implicated in the evolution of this codon reassignment. However, such a sense codon reassignment might also require a change in the specificity of the cognate aminoacyl tRNA-synthetase, in this case the ser-tRNA synthetase. Here we describe the cloning and sequence analysis of the C. albicans seryl aminoacyl-tRNA synthetase (CaSerRS) gene (CaSES1). The predicted CaSerRS sequence shows a significant level of amino acid identity to SerRs from other organisms and fully complements a S. cerevisiae SerRS null strain without any apparent defect in growth rate. This suggests that the SerRS recognizes and charges S. cerevisiae ser-tRNAs with similar efficiency to that of the S. cerevisiae SerRS. Using an antibody raised against CaSerRS, we also demonstrate the presence of SerRS in a range of Candida spp. showing CUG codon reassignment. We conclude that the key element in CUG reassigment in Candida spp. is the tRNA that decodes the CUG codon rather than a SerRS structural change. The nucleotide sequence of the CaSES1 gene has been deposited at GenBank under Accession No. AF290915.

Equilibrium unfolding of Bombyx mori glycyl-tRNA synthetase

Unfolding of Bombyx mori glycyl-tRNA synthetase was examined by multiple spectroscopic techniques. Tryptophan fluorescence of wild type enzyme and an N-terminally truncated form (N55) increased at low concentrations of urea or guanidine-HCl followed by a reduction in intensity at intermediate denaturant concentrations; a transition at higher denaturant was detected as decreased fluorescence intensity and a red-shifted emission. Solute quenching of fluorescence indicated that tryptophans become progressively solvent-exposed during unfolding. Wild type enzyme had stronger negative CD bands between 220 and 230 nm than the mutant, indicative of greater alpha-helical content. Urea or guanidine-HCl caused a reduction in ellipticity at 222 nm at low denaturant concentration with the wild type enzyme, a transition that is absent in the mutant; both enzymes exhibited a cooperative transition at higher denaturant concentrations. Both enzymes dissociate to monomers in 1.5 m urea. Unfolding of wild type enzyme is described by a multistate unfolding and a parallel two state unfolding; the two-state component is absent in the mutant. Changes in spectral properties associated with unfolding were largely reversible after dilution to low denaturant. Unfolding of glycyl-tRNA synthetase is complex with a native state, a native-like monomer, partially unfolded states, and the unfolded state.

The phylogeny of proteobacteria: relationships to other eubacterial phyla and eukaryotes

The evolutionary relationships of proteobacteria, which comprise the largest and phenotypically most diverse division among prokaryotes, are examined based on the analyses of available molecular sequence data. Sequence alignments of different proteins have led to the identification of numerous conserved inserts and deletions (referred to as signature sequences), which either are unique characteristics of various proteobacterial species or are shared by only members from certain subdivisions of proteobacteria. These signature sequences provide molecular means to define the proteobacterial phyla and their various subdivisions and to understand their evolutionary relationships to the other groups of eubacteria as well as the eukaryotes. Based on signature sequences that are present in different proteins it is now possible to infer that the various eubacterial phyla evolved from a common ancestor in the following order: low-G+C Gram-positive-->high-G+C Gram-positive-->Deinococcus-Thermus (green nonsulfur bacteria)-->cyanobacteria-->Spirochetes-->Chlamydia-Cytophaga-Aquifex -green sulfur bacteria-->Proteobacteria-1 (epsilon and delta)-->Proteobacteria-2 (alpha)-->Proteobacteria-3 (beta)-->Proteobacteria-4 (gamma). An unexpected but important aspect of the relationship deduced here is that the main eubacterial phyla are related to each other linearly rather than in a tree-like manner, suggesting that the major evolutionary changes within Bacteria have taken place in a directional manner. The identified signatures permit placement of prokaryotes into different groups/divisions and could be used for determinative purposes. These signatures generally support the origin of mitochondria from an alpha-proteobacterium and provide evidence that the nuclear cytosolic homologs of many genes are also derived from proteobacteria.

A single base substitution in the variable pocket of yeast tRNA(Arg) eliminates species-specific aminoacylation

Early biochemical data showed that aminoacyl-tRNA synthetases often displayed species-specific recognition of tRNA. We compared the ability of purified Saccharomyces cerevisiae and Escherichia coli arginyl-tRNA synthetases to aminoacylate native and transcribed yeast tRNA(Arg) as well as E. coli tRNA(Arg). The kinetic data revealed that yeast ArgRS could charge E. coli tRNA(Arg), but at a lower efficiency than it charged either the transcribed or native yeast tRNA(Arg). E. coli ArgRS can acylate only its cognate E. coli tRNA. Strikingly, a single base change from C to A at position 20 in yeast tRNA(3)(Arg) altered the species specificity. The transcript of yeast tRNA(3)(Arg)CA20 mutant was aminoacylated by E. coli ArgRS with a 10(6) increase in k(cat)/K(m) over that for aminoacylation of yeast tRNA(3)(Arg) transcript. This indicates that A20 is not only an important identity of E. coli tRNA(Arg), but is also the key to altering species-specific aminoacylation of yeast tRNA(Arg).

Evolution of the genetic code

Comparative path lengths in amino acid biosynthesis and other molecular indicators of the timing of codon assignment were examined to reconstruct the main stages of code evolution. The codon tree obtained was rooted in the 4 N-fixing amino acids (Asp, Glu, Asn, Gln) and 16 triplets of the NAN set. This small, locally phased (commaless) code evidently arose from ambiguous translation on a poly(A) collector strand, in a surface reaction network. Copolymerisation of these amino acids yields polyanionic peptide chains, which could anchor uncharged amide residues to a positively charged mineral surface. From RNA virus structure and replication in vitro, the first genes seemed to be RNA segments spliced into tRNA. Expansion of the code reduced the risk of mutation to an unreadable codon. This step was conditional on initiation at the 5'-codon of a translated sequence. Incorporation of increasingly hydrophobic amino acids accompanied expansion. As codons of the NUN set were assigned most slowly, they received the most nonpolar amino acids. The origin of ferredoxin and Gln synthetase was traced to mid-expansion phase. Surface metabolism ceased by the end of code expansion, as cells bounded by a proteo-phospholipid membrane, with a protoATPase, had emerged. Incorporation of positively charged and aromatic amino acids followed. They entered the post-expansion code by codon capture. Synthesis of efficient enzymes with acid-base catalysis was then possible. Both types of aminoacyl-tRNA synthetases were attributed to this stage. tRNA sequence diversity and error rates in RNA replication indicate the code evolved within 20 million yr in the preIsuan era. These findings on the genetic code provide empirical evidence, from a contemporaneous source, that a surface reaction network, centred on C-fixing autocatalytic cycles, rapidly led to cellular life on Earth.

tRNA recognition and evolution of determinants in seryl-tRNA synthesis

We have analyzed the evolution of recognition of tRNAsSerby seryl-tRNA synthetases, and compared it to other type 2 tRNAs, which contain a long extra arm. In Eubacteria and chloroplasts this type of tRNA is restricted to three families: tRNALeu, tRNASer and tRNATyr. tRNALeuand tRNASer also carry a long extra arm in Archaea, Eukarya and all organelles with the exception of animal mitochondria. In contrast, the long extra arm of tRNATyr is far less conserved: it was drastically shortened after the separation of Archaea and Eukarya from Eubacteria, and it is also truncated in animal mitochondria. The high degree of phylo-genetic divergence in the length of tRNA variable arms, which are recognized by both class I and class II aminoacyl-tRNA synthetases, makes type 2 tRNA recognition an ideal system with which to study how tRNA discrimination may have evolved in tandem with the evolution of other components of the translation machinery.

Aminoacylation of tRNA in the evolution of an aminoacyl-tRNA synthetase

Aminoacyl-tRNA synthetases catalyze aminoacylation of tRNAs by joining an amino acid to its cognate tRNA. The selection of the cognate tRNA is jointly determined by separate structural domains that examine different regions of the tRNA. The cysteine-tRNA synthetase of Escherichia coli has domains that select for tRNAs containing U73, the GCA anticodon, and a specific tertiary structure at the corner of the tRNA L shape. The E. coli enzyme does not efficiently recognize the yeast or human tRNACys, indicating the evolution of determinants for tRNA aminoacylation from E. coli to yeast to human and the coevolution of synthetase domains that interact with these determinants. By successively modifying the yeast and human tRNACys to ones that are efficiently aminoacylated by the E. coli enzyme, we have identified determinants of the tRNA that are important for aminoacylation but that have diverged in the course of evolution. These determinants provide clues to the divergence of synthetase domains. We propose that the domain for selecting U73 is conserved in evolution. In contrast, we propose that the domain for selecting the corner of the tRNA L shape diverged early, after the separation between E. coli and yeast, while that for selecting the GCA-containing anticodon loop diverged late, after the separation between yeast and human.

Molecular biology of hyperthermophilic Archaea

The sequences of a number of archaeal genomes have recently been completed, and many more are expected shortly. Consequently, the research of Archaea in general and hyperthermophiles in particular has entered a new phase, with many exciting discoveries to be expected. The wealth of sequence information has already led, and will continue to lead to the identification of many enzymes with unique properties, some of which have potential for industrial applications. Subsequent functional genomics will help reveal fundamental matters such as details concerning the genetic, biochemical and physiological adaptation of extremophiles, and hence give insight into their genomic evolution, polypeptide structure-function relations, and metabolic regulation. In order to optimally exploit many unique features that are now emerging, the development of genetic systems for hyperthermophilic Archaea is an absolute requirement. Such systems would allow the application of this class of Archaea as so-called "cell factories": (i) expression of certain archaeal enzymes for which no suitable conventional (mesophilic bacterial or eukaryal) systems are available, (ii) selection for thermostable variants of potentially interesting enzymes from mesophilic origin, and (iii) the development of in vivo production systems by metabolic engineering. An overview is given of recent insight in the molecular biology of hyperthermophilic Archaea, as well as of a number of promising developments that should result in the generation of suitable genetic systems in the near future.

Genetic code in evolution: switching species-specific aminoacylation with a peptide transplant

The genetic code is established in aminoacylation reactions whereby amino acids are joined to tRNAs bearing the anticodons of the genetic code. Paradoxically, while the code is universal there are many examples of species-specific aminoacylations, where a tRNA from one taxonomic domain cannot be acylated by a synthetase from another. Here we consider an example where a human, but not a bacterial, tRNA synthetase charges its cognate eukaryotic tRNA and where the bacterial, but not the human, enzyme charges the cognate bacterial tRNA. While the bacterial enzyme has less than 10% sequence identity with the human enzyme, transplantation of a 39 amino acid peptide from the human into the bacterial enzyme enabled the latter to charge its eukaryotic tRNA counterpart in vitro and in vivo. Conversely, substitution of the corresponding peptide of the bacterial enzyme for that of the human enabled the human enzyme to charge bacterial tRNA. This peptide element discriminates a base pair difference in the respective tRNA acceptor stems. Thus, functionally important co-adaptations of a synthetase to its tRNA act as small modular units that can be moved across taxonomic domains and thereby preserve the universality of the code.

Archaea and the prokaryote-to-eukaryote transition

Since the late 1970s, determining the phylogenetic relationships among the contemporary domains of life, the Archaea (archaebacteria), Bacteria (eubacteria), and Eucarya (eukaryotes), has been central to the study of early cellular evolution. The two salient issues surrounding the universal tree of life are whether all three domains are monophyletic (i.e., all equivalent in taxanomic rank) and where the root of the universal tree lies. Evaluation of the status of the Archaea has become key to answering these questions. This review considers our cumulative knowledge about the Archaea in relationship to the Bacteria and Eucarya. Particular attention is paid to the recent use of molecular phylogenetic approaches to reconstructing the tree of life. In this regard, the phylogenetic analyses of more than 60 proteins are reviewed and presented in the context of their participation in major biochemical pathways. Although many gene trees are incongruent, the majority do suggest a sisterhood between Archaea and Eucarya. Altering this general pattern of gene evolution are two kinds of potential interdomain gene transferrals. One horizontal gene exchange might have involved the gram-positive Bacteria and the Archaea, while the other might have occurred between proteobacteria and eukaryotes and might have been mediated by endosymbiosis.

Molecular divergence of lysozymes and alpha-lactalbumin

The vast number of proteins that sustain the currently living organisms have been generated from a relatively small number of ancestral genes that has involved a variety of processes. Lysozyme is an ancient protein whose origin goes back an estimated 400 to 600 million years. This protein was originally a bacteriolytic defensive agent and has been adapted to serve a digestive function on at least two occasions, separated by nearly 40 million years. The origins of the related goose type and T4 phage lysozyme that are distinct from the more common C type are obscure. They share no discernable amino acid sequence identity and yet they possess common secondary and tertiary structures. Lysozyme C gene also gave rise, after gene duplication 300 to 400 million years ago, to a gene that currently codes for alpha-lactalbumin, a protein expressed only in the lactating mammary gland of all but a few species of mammals. It is required for the synthesis of lactose, the sugar secreted in milk. alpha-Lactalbumin shares only 40% identity in amino acid sequence with lysozyme C, but it has a closer spatial structure and gene organization. Although structurally similar, functionally they are quite distinct. Specific amino acid substitutions in alpha-lactalbumin account for the loss of the enzyme activity of lysozyme and the acquisition of the features necessary for its role in lactose synthesis. Evolutionary implications are as yet unclear but are being unraveled in many laboratories.

Vestiges of early molecular processes leading to the genetic code

We compare predictions from a proposed model for the origin of the genetic code (J. Theor. Biol (1993) 164, 291-305) with existing information on the base content of codons and abundance of amino acid in different organisms. A comparison is also made between the three groups of amino acids suggested by the model and the two classes of aminoacetyl-tRNA synthetases. The observed agreements tend to support the model.

Evidence for the early divergence of tryptophanyl- and tyrosyl-tRNA synthetases

Each amino acid is attached to its cognate tRNA by a distinct aminoacyl-tRNA synthetase (aaRS). The conventional evolutionary view is that the modern complement of synthetases existed prior to the divergence of eubacteria and eukaryotes. Thus comparisons of prokaryotic and eukaryotic aminoacyl-tRNA synthetases of the same type (charging specificity) should show greater sequence similarities than comparisons between synthetases of different types-and this is almost always so. However, a recent study [Ribas de Pouplana L, Furgier M, Quinn CL, Schimmel P (1996) Proc Natl Acad Sci USA 93:166-170] suggested that tryptophanyl- (TrpRS) and tyrosyl-tRNA (TyrRS) synthetases of the Eucarya (eukaryotes) are more similar to each other than either is to counterparts in the Bacteria (eubacteria). Here, we reexamine the evolutionary relationships of TyrRS and TrpRS using a broader range of taxa, including new sequence data from the Archaea (archaebacteria) as well as species of Eucarya and Bacteria. Our results differ from those of Ribas de Pouplana et al.: All phylogenetic methods support the separate monophyly of TrpRS and TyrRS. We attribute this result to the inclusion of the archaeal data which might serve to reduce long branch effects possibly associated with eukaryotic TrpRS and TyrRS sequences. Furthermore, reciprocally rooted phylogenies of TrpRS and TyrRS sequences confirm the closer evolutionary relationship of Archaea to eukaryotes by placing the root of the universal tree in the Bacteria.

Alanyl-tRNA synthetase gene of the extreme acidophilic chemolithoautotrophic Thiobacillus ferrooxidans is highly homologous to alaS genes from all living kingdoms but cannot be transcribed from its promoter in Escherichia coli

The alaS gene of Thiobacillus ferrooxidans has been cloned and sequenced and its expression in Escherichia coli and T. ferrooxidans analysed. The same genomic organization to that in E. coli (recA-recX-alaS) has been found in T. ferrooxidans. The recA and alaS genes cannot be transcribed from their own promoters in E. coli. In addition to the well-known homology at the protein level between AlaS proteins from various organisms, a strong homology was found between all the known alaS genes from bacteria, archaea and eucarya. Two regions, one of which corresponds to the catalytic core, are particularly well-conserved at the nucleotide sequence level, a possible indication of strong constraints during evolution on these parts of the genes.

General structure/function properties of microbial methionyl-tRNA synthetases

Alignment of the sequences of methionyl-tRNA synthetases from various microbial sources shows low levels of identities. However, sequence identities are clustered in a limited number of sites, most of which contain peptide patterns known to support the activity of the Escherichia coli enzyme. In the present study, site-directed mutagenesis was used to probe the role of these conserved residues in the case of the Bacillus stearothermophilus methionyl-tRNA synthetase. The B. stearothermophilus enzyme was chosen in this study because it can be produced as an active truncated monomeric form, similar to the monomeric derivative of E. coli methionyl-tRNA synthetase produced by mild proteolysis. The two core enzyme molecules share only 27% identical residues. The results allowed the identification of the binding sites for ATP, methionine and tRNA, as well as that responsible for the tight binding of the zinc ion to the enzyme. It is concluded that the thermostable synthetase adopts a three-dimensional folding very similar to that of the E. coli one. Therefore, the two methionyl-tRNA synthetase sequences, although significantly different, maintain a common scaffold with the functionally important residues exposed at constant positions. Sequence alignments suggest that the above conclusion can be generalized to the known methionyl-tRNA synthetases from various sources.

Aminoacylation identity switch of turnip yellow mosaic virus RNA from valine to methionine results in an infectious virus

The turnip yellow mosaic virus genomic RNA terminates at its 3' end in a tRNA-like structure that is capable of specific valylation. By directed mutation, the aminoacylation specificity has been switched from valine to methionine, a novel specificity for viral tRNA-like structures. The switch to methionine specificity, assayed in vitro under physiological buffer conditions with wheat germ methionyl-tRNA synthetase, required mutation of the anticodon loop and the acceptor stem pseudoknot. The resultant methionylatable genomes are infectious and stable in plants, but genomes that lack strong methionine acceptance (as previously shown with regard to valine acceptance) replicate poorly. The results indicate that amplification of turnip yellow mosaic virus RNA requires aminoacylation, but that neither the natural (valine) specificity nor interaction specifically with valyl-tRNA synthetase is crucial.

Widespread use of the glu-tRNAGln transamidation pathway among bacteria. A member of the alpha purple bacteria lacks glutaminyl-trna synthetase

The expression of the Rhizobium meliloti glutamyl-tRNA synthetase gene in Escherichia coli under the control of a trc promoter results in a toxic effect upon isopropyl-beta-D-thiogalactopyranoside induction, which is probably caused by a misacylation activity. To further investigate this unexpected result, we looked at the pathway of Gln-tRNAGln formation in R. meliloti. No glutaminyl-tRNA synthetase activity has been found in R. meliloti crude extract, but we detected a specific aminotransferase activity that changes Glu-tRNAGln to Gln-tRNAGln. Our results show that R. meliloti, a member of the alpha-subdivision of the purple bacteria, is the first Gram-negative bacteria reported to use a transamidation pathway for Gln-tRNAGln synthesis. A phylogenetic analysis of the contemporary glutamyl-tRNA synthetase and glutaminyl-tRNA synthetase amino acid sequences reveals that a close evolutionary relationship exists between R. meliloti and yeast mitochondrial glutamyl-tRNA synthetases, which is consistent with an origin of mitochondria in the alpha-subdivision of Gram-negative purple bacteria. A 256-amino acid open reading frame closely related to bacterial glutamyl-tRNA synthetases, which probably originates from a glutamyl-tRNA synthetase gene duplication, was found in the 4-min region of the E. coli chromosome. We suggest that this open reading frame is a relic of an ancient transamidation pathway that occurred in an E. coli ancestor before the horizontal transfer of a eukaryotic glutaminyl-tRNA synthetase (Lamour, V., Quevillon, S., Diriong, S., N'Guyen, V. C., Lipinski, M., and Mirande, M.(1994) Proc. Natl. Acad. Sci. U. S. A. 91, 8670-8674) and that it favored its stable acquisition. From these observations, a revisited model for the evolution of the contemporary glutamyl-tRNA synthetases and glutaminyl-tRNA synthetases that differs from the generally accepted model for the evolution of aminoacyl-tRNA synthetases is proposed.