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

Non-canonical functions of human cytoplasmic tyrosyl-, tryptophanyl- and other aminoacyl-tRNA synthetases

Aminoacyl-tRNA synthetases catalyze the aminoacylation of their cognate tRNAs. Here we review the accumulated knowledge of non-canonical functions of human cytoplasmic aminoacyl-tRNA synthetases, especially tyrosyl- (TyrRS) and tryptophanyl-tRNA synthetase (TrpRS). Human TyrRS and TrpRS have an extra domain. Two distinct cytokines, i.e., the core catalytic "mini TyrRS" and the extra C-domain, are generated from human TyrRS by proteolytic cleavage. Moreover, the core catalytic domains of human TyrRS and TrpRS function as angiogenic and angiostatic factors, respectively, whereas the full-length forms are inactive for this function. It is also known that many synthetases change their localization in response to a specific signal and subsequently exhibit alternative functions. Furthermore, some synthetases function as sensors for amino acids by changing their protein interactions in an amino acid-dependent manner. Further studies will be necessary to elucidate regulatory mechanisms of non-canonical functions of aminoacyl-tRNA synthetases in particular, by analyzing the effect of their post-translational modifications.

Tryptophanyl-tRNA synthetase mediates high-affinity tryptophan uptake into human cells

The tryptophan (Trp) transport system has a high affinity and selectivity toward Trp, and has been reported to exist in both human and mouse macrophages. Although this system is highly expressed in interferon-γ (IFN-γ)-treated cells and indoleamine 2,3-dioxygenase 1 (IDO1)-expressing cells, its identity remains incompletely understood. Tryptophanyl-tRNA synthetase (TrpRS) is also highly expressed in IFN-γ-treated cells and also has high affinity and selectivity for Trp. Here, we investigated the effects of human TrpRS expression on Trp uptake into IFN-γ-treated human THP-1 monocytes or HeLa cells. Inhibition of human TrpRS expression by TrpRS-specific siRNAs decreased and overexpression of TrpRS increased Trp uptake into the cells. Of note, the TrpRS-mediated uptake system had more than hundred-fold higher affinity for Trp than the known System L amino acid transporter, promoted uptake of low Trp concentrations, and had very high Trp selectivity. Moreover, site-directed mutagenesis experiments indicated that Trp- and ATP-binding sites, but not tRNA-binding sites, in TrpRS are essential for TrpRS-mediated Trp uptake into the human cells. We further demonstrate that the addition of purified TrpRS to cell culture medium increases Trp uptake into cells. Taken together, our results reveal that TrpRS plays an important role in high-affinity Trp uptake into human cells.

Ancestral Reconstruction of a Pre-LUCA Aminoacyl-tRNA Synthetase Ancestor Supports the Late Addition of Trp to the Genetic Code

The genetic code was likely complete in its current form by the time of the last universal common ancestor (LUCA). Several scenarios have been proposed for explaining the code's pre-LUCA emergence and expansion, and the relative order of the appearance of amino acids used in translation. One co-evolutionary model of genetic code expansion proposes that at least some amino acids were added to the code by the ancient divergence of aminoacyl-tRNA synthetase (aaRS) families. Of all the amino acids used within the genetic code, Trp is most frequently claimed as a relatively recent addition. We observe that, since TrpRS and TyrRS are paralogous protein families retaining significant sequence similarity, the inferred sequence composition of their ancestor can be used to evaluate this co-evolutionary model of genetic code expansion. We show that ancestral sequence reconstructions of the pre-LUCA paralog ancestor of TyrRS and TrpRS have several sites containing Tyr, yet a complete absence of sites containing Trp. This is consistent with the paralog ancestor being specific for the utilization of Tyr, with Trp being a subsequent addition to the genetic code facilitated by a process of aaRS divergence and neofunctionalization. Only after this divergence could Trp be specifically encoded and incorporated into proteins, including the TyrRS and TrpRS descendant lineages themselves. This early absence of Trp is observed under both homogeneous and non-homogeneous models of ancestral sequence reconstruction. Simulations support that this observed absence of Trp is unlikely to be due to chance or model bias. These results support that the final stages of genetic code evolution occurred well within the "protein world," and that the presence-absence of Trp within conserved sites of ancient protein domains is a likely measure of their relative antiquity, permitting the relative timing of extremely early events within protein evolution before LUCA.

VEGF, not VEGFR2, is associated with the angiogenesis effect of mini-TyrRS/mini-TrpRS in human umbilical vein endothelial cells in hypoxia

The purpose of this study was to determine the relationship between VEGF and mini-TyrRS/mini-TrpRS in angiogenesis in hypoxic culture and to begin to comprehend their mechanism in angiogenesis. We designed a VEGF gene silencing assay by using lentivirus vectors, and then western blotting was used to determine the protein expression of VEGF, VEGFR2 and pVEGFR2 in three groups in hypoxic culture at 3, 6, 12, or 24 h: (1) untransfected human umbilical vein endothelial cells (HUVECs) (Control); (2) pGCSIL-GFP lentivirus vector-transduced HUVECs (Mock); and (3) pGCSIL-shVEGF lentivirus vector-transduced HUVECs (Experimental). We also detected the effects of mini-TyrRS/mini-TrpRS peptides on HUVEC proliferation, migration and tube formation after lentivirus vector transfection and VEGFR2 antibody injection. The results indicated that expression of the mini-TyrRS protein was increased, whereas that of mini-TrpRS was specifically decreased in hypoxic culture both in control and mock groups. However, this trend in protein levels of mini-TyrRS and mini-TrpRS was lost in the experimental group after transduction with the pGCSIL-shVEGF lentivirus vector. The protein expression of VEGF was increased in hypoxic culture both in control and mock groups. After transduction with the pGCSIL-shVEGF lentivirus vector, the protein level of VEGF was noticeably decreased in the experimental group; however, for VEGFR2, the results showed no significant difference in VEGFR2 protein expression in any of the groups. For pVEGFR2, we found a distinct trend from that seen with VEGF. The protein expression of pVEGFR2 was sharply increased in hypoxic culture in the three groups. The addition of mini-TyrRS significantly promoted proliferation, migration and tube formation of HUVECs, while mini-TrpRS inhibited these processes in both control and mock groups in hypoxic culture. However, these effects disappeared after transduction with the pGCSIL-shVEGF lentivirus vector in the experimental group, but no significant difference was observed after VEGFR2 antibody injection. The protein expression of VEGF is similar to that of mini-TyrRS in hypoxic culture and plays an important role in the mini-TyrRS/mini-TrpRS-stimulated proliferation, migration and tube formation of HUVECs in hypoxia. These results also suggest that the change in mini-TyrRS and mini-TrpRS expression in hypoxic culture is not related to VEGFR2 and that some other possible mechanisms, are involved in the phosphorylation of VEGFR2.

Statistical evaluation of the Rodin-Ohno hypothesis: sense/antisense coding of ancestral class I and II aminoacyl-tRNA synthetases

We tested the idea that ancestral class I and II aminoacyl-tRNA synthetases arose on opposite strands of the same gene. We assembled excerpted 94-residue Urgenes for class I tryptophanyl-tRNA synthetase (TrpRS) and class II Histidyl-tRNA synthetase (HisRS) from a diverse group of species, by identifying and catenating three blocks coding for secondary structures that position the most highly conserved, active-site residues. The codon middle-base pairing frequency was 0.35 ± 0.0002 in all-by-all sense/antisense alignments for 211 TrpRS and 207 HisRS sequences, compared with frequencies between 0.22 ± 0.0009 and 0.27 ± 0.0005 for eight different representations of the null hypothesis. Clustering algorithms demonstrate further that profiles of middle-base pairing in the synthetase antisense alignments are correlated along the sequences from one species-pair to another, whereas this is not the case for similar operations on sets representing the null hypothesis. Most probable reconstructed sequences for ancestral nodes of maximum likelihood trees show that middle-base pairing frequency increases to approximately 0.42 ± 0.002 as bacterial trees approach their roots; ancestral nodes from trees including archaeal sequences show a less pronounced increase. Thus, contemporary and reconstructed sequences all validate important bioinformatic predictions based on descent from opposite strands of the same ancestral gene. They further provide novel evidence for the hypothesis that bacteria lie closer than archaea to the origin of translation. Moreover, the inverse polarity of genetic coding, together with a priori α-helix propensities suggest that in-frame coding on opposite strands leads to similar secondary structures with opposite polarity, as observed in TrpRS and HisRS crystal structures.

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.

The double-length tyrosyl-tRNA synthetase from the eukaryote Leishmania major forms an intrinsically asymmetric pseudo-dimer

The single tyrosyl-tRNA synthetase (TyrRS) gene in trypanosomatid genomes codes for a protein that is twice the length of TyrRS from virtually all other organisms. Each half of the double-length TyrRS contains a catalytic domain and an anticodon-binding domain; however, the two halves retain only 17% sequence identity to each other. The structural and functional consequences of this duplication and divergence are unclear. TyrRS normally forms a homodimer in which the active site of one monomer pairs with the anticodon-binding domain from the other. However, crystal structures of Leishmania major TyrRS show that, instead, the two halves of a single molecule form a pseudo-dimer resembling the canonical TyrRS dimer. Curiously, the C-terminal copy of the catalytic domain has lost the catalytically important HIGH and KMSKS motifs characteristic of class I aminoacyl-tRNA synthetases. Thus, the pseudo-dimer contains only one functional active site (contributed by the N-terminal half) and only one functional anticodon recognition site (contributed by the C-terminal half). Despite biochemical evidence for negative cooperativity between the two active sites of the usual TyrRS homodimer, previous structures have captured a crystallographically-imposed symmetric state. As the L. major TyrRS pseudo-dimer is inherently asymmetric, conformational variations observed near the active site may be relevant to understanding how the state of a single active site is communicated across the dimer interface. Furthermore, substantial differences between trypanosomal TyrRS and human homologs are promising for the design of inhibitors that selectively target the parasite enzyme.

Deep phylogeny--how a tree can help characterize early life on Earth

The Darwinian concept of biological evolution assumes that life on Earth shares a common ancestor. The diversification of this common ancestor through speciation events and vertical transmission of genetic material implies that the classification of life can be illustrated in a tree-like manner, commonly referred to as the Tree of Life. This article describes features of the Tree of Life, such as how the tree has been both pruned and become bushier throughout the past century as our knowledge of biology has expanded. We present current views that the classification of life may be best illustrated as a ring or even a coral with tree-like characteristics. This article also discusses how the organization of the Tree of Life offers clues about ancient life on Earth. In particular, we focus on the environmental conditions and temperature history of Precambrian life and show how chemical, biological, and geological data can converge to better understand this history."You know, a tree is a tree. How many more do you need to look at?"--Ronald Reagan (Governor of California), quoted in the Sacramento Bee, opposing expansion of Redwood National Park, March 3, 1966.

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.

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

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

The tree of life might be rooted in the branch leading to Nanoarchaeota

It is suggested that the tree of life might be rooted in the domain of the Archaea, in the branch leading to the phylum of Nanoarchaeota. This hypothesis seems to be corroborated by the uniqueness and ancestrality of some traits possessed by Nanoarchaeum equitans, such as split genes separately codifying for the 5' and 3' halves of the tRNA molecule. These half genes are the oldest ancestral form of gene we have ever seen. This, along with the absence of operons from the genome of N. equitans, would seem to indicate that this genome is a molecular fossil of the evolutionary stage which the ancestral genomes had reached when the first lines of divergence were established. Moreover, the late appearance of DNA coinciding with the rooting of the universal phylogenetic tree would make the genome of N. equitans a witness to this fundamental event.

Horizontal Gene Transfer in Aminoacyl-tRNA Synthetases Including Leucine-Specific Subtypes

Aminoacyl-tRNA synthetases catalyze a fundamental reaction for the flow of genetic information from RNA to protein. Their presence in all organisms known today highlights their important role in the early evolution of life. We investigated the evolutionary history of aminoacyl-tRNA synthetases on the basis of sequence data from more than 200 Archaea, Bacteria, and Eukaryota. Phylogenetic profiles are in agreement with previous observations that many genes for aminoacyl-tRNA synthetases were transferred horizontally between species from all domains of life. We extended these findings by a detailed analysis of the history of leucyl-tRNA synthetases. Thereby, we identified a previously undetected case of horizontal gene transfer from Bacteria to Archaea based on phylogenetic profiles, trees, and networks. This means that, finally, the last subfamily of aminoacyl-tRNA synthetases has lost its exceptional position as the sole subfamily that is devoid of horizontal gene transfer. Furthermore, the leucyl-tRNA synthetase phylogenetic tree suggests a dichotomy of the archaeal/eukaryotic-cytosolic and bacterial/eukaryotic-mitochondrial proteins. We argue that the traditional division of life into Prokaryota (non-chimeric) and Eukaryota (chimeric) is favorable compared to Woese's trichotomy into Archaea/Bacteria/Eukaryota.

Structure of human tryptophanyl-tRNA synthetase in complex with tRNATrp reveals the molecular basis of tRNA recognition and specificity

Aminoacyl-tRNA synthetases (aaRSs) are a family of enzymes responsible for the covalent link of amino acids to their cognate tRNAs. The selectivity and species-specificity in the recognitions of both amino acid and tRNA by aaRSs play a vital role in maintaining the fidelity of protein synthesis. We report here the first crystal structure of human tryptophanyl-tRNA synthetase (hTrpRS) in complex with tRNA^Trp and Trp which, together with biochemical data, reveals the molecular basis of a novel tRNA binding and recognition mechanism. hTrpRS recognizes the tRNA acceptor arm from the major groove; however, the 3′ end CCA of the tRNA makes a sharp turn to bind at the active site with a deformed conformation. The discriminator base A73 is specifically recognized by an α-helix of the unique N-terminal domain and the anticodon loop by an α-helix insertion of the C-terminal domain. The N-terminal domain appears to be involved in Trp activation, but not essential for tRNA binding and acylation. Structural and sequence comparisons suggest that this novel tRNA binding and recognition mechanism is very likely shared by other archaeal and eukaryotic TrpRSs, but not by bacterial TrpRSs. Our findings provide insights into the molecular basis of tRNA specificity and species-specificity.

Domain-based small molecule binding site annotation

BACKGROUND

Accurate small molecule binding site information for a protein can facilitate studies in drug docking, drug discovery and function prediction, but small molecule binding site protein sequence annotation is sparse. The Small Molecule Interaction Database (SMID), a database of protein domain-small molecule interactions, was created using structural data from the Protein Data Bank (PDB). More importantly it provides a means to predict small molecule binding sites on proteins with a known or unknown structure and unlike prior approaches, removes large numbers of false positive hits arising from transitive alignment errors, non-biologically significant small molecules and crystallographic conditions that overpredict ion binding sites.

DESCRIPTION

Using a set of co-crystallized protein-small molecule structures as a starting point, SMID interactions were generated by identifying protein domains that bind to small molecules, using NCBI's Reverse Position Specific BLAST (RPS-BLAST) algorithm. SMID records are available for viewing at http://smid.blueprint.org. The SMID-BLAST tool provides accurate transitive annotation of small-molecule binding sites for proteins not found in the PDB. Given a protein sequence, SMID-BLAST identifies domains using RPS-BLAST and then lists potential small molecule ligands based on SMID records, as well as their aligned binding sites. A heuristic ligand score is calculated based on E-value, ligand residue identity and domain entropy to assign a level of confidence to hits found. SMID-BLAST predictions were validated against a set of 793 experimental small molecule interactions from the PDB, of which 472 (60%) of predicted interactions identically matched the experimental small molecule and of these, 344 had greater than 80% of the binding site residues correctly identified. Further, we estimate that 45% of predictions which were not observed in the PDB validation set may be true positives.

CONCLUSION

By focusing on protein domain-small molecule interactions, SMID is able to cluster similar interactions and detect subtle binding patterns that would not otherwise be obvious. Using SMID-BLAST, small molecule targets can be predicted for any protein sequence, with the only limitation being that the small molecule must exist in the PDB. Validation results and specific examples within illustrate that SMID-BLAST has a high degree of accuracy in terms of predicting both the small molecule ligand and binding site residue positions for a query protein.

Ancient gene duplications and the root(s) of the tree of life

Tracing organismal histories on the timescale of the tree of life remains one of the challenging tasks in evolutionary biology. The hotly debated questions include the evolutionary relationship between the three domains of life (e.g., which of the three domains are sister domains, are the domains para-, poly-, or monophyletic) and the location of the root within the universal tree of life. For the latter, many different points of view have been considered but so far no consensus has been reached. The only widely accepted rationale to root the universal tree of life is to use anciently duplicated paralogous genes that are present in all three domains of life. To date only few anciently duplicated gene families useful for phylogenetic reconstruction have been identified. Here we present results from a systematic search for ancient gene duplications using twelve representative, completely sequenced, archaeal and bacterial genomes. Phylogenetic analyses of identified cases show that the majority of datasets support a root between Archaea and Bacteria; however, some datasets support alternative hypotheses, and all of them suffer from a lack of strong phylogenetic signal. The results are discussed with respect to the impact of horizontal gene transfer on the ability to reconstruct organismal evolution. The exchange of genetic information between divergent organisms gives rise to mosaic genomes, where different genes in a genome have different histories. Simulations show that even low rates of horizontal gene transfer dramatically complicate the reconstruction of organismal evolution, and that the different most recent common molecular ancestors likely existed at different times and in different lineages.

Evolution of the tRNATyr/TyrRS aminoacylation systems

The tRNA identity rules ensuring fidelity of translation are globally conserved throughout evolution except for tyrosyl-tRNA synthetases (TyrRSs) that display species-specific tRNA recognition. This discrimination originates from the presence of a conserved identity pair, G1-C72, located at the top of the acceptor stem of tRNA(Tyr) from eubacteria that is invariably replaced by an unusual C1-G72 pair in archaeal and eubacterial tRNA(Tyr). In addition to the key role of pair 1-72 in tyrosylation, discriminator base A73, the anticodon triplet and the large variable region (present in eubacterial tRNA(Tyr) but not found in eukaryal tRNA(Tyr)) contribute to tyrosylation with variable strengths. Crystallographic structures of two tRNA(Tyr)/TyrRS complexes revealed different interaction modes in accordance with the phylum-specificity. Recent functional studies on the human mitochondrial tRNA(Tyr)/TyrRS system indicates strong deviations from the canonical tyrosylation rules. These differences are discussed in the light of the present knowledge on TyrRSs.

The presence of a haloarchaeal type tyrosyl-tRNA synthetase marks the opisthokonts as monophyletic

Lateral gene transfer plays an important role in the evolution of life. Events of ancient gene transfer can transmit genetic novelties to descendent lineages and subsequently shape their genetic systems. We here present the analyses of the gene encoding tyrosyl-tRNA synthetase (tyrRS), which reveal two eukaryotic tyrRS lineages, one including the opisthokonts and the other the remaining eukaryotes. The different origins of tyrRS lineages between the opisthokonts and the remaining eukaryotes indicate a likely case of ancient lateral gene transfer of tyrRS from an archaeon to the opisthokonts, which lends further support for the monophyly of the latter group. Ancient paralogy followed by differential gene loss is an alternative, albeit less parsimonious explanation for the distribution of the two eukaryotic tyrRS types. In either case, the presence of a haloarchaeal tyrRS type in the opisthokonts marks this group as monophyletic. This finding also points to the potential utility of ancient gene transfer events as molecular markers for major organismal lineages.

Variant of human enzyme sequesters reactive intermediate

In cellular environments, coupled hydrolytic reactions are used to force efficient product formation in enzyme-catalyzed reactions. In the first step of protein synthesis, aminoacyl-tRNA synthetases react with amino acid and ATP to form an enzyme-bound adenylate that, in the next step, reacts with tRNA to form aminoacyl-tRNA. The reaction liberates pyrophosphate (PP(i)) which, in turn, can be hydrolyzed by pyrophosphatase to drive efficient aminoacylation. A potential polymorphic variant of human tryptophanyl-tRNA synthetase is shown here to sequester tryptophanyl adenylate. The bound adenylate does not react efficiently with the liberated PP(i) that normally competes with tRNA to resynthesize ATP and free amino acid. Structural analysis of this variant showed that residues needed for binding ATP phosphates and thus PP(i) were reoriented from their conformations in the structure of the more common sequence variant. Significantly, the reorientation does not affect reaction with tRNA, so that efficient aminoacylation is achieved.

A novel firmicute protein family related to the actinobacterial resuscitation-promoting factors by non-orthologous domain displacement

BACKGROUND

In Micrococcus luteus growth and resuscitation from starvation-induced dormancy is controlled by the production of a secreted growth factor. This autocrine resuscitation-promoting factor (Rpf) is the founder member of a family of proteins found throughout and confined to the actinobacteria (high G + C Gram-positive bacteria). The aim of this work was to search for and characterise a cognate gene family in the firmicutes (low G + C Gram-positive bacteria) and obtain information about how they may control bacterial growth and resuscitation.

RESULTS

In silico analysis of the accessory domains of the Rpf proteins permitted their classification into several subfamilies. The RpfB subfamily is related to a group of firmicute proteins of unknown function, represented by YabE of Bacillus subtilis. The actinobacterial RpfB and firmicute YabE proteins have very similar domain structures and genomic contexts, except that in YabE, the actinobacterial Rpf domain is replaced by another domain, which we have called Sps. Although totally unrelated in both sequence and secondary structure, the Rpf and Sps domains fulfil the same function. We propose that these proteins have undergone "non-orthologous domain displacement", a phenomenon akin to "non-orthologous gene displacement" that has been described previously. Proteins containing the Sps domain are widely distributed throughout the firmicutes and they too fall into a number of distinct subfamilies. Comparative analysis of the accessory domains in the Rpf and Sps proteins, together with their weak similarity to lytic transglycosylases, provide clear evidence that they are muralytic enzymes.

CONCLUSIONS

The results indicate that the firmicute Sps proteins and the actinobacterial Rpf proteins are cognate and that they control bacterial culturability via enzymatic modification of the bacterial cell envelope.

DbW: automatic update of a functional family-specific multiple alignment

MOTIVATION

Recent advances in gene sequencing have provided complete sequence information for a number of genomes and as a result the amount of data in the sequence databases is growing at an exponential rate. We introduce here a new program, DbW, to automate the update of a functional family-specific multiple alignment that tries to include relevant sequences. The program is based on the use of different sources of information: sequences and annotations in databases.

RESULTS

The advantages of DbW are demonstrated using the 20 families of aminoacyl-tRNA synthetases, where DbW detects a maximum of homologous sequences in the Swiss-Prot and SPTREMBL databases. The global specificity of DbW in this test is 98.4% (1.6% of the sequences included in the alignment did not belong to the family according to their function), and the global sensitivity of DbW is estimated to be 95.2%. Thus, DbW provides a reliable basis for the many applications that rely on accurate multiple alignments, e.g. functional residue identification, 2D/3D structure prediction or homology modeling.

AVAILABILITY

The DbW software is available for download at ftp://ftp-igbmc.u-strasbg.fr/pub/DbW/DbW.tar and online at http://titus.u-strasbg.fr/DbW CONTACT: prigent@igbmc.u-strasbg.fr.

On the evolution of structure in aminoacyl-tRNA synthetases

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

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

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

Comparison of the catalytic roles played by the KMSKS motif in the human and Bacillus stearothermophilus trosyl-tRNA synthetases

The Class I aminoacyl-tRNA synthetases are characterized by two signature sequence motifs, "HIGH" and "KMSKS." In Bacillus stearothermophilus tyrosyl-tRNA synthetase, the KMSKS motif (230KFGKT234) has been shown to stabilize the transition state for tyrosine activation through interactions with the pyrophosphate moiety of ATP. In most eukaryotic tyrosyl-tRNA synthetases, the second lysine in the KMSKS motif is replaced by a serine or an alanine residue. Recent kinetic studies indicate that potassium functionally compensates for the absence of the second lysine in the human tyrosyl-tRNA synthetase (222KKSSS226). In this paper, site-directed mutagenesis and pre-steady state kinetics are used to determine the roles that serines 224, 225, and 226 play in catalysis of the tyrosine activation reaction. In addition, the catalytic role played by a downstream lysine conserved in eukaryotic tyrosyl-tRNA synthetases, Lys-231, is investigated. Replacing Ser-224 and Ser-226 with alanine decreases the forward rate constant 7.5- and 60-fold, respectively. In contrast, replacing either Ser-225 or Lys-231 with alanine has no effect on the catalytic activity of the enzyme. These results are consistent with the hypothesis that the KMSSS sequence in human tyrosyl-tRNA synthetase stabilizes the transition state for the tyrosine activation reaction by interacting with the pyrophosphate moiety of ATP. In addition, although they play similar roles in catalysis, the overall contribution of the KMSKS motif to catalysis appears to be significantly less in human tyrosyl-tRNA synthetase than it is in the B. stearothermophilus enzyme.

Ancient phylogenetic relationships

Traditional views on deep evolutionary events have been seriously challenged over the last few years, following the identification of major pitfalls affecting molecular phylogeny reconstruction. Here we describe the principally encountered artifacts, notably long branch attraction, and their causes (i.e., difference in evolutionary rates, mutational saturation, compositional biases). Additional difficulties due to phenomena of biological nature (i.e., lateral gene transfer, recombination, hidden paralogy) are also discussed. Moreover, contrary to common beliefs, we show that the use of rare genomic events can also be misleading and should be treated with the same caution as standard molecular phylogeny. The universal tree of life, as described in most textbooks, is partly affected by tree reconstruction artifacts, e.g. (i) the bacterial rooting of the universal tree of life; (ii) the early emergence of amitochondriate lineages in eukaryotic phylogenies; and (iii) the position of hyperthermophilic taxa in bacterial phylogenies. We present an alternative view of this tree, based on recent evidence obtained from reanalyses of ancient data sets and from novel analyses of large combination of genes.

Recent emergence of the modern genetic code: a proposal

This article proposes that the genetic code was not fully formed before the divergence of life into three kingdoms. Rather, at least arginine and tryptophan evolved after the diversification of archaea, bacteria and eukaryotes, and were spread by horizontal gene transfer. Evidence for this hypothesis is based on data suggesting that enzymes for biosynthesis of arginine and tryptophan, and for arginine tRNA ligase, have shorter divergence times than the underlying lineages. Also, many of these genes display "star" phylogenies. This proposal is an extension of the idea that the genetic code was unified because of the evolutionary pressure from horizontal gene transfer. These considerations further undermine the need to postulate the existence of a "last common ancestor"; a simpler model would be that multiple lineages gave rise to life today.

Gene family evolution and homology: genomics meets phylogenetics

With the advent of high-throughput DNA sequencing and whole-genome analysis, it has become clear that the coding portions of the genome are organized hierarchically in gene families and superfamilies. Because the hierarchy of genes, like that of living organisms, reflects an ancient and continuing process of gene duplication and divergence, many of the conceptual and analytical tools used in phylogenetic systematics can and should be used in comparative genomics. Phylogenetic principles and techniques for assessing homology, inferring relationships among genes, and reconstructing evolutionary events provide a powerful way to interpret the ever increasing body of sequence data. In this review, we outline the application of phylogenetic approaches to comparative genomics, beginning with the inference of phylogeny and the assessment of gene orthology and paralogy. We also show how the phylogenetic approach makes possible novel kinds of comparative analysis, including detection of domain shuffling and lateral gene transfer, reconstruction of the evolutionary diversification of gene families, tracing of evolutionary change in protein function at the amino acid level, and prediction of structure-function relationships. A marriage of the principles of phylogenetic systematics with the copious data generated by genomics promises unprecedented insights into the nature of biological organization and the historical processes that created it.

A human aminoacyl-tRNA synthetase as a regulator of angiogenesis

Aminoacyl-tRNA synthetases catalyze the first step of protein synthesis. It was shown recently that human tyrosyl-tRNA synthetase (TyrRS) can be split into two fragments having distinct cytokine activities, thereby linking protein synthesis to cytokine signaling pathways. Tryptophanyl-tRNA synthetase (TrpRS) is a close homologue of TyrRS. A natural fragment, herein designated as mini TrpRS, was shown by others to be produced by alternative splicing. Production of this fragment is reported to be stimulated by IFN-gamma, a cytokine that also stimulates production of angiostatic factors. Mini TrpRS is shown here to be angiostatic in a mammalian cell culture system, the chicken embryo, and two independent angiogenesis assays in the mouse. The full-length enzyme is inactive in the same assays. Thus, protein synthesis may be linked to the regulation of angiogenesis by a natural fragment of TrpRS.

Ancient adaptation of the active site of tryptophanyl-tRNA synthetase for tryptophan binding

The amino acid binding domains of the tryptophanyl (TrpRS)- and tyrosyl-tRNA synthetases (TyrRS) of Bacillus stearothermophilus are highly homologous. These similarities suggest that conserved residues in TrpRS may be responsible for both determining tryptophan recognition and discrimination against tyrosine. This was investigated by the systematic mutation of TrpRS residues based upon the identity of homologous positions in TyrRS. Of the four residues which interact directly with the aromatic side chain of tryptophan (Phe5, Met129, Asp132, and Val141) replacements of Asp132 led to significant changes in the catalytic efficiency of Trp aminoacylation (200-1250-fold reduction in k(cat)/K(M)) and substitution of Val141 by the larger Glu side chain reduced k(cat)/K(M) by 300-fold. Mutation of Pro127, which determines the position of active-site residues, did not significantly effect Trp binding. Of the mutants tested, D132N TrpRS also showed a significant reduction in discrimination against Tyr, with Tyr acting as a competitive inhibitor but not a substrate. The analogous residue in B. stearothermophilusTyrRS (Asp176) has also been implicated as a determinant of amino acid specificity in earlier studies [de Prat Gay, G., Duckworth, H. W., and Fersht, A. R. (1993) FEBS Lett. 318, 167-171]. This striking similarity in the function of a highly conserved residue found in both TrpRS and TyrRS provides mechanistic support for a common origin of the two enzymes.

Aminoacyl-tRNA synthetases, the genetic code, and the evolutionary process

The aminoacyl-tRNA synthetases (AARSs) and their relationship to the genetic code are examined from the evolutionary perspective. Despite a loose correlation between codon assignments and AARS evolutionary relationships, the code is far too highly structured to have been ordered merely through the evolutionary wanderings of these enzymes. Nevertheless, the AARSs are very informative about the evolutionary process. Examination of the phylogenetic trees for each of the AARSs reveals the following. (i) Their evolutionary relationships mostly conform to established organismal phylogeny: a strong distinction exists between bacterial- and archaeal-type AARSs. (ii) Although the evolutionary profiles of the individual AARSs might be expected to be similar in general respects, they are not. It is argued that these differences in profiles reflect the stages in the evolutionary process when the taxonomic distributions of the individual AARSs became fixed, not the nature of the individual enzymes. (iii) Horizontal transfer of AARS genes between Bacteria and Archaea is asymmetric: transfer of archaeal AARSs to the Bacteria is more prevalent than the reverse, which is seen only for the "gemini group. " (iv) The most far-ranging transfers of AARS genes have tended to occur in the distant evolutionary past, before or during formation of the primary organismal domains. These findings are also used to refine the theory that at the evolutionary stage represented by the root of the universal phylogenetic tree, cells were far more primitive than their modern counterparts and thus exchanged genetic material in far less restricted ways, in effect evolving in a communal sense.

Evolution of the genetic code

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

Archaeal genomics

Four euryarchaeal genomes have been completely sequenced and are publicly available: Methanococcus jannaschii, Methanobacterium thermoautotrophicum, Pyrococcus horikoshii and Archaeoglobus fulgidus. Four more genome sequences, two crenarchaeal and two pyrococci, will soon be released. In addition, seven more archaeal genome sequencing projects are under way, including two halophiles, two Thermoplasma, and a methanogen. These projects cover all branches of the archaeal domain and will lead to new insights into archaeal metabolism, DNA processing, and evolutionary relationships with the Bacteria and Eukarya.

Highly differentiated motifs responsible for two cytokine activities of a split human tRNA synthetase

While native human tyrosyl-tRNA synthetase (TyrRS) is inactive as a cell-signaling molecule, it can be split into two distinct cytokines. The enzyme is secreted under apoptotic conditions in culture where it is cleaved into an N-terminal fragment that harbors the catalytic site and into a C-domain fragment found only in the mammalian enzymes. The N-terminal fragment is an interleukin-8 (IL-8)-like cytokine, whereas the released C-domain is an endothelial-monocyte-activating polypeptide II (EMAP II)-like cytokine. Although the IL-8-like activity of the N-fragment depends on an ELR motif found in alpha-chemokines and conserved among mammalian TyrRSs, here we show that a similar (NYR) motif in the context of a lower eukaryote TyrRS does not confer the IL8-like activity. We also show that a heptapeptide from the C-domain has EMAP II-like chemotaxis activity for mononuclear phagocytes and polymorphonuclear leukocytes. Eukaryote proteins other than human TyrRS that have EMAP II-like domains have variants of the heptapeptide motif. Peptides based on these sequences are inactive as cytokines. Thus, the cytokine activities of split human TyrRS depend on highly differentiated motifs that are idiosyncratic to the mammalian system.

Evolution of Aminoacyl-tRNA Synthetases—Analysis of Unique Domain Architectures and Phylogenetic Trees Reveals a Complex History of Horizontal Gene Transfer Events

Phylogenetic analysis of aminoacyl-tRNA synthetases (aaRSs) of all 20 specificities from completely sequenced bacterial, archaeal, and eukaryotic genomes reveals a complex evolutionary picture. Detailed examination of the domain architecture of aaRSs using sequence profile searches delineated a network of partially conserved domains that is even more elaborate than previously suspected. Several unexpected evolutionary connections were identified, including the apparent origin of the β-subunit of bacterial GlyRS from the HD superfamily of hydrolases, a domain shared by bacterial AspRS and the B subunit of archaeal glutamyl-tRNA amidotransferases, and another previously undetected domain that is conserved in a subset of ThrRS, guanosine polyphosphate hydrolases and synthetases, and a family of GTPases. Comparison of domain architectures and multiple alignments resulted in the delineation of synapomorphies—shared derived characters, such as extra domains or inserts—for most of the aaRSs specificities. These synapomorphies partition sets of aaRSs with the same specificity into two or more distinct and apparently monophyletic groups. In conjunction with cluster analysis and a modification of the midpoint-rooting procedure, this partitioning was used to infer the likely root position in phylogenetic trees. The topologies of the resulting rooted trees for most of the aaRSs specificities are compatible with the evolutionary “standard model” whereby the earliest radiation event separated bacteria from the common ancestor of archaea and eukaryotes as opposed to the two other possible evolutionary scenarios for the three major divisions of life. For almost all aaRSs specificities, however, this simple scheme is confounded by displacement of some of the bacterial aaRSs by their eukaryotic or, less frequently, archaeal counterparts. Displacement of ancestral eukaryotic aaRS genes by bacterial ones, presumably of mitochondrial origin, was observed for three aaRSs. In contrast, there was no convincing evidence of displacement of archaeal aaRSs by bacterial ones. Displacement of aaRS genes by eukaryotic counterparts is most common among parasitic and symbiotic bacteria, particularly the spirochaetes, in which 10 of the 19 aaRSs seem to have been displaced by the respective eukaryotic genes and two by the archaeal counterpart. Unlike the primary radiation events between the three main divisions of life, that were readily traceable through the phylogenetic analysis of aaRSs, no consistent large-scale bacterial phylogeny could be established. In part, this may be due to additional gene displacement events among bacterial lineages. Argument is presented that, although lineage-specific gene loss might have contributed to the evolution of some of the aaRSs, this is not a viable alternative to horizontal gene transfer as the principal evolutionary phenomenon in this gene class.

[Complete multiple alignments of all aaRSs from complete genomes as well as the alignments of conserved regions used for phylogenetic tree construction are available at ftp://ncbi.nlm.nih.gov/pub/koonin/aaRS]

tRNA recognition and evolution of determinants in seryl-tRNA synthesis

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

Evolutionary anomalies among the aminoacyl-tRNA synthetases

Unexpected relationships among the various aminoacyl-tRNA synthetases continue to be uncovered. The question arises - is this mainly the result of promiscuous exchange, or is the confusion really a reflection of the differential loss of past duplications? Phylogenetic analysis may yet provide the answer.

Archaea and the prokaryote-to-eukaryote transition

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

History assignment: when was the mitochondrion founded?

The near simultaneous radiation of the major eukaryotic evolutionary assemblages-plants, animals, fungi, and at least three other complex protist assemblages worthy of 'kingdom level' status-was preceded by the divergence of many independent protist lineages. The earliest branches are represented by organisms that do not contain mitochondria or plastids, suggesting that the primitive eukaryotic state did not include these organelles. New information about nuclear-coded proteins that localize in the mitochondrion, however, suggests that the ancestral symbionts for mitochondria were present in the first eukaryotes. Phylogenetic support for this hypothesis is persuasive but it is not possible to account for the relative times of divergence for mitochondria and their ancestral symbionts relative to eukaryotic branching patterns inferred from nuclear genes.