Structure of an archaeal non-discriminating glutamyl-tRNA synthetase: a missing link in the evolution of Gln-tRNAGln formation

Evolution and variation in amide aminoacyl-tRNA synthesis

The amide proteogenic amino acids, asparagine and glutamine, are two of the twenty amino acids used in translation by all known life. The aminoacyl-tRNA synthetases for asparagine and glutamine, asparaginyl-tRNA synthetase and glutaminyl tRNA synthetase, evolved after the split in the last universal common ancestor of modern organisms. Before that split, life used two-step indirect pathways to synthesize asparagine and glutamine on their cognate tRNAs to form the aminoacyl-tRNA used in translation. These two-step pathways were retained throughout much of the bacterial and archaeal domains of life and eukaryotic organelles. The indirect routes use non-discriminating aminoacyl-tRNA synthetases (non-discriminating aspartyl-tRNA synthetase and non-discriminating glutamyl-tRNA synthetase) to misaminoacylate the tRNA. The misaminoacylated tRNA formed is then transamidated into the amide aminoacyl-tRNA used in protein synthesis by tRNA-dependent amidotransferases (GatCAB and GatDE). The enzymes and tRNAs involved assemble into complexes known as transamidosomes to help maintain translational fidelity. These pathways have evolved to meet the varied cellular needs across a diverse set of organisms, leading to significant variation. In certain bacteria, the indirect pathways may provide a means to adapt to cellular stress by reducing the fidelity of protein synthesis. The retention of these indirect pathways versus acquisition of asparaginyl-tRNA synthetase and glutaminyl tRNA synthetase in lineages likely involves a complex interplay of the competing uses of glutamine and asparagine beyond translation, energetic costs, co-evolution between enzymes and tRNA, and involvement in stress response that await further investigation.

Enzymic recognition of amino acids drove the evolution of primordial genetic codes

How genetic information gained its exquisite control over chemical processes needed to build living cells remains an enigma. Today, the aminoacyl-tRNA synthetases (AARS) execute the genetic codes in all living systems. But how did the AARS that emerged over three billion years ago as low-specificity, protozymic forms then spawn the full range of highly-specific enzymes that distinguish between 22 diverse amino acids? A phylogenetic reconstruction of extant AARS genes, enhanced by analysing modular acquisitions, reveals six AARS with distinct bacterial, archaeal, eukaryotic, or organellar clades, resulting in a total of 36 families of AARS catalytic domains. Small structural modules that differentiate one AARS family from another played pivotal roles in discriminating between amino acid side chains, thereby expanding the genetic code and refining its precision. The resulting model shows a tendency for less elaborate enzymes, with simpler catalytic domains, to activate amino acids that were not synthesised until later in the evolution of the code. The most probable evolutionary route for an emergent amino acid type to establish a place in the code was by recruiting older, less specific AARS, rather than adapting contemporary lineages. This process, retrofunctionalisation, differs from previously described mechanisms through which amino acids would enter the code.

Signatures of tRNAGlx -specificity in proteobacterial glutamyl-tRNA synthetases

The canonical function of glutamyl-tRNA synthetase (GluRS) is to glutamylate tRNAGlu . Yet not all bacterial GluRSs glutamylate tRNAGlu ; many glutamylate both tRNAGlu and tRNAGln , while some glutamylate only tRNAGln and not the cognate substrate tRNAGlu . Understanding the basis of the unique specificity of tRNAGlx is important. Mutational studies have hinted at hotspot residues, both on tRNAGlx and GluRS, which play crucial roles in tRNAGlx -specificity. However, its underlying structural basis remains unexplored. The majority of biochemical studies related to tRNAGlx -specificity have been performed on GluRS from Escherichia coli and other proteobacterial species. However, since the early crystal structures of GluRS and tRNAGlu -bound GluRS were from non-proteobacterial species (Thermus thermophilus), proteobacterial biochemical data have often been interpreted in the context of non-proteobacterial GluRS structures. Marked differences between proteobacterial and non-proteobacterial GluRSs have been demonstrated; therefore, it is important to understand tRNAGlx -specificity vis-a-vis proteobacterial GluRS structures. To this end, we solved the crystal structure of a double mutant GluRS from E. coli. Using the solved structure and several other currently available proteo- and non-proteobacterial GluRS crystal structures, we probed the structural basis of the tRNAGlx -specificity of bacterial GluRSs. Specifically, our analyses suggest a unique role played by the tRNAGlx D-helix contacting loop of GluRS in the modulation of tRNAGln -specificity. While earlier studies have identified functional hotspots on tRNAGlx that control the tRNAGlx -specificity of GluRS, this is the first report of complementary signatures of tRNAGlx -specificity in GluRS.

Regulation of BRCA1 stability through the tandem UBX domains of isoleucyl-tRNA synthetase 1

Aminoacyl-tRNA synthetases (ARSs) have evolved to acquire various additional domains. These domains allow ARSs to communicate with other cellular proteins in order to promote non-translational functions. Vertebrate cytoplasmic isoleucyl-tRNA synthetases (IARS1s) have an uncharacterized unique domain, UNE-I. Here, we present the crystal structure of the chicken IARS1 UNE-I complexed with glutamyl-tRNA synthetase 1 (EARS1). UNE-I consists of tandem ubiquitin regulatory X (UBX) domains that interact with a distinct hairpin loop on EARS1 and protect its neighboring proteins in the multi-synthetase complex from degradation. Phosphomimetic mutation of the two serine residues in the hairpin loop releases IARS1 from the complex. IARS1 interacts with BRCA1 in the nucleus, regulates its stability by inhibiting ubiquitylation via the UBX domains, and controls DNA repair function.

Remitting and exacerbating white matter lesions in leukoencephalopathy with thalamus and brainstem involvement and high lactate

BACKGROUND

Leukoencephalopathy with thalamus and brainstem involvement and high lactate (LTBL) is a hereditary disorder caused by biallelic variants in the EARS2 gene. Patients exhibit developmental delay, hypotonia, and hyperreflexia. Brain magnetic resonance imaging (MRI) reveals T2-hyperintensities in the deep white matter, thalamus, and brainstem, which generally stabilize over time. Herein, we report a case of LTBL, showing remitting and exacerbating white matter lesions.

CASE DESCRIPTION

A non-consanguineous Japanese boy exhibited unsteady head control with prominent hypotonia, with no family history of neurological diseases. Brain MRI at one year of age revealed extensive T2-hyperintensities on the cerebral white matter, cerebellum, thalamus, basal ganglia, pons, and medulla oblongata. Magnetic resonance spectroscopy of the lesions showed lactate and myoinositol peaks. Whole-exome sequencing yielded novel compound heterozygous EARS2 variants of c.164G>T, p.Arg55Leu and c.484C>T, p.Arg162Trp. Interestingly, the lesions were reduced at three years of age, and new lesions emerged at eight years of age. At 10 years of age, the lesions were changed in the corpus callosum, deep cerebral white matter, and cerebellum, without physical exacerbation. The lesions improved one year later.

CONCLUSION

We present the first case with remitting and exacerbating brain lesions in LTBL. EARS2 could relate to selective and specific brain regions and age dependency. Although the exact role of EARS2 remains unknown, the remitting and exacerbating imaging changes may be a clue in elucidating a novel EARS2 function in LTBL.

On the Track of the Missing tRNA Genes: A Source of Non-Canonical Functions?

Cellular tRNAs appear today as a diverse population of informative macromolecules with conserved general elements ensuring essential common functions and different and distinctive features securing specific interactions and activities. Their differential expression and the variety of post-transcriptional modifications they are subject to, lead to the existence of complex repertoires of tRNA populations adjusted to defined cellular states. Despite the tRNA-coding genes redundancy in prokaryote and eukaryote genomes, it is surprising to note the absence of genes coding specific translational-active isoacceptors throughout the phylogeny. Through the analysis of different releases of tRNA databases, this review aims to provide a general summary about those “missing tRNA genes.” This absence refers to both tRNAs that are not encoded in the genome, as well as others that show critical sequence variations that would prevent their activity as canonical translation adaptor molecules. Notably, while a group of genes are universally missing, others are absent in particular kingdoms. Functional information available allows to hypothesize that the exclusion of isodecoding molecules would be linked to: 1) reduce ambiguities of signals that define the specificity of the interactions in which the tRNAs are involved; 2) ensure the adaptation of the translational apparatus to the cellular state; 3) divert particular tRNA variants from ribosomal protein synthesis to other cellular functions. This leads to consider the “missing tRNA genes” as a source of putative non-canonical tRNA functions and to broaden the concept of adapter molecules in ribosomal-dependent protein synthesis.

Evolution of Life on Earth: tRNA, Aminoacyl-tRNA Synthetases and the Genetic Code

Life on Earth and the genetic code evolved around tRNA and the tRNA anticodon. We posit that the genetic code initially evolved to synthesize polyglycine as a cross-linking agent to stabilize protocells. We posit that the initial amino acids to enter the code occupied larger sectors of the code that were then invaded by incoming amino acids. Displacements of amino acids follow selection rules. The code sectored from a glycine code to a four amino acid code to an eight amino acid code to an ~16 amino acid code to the standard 20 amino acid code with stops. The proposed patterns of code sectoring are now most apparent from patterns of aminoacyl-tRNA synthetase evolution. The Elongation Factor-Tu GTPase anticodon-codon latch that checks the accuracy of translation appears to have evolved at about the eight amino acid to ~16 amino acid stage. Before evolution of the EF-Tu latch, we posit that both the 1st and 3rd anticodon positions were wobble positions. The genetic code evolved via tRNA charging errors and via enzymatic modifications of amino acids joined to tRNAs, followed by tRNA and aminoacyl-tRNA synthetase differentiation. Fidelity mechanisms froze the code by inhibiting further innovation.

PhyloToL: A Taxon/Gene-Rich Phylogenomic Pipeline to Explore Genome Evolution of Diverse Eukaryotes

Estimating multiple sequence alignments (MSAs) and inferring phylogenies are essential for many aspects of comparative biology. Yet, many bioinformatics tools for such analyses have focused on specific clades, with greatest attention paid to plants, animals, and fungi. The rapid increase in high-throughput sequencing (HTS) data from diverse lineages now provides opportunities to estimate evolutionary relationships and gene family evolution across the eukaryotic tree of life. At the same time, these types of data are known to be error-prone (e.g., substitutions, contamination). To address these opportunities and challenges, we have refined a phylogenomic pipeline, now named PhyloToL, to allow easy incorporation of data from HTS studies, to automate production of both MSAs and gene trees, and to identify and remove contaminants. PhyloToL is designed for phylogenomic analyses of diverse lineages across the tree of life (i.e., at scales of >100 My). We demonstrate the power of PhyloToL by assessing stop codon usage in Ciliophora, identifying contamination in a taxon- and gene-rich database and exploring the evolutionary history of chromosomes in the kinetoplastid parasite Trypanosoma brucei, the causative agent of African sleeping sickness. Benchmarking PhyloToL’s homology assessment against that of OrthoMCL and a published paper on superfamilies of bacterial and eukaryotic organellar outer membrane pore-forming proteins demonstrates the power of our approach for determining gene family membership and inferring gene trees. PhyloToL is highly flexible and allows users to easily explore HTS data, test hypotheses about phylogeny and gene family evolution and combine outputs with third-party tools (e.g., PhyloChromoMap, iGTP).

Aminoacyl-tRNA synthetases: Structure, function, and drug discovery

Aminoacyl-tRNA synthetases (AARSs) are the enzymes that catalyze the aminoacylation reaction by covalently linking an amino acid to its cognate tRNA in the first step of protein translation. Beyond this classical function, these enzymes are also known to have a role in several metabolic and signaling pathways that are important for cell viability. Study of these enzymes is of great interest to the researchers due to its pivotal role in the growth and survival of an organism. Further, unfolding the interesting structural and functional aspects of these enzymes in the last few years has qualified them as a potential drug target against various diseases. Here we review the classification, function, and the conserved as well the appended structural architecture of these enzymes in detail, including its association with multi-synthetase complexes. We also considered their role in human diseases in terms of mutations and autoantibodies against AARSs. Finally, we have discussed the available inhibitors against AARSs. This review offers comprehensive information on AARSs under a single canopy that would be a good inventory for researchers working in this area.

Nuclear codon reassignments in the genomics era and mechanisms behind their evolution

The canonical genetic code ubiquitously translates nucleotide into peptide sequence with several alterations known in viruses, bacteria, mitochondria, plastids, and single-celled eukaryotes. A new hypothesis to explain genetic code changes, termed tRNA loss driven codon reassignment, has been proposed recently when the polyphyly of the yeast codon reassignment events has been uncovered. According to this hypothesis, the driving force for genetic code changes are tRNA or translation termination factor loss-of-function mutations or loss-of-gene events. The free codon can subsequently be captured by all tRNAs that have an appropriately mutated anticodon and are efficiently charged. Thus, codon capture most likely happens by near-cognate tRNAs and tRNAs whose anticodons are not part of the recognition sites of the respective aminoacyl-tRNA-synthetases. This hypothesis comprehensively explains the CTG codon translation as alanine in Pachysolen yeast together with the long known translation of the same codon as serine in Candida albicans and related species, and can also be applied to most other known reassignments.

A water-mediated and substrate-assisted aminoacylation mechanism in the discriminating aminoacyl-tRNA synthetase GlnRS and non-discriminating GluRS

Computational methods have been used to elucidate key differences and similarities between the distinct aminoacyl-tRNA synthetases (aaRS) GlnRS and non-discriminating-GluRS.

Quest for Ancestors of Eukaryal Cells Based on Phylogenetic Analyses of Aminoacyl-tRNA Synthetases

The three-domain phylogenetic system of life has been challenged, particularly with regard to the position of Eukarya. The recent increase of known genome sequences has allowed phylogenetic analyses of all extant organisms using concatenated sequence alignment of universally conserved genes; these data supported the two-domain hypothesis, which place eukaryal species as ingroups of the Domain Archaea. However, the origin of Eukarya is complicated: the closest archaeal species to Eukarya differs in single-gene phylogenetic analyses depending on the genes. In this report, we performed molecular phylogenetic analyses of 23 aminoacyl-tRNA synthetases (ARS). Cytoplasmic ARSs in 12 trees showed a monophyletic Eukaryotic branch. One ARS originated from TACK superphylum. One ARS originated from Euryarchaeota and three originated from DPANN superphylum. Four ARSs originated from different bacterial species. The other 8 cytoplasmic ARSs were split into two or three groups in respective trees, which suggested that the cytoplasmic ARSs were replaced by secondary ARSs, and the original ARSs have been lost during evolution of Eukarya. In these trees, one original cytoplasmic ARS was derived from Euryarchaeota and three were derived from DPANN superphylum. Our results strongly support the two-domain hypothesis. We discovered that rampant-independent lateral gene transfers from several archaeal species of DPANN superphylum have contributed to the formation of Eukaryal cells. Based on our phylogenetic analyses, we proposed a model for the establishment of Eukarya.

A novel nuclear genetic code alteration in yeasts and the evolution of codon reassignment in eukaryotes

The genetic code is the cellular translation table for the conversion of nucleotide sequences into amino acid sequences. Changes to the meaning of sense codons would introduce errors into almost every translated message and are expected to be highly detrimental. However, reassignment of single or multiple codons in mitochondria and nuclear genomes, although extremely rare, demonstrates that the code can evolve. Several models for the mechanism of alteration of nuclear genetic codes have been proposed (including “codon capture,” “genome streamlining,” and “ambiguous intermediate” theories), but with little resolution. Here, we report a novel sense codon reassignment in Pachysolen tannophilus, a yeast related to the Pichiaceae. By generating proteomics data and using tRNA sequence comparisons, we show that Pachysolen translates CUG codons as alanine and not as the more usual leucine. The Pachysolen tRNA_CAG is an anticodon-mutated tRNA^Ala containing all major alanine tRNA recognition sites. The polyphyly of the CUG-decoding tRNAs in yeasts is best explained by a tRNA loss driven codon reassignment mechanism. Loss of the CUG-tRNA in the ancient yeast is followed by gradual decrease of respective codons and subsequent codon capture by tRNAs whose anticodon is not part of the aminoacyl-tRNA synthetase recognition region. Our hypothesis applies to all nuclear genetic code alterations and provides several testable predictions. We anticipate more codon reassignments to be uncovered in existing and upcoming genome projects.

The crystal structure of human GlnRS provides basis for the development of neurological disorders

Cytosolic glutaminyl-tRNA synthetase (GlnRS) is the singular enzyme responsible for translation of glutamine codons. Compound heterozygous mutations in GlnRS cause severe brain disorders by a poorly understood mechanism. Herein, we present crystal structures of the wild type and two pathological mutants of human GlnRS, which reveal, for the first time, the domain organization of the intact enzyme and the structure of the functionally important N-terminal domain (NTD). Pathological mutations mapping in the NTD alter the domain structure, and decrease catalytic activity and stability of GlnRS, whereas missense mutations in the catalytic domain induce misfolding of the enzyme. Our results suggest that the reduced catalytic efficiency and a propensity of GlnRS mutants to misfold trigger the disease development. This report broadens the spectrum of brain pathologies elicited by protein misfolding and provides a paradigm for understanding the role of mutations in aminoacyl-tRNA synthetases in neurological diseases.

Assembly of Multi-tRNA Synthetase Complex via Heterotetrameric Glutathione Transferase-homology Domains

Many multicomponent protein complexes mediating diverse cellular processes are assembled through scaffolds with specialized protein interaction modules. The multi-tRNA synthetase complex (MSC), consisting of nine different aminoacyl-tRNA synthetases and three non-enzymatic factors (AIMP1-3), serves as a hub for many signaling pathways in addition to its role in protein synthesis. However, the assembly process and structural arrangement of the MSC components are not well understood. Here we show the heterotetrameric complex structure of the glutathione transferase (GST) domains shared among the four MSC components, methionyl-tRNA synthetase (MRS), glutaminyl-prolyl-tRNA synthetase (EPRS), AIMP2 and AIMP3. The MRS-AIMP3 and EPRS-AIMP2 using interface 1 are bridged via interface 2 of AIMP3 and EPRS to generate a unique linear complex of MRS-AIMP3:EPRS-AIMP2 at the molar ratio of (1:1):(1:1). Interestingly, the affinity at interface 2 of AIMP3:EPRS can be varied depending on the occupancy of interface 1, suggesting the dynamic nature of the linear GST tetramer. The four components are optimally arranged for maximal accommodation of additional domains and proteins. These characteristics suggest the GST tetramer as a unique and dynamic structural platform from which the MSC components are assembled. Considering prevalence of the GST-like domains, this tetramer can also provide a tool for the communication of the MSC with other GST-containing cellular factors.

Dispensability of zinc and the putative zinc-binding domain in bacterial glutamyl-tRNA synthetase

The putative zinc-binding domain (pZBD) in Escherichia coli glutamyl-tRNA synthetase (GluRS) is known to correctly position the tRNA acceptor arm and modulate the amino acid-binding site. However, its functional role in other bacterial species is not clear since many bacterial GluRSs lack a zinc-binding motif in the pZBD. From experimental studies on pZBD-swapped E. coli GluRS, with Thermosynechoccus elongatus GluRS, Burkholderia thailandensis GluRS and E. coli glutamyl-queuosine-tRNA^Asp synthetase (Glu-Q-RS), we show that E. coli GluRS, containing the zinc-free pZBD of B. thailandensis, is as functional as the zinc-bound wild-type E. coli GluRS, whereas the other constructs, all zinc-bound, show impaired function. A pZBD-tinkered version of E. coli GluRS that still retained Zn-binding capacity, also showed reduced activity. This suggests that zinc is not essential for the pZBD to be functional. From extensive structural and sequence analyses from whole genome database of bacterial GluRS, we further show that in addition to many bacterial GluRS lacking a zinc-binding motif, the pZBD is actually deleted in some bacteria, all containing either glutaminyl-tRNA synthetase (GlnRS) or a second copy of GluRS (GluRS2). Correlation between the absence of pZBD and the occurrence of glutamine amidotransferase CAB (GatCAB) in the genome suggests that the primordial role of the pZBD was to facilitate transamidation of misacylated Glu-tRNA^Gln via interaction with GatCAB, whereas its role in tRNA^Glu interaction may be a consequence of the presence of pZBD.

Zinc is functionally important in glutamylation of tRNA^Glu in Escherichia coli, yet, it is absent from many bacterial glutamyl-tRNA synthetases (GluRSs). We demonstrate and rationalize why zinc or the putative zinc-binding domain (pZBD) is not indispensable in all bacterial GluRSs.

Expanded microbial genome coverage and improved protein family annotation in the COG database

Microbial genome sequencing projects produce numerous sequences of deduced proteins, only a small fraction of which have been or will ever be studied experimentally. This leaves sequence analysis as the only feasible way to annotate these proteins and assign to them tentative functions. The Clusters of Orthologous Groups of proteins (COGs) database (http://www.ncbi.nlm.nih.gov/COG/), first created in 1997, has been a popular tool for functional annotation. Its success was largely based on (i) its reliance on complete microbial genomes, which allowed reliable assignment of orthologs and paralogs for most genes; (ii) orthology-based approach, which used the function(s) of the characterized member(s) of the protein family (COG) to assign function(s) to the entire set of carefully identified orthologs and describe the range of potential functions when there were more than one; and (iii) careful manual curation of the annotation of the COGs, aimed at detailed prediction of the biological function(s) for each COG while avoiding annotation errors and overprediction. Here we present an update of the COGs, the first since 2003, and a comprehensive revision of the COG annotations and expansion of the genome coverage to include representative complete genomes from all bacterial and archaeal lineages down to the genus level. This re-analysis of the COGs shows that the original COG assignments had an error rate below 0.5% and allows an assessment of the progress in functional genomics in the past 12 years. During this time, functions of many previously uncharacterized COGs have been elucidated and tentative functional assignments of many COGs have been validated, either by targeted experiments or through the use of high-throughput methods. A particularly important development is the assignment of functions to several widespread, conserved proteins many of which turned out to participate in translation, in particular rRNA maturation and tRNA modification. The new version of the COGs is expected to become an important tool for microbial genomics.

Gene rearrangements in gekkonid mitochondrial genomes with shuffling, loss, and reassignment of tRNA genes

BACKGROUND

Vertebrate mitochondrial genomes (mitogenomes) are 16-18 kbp double-stranded circular DNAs that encode a set of 37 genes. The arrangement of these genes and the major noncoding region is relatively conserved through evolution although gene rearrangements have been described for diverse lineages. The tandem duplication-random loss model has been invoked to explain the mechanisms of most mitochondrial gene rearrangements. Previously reported mitogenomic sequences for geckos rarely included gene rearrangements, which we explore in the present study.

RESULTS

We determined seven new mitogenomic sequences from Gekkonidae using a high-throughput sequencing method. The Tropiocolotes tripolitanus mitogenome involves a tandem duplication of the gene block: tRNAArg, NADH dehydrogenase subunit 4L, and NADH dehydrogenase subunit 4. One of the duplicate copies for each protein-coding gene may be pseudogenized. A duplicate copy of the tRNAArg gene appears to have been converted to a tRNAGln gene by a C to T base substitution at the second anticodon position, although this gene may not be fully functional in protein synthesis. The Stenodactylus petrii mitogenome includes several tandem duplications of tRNALeu genes, as well as a translocation of the tRNAAla gene and a putative origin of light-strand replication within a tRNA gene cluster. Finally, the Uroplatus fimbriatus and U. ebenaui mitogenomes feature the apparent loss of the tRNAGlu gene from its original position. Uroplatus fimbriatus appears to retain a translocated tRNAGlu gene adjacent to the 5' end of the major noncoding region.

CONCLUSIONS

The present study describes several new mitochondrial gene rearrangements from Gekkonidae. The loss and reassignment of tRNA genes is not very common in vertebrate mitogenomes and our findings raise new questions as to how missing tRNAs are supplied and if the reassigned tRNA gene is fully functional. These new examples of mitochondrial gene rearrangements in geckos should broaden our understanding of the evolution of mitochondrial gene arrangements.

Coevolution of specificity determinants in eukaryotic glutamyl- and glutaminyl-tRNA synthetases

The glutaminyl-tRNA synthetase (GlnRS) enzyme, which pairs glutamine with tRNA(Gln) for protein synthesis, evolved by gene duplication in early eukaryotes from a nondiscriminating glutamyl-tRNA synthetase (GluRS) that aminoacylates both tRNA(Gln) and tRNA(Glu) with glutamate. This ancient GluRS also separately differentiated to exclude tRNA(Gln) as a substrate, and the resulting discriminating GluRS and GlnRS further acquired additional protein domains assisting function in cis (the GlnRS N-terminal Yqey domain) or in trans (the Arc1p protein associating with GluRS). These added domains are absent in contemporary bacterial GlnRS and GluRS. Here, using Saccharomyces cerevisiae enzymes as models, we find that the eukaryote-specific protein domains substantially influence amino acid binding, tRNA binding and aminoacylation efficiency, but they play no role in either specific nucleotide readout or discrimination against noncognate tRNA. Eukaryotic tRNA(Gln) and tRNA(Glu) recognition determinants are found in equivalent positions and are mutually exclusive to a significant degree, with key nucleotides located adjacent to portions of the protein structure that differentiated during the evolution of archaeal nondiscriminating GluRS to GlnRS. These findings provide important corroboration for the evolutionary model and suggest that the added eukaryotic domains arose in response to distinctive selective pressures associated with the greater complexity of the eukaryotic translational apparatus. We also find that the affinity of GluRS for glutamate is significantly increased when Arc1p is not associated with the enzyme. This is consistent with the lower concentration of intracellular glutamate and the dissociation of the Arc1p:GluRS complex upon the diauxic shift to respiratory conditions.

Evolutionary insights about bacterial GlxRS from whole genome analyses: is GluRS2 a chimera?

Background

Evolutionary histories of glutamyl-tRNA synthetase (GluRS) and glutaminyl-tRNA synthetase (GlnRS) in bacteria are convoluted. After the divergence of eubacteria and eukarya, bacterial GluRS glutamylated both tRNA^Gln and tRNA^Glu until GlnRS appeared by horizontal gene transfer (HGT) from eukaryotes or a duplicate copy of GluRS (GluRS2) that only glutamylates tRNA^Gln appeared. The current understanding is based on limited sequence data and not always compatible with available experimental results. In particular, the origin of GluRS2 is poorly understood.

Results

A large database of bacterial GluRS, GlnRS, tRNA^Gln and the trimeric aminoacyl-tRNA-dependent amidotransferase (gatCAB), constructed from whole genomes by functionally annotating and classifying these enzymes according to their mutual presence and absence in the genome, was analyzed. Phylogenetic analyses showed that the catalytic and the anticodon-binding domains of functional GluRS2 (as in Helicobacter pylori) were independently acquired from evolutionarily distant hosts by HGT. Non-functional GluRS2 (as in Thermotoga maritima), on the other hand, was found to contain an anticodon-binding domain appended to a gene-duplicated catalytic domain. Several genomes were found to possess both GluRS2 and GlnRS, even though they share the common function of aminoacylating tRNA^Gln. GlnRS was widely distributed among bacterial phyla and although phylogenetic analyses confirmed the origin of most bacterial GlnRS to be through a single HGT from eukarya, many GlnRS sequences also appeared with evolutionarily distant phyla in phylogenetic tree. A GlnRS pseudogene could be identified in Sorangium cellulosum.

Conclusions

Our analysis broadens the current understanding of bacterial GlxRS evolution and highlights the idiosyncratic evolution of GluRS2. Specifically we show that: i) GluRS2 is a chimera of mismatching catalytic and anticodon-binding domains, ii) the appearance of GlnRS and GluRS2 in a single bacterial genome indicating that the evolutionary histories of the two enzymes are distinct, iii) GlnRS is more widespread in bacteria than is believed, iv) bacterial GlnRS appeared both by HGT from eukarya and intra-bacterial HGT, v) presence of GlnRS pseudogene shows that many bacteria could not retain the newly acquired eukaryal GlnRS. The functional annotation of GluRS, without recourse to experiments, performed in this work, demonstrates the inherent and unique advantages of using whole genome over isolated sequence databases.

The structure of yeast glutaminyl-tRNA synthetase and modeling of its interaction with tRNA

Eukaryotic glutaminyl-tRNA synthetase (GlnRS) contains an appended N-terminal domain (NTD) whose precise function is unknown. Although GlnRS structures from two prokaryotic species are known, no eukaryotic GlnRS structure has been reported. Here we present the first crystallographic structure of yeast GlnRS, finding that the structure of the C-terminal domain is highly similar to Escherichia coli GlnRS but that 214 residues, including the NTD, are crystallographically disordered. We present a model of the full-length enzyme in solution, using the structures of the C-terminal domain, and the isolated NTD, with small-angle X-ray scattering data of the full-length molecule. We proceed to model the enzyme bound to tRNA, using the crystallographic structures of GatCAB and GlnRS-tRNA complex from bacteria. We contrast the tRNA-bound model with the tRNA-free solution state and perform molecular dynamics on the full-length GlnRS-tRNA complex, which suggests that tRNA binding involves the motion of a conserved hinge in the NTD.

Emergence and evolution

The aminoacyl-tRNA synthetases (aaRSs) are essential components of the protein synthesis machinery responsible for defining the genetic code by pairing the correct amino acids to their cognate tRNAs. The aaRSs are an ancient enzyme family believed to have origins that may predate the last common ancestor and as such they provide insights into the evolution and development of the extant genetic code. Although the aaRSs have long been viewed as a highly conserved group of enzymes, findings within the last couple of decades have started to demonstrate how diverse and versatile these enzymes really are. Beyond their central role in translation, aaRSs and their numerous homologs have evolved a wide array of alternative functions both inside and outside translation. Current understanding of the emergence of the aaRSs, and their subsequent evolution into a functionally diverse enzyme family, are discussed in this chapter.

Structural diversity and protein engineering of the aminoacyl-tRNA synthetases

Aminoacyl-tRNA synthetases (aaRS) are the enzymes that ensure faithful transmission of genetic information in all living cells, and are central to the developing technologies for expanding the capacity of the translation apparatus to incorporate nonstandard amino acids into proteins in vivo. The 24 known aaRS families are divided into two classes that exhibit functional evolutionary convergence. Each class features an active site domain with a common fold that binds ATP, the amino acid, and the 3'-terminus of tRNA, embellished by idiosyncratic further domains that bind distal portions of the tRNA and enhance specificity. Fidelity in the expression of the genetic code requires that the aaRS be selective for both amino acids and tRNAs, a substantial challenge given the presence of structurally very similar noncognate substrates of both types. Here we comprehensively review central themes concerning the architectures of the protein structures and the remarkable dual-substrate selectivities, with a view toward discerning the most important issues that still substantially limit our capacity for rational protein engineering. A suggested general approach to rational design is presented, which should yield insight into the identities of the protein-RNA motifs at the heart of the genetic code, while also offering a basis for improving the catalytic properties of engineered tRNA synthetases emerging from genetic selections.

Rational design and directed evolution of a bacterial-type glutaminyl-tRNA synthetase precursor

Protein biosynthesis requires aminoacyl-transfer RNA (tRNA) synthetases to provide aminoacyl-tRNA substrates for the ribosome. Most bacteria and all archaea lack a glutaminyl-tRNA synthetase (GlnRS); instead, Gln-tRNA^Gln is produced via an indirect pathway: a glutamyl-tRNA synthetase (GluRS) first attaches glutamate (Glu) to tRNA^Gln, and an amidotransferase converts Glu-tRNA^Gln to Gln-tRNA^Gln. The human pathogen Helicobacter pylori encodes two GluRS enzymes, with GluRS2 specifically aminoacylating Glu onto tRNA^Gln. It was proposed that GluRS2 is evolving into a bacterial-type GlnRS. Herein, we have combined rational design and directed evolution approaches to test this hypothesis. We show that, in contrast to wild-type (WT) GlnRS2, an engineered enzyme variant (M110) with seven amino acid changes is able to rescue growth of the temperature-sensitive Escherichia coli glnS strain UT172 at its non-permissive temperature. In vitro kinetic analyses reveal that WT GluRS2 selectively acylates Glu over Gln, whereas M110 acylates Gln 4-fold more efficiently than Glu. In addition, M110 hydrolyzes adenosine triphosphate 2.5-fold faster in the presence of Glu than Gln, suggesting that an editing activity has evolved in this variant to discriminate against Glu. These data imply that GluRS2 is a few steps away from evolving into a GlnRS and provides a paradigm for studying aminoacyl-tRNA synthetase evolution using directed engineering approaches.

A functional loop spanning distant domains of glutaminyl-tRNA synthetase also stabilizes a molten globule state

Molten globule and other disordered states of proteins are now known to play important roles in many cellular processes. From equilibrium unfolding studies of two paralogous proteins and their variants, glutaminyl-tRNA synthetase (GlnRS) and two of its variants [glutamyl-tRNA synthetase (GluRS) and its isolated domains, and a GluRS-GlnRS chimera], we demonstrate that only GlnRS forms a molten globule-like intermediate at low urea concentrations. We demonstrated that a loop in the GlnRS C-terminal anticodon binding domain that promotes communication with the N-terminal domain and indirectly modulates amino acid binding is also responsible for stabilization of the molten globule state. This loop was inserted into GluRS in the eukaryotic branch after the archaea-eukarya split, right around the time when GlnRS evolved. Because of the structural and functional importance of the loop, it is proposed that the insertion of the loop into a putative ancestral GluRS in eukaryotes produced a catalytically active molten globule state. Because of their enhanced dynamic nature, catalytically active molten globules are likely to possess broad substrate specificity. It is further proposed that the putative broader substrate specificity allowed the catalytically active molten globule to accept glutamine in addition to glutamic acid, leading to the evolution of GlnRS.

Riboswitch (T-box)-mediated control of tRNA-dependent amidation in Clostridium acetobutylicum rationalizes gene and pathway redundancy for asparagine and asparaginyl-trnaasn synthesis

Analysis of the Gram-positive Clostridium acetobutylicum genome reveals an inexplicable level of redundancy for the genes putatively involved in asparagine (Asn) and Asn-tRNA(Asn) synthesis. Besides a duplicated set of gatCAB tRNA-dependent amidotransferase genes, there is a triplication of aspartyl-tRNA synthetase genes and a duplication of asparagine synthetase B genes. This genomic landscape leads to the suspicion of the incoherent simultaneous use of the direct and indirect pathways of Asn and Asn-tRNA(Asn) formation. Through a combination of biochemical and genetic approaches, we show that C. acetobutylicum forms Asn and Asn-tRNA(Asn) by tRNA-dependent amidation. We demonstrate that an entire transamidation pathway composed of aspartyl-tRNA synthetase and one set of GatCAB genes is organized as an operon under the control of a tRNA(Asn)-dependent T-box riboswitch. Finally, our results suggest that this exceptional gene redundancy might be interconnected to control tRNA-dependent Asn synthesis, which in turn might be involved in controlling the metabolic switch from acidogenesis to solventogenesis in C. acetobutylicum.

Structural conservation of an ancient tRNA sensor in eukaryotic glutaminyl-tRNA synthetase

In all organisms, aminoacyl tRNA synthetases covalently attach amino acids to their cognate tRNAs. Many eukaryotic tRNA synthetases have acquired appended domains, whose origin, structure and function are poorly understood. The N-terminal appended domain (NTD) of glutaminyl-tRNA synthetase (GlnRS) is intriguing since GlnRS is primarily a eukaryotic enzyme, whereas in other kingdoms Gln-tRNA^Gln is primarily synthesized by first forming Glu-tRNA^Gln, followed by conversion to Gln-tRNA^Gln by a tRNA-dependent amidotransferase. We report a functional and structural analysis of the NTD of Saccharomyces cerevisiae GlnRS, Gln4. Yeast mutants lacking the NTD exhibit growth defects, and Gln4 lacking the NTD has reduced complementarity for tRNA^Gln and glutamine. The 187-amino acid Gln4 NTD, crystallized and solved at 2.3 Å resolution, consists of two subdomains, each exhibiting an extraordinary structural resemblance to adjacent tRNA specificity-determining domains in the GatB subunit of the GatCAB amidotransferase, which forms Gln-tRNA^Gln. These subdomains are connected by an apparent hinge comprised of conserved residues. Mutation of these amino acids produces Gln4 variants with reduced affinity for tRNA^Gln, consistent with a hinge-closing mechanism proposed for GatB recognition of tRNA. Our results suggest a possible origin and function of the NTD that would link the phylogenetically diverse mechanisms of Gln-tRNA^Gln synthesis.

Rational design of an evolutionary precursor of glutaminyl-tRNA synthetase

The specificity of most aminoacyl-tRNA synthetases for an amino acid and cognate tRNA pair evolved before the divergence of the three domains of life. Glutaminyl-tRNA synthetase (GlnRS) evolved later and is derived from the archaeal-type nondiscriminating glutamyl-tRNA synthetase (GluRS), an enzyme with relaxed tRNA specificity capable of forming both Glu-tRNA(Glu) and Glu-tRNA(Gln). The archaea lack GlnRS and use a specialized amidotransferase to convert Glu-tRNA(Gln) to Gln-tRNA(Gln) needed for protein synthesis. We show that the Methanothermobacter thermautotrophicus GluRS is active toward tRNA(Glu) and the two tRNA(Gln) isoacceptors the organism encodes, but with a significant catalytic preference for tRNA(Gln2)(CUG). The less active tRNA(Gln1)(UUG) responds to the less common CAA codon for Gln. From a biochemical characterization of M. thermautotrophicus GluRS variants, we found that the evolution of tRNA specificity in GlnRS could be recapitulated by converting the M. thermautotrophicus GluRS to a tRNA(Gln) specific enzyme, solely through the addition of an acceptor stem loop present in bacterial GlnRS. One designed GluRS variant is also highly specific for the tRNA(Gln2)(CUG) isoacceptor, which responds to the CAG codon, and shows no activity toward tRNA(Gln1)(UUG). Because it is now possible to eliminate particular codons from the genome of Escherichia coli, additional codons will become available for genetic code engineering. Isoacceptor-specific aminoacyl-tRNA synthetases will enable the reassignment of more open codons while preserving accurate encoding of the 20 canonical amino acids.

Archaeal 3'-phosphate RNA splicing ligase characterization identifies the missing component in tRNA maturation

Intron removal from tRNA precursors involves cleavage by a tRNA splicing endonuclease to yield tRNA 3'-halves beginning with a 5'-hydroxyl, and 5'-halves ending in a 2',3'-cyclic phosphate. A tRNA ligase then incorporates this phosphate into the internucleotide bond that joins the two halves. Although this 3'-P RNA splicing ligase activity was detected almost three decades ago in extracts from animal and later archaeal cells, the protein responsible was not yet identified. Here we report the purification of this ligase from Methanopyrus kandleri cells, and its assignment to the still uncharacterized RtcB protein family. Studies with recombinant Pyrobaculum aerophilum RtcB showed that the enzyme is able to join spliced tRNA halves to mature-sized tRNAs where the joining phosphodiester linkage contains the phosphate originally present in the 2',3'-cyclic phosphate. The data confirm RtcB as the archaeal RNA 3'-P ligase. Structural genomics efforts previously yielded a crystal structure of the Pyrococcus horikoshii RtcB protein containing a new protein fold and a conserved putative Zn(2+) binding cleft. This structure guided our mutational analysis of the P. aerophilum enzyme. Mutations of highly conserved residues in the cleft (C100A, H205A, H236A) rendered the enzyme inactive suggesting these residues to be part of the active site of the P. aerophilum ligase. There is no significant sequence similarity between the active sites of P. aerophilum ligase and that of T4 RNA ligase, nor ligases from plants and fungi. RtcB sequence conservation in archaea and in eukaryotes implicates eukaryotic RtcB as the long-sought animal 3'-P RNA ligase.

Synthesis of Glu-tRNA(Gln) by engineered and natural aminoacyl-tRNA synthetases

A protein engineering approach to delineating which distinct elements of homologous tRNA synthetase architectures are responsible for divergent RNA-amino acid pairing specificities is described. Previously, we constructed a hybrid enzyme in which 23 amino acids from the catalytic domain of Escherichia coli glutaminyl-tRNA synthetase (GlnRS) were replaced with the corresponding residues of human glutamyl-tRNA synthetase (GluRS). The engineered hybrid (GlnRS S1/L1/L2) synthesizes Glu-tRNA(Gln) more than 10(4)-fold more efficiently than GlnRS. Detailed comparison of kinetic parameters between GlnRS S1/L1/L2 and the naturally occurring Methanothermobacter thermautotrophicus GluRS(ND), which is also capable of Glu-tRNA(Gln) synthesis, now shows that both k(cat) and K(m) for glutamate are recapitulated in the engineered enzyme, but that K(m) for tRNA is 200-fold higher. Thus, the simultaneous optimization of paired amino acid and tRNA binding sites found in a naturally occurring enzyme is not recapitulated in a hybrid that is successfully engineered for amino acid complementarity. We infer that the GlnRS architecture has differentiated to match only cognate amino acid-RNA pairs, and that the substrate selection functions do not operate independently of each other. Design and characterization of four additional hybrids identify further residues involved in improving complementarity for glutamate and in communicating between amino acid and tRNA binding sites. The robust catalytic function demonstrated in this engineered system offers a novel platform for exploring the stereochemical origins of coding as a property of the ancient Rossmann fold.