Evolution of the Glx-tRNA synthetase family: the glutaminyl enzyme as a case of horizontal gene transfer

Order of amino acid recruitment into the genetic code resolved by Last Universal Common Ancestor's protein domains

The current "consensus" order in which amino acids were added to the genetic code is based on potentially biased criteria such as absence of sulfur-containing amino acids from the Urey-Miller experiment which lacked sulfur. Even if inferred perfectly, abiotic abundance might not reflect abundance in the organisms in which the genetic code evolved. Here, we instead exploit the fact that proteins that emerged prior to the genetic code's completion are likely enriched in early amino acids and depleted in late amino acids. We identify the most ancient protein-coding sequences born prior to the archaeal-bacterial split. Amino acid usage in protein sequences whose ancestors date back to a single homolog in the Last Universal Common Ancestor (LUCA) largely matches the consensus order. However, our findings indicate that metal-binding (cysteine and histidine) and sulfur-containing (cysteine and methionine) amino acids were added to the genetic code much earlier than previously thought. Surprisingly, even more ancient protein sequences - those that had already diversified into multiple distinct copies in LUCA - show a different pattern to single copy LUCA sequences: significantly less depleted in the late amino acids tryptophan and tyrosine, and enriched rather than depleted in phenylalanine. This is compatible with at least some of these sequences predating the current genetic code. Their distinct enrichment patterns thus provide hints about earlier, alternative genetic codes.

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.

Microbial driven genetic variation in holobionts

Genetic variation in holobionts, (host and microbiome), occurring by changes in both host and microbiome genomes, can be observed from two perspectives: observable variations and the processes that bring about the variation. The observable includes the enormous genetic diversity of prokaryotes, which gave rise to eukaryotic organisms. Holobionts then evolved a rich microbiome with a stable core containing essential genes, less so common taxa, and a more diverse non-core enabling considerable genetic variation. The result being that, the human gut microbiome, for example, contains 1,000 times more unique genes than are present in the human genome. Microbial driven genetic variation processes in holobionts include: (1) Acquisition of novel microbes from the environment, which bring in multiple genes in one step, (2) amplification/reduction of certain microbes in the microbiome, that contribute to holobiont` s adaptation to changing conditions, (3) horizontal gene transfer between microbes and between microbes and host, (4) mutation, which plays an important role in optimizing interactions between different microbiota and between microbiota and host. We suggest that invertebrates and plants, where microbes can live intracellularly, have a greater chance of genetic exchange between microbiota and host, thus a greater chance of vertical transmission and a greater effect of microbiome on evolution of host than vertebrates. However, even in vertebrates the microbiome can aid in environmental fluctuations by amplification/reduction and by acquisition of novel microorganisms.

Did Amino Acid Side Chain Reactivity Dictate the Composition and Timing of Aminoacyl-tRNA Synthetase Evolution?

The twenty amino acids in the standard genetic code were fixed prior to the last universal common ancestor (LUCA). Factors that guided this selection included establishment of pathways for their metabolic synthesis and the concomitant fixation of substrate specificities in the emerging aminoacyl-tRNA synthetases (aaRSs). In this conceptual paper, we propose that the chemical reactivity of some amino acid side chains (e.g., lysine, cysteine, homocysteine, ornithine, homoserine, and selenocysteine) delayed or prohibited the emergence of the corresponding aaRSs and helped define the amino acids in the standard genetic code. We also consider the possibility that amino acid chemistry delayed the emergence of the glutaminyl- and asparaginyl-tRNA synthetases, neither of which are ubiquitous in extant organisms. We argue that fundamental chemical principles played critical roles in fixation of some aspects of the genetic code pre- and post-LUCA.

Aminoacyl-tRNA synthetases

The aminoacyl-tRNA synthetases are an essential and universally distributed family of enzymes that plays a critical role in protein synthesis, pairing tRNAs with their cognate amino acids for decoding mRNAs according to the genetic code. Synthetases help to ensure accurate translation of the genetic code by using both highly accurate cognate substrate recognition and stringent proofreading of noncognate products. While alterations in the quality control mechanisms of synthetases are generally detrimental to cellular viability, recent studies suggest that in some instances such changes facilitate adaption to stress conditions. Beyond their central role in translation, synthetases are also emerging as key players in an increasing number of other cellular processes, with far-reaching consequences in health and disease. The biochemical versatility of the synthetases has also proven pivotal in efforts to expand the genetic code, further emphasizing the wide-ranging roles of the aminoacyl-tRNA synthetase family in synthetic and natural biology.

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.

The phylogenetic distribution of the glutaminyl-tRNA synthetase and Glu-tRNAGln amidotransferase in the fundamental lineages would imply that the ancestor of archaea, that of eukaryotes and LUCA were progenotes

The function of the glutaminyl-tRNA synthetase and Glu-tRNA^Gln amidotransferase might be related to the origin of the genetic code because, for example, glutaminyl-tRNA synthetase catalyses the fundamental reaction that makes the genetic code. If the evolutionary stage of the origin of these two enzymes could be unambiguously identified, then the genetic code should still have been originating at that particular evolutionary stage because the fundamental reaction that makes the code itself was still evidently evolving. This would result in that particular evolutionary moment being attributed to the evolutionary stage of the progenote because it would have a relationship between the genotype and the phenotype not yet fully realized because the genetic code was precisely still originating. I then analyzed the distribution of the glutaminyl-tRNA synthetase and Glu-tRNA^Gln aminodotrasferase in the main phyletic lineages. Since in some cases the origin of these two enzymes can be related to the evolutionary stages of ancestors of archaea and eukaryotes, this would indicate these ancestors as progenotes because at that evolutionary moment the genetic code was evidently still evolving, thus realizing the definition of progenote. The conclusion that the ancestor of archaea and that of eukaryotes were progenotes would imply that even the last universal common ancestor (LUCA) was a progenote because it appeared, on the tree of life, temporally before these ancestors.

Universal Codons with Enrichment from GC to AU Nucleotide Composition Reveal a Chronological Assignment from Early to Late Along with LUCA Formation

The emergence of a primitive genetic code should be considered the most essential event during the origin of life. Almost a complete set of codons (as we know them) should have been established relatively early during the evolution of the last universal common ancestor (LUCA) from which all known organisms descended. Many hypotheses have been proposed to explain the driving forces and chronology of the evolution of the genetic code; however, none is commonly accepted. In the current paper, we explore the features of the genetic code that, in our view, reflect the mechanism and the chronological order of the origin of the genetic code. Our hypothesis postulates that the primordial RNA was mostly GC-rich, and this bias was reflected in the order of amino acid codon assignment. If we arrange the codons and their corresponding amino acids from GC-rich to AU-rich, we find that: 1. The amino acids encoded by GC-rich codons (Ala, Gly, Arg, and Pro) are those that contribute the most to the interactions with RNA (if incorporated into short peptides). 2. This order correlates with the addition of novel functions necessary for the evolution from simple to longer folded peptides. 3. The overlay of aminoacyl-tRNA synthetases (aaRS) to the amino acid order produces a distinctive zonal distribution for class I and class II suggesting an interdependent origin. These correlations could be explained by the active role of the bridge peptide (BP), which we proposed earlier in the evolution of the genetic code.

How evolution shapes enzyme selectivity - lessons from aminoacyl-tRNA synthetases and other amino acid utilizing enzymes

Aminoacyl-tRNA synthetases (AARSs) charge tRNA with their cognate amino acids. Many other enzymes use amino acids as substrates, yet discrimination against noncognate amino acids that threaten the accuracy of protein translation is a hallmark of AARSs. Comparing AARSs to these other enzymes allowed us to recognize patterns in molecular recognition and strategies used by evolution for exercising selectivity. Overall, AARSs are 2-3 orders of magnitude more selective than most other amino acid utilizing enzymes. AARSs also reveal the physicochemical limits of molecular discrimination. For example, amino acids smaller by a single methyl moiety present a discrimination ceiling of ~200, while larger ones can be discriminated by up to 10⁵ -fold. In contrast, substrates larger by a hydroxyl group challenge AARS selectivity, due to promiscuous H-bonding with polar active site groups. This 'hydroxyl paradox' is resolved by editing. Indeed, when the physicochemical discrimination limits are reached, post-transfer editing - hydrolysis of tRNAs charged with noncognate amino acids, evolved. The editing site often selectively recognizes the edited noncognate substrate using the very same feature that the synthetic site could not efficiently discriminate against. Finally, the comparison to other enzymes also reveals that the selectivity of AARSs is an explicitly evolved trait, showing some clear examples of how selection acted not only to optimize catalytic efficiency with the target substrate, but also to abolish activity with noncognate threat substrates ('negative selection').

The complex evolutionary history of aminoacyl-tRNA synthetases

Aminoacyl-tRNA synthetases (AARSs) are a superfamily of enzymes responsible for the faithful translation of the genetic code and have lately become a prominent target for synthetic biologists. Our large-scale analysis of >2500 prokaryotic genomes reveals the complex evolutionary history of these enzymes and their paralogs, in which horizontal gene transfer played an important role. These results show that a widespread belief in the evolutionary stability of this superfamily is misconceived. Although AlaRS, GlyRS, LeuRS, IleRS, ValRS are the most stable members of the family, GluRS, LysRS and CysRS often have paralogs, whereas AsnRS, GlnRS, PylRS and SepRS are often absent from many genomes. In the course of this analysis, highly conserved protein motifs and domains within each of the AARS loci were identified and used to build a web-based computational tool for the genome-wide detection of AARS coding sequences. This is based on hidden Markov models (HMMs) and is available together with a cognate database that may be used for specific analyses. The bioinformatics tools that we have developed may also help to identify new antibiotic agents and targets using these essential enzymes. These tools also may help to identify organisms with alternative pathways that are involved in maintaining the fidelity of the genetic code.

Difficulty in obtaining the complete mRNA coding sequence at 5' region (5' end mRNA artifact): Causes, consequences in biology and medicine and possible solutions for obtaining the actual amino acid sequence of proteins (Review)

The known difficulty in obtaining the actual full length, complete sequence of a messenger RNA (mRNA) may lead to the erroneous determination of its coding sequence at the 5' region (5' end mRNA artifact), and consequently to the wrong assignment of the translation start codon, leading to the inaccurate prediction of the encoded polypeptide at its amino terminus. Among the known human genes whose study was affected by this artifact, we can include disco interacting protein 2 homolog A (DIP2A; KIAA0184), Down syndrome critical region 1 (DSCR1), SON DNA binding protein (SON), trefoil factor 3 (TFF3) and URB1 ribosome biogenesis 1 homolog (URB1; KIAA0539) on chromosome 21, as well as receptor for activated C kinase 1 (RACK1, also known as GNB2L1), glutaminyl‑tRNA synthetase (QARS) and tyrosyl-DNA phosphodiesterase 2 (TDP2) along with another 474 loci, including interleukin 16 (IL16). In this review, we discuss the causes of this issue, its quantitative incidence in biomedical research, the consequences in biology and medicine, and the possible solutions for obtaining the actual amino acid sequence of proteins in the post-genomics era.

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.

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.

Evolutionary Ancestry of Eukaryotic Protein Kinases and Choline Kinases

The reversible phosphorylation of proteins catalyzed by protein kinases in eukaryotes supports an important role for eukaryotic protein kinases (ePKs) in the emergence of nucleated cells in the third superkingdom of life. Choline kinases (ChKs) could also be critical in the early evolution of eukaryotes, because of their function in the biosynthesis of phosphatidylcholine, which is unique to eukaryotic membranes. However, the genomic origins of ePKs and ChKs are unclear. The high degeneracy of protein sequences and broad expansion of ePK families have made this fundamental question difficult to answer. In this study, we identified two class-I aminoacyl-tRNA synthetases with high similarities to consensus amino acid sequences of human protein-serine/threonine kinases. Comparisons of primary and tertiary structures supported that ePKs and ChKs evolved from a common ancestor related to glutaminyl aminoacyl-tRNA synthetases, which may have been one of the key factors in the successful of emergence of ancient eukaryotic cells from bacterial colonies.

Plasmodium Apicoplast Gln-tRNAGln Biosynthesis Utilizes a Unique GatAB Amidotransferase Essential for Erythrocytic Stage Parasites

The malaria parasite Plasmodium falciparum apicoplast indirect aminoacylation pathway utilizes a non-discriminating glutamyl-tRNA synthetase to synthesize Glu-tRNA(Gln) and a glutaminyl-tRNA amidotransferase to convert Glu-tRNA(Gln) to Gln-tRNA(Gln). Here, we show that Plasmodium falciparum and other apicomplexans possess a unique heterodimeric glutamyl-tRNA amidotransferase consisting of GatA and GatB subunits (GatAB). We localized the P. falciparum GatA and GatB subunits to the apicoplast in blood stage parasites and demonstrated that recombinant GatAB converts Glu-tRNA(Gln) to Gln-tRNA(Gln) in vitro. We demonstrate that the apicoplast GatAB-catalyzed reaction is essential to the parasite blood stages because we could not delete the Plasmodium berghei gene encoding GatA in blood stage parasites in vivo. A phylogenetic analysis placed the split between Plasmodium GatB, archaeal GatE, and bacterial GatB prior to the phylogenetic divide between bacteria and archaea. Moreover, Plasmodium GatA also appears to have emerged prior to the bacterial-archaeal phylogenetic divide. Thus, although GatAB is found in Plasmodium, it emerged prior to the phylogenetic separation of archaea and bacteria.

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.

Preliminary X-ray crystallographic analysis of an engineered glutamyl-tRNA synthetase from Escherichia coli

The nature of interaction between glutamyl-tRNA synthetase (GluRS) and its tRNA substrate is unique in bacteria in that many bacterial GluRS are capable of recognizing two tRNA substrates: tRNAGlu and tRNAGln. To properly understand this distinctive GluRS-tRNA interaction it is important to pursue detailed structure-function studies; however, because of the fact that tRNA-GluRS interaction in bacteria is also associated with phylum-specific idiosyncrasies, the structure-function correlation studies must also be phylum-specific. GluRS from Thermus thermophilus and Escherichia coli, which belong to evolutionarily distant phyla, are the biochemically best characterized. Of these, only the structure of T. thermophilus GluRS is available. To fully unravel the subtleties of tRNAGlu-GluRS interaction in E. coli, a model bacterium that can also be pathogenic, determination of the E. coli GluRS structure is essential. However, previous attempts have failed to crystallize E. coli GluRS. By mapping crystal contacts of a homologous GluRS onto the E. coli GluRS sequence, two surface residues were identified that might have been hindering crystallization attempts. Accordingly, these two residues were mutated and crystallization of the double mutant was attempted. Here, the design, expression, purification and crystallization of an engineered E. coli GluRS in which two surface residues were mutated to optimize crystal contacts are reported.

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.

Trans-kingdom rescue of Gln-tRNAGln synthesis in yeast cytoplasm and mitochondria

Aminoacylation of transfer RNA^Gln (tRNA^Gln) is performed by distinct mechanisms in different kingdoms and represents the most diverged route of aminoacyl-tRNA synthesis found in nature. In Saccharomyces cerevisiae, cytosolic Gln-tRNA^Gln is generated by direct glutaminylation of tRNA^Gln by glutaminyl-tRNA synthetase (GlnRS), whereas mitochondrial Gln-tRNA^Gln is formed by an indirect pathway involving charging by a non-discriminating glutamyl-tRNA synthetase and the subsequent transamidation by a specific Glu-tRNA^Gln amidotransferase. Previous studies showed that fusion of a yeast non-specific tRNA-binding cofactor, Arc1p, to Escherichia coli GlnRS enables the bacterial enzyme to substitute for its yeast homologue in vivo. We report herein that the same fusion enzyme, upon being imported into mitochondria, substituted the indirect pathway for Gln-tRNA^Gln synthesis as well, despite significant differences in the identity determinants of E. coli and yeast cytosolic and mitochondrial tRNA^Gln isoacceptors. Fusion of Arc1p to the bacterial enzyme significantly enhanced its aminoacylation activity towards yeast tRNA^Gln isoacceptors in vitro. Our study provides a mechanism by which trans-kingdom rescue of distinct pathways of Gln-tRNA^Gln synthesis can be conferred by a single enzyme.

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.

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.

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.

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

The molecular basis of the genetic code relies on the specific ligation of amino acids to their cognate tRNA molecules. However, two pathways exist for the formation of Gln-tRNA^Gln. The evolutionarily older indirect route utilizes a non-discriminating glutamyl-tRNA synthetase (ND-GluRS) that can form both Glu-tRNA^Glu and Glu-tRNA^Gln. The Glu-tRNA^Gln is then converted to Gln-tRNA^Gln by an amidotransferase. Since the well-characterized bacterial ND-GluRS enzymes recognize tRNA^Glu and tRNA^Gln with an unrelated α-helical cage domain in contrast to the β-barrel anticodon-binding domain in archaeal and eukaryotic GluRSs, the mode of tRNA^Glu/tRNA^Gln discrimination in archaea and eukaryotes was unknown. Here, we present the crystal structure of the Methanothermobacter thermautotrophicus ND-GluRS, which is the evolutionary predecessor of both the glutaminyl-tRNA synthetase (GlnRS) and the eukaryotic discriminating GluRS. Comparison with the previously solved structure of the Escherichia coli GlnRS-tRNA^Gln complex reveals the structural determinants responsible for specific tRNA^Gln recognition by GlnRS compared to promiscuous recognition of both tRNAs by the ND-GluRS. The structure also shows the amino acid recognition pocket of GluRS is more variable than that found in GlnRS. Phylogenetic analysis is used to reconstruct the key events in the evolution from indirect to direct genetic encoding of glutamine.

Exit strategies for charged tRNA from GluRS

For several class I aminoacyl-tRNA synthetases (aaRSs), the rate-determining step in aminoacylation is the dissociation of charged tRNA from the enzyme. In this study, the following factors affecting the release of the charged tRNA from aaRSs are computationally explored: the protonation states of amino acids and substrates present in the active site, and the presence and the absence of AMP and elongation factor Tu. Through molecular modeling, internal pK(a) calculations, and molecular dynamics simulations, distinct, mechanistically relevant post-transfer states with charged tRNA bound to glutamyl-tRNA synthetase from Thermus thermophilus (Glu-tRNA(Glu)) are considered. The behavior of these nonequilibrium states is characterized as a function of time using dynamical network analysis, local energetics, and changes in free energies to estimate transitions that occur during the release of the tRNA. The hundreds of nanoseconds of simulation time reveal system characteristics that are consistent with recent experimental studies. Energetic and network results support the previously proposed mechanism in which the transfer of amino acid to tRNA is accompanied by the protonation of AMP to H-AMP. Subsequent migration of proton to water reduces the stability of the complex and loosens the interface both in the presence and in the absence of AMP. The subsequent undocking of AMP or tRNA then proceeds along thermodynamically competitive pathways. Release of the tRNA acceptor stem is further accelerated by the deprotonation of the alpha-ammonium group on the charging amino acid. The proposed general base is Glu41, a residue binding the alpha-ammonium group that is conserved in both structure and sequence across nearly all class I aaRSs. This universal handle is predicted through pK(a) calculations to be part of a proton relay system for destabilizing the bound charging amino acid following aminoacylation. Addition of elongation factor Tu to the aaRS.tRNA complex stimulates the dissociation of the tRNA core and the tRNA acceptor stem.

Crystal structure of Pyrococcus horikoshii tryptophanyl-tRNA synthetase and structure-based phylogenetic analysis suggest an archaeal origin of tryptophanyl-tRNA synthetase

The ancient and ubiquitous aminoacyl-tRNA synthetases constitute a valuable model system for studying early evolutionary events. So far, the evolutionary relationship of tryptophanyl- and tyrosyl-tRNA synthetase (TrpRS and TyrRS) remains controversial. As TrpRS and TyrRS share low sequence homology but high structural similarity, a structure-based method would be advantageous for phylogenetic analysis of the enzymes. Here, we present the first crystal structure of an archaeal TrpRS, the structure of Pyrococcus horikoshii TrpRS (pTrpRS) in complex with tryptophanyl-5′ AMP (TrpAMP) at 3.0 Å resolution which demonstrates more similarities to its eukaryotic counterparts. With the pTrpRS structure, we perform a more complete structure-based phylogenetic study of TrpRS and TyrRS, which for the first time includes representatives from all three domains of life. Individually, each enzyme shows a similar evolutionary profile as observed in the sequence-based phylogenetic studies. However, TyrRSs from Archaea/Eucarya cluster with TrpRSs rather than their bacterial counterparts, and the root of TrpRS locates in the archaeal branch of TyrRS, indicating the archaeal origin of TrpRS. Moreover, the short distance between TrpRS and archaeal TyrRS and that between bacterial and archaeal TrpRS, together with the wide distribution of TrpRS, suggest that the emergence of TrpRS and subsequent acquisition by Bacteria occurred at early stages of evolution.

Recognition of tRNAGln by Helicobacter pylori GluRS2--a tRNAGln-specific glutamyl-tRNA synthetase

Accurate aminoacylation of tRNAs by the aminoacyl-tRNA synthetases (aaRSs) plays a critical role in protein translation. However, some of the aaRSs are missing in many microorganisms. Helicobacter pylori does not have a glutaminyl-tRNA synthetase (GlnRS) but has two divergent glutamyl-tRNA synthetases: GluRS1 and GluRS2. Like a canonical GluRS, GluRS1 aminoacylates tRNA^Glu1 and tRNA^Glu2. In contrast, GluRS2 only misacylates tRNA^Gln to form Glu-tRNA^Gln. It is not clear how GluRS2 achieves specific recognition of tRNA^Gln while rejecting the two H. pylori tRNA^Glu isoacceptors. Here, we show that GluRS2 recognizes major identity elements clustered in the tRNA^Gln acceptor stem. Mutations in the tRNA anticodon or at the discriminator base had little to no impact on enzyme specificity and activity.

The role of the catalytic domain of E. coli GluRS in tRNAGln discrimination

Discrimination of tRNA(Gln) is an integral function of several bacterial glutamyl-tRNA synthetases (GluRS). The origin of the discrimination is thought to arise from unfavorable interactions between tRNA(Gln) and the anticodon-binding domain of GluRS. From experiments on an anticodon-binding domain truncated Escherichia coli (E. coli) GluRS (catalytic domain) and a chimeric protein, constructed from the catalytic domain of E. coli GluRS and the anticodon-binding domain of E. coli glutaminyl-tRNA synthetase (GlnRS), we show that both proteins discriminate against E. coli tRNA(Gln). Our results demonstrate that in addition to the anticodon-binding domain, tRNA(Gln) discriminatory elements may be present in the catalytic domain in E. coli GluRS as well.

A chimaeric glutamyl:glutaminyl-tRNA synthetase: implications for evolution

aaRSs (aminoacyl-tRNA synthetases) are multi-domain proteins that have evolved by domain acquisition. The anti-codon binding domain was added to the more ancient catalytic domain during aaRS evolution. Unlike in eukaryotes, the anti-codon binding domains of GluRS (glutamyl-tRNA synthetase) and GlnRS (glutaminyl-tRNA synthetase) in bacteria are structurally distinct. This originates from the unique evolutionary history of GlnRSs. Starting from the catalytic domain, eukaryotic GluRS evolved by acquiring the archaea/eukaryote-specific anti-codon binding domain after branching away from the eubacteria family. Subsequently, eukaryotic GlnRS evolved from GluRS by gene duplication and horizontally transferred to bacteria. In order to study the properties of the putative ancestral GluRS in eukaryotes, formed immediately after acquiring the anti-codon binding domain, we have designed and constructed a chimaeric protein, cGluGlnRS, consisting of the catalytic domain, Ec GluRS (Escherichia coli GluRS), and the anti-codon binding domain of EcGlnRS (E. coli GlnRS). In contrast to the isolated EcN-GluRS, cGluGlnRS showed detectable activity of glutamylation of E. coli tRNAglu and was capable of complementing an E. coli ts (temperature-sensitive)-GluRS strain at non-permissive temperatures. Both cGluGlnRS and EcN-GluRS were found to bind E. coli tRNAglu with native EcGluRS-like affinity, suggesting that the anticodon-binding domain in cGluGlnRS enhances kcat for glutamylation. This was further confirmed from similar experiments with a chimaera between EcN-GluRS and the substrate-binding domain of EcDnaK (E. coli DnaK). We also show that an extended loop, present in the anticodon-binding domains of GlnRSs, is absent in archaeal GluRS, suggesting that the loop was a later addition, generating additional anti-codon discrimination capability in GlnRS as it evolved from GluRS in eukaryotes.

Amino acid modifications on tRNA

The accurate formation of cognate aminoacyl-transfer RNAs (aa-tRNAs) is essential for the fidelity of translation. Most amino acids are esterified onto their cognate tRNA isoacceptors directly by aa-tRNA synthetases. However, in the case of four amino acids (Gln, Asn, Cys and Sec), aminoacyl-tRNAs are made through indirect pathways in many organisms across all three domains of life. The process begins with the charging of noncognate amino acids to tRNAs by a specialized synthetase in the case of Cys-tRNA(Cys) formation or by synthetases with relaxed specificity, such as the non-discriminating glutamyl-tRNA, non-discriminating aspartyl-tRNA and seryl-tRNA synthetases. The resulting misacylated tRNAs are then converted to cognate pairs through transformation of the amino acids on the tRNA, which is catalyzed by a group of tRNA-dependent modifying enzymes, such as tRNA-dependent amidotransferases, Sep-tRNA:Cys-tRNA synthase, O-phosphoseryl-tRNA kinase and Sep-tRNA:Sec-tRNA synthase. The majority of these indirect pathways are widely spread in all domains of life and thought to be part of the evolutionary process.

From one amino acid to another: tRNA-dependent amino acid biosynthesis

Aminoacyl-tRNAs (aa-tRNAs) are the essential substrates for translation. Most aa-tRNAs are formed by direct aminoacylation of tRNA catalyzed by aminoacyl-tRNA synthetases. However, a smaller number of aa-tRNAs (Asn-tRNA, Gln-tRNA, Cys-tRNA and Sec-tRNA) are made by synthesizing the amino acid on the tRNA by first attaching a non-cognate amino acid to the tRNA, which is then converted to the cognate one catalyzed by tRNA-dependent modifying enzymes. Asn-tRNA or Gln-tRNA formation in most prokaryotes requires amidation of Asp-tRNA or Glu-tRNA by amidotransferases that couple an amidase or an asparaginase to liberate ammonia with a tRNA-dependent kinase. Both archaeal and eukaryotic Sec-tRNA biosynthesis and Cys-tRNA synthesis in methanogens require O-phosophoseryl-tRNA formation. For tRNA-dependent Cys biosynthesis, O-phosphoseryl-tRNA synthetase directly attaches the amino acid to the tRNA which is then converted to Cys by Sep-tRNA: Cys-tRNA synthase. In Sec-tRNA synthesis, O-phosphoseryl-tRNA kinase phosphorylates Ser-tRNA to form the intermediate which is then modified to Sec-tRNA by Sep-tRNA:Sec-tRNA synthase. Complex formation between enzymes in the same pathway may protect the fidelity of protein synthesis. How these tRNA-dependent amino acid biosynthetic routes are integrated into overall metabolism may explain why they are still retained in so many organisms.

On the evolution of the tRNA-dependent amidotransferases, GatCAB and GatDE

Glutaminyl-tRNA synthetase and asparaginyl-tRNA synthetase evolved from glutamyl-tRNA synthetase and aspartyl-tRNA synthetase, respectively, after the split in the last universal communal ancestor (LUCA). Glutaminyl-tRNA(Gln) and asparaginyl-tRNA(Asn) were likely formed in LUCA by amidation of the mischarged species, glutamyl-tRNA(Gln) and aspartyl-tRNA(Asn), by tRNA-dependent amidotransferases, as is still the case in most bacteria and all known archaea. The amidotransferase GatCAB is found in both domains of life, while the heterodimeric amidotransferase GatDE is found only in Archaea. The GatB and GatE subunits belong to a unique protein family that includes Pet112 that is encoded in the nuclear genomes of numerous eukaryotes. GatE was thought to have evolved from GatB after the emergence of the modern lines of decent. Our phylogenetic analysis though places the split between GatE and GatB, prior to the phylogenetic divide between Bacteria and Archaea, and Pet112 to be of mitochondrial origin. In addition, GatD appears to have emerged prior to the bacterial-archaeal phylogenetic divide. Thus, while GatDE is an archaeal signature protein, it likely was present in LUCA together with GatCAB. Archaea retained both amidotransferases, while Bacteria emerged with only GatCAB. The presence of GatDE has favored a unique archaeal tRNA(Gln) that may be preventing the acquisition of glutaminyl-tRNA synthetase in Archaea. Archaeal GatCAB, on the other hand, has not favored a distinct tRNA(Asn), suggesting that tRNA(Asn) recognition is not a major barrier to the retention of asparaginyl-tRNA synthetase in many Archaea.

Deinococcus glutaminyl-tRNA synthetase is a chimer between proteins from an ancient and the modern pathways of aminoacyl-tRNA formation

Glutaminyl-tRNA synthetase from Deinococcus radiodurans possesses a C-terminal extension of 215 residues appending the anticodon-binding domain. This domain constitutes a paralog of the Yqey protein present in various organisms and part of it is present in the C-terminal end of the GatB subunit of GatCAB, a partner of the indirect pathway of Gln-tRNA^Gln formation. To analyze the peculiarities of the structure–function relationship of this GlnRS related to the Yqey domain, a structure of the protein was solved from crystals diffracting at 2.3 Å and a docking model of the synthetase complexed to tRNA^Gln constructed. The comparison of the modeled complex with the structure of the E. coli complex reveals that all residues of E. coli GlnRS contacting tRNA^Gln are conserved in D. radiodurans GlnRS, leaving the functional role of the Yqey domain puzzling. Kinetic investigations and tRNA-binding experiments of full length and Yqey-truncated GlnRSs reveal that the Yqey domain is involved in tRNA^Gln recognition. They demonstrate that Yqey plays the role of an affinity-enhancer of GlnRS for tRNA^Gln acting only in cis. However, the presence of Yqey in free state in organisms lacking GlnRS, suggests that this domain may exert additional cellular functions.

Novel tRNA aminoacylation mechanisms

In nature, ribosomally synthesized proteins can contain at least 22 different amino acids: the 20 common amino acids as well as selenocysteine and pyrrolysine. Each of these amino acids is inserted into proteins codon-specifically via an aminoacyl-transfer RNA (aa-tRNA). In most cases, these aa-tRNAs are biosynthesized directly by a set of highly specific and accurate aminoacyl-tRNA synthetases (aaRSs). However, in some cases aaRSs with relaxed or novel substrate specificities cooperate with other enzymes to generate specific canonical and non-canonical aminoacyl-tRNAs.

Steroids, triterpenoids and molecular oxygen

There is a close connection between modern-day biosynthesis of particular triterpenoid biomarkers and presence of molecular oxygen in the environment. Thus, the detection of steroid and triterpenoid hydrocarbons far back in Earth history has been used to infer the antiquity of oxygenic photosynthesis. This prompts the question: were these compounds produced similarly in the past? In this paper, we address this question with a review of the current state of knowledge surrounding the oxygen requirement for steroid biosynthesis and phylogenetic patterns in the distribution of steroid and triterpenoid biosynthetic pathways. The hopanoid and steroid biosynthetic pathways are very highly conserved within the bacterial and eukaryotic domains, respectively. Bacteriohopanepolyols are produced by a wide range of bacteria, and are methylated in significant abundance at the C2 position by oxygen-producing cyanobacteria. On the other hand, sterol biosynthesis is sparsely distributed in distantly related bacterial taxa and the pathways do not produce the wide range of products that characterize eukaryotes. In particular, evidence for sterol biosynthesis by cyanobacteria appears flawed. Our experiments show that cyanobacterial cultures are easily contaminated by sterol-producing rust fungi, which can be eliminated by treatment with cycloheximide affording sterol-free samples. Sterols are ubiquitous features of eukaryotic membranes, and it appears likely that the initial steps in sterol biosynthesis were present in their modern form in the last common ancestor of eukaryotes. Eleven molecules of O2 are required by four enzymes to produce one molecule of cholesterol. Thermodynamic arguments, optimization of function and parsimony all indicate that an ancestral anaerobic pathway is highly unlikely. The known geological record of molecular fossils, especially steranes and triterpanes, is notable for the limited number of structural motifs that have been observed. With a few exceptions, the carbon skeletons are the same as those found in the lipids of extant organisms and no demonstrably extinct structures have been reported. Furthermore, their patterns of occurrence over billion year time-scales correlate strongly with environments of deposition. Accordingly, biomarkers are excellent indicators of environmental conditions even though the taxonomic affinities of all biomarkers cannot be precisely specified. Biomarkers are ultimately tied to biochemicals with very specific functional properties, and interpretations of the biomarker record will benefit from increased understanding of the biological roles of geologically durable molecules.

Structural basis of RNA-dependent recruitment of glutamine to the genetic code

Glutaminyl-transfer RNA (Gln-tRNA(Gln)) in archaea is synthesized in a pretranslational amidation of misacylated Glu-tRNA(Gln) by the heterodimeric Glu-tRNA(Gln) amidotransferase GatDE. Here we report the crystal structure of the Methanothermobacter thermautotrophicus GatDE complexed to tRNA(Gln) at 3.15 angstroms resolution. Biochemical analysis of GatDE and of tRNA(Gln) mutants characterized the catalytic centers for the enzyme's three reactions (glutaminase, kinase, and amidotransferase activity). A 40 angstrom-long channel for ammonia transport connects the active sites in GatD and GatE. tRNA(Gln) recognition by indirect readout based on shape complementarity of the D loop suggests an early anticodon-independent RNA-based mechanism for adding glutamine to the genetic code.

Crystal structure of a non-discriminating glutamyl-tRNA synthetase

Error-free protein biosynthesis is dependent on the reliable charging of each tRNA with its cognate amino acid. Many bacteria, however, lack a glutaminyl-tRNA synthetase. In these organisms, tRNA(Gln) is initially mischarged with glutamate by a non-discriminating glutamyl-tRNA synthetase (ND-GluRS). This enzyme thus charges both tRNA(Glu) and tRNA(Gln) with glutamate. Discriminating GluRS (D-GluRS), found in some bacteria and all eukaryotes, exclusively generates Glu-tRNA(Glu). Here we present the first crystal structure of a non-discriminating GluRS from Thermosynechococcus elongatus (ND-GluRS(Tel)) in complex with glutamate at a resolution of 2.45 A. Structurally, the enzyme shares the overall architecture of the discriminating GluRS from Thermus thermophilus (D-GluRS(Tth)). We confirm experimentally that GluRS(Tel) is non-discriminating and present kinetic parameters for synthesis of Glu-tRNA(Glu) and of Glu-tRNA(Gln). Anticodons of tRNA(Glu) (34C/UUC36) and tRNA(Gln) (34C/UUG36) differ only in base 36. The pyrimidine base of C36 is specifically recognized in D-GluRS(Tth) by the residue Arg358. In ND-GluRS(Tel) this arginine residue is replaced by glycine (Gly366) presumably allowing both cytosine and the bulkier purine base G36 of tRNA(Gln) to be tolerated. Most other ND-GluRS share this structural feature, leading to relaxed substrate specificity.

Structural basis of the water-assisted asparagine recognition by asparaginyl-tRNA synthetase

Asparaginyl-tRNA synthetase (AsnRS) is a member of the class-II aminoacyl-tRNA synthetases, and is responsible for catalyzing the specific aminoacylation of tRNA(Asn) with asparagine. Here, the crystal structure of AsnRS from Pyrococcus horikoshii, complexed with asparaginyl-adenylate (Asn-AMP), was determined at 1.45 A resolution, and those of free AsnRS and AsnRS complexed with an Asn-AMP analog (Asn-SA) were solved at 1.98 and 1.80 A resolutions, respectively. All of the crystal structures have many solvent molecules, which form a network of hydrogen-bonding interactions that surrounds the entire AsnRS molecule. In the AsnRS/Asn-AMP complex (or the AsnRS/Asn-SA), one side of the bound Asn-AMP (or Asn-SA) is completely covered by the solvent molecules, which complement the binding site. In particular, two of these water molecules were found to interact directly with the asparagine amide and carbonyl groups, respectively, and to contribute to the formation of a pocket highly complementary to the asparagine side-chain. Thus, these two water molecules appear to play a key role in the strict recognition of asparagine and the discrimination against aspartic acid by the AsnRS. This water-assisted asparagine recognition by the AsnRS strikingly contrasts with the fact that the aspartic acid recognition by the closely related aspartyl-tRNA synthetase is achieved exclusively through extensive interactions with protein amino acid residues. Furthermore, based on a docking model of AsnRS and tRNA, a single arginine residue (Arg83) in the AsnRS was postulated to be involved in the recognition of the third position of the tRNA(Asn) anticodon (U36). We performed a mutational analysis of this particular arginine residue, and confirmed its significance in the tRNA recognition.

COCO-CL: hierarchical clustering of homology relations based on evolutionary correlations

MOTIVATION

Determining orthology relations among genes across multiple genomes is an important problem in the post-genomic era. Identifying orthologous genes can not only help predict functional annotations for newly sequenced or poorly characterized genomes, but can also help predict new protein-protein interactions. Unfortunately, determining orthology relation through computational methods is not straightforward due to the presence of paralogs. Traditional approaches have relied on pairwise sequence comparisons to construct graphs, which were then partitioned into putative clusters of orthologous groups. These methods do not attempt to preserve the non-transitivity and hierarchic nature of the orthology relation.

RESULTS

We propose a new method, COCO-CL, for hierarchical clustering of homology relations and identification of orthologous groups of genes. Unlike previous approaches, which are based on pairwise sequence comparisons, our method explores the correlation of evolutionary histories of individual genes in a more global context. COCO-CL can be used as a semi-independent method to delineate the orthology/paralogy relation for a refined set of homologous proteins obtained using a less-conservative clustering approach, or as a refiner that removes putative out-paralogs from clusters computed using a more inclusive approach. We analyze our clustering results manually, with support from literature and functional annotations. Since our orthology determination procedure does not employ a species tree to infer duplication events, it can be used in situations when the species tree is unknown or uncertain.

CONTACT

jothi@mail.nih.gov, przytyck@mail.nih.gov

SUPPLEMENTARY INFORMATION

Supplementary materials are available at Bioinformatics online.

Reduced stress tolerance of glutamine-deprived human monocytic cells is associated with selective down-regulation of Hsp70 by decreased mRNA stability

In critically ill patients, clinicians observe a reverse correlation of survival and a decreased plasma concentration of the most abundant free amino acid, glutamine (Gln). However, in this context, the role of Gln remains largely elusive. Gln is used as an energy substrate by monocytes. Gln deprivation of these cells results in an increased susceptibility to cell stress and apoptosis, as well as in a reduced responsiveness to pro-inflammatory stimuli. We performed a systematic study to elucidate the molecular mechanism by which Gln depletion affects the heat stress response of the monocytic cell line U937. Proteomic analysis revealed that Gln depletion was associated with specific changes in the protein expression pattern. However, the overall level of tRNA-bound Gln remained unaffected. The stress protein heat shock protein (Hsp) 70 showed the highest reduction in protein synthesis. This was due to enhanced mRNA decay during Gln starvation while the transcriptional and the translational control of Hsp70 expression remained unchanged. A physiological Gln concentration and above was found to be necessary for maximum Hsp70 accumulation upon heat shock. Thus, the study shows a specific link between Gln metabolism and the regulation of heat shock proteins.

Glu-Q-tRNAAsp synthetase coded by the yadB gene, a new paralog of aminoacyl-tRNA synthetase that glutamylates tRNAAsp anticodon

Analysis of the completed genome sequences revealed presence in various bacteria of an open reading frame (ORF) encoding a polypeptide chain presenting important similarities with the catalytic domain of glutamyl-tRNA synthetases but deprived of the C-terminal anticodon-binding domain. This paralog of glutamyl-tRNA synthetases, the YadB protein, activates glutamate in the absence of tRNA and transfers the activated glutamate not on tRNA(Glu) but instead on tRNA(Asp). It has been shown that tRNA(Asp) is able to accept two amino acids: aspartate charged by aspartyl-tRNA synthetase and glutamate charged by YadB. The functional properties of YadB contrast with those of the canonical glutamyl-tRNA synthetases, which activate Glu only in presence of the cognate tRNA before aminoacylation of the 3'-end of tRNA. Biochemical approaches and mass spectrometry investigations revealed that YadB transfers the activated glutamate on the cyclopenthene-diol ring of the modified nucleoside queuosine posttranscriptionally inserted at the wobble position of the anticodon-loop to form glutamyl-queuosine. Unstability of the ester bond between the glutamate residue and the cyclopenthene-diol (half-life 7.5 min) explains why until now this modification escaped detection. Among Escherichia coli tRNAs containing queuosine in the wobble position, only tRNA(Asp) is substrate of YadB. Sequence comparison reveals a structural mimicry between the anticodon-stem and loop of tRNA(Asp) and the amino acid acceptor-stem of tRNA(Glu). YadB, renamed glutamyl-Q-tRNA(Asp) synthetase, constitutes the first enzyme structurally related to aminoacyl-tRNA synthetases which catalyzes a hypermodification in tRNA, and whose function seems to be conserved among prokaryotes. The discovery of glutamyl-Q-tRNA(Asp) synthetase breaks down the current paradigm according to which the catalytic domain of aminoacyl-tRNA synthetases recognizes the amino acid acceptor-stem of tRNA and aminoacylates the 3'-terminal ribose. The evolutionary significance of the existence of an aminoacyl-tRNA synthetase paralog dedicated to the hypermodification of a tRNA anticodon will be discussed.

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

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

Saccharomyces cerevisiae imports the cytosolic pathway for Gln-tRNA synthesis into the mitochondrion

Aminoacyl-tRNA (aa-tRNA) formation, an essential process in protein biosynthesis, is generally achieved by direct attachment of an amino acid to tRNA by the aa-tRNA synthetases. An exception is Gln-tRNA synthesis, which in eukaryotes is catalyzed by glutaminyl-tRNA synthetase (GlnRS), while most bacteria, archaea, and chloroplasts employ the transamidation pathway, in which a tRNA-dependent glutamate modification generates Gln-tRNA. Mitochondrial protein synthesis is carried out normally by mitochondrial enzymes and organelle-encoded tRNAs that are different from their cytoplasmic counterparts. Early work suggested that mitochondria use the transamidation pathway for Gln-tRNA formation. We found no biochemical support for this in Saccharomyces cerevisiae mitochondria, but demonstrated the presence of the cytoplasmic GlnRS in the organelle and its involvement in mitochondrial Gln-tRNA synthesis. In addition, we showed in vivo localization of cytoplasmic tRNAGln in mitochondria and demonstrated its role in mitochondrial translation. We furthermore reconstituted in vitro cytoplasmic tRNAGln import into mitochondria by a novel mechanism. This tRNA import mechanism expands our knowledge of RNA trafficking in the eukaryotic cell. These findings change our view of the evolution of organellar protein synthesis.

Divergent anticodon recognition in contrasting glutamyl-tRNA synthetases

The pathogenic bacterium Helicobacter pylori utilizes two essential glutamyl-tRNA synthetases (GluRS1 and GluRS2). These two enzymes are closely related in evolution and yet they aminoacylate contrasting tRNAs. GluRS1 is a canonical discriminating GluRS (D-GluRS) that biosynthesizes Glu-tRNA(Glu) and cannot make Glu-tRNA(Gln). In contrast, GluRS2 is non-canonical as it is only essential for the production of misacylated Glu-tRNA(Gln). The co-existence and evident divergence of these two enzymes was capitalized upon to directly examine how GluRS2 acquired tRNA(Gln) specificity. One key feature that distinguishes tRNA(Glu) from tRNA(Gln) is the third position in the anticodon of each tRNA (C36 versus G36, respectively). By comparing sequence alignments of different GluRSs, including GluRS1s and GluRS2s, to the crystal structure of the Thermus thermophilus D-GluRS:tRNA(Glu) complex, a divergent pattern of conservation in enzymes that aminoacylate tRNA(Glu)versus those specific for tRNA(Gln) emerged and was experimentally validated. In particular, when an arginine conserved in discriminating GluRSs and GluRS1s was inserted into Hp GluRS2 (Glu334Arg GluRS2), the catalytic efficiency of the mutant enzyme (k(cat)/K(Mapp)) was reduced by approximately one order of magnitude towards tRNA(Gln). However, this mutation did not introduce activity towards tRNA(Glu). In contrast, disruption of a glycine that is conserved in all GluRS2s but not in other GluRSs (Gly417Thr GluRS2) generated a mutant GluRS2 with weak activity towards tRNA(Glu1). Synergy between these two mutations was observed in the double mutant (Glu334Arg/Gly417Thr GluRS2), which specifically and more robustly aminoacylates tRNA(Glu1) instead of tRNA(Gln). As GluRS1 and GluRS2 are related by an apparent gene duplication event, these results demonstrate that we can experimentally map critical evolutionary events in the emergence of new tRNA specificities.