Coexistence of bacterial leucyl-tRNA synthetases with archaeal tRNA binding domains that distinguish tRNA(Leu) in the archaeal mode

EcoGenoRisk: Developing a computational ecological risk assessment tool for synthetic biology

The expanding field of synthetic biology (synbio) supports new opportunities in the design of targeted bioproducts or modified microorganisms. However, this rapid development of synbio products raises concerns surrounding the potential risks of modified microorganisms contaminating unintended environments. These potential invasion risks require new bioinformatic tools to inform the design phase. EcoGenoRisk is a newly constructed computational risk assessment tool for invasiveness that aims to predict where synbio microorganisms may establish a population by screening for habitats of genetically similar microorganisms. The first module of the tool identifies genetically similar microorganisms and potential ecological relationships such as competition, mutualism, and inhibition. In total, 520 archaeal and 32,828 bacterial complete assembly genomes were analyzed to test the specificity and accuracy of the tool as well as to characterize the enzymatic profiles of different taxonomic lineages. Additionally, ecological relationships were analyzed to determine which would result in the greatest potential overlap between shared functional profiles. Notably, competition displayed the significantly highest overlap of shared functions between compared genomes. Overall, EcoGenoRisk is a flexible software pipeline that assists environmental risk assessors to query large databases of known microorganisms and prioritize follow-up bench scale studies.

Diversification of aminoacyl-tRNA synthetase activities via genomic duplication

Intricate evolutionary events enabled the emergence of the full set of aminoacyl-tRNA synthetase (aaRS) families that define the genetic code. The diversification of aaRSs has continued in organisms from all domains of life, yielding aaRSs with unique characteristics as well as aaRS-like proteins with innovative functions outside translation. Recent bioinformatic analyses have revealed the extensive occurrence and phylogenetic diversity of aaRS gene duplication involving every synthetase family. However, only a fraction of these duplicated genes has been characterized, leaving many with biological functions yet to be discovered. Here we discuss how genomic duplication is associated with the occurrence of novel aaRSs and aaRS-like proteins that provide adaptive advantages to their hosts. We illustrate the variety of activities that have evolved from the primordial aaRS catalytic sites. This precedent underscores the need to investigate currently unexplored aaRS genomic duplications as they may hold a key to the discovery of exciting biological processes, new drug targets, important bioactive molecules, and tools for synthetic biology applications.

Distinct pathogenic mechanisms of various RARS1 mutations in Pelizaeus-Merzbacher-like disease

Mutations of the genes encoding aminoacyl-tRNA synthetases are highly associated with various central nervous system disorders. Recurrent mutations, including c.5A>G, p.D2G; c.1367C>T, p.S456L; c.1535G>A, p.R512Q and c.1846_1847del, p. Y616Lfs*6 of RARS1 gene, which encodes two forms of human cytoplasmic arginyl-tRNA synthetase (hArgRS), are linked to Pelizaeus-Merzbacher-like disease (PMLD) with unclear pathogenesis. Among these mutations, c.5A>G is the most extensively reported mutation, leading to a p.D2G mutation in the N-terminal extension of the long-form hArgRS. Here, we showed the detrimental effects of R512Q substitution and ΔC mutations on the structure and function of hArgRS, while the most frequent mutation c.5A>G, p.D2G acted in a different manner without impairing hArgRS activity. The nucleotide substitution c.5A>G reduced translation of hArgRS mRNA, and an upstream open reading frame contributed to the suppressed translation of the downstream main ORF. Taken together, our results elucidated distinct pathogenic mechanisms of various RARS1 mutations in PMLD.

Duplication of leucyl-tRNA synthetase in an archaeal extremophile may play a role in adaptation to variable environmental conditions

Aminoacyl-tRNA synthetases (aaRSs) are ancient enzymes that play a fundamental role in protein synthesis. They catalyze the esterification of specific amino acids to the 3'-end of their cognate tRNAs and therefore play a pivotal role in protein synthesis. Although previous studies suggest that aaRS-dependent errors in protein synthesis can be beneficial to some microbial species, evidence that reduced aaRS fidelity can be adaptive is limited. Using bioinformatics analyses, we identified two distinct leucyl-tRNA synthetase (LeuRS) genes within all genomes of the archaeal family Sulfolobaceae. Remarkably, one copy, designated LeuRS-I, had key amino acid substitutions within its editing domain that would be expected to disrupt hydrolytic editing of mischarged tRNA^Leu and to result in variation within the proteome of these extremophiles. We found that another copy, LeuRS-F, contains canonical active sites for aminoacylation and editing. Biochemical and genetic analyses of the paralogs within Sulfolobus islandicus supported the hypothesis that LeuRS-F, but not LeuRS-I, functions as an essential tRNA synthetase that accurately charges leucine to tRNA^Leu for protein translation. Although LeuRS-I was not essential, its expression clearly supported optimal S. islandicus growth. We conclude that LeuRS-I may have evolved to confer a selective advantage under the extreme and fluctuating environmental conditions characteristic of the volcanic hot springs in which these archaeal extremophiles reside.

LeuRS can leucylate type I and type II tRNALeus in Streptomyces coelicolor

Transfer RNAs (tRNAs) are divided into two types, type I with a short variable loop and type II with a long variable loop. Aminoacylation of type I or type II tRNA^Leu is catalyzed by their cognate leucyl-tRNA synthetases (LeuRSs). However, in Streptomyces coelicolor, there are two types of tRNA^Leu and only one LeuRS (ScoLeuRS). We found that the enzyme could leucylate both types of ScotRNA^Leu, and had a higher catalytic efficiency for type II ScotRNA^Leu(UAA) than for type I ScotRNA^Leu(CAA). The results from tRNA and enzyme mutagenesis showed that ScoLeuRS did not interact with the canonical discriminator A73. The number of nucleotides, rather than the type of base of the variable loop in the two types of ScotRNA^Leus, was determined as important for aminoacylation. In vitro and in vivo assays showed that the tertiary structure formed by the D-loop and TψC-loop is more important for ScotRNA^Leu(UAA)_. We showed that the leucine-specific domain (LSD) of ScoLeuRS could help LeuRS, which originally only leucylates type II tRNA^Leu, to aminoacylate type I ScotRNA^Leu(CAA) and identified the crucial amino acid residues at the C-terminus of the LSD to recognize type I ScotRNA^Leu(CAA). Overall, our findings identified a rare recognition mechanism of LeuRS to tRNA^Leu.

Recent development of leucyl-tRNA synthetase inhibitors as antimicrobial agents

Aminoacyl-tRNA synthetases (aaRSs) widely exist in organisms and mediate protein synthesis. Inhibiting these synthetases can lead to the termination of protein synthesis and subsequently achieve antibacterial and antiparasitic purposes. Moreover, the structures of aaRSs found in eukaryotes have considerable structural differences compared to those in prokaryotes, based on which it is possible to develop highly selective inhibitors. Leucyl-tRNA synthetase (LeuRS) with unique synthesis and editing sites is one of 20 kinds of aaRSs. Many inhibitors targeting LeuRS have been designed and synthesized, some of which have entered clinical use. For example, the benzoxaborole compound AN2690 has been approved by the FDA for the treatment of onychomycosis. AN3365 is suspended in the phase II clinical trial due to the rapid development of AN3365 resistance, but it may be used in combination with other antibiotics. The aaRSs, especially LeuRS, are being considered as targets of new potential anti-infective drugs for the treatment of not only bacterial or fungal infections but also infections by trypanosomes and malaria parasites. This review mainly describes the development of LeuRS inhibitors, focusing on their mechanisms of action, structure-activity relationships (SARs), and in vitro and in vivo activities.

A threonyl-tRNA synthetase-like protein has tRNA aminoacylation and editing activities

TARS and TARS2 encode cytoplasmic and mitochondrial threonyl-tRNA synthetases (ThrRSs) in mammals, respectively. Interestingly, in higher eukaryotes, a third gene, TARSL2, encodes a ThrRS-like protein (ThrRS-L), which is highly homologous to cytoplasmic ThrRS but with a different N-terminal extension (N-extension). Whether ThrRS-L has canonical functions is unknown. In this work, we studied the organ expression pattern, cellular localization, canonical aminoacylation and editing activities of mouse ThrRS-L (mThrRS-L). Tarsl2 is ubiquitously but unevenly expressed in mouse tissues. Different from mouse cytoplasmic ThrRS (mThrRS), mThrRS-L is located in both the cytoplasm and nucleus; the nuclear distribution is mediated via a nuclear localization sequence at its C-terminus. Native mThrRS-L enriched from HEK293T cells was active in aminoacylation and editing. To investigate the in vitro catalytic properties of mThrRS-L accurately, we replaced the N-extension of mThrRS-L with that of mThrRS. The chimeric protein (mThrRS-L-N^T) has amino acid activation, aminoacylation and editing activities. We compared the activities and cross-species tRNA recognition between mThrRS-L-N^T and mThrRS. Despite having a similar aminoacylation activity, mThrRS-L-N^T and mThrRS exhibit differences in tRNA recognition and editing capacity. Our results provided the first analysis of the aminoacylation and editing activities of ThrRS-L, and improved our understanding of Tarsl2.

A Rice Glutamyl-tRNA Synthetase Modulates Early Anther Cell Division and Patterning

Aminoacyl-tRNA synthetases (aaRSs) have housekeeping roles in protein synthesis, but little is known about how these aaRSs are involved in organ development. Here, we report that a rice (Oryza sativa) glutamyl-tRNA synthetase (OsERS1) maintains proper somatic cell organization and limits the overproliferation of male germ cells during early anther development. The expression of OsERS1 is specifically detectable in meristematic layer 2-derived cells of the early anther, and osers1 anthers exhibit overproliferation and disorganization of layer 2-derived cells, producing fused lobes and extra germ cells in early anthers. The conserved biochemical function of OsERS1 in ligating glutamate to tRNA^Glu is enhanced by its cofactor aaRS OsARC. Furthermore, metabolomics profiling revealed that OsERS1 is an important node for multiple metabolic pathways, indicated by the accumulation of amino acids and tricarboxylic acid cycle components in osers1 anthers. Notably, the anther defects of the osers1 mutant are causally associated with the abnormal accumulation of hydrogen peroxide, which can reconstitute the osers1 phenotype when applied to wild-type anthers. Collectively, these findings demonstrate how aaRSs affect male organ development in plants, likely through protein synthesis, metabolic homeostasis, and redox status.

Self-protective responses to norvaline-induced stress in a leucyl-tRNA synthetase editing-deficient yeast strain

The editing function of aminoacyl-tRNA synthetases (aaRSs) is indispensible for formation of the correct aminoacyl-tRNAs. Editing deficiency may lead to growth inhibition and the pathogenesis of various diseases. Herein, we confirmed that norvaline (Nva) but not isoleucine or valine is the major threat to the editing function of Saccharomyces cerevisiae leucyl-tRNA synthetase (ScLeuRS), both in vitro and in vivo. Nva could be misincorporated into the proteome of the LeuRS editing-deficient yeast strain (D419A/ScΔleuS), potentially resulting in dysfunctional protein folding and growth delay. Furthermore, the exploration of the Nva-induced intracellular stress response mechanism in D419A/ScΔleuS revealed that Hsp70 chaperones were markedly upregulated in response to the potential protein misfolding. Additionally, proline (Pro), glutamate (Glu) and glutamine (Gln), which may accumulate due to the conversion of Nva, collectively contributed to the reduction of reactive oxygen species (ROS) levels in Nva-treated D419A/ScΔleuS cells. In conclusion, our study highlights the significance of the editing function of LeuRS and provides clues for understanding the intracellular stress protective mechanisms that are triggered in aaRS editing-deficient organisms.

Aminoacetylation Reaction Catalyzed by Leucyl-tRNA Synthetase Operates via a Self-Assisted Mechanism Using a Conserved Residue and the Aminoacyl Substrate

Leucyl-tRNA synthetase catalyzes attachment of leucine amino acid to its cognate tRNA. During the second, aminoacetylation, step of the reaction, the leucyl moiety is transferred from leucyl-adenylate to the terminal A76 adenosine of tRNA. In this work, we have investigated the aminoacetylation step catalyzed by leucyl-tRNA synthase, using ab initio quantum chemical/molecular mechanical hybrid potentials in conjunction with reaction-path-location algorithms and molecular dynamics free energy simulations. We have modeled reaction mechanisms arising from both crystallographic studies and computational work. We invoke various groups as potential proton acceptors-namely, the phosphate and leucyl amino groups of leucyl-adenylate, the A76 base of tRNA, and the Asp80 and Glu532 residues of the protein-and consider both metal-assisted and metal-free reactions. Free energy calculations indicate that both the phosphate group of leucyl adenylate and Glu532 are not strong bases. This agrees with the results of the quantum chemical/molecular mechanical reaction path calculations which give high free energy barriers for the studied pathways involving these groups. A self-assisted mechanism with the leucyl amino group and Asp80 as proton acceptors is the most likely. Furthermore, in this mechanism the presence of a metal ion coordinated by the phosphate group and Glu532 strongly activates the reaction.

C-terminal Domain of Leucyl-tRNA Synthetase from Pathogenic Candida albicans Recognizes both tRNASer and tRNALeu

Leucyl-tRNA synthetase (LeuRS) is a multidomain enzyme that catalyzes Leu-tRNA(Leu) formation and is classified into bacterial and archaeal/eukaryotic types with significant diversity in the C-terminal domain (CTD). CTDs of both bacterial and archaeal LeuRSs have been reported to recognize tRNA(Leu) through different modes of interaction. In the human pathogen Candida albicans, the cytoplasmic LeuRS (CaLeuRS) is distinguished by its capacity to recognize a uniquely evolved chimeric tRNA(Ser) (CatRNA(Ser)(CAG)) in addition to its cognate CatRNA(Leu), leading to CUG codon reassignment. Our previous study showed that eukaryotic but not archaeal LeuRSs recognize this peculiar tRNA(Ser), suggesting the significance of their highly divergent CTDs in tRNA(Ser) recognition. The results of this study provided the first evidence of the indispensable function of the CTD of eukaryotic LeuRS in recognizing non-cognate CatRNA(Ser) and cognate CatRNA(Leu). Three lysine residues were identified as involved in mediating enzyme-tRNA interaction in the leucylation process: mutation of all three sites totally ablated the leucylation activity. The importance of the three lysine residues was further verified by gel mobility shift assays and complementation of a yeast leuS gene knock-out strain.

Degenerate connective polypeptide 1 (CP1) domain from human mitochondrial leucyl-tRNA synthetase

The connective polypeptide 1 (CP1) editing domain of leucyl-tRNA synthetase (LeuRS) from various species either harbors a conserved active site to exclude tRNA mis-charging with noncognate amino acids or is evolutionarily truncated or lost because there is no requirement for high translational fidelity. However, human mitochondrial LeuRS (hmtLeuRS) contains a full-length but degenerate CP1 domain that has mutations in some residues important for post-transfer editing. The significance of such an inactive CP1 domain and a translational accuracy mechanism with different noncognate amino acids are not completely understood. Here, we identified the essential role of the evolutionarily divergent CP1 domain in facilitating hmtLeuRS's catalytic efficiency and endowing enzyme with resistance to AN2690, a broad-spectrum drug acting on LeuRSs. In addition, the canonical core of hmtLeuRS is not stringent for noncognate norvaline (Nva) and valine (Val). hmtLeuRS has a very weak tRNA-independent pre-transfer editing activity for Nva, which is insufficient to remove mis-activated Nva. Moreover, hmtLeuRS chimeras fused with a functional CP1 domain from LeuRSs of other species, regardless of origin, showed restored post-transfer editing activity and acquired fidelity during aminoacylation. This work offers a novel perspective on the role of the CP1 domain in optimizing aminoacylation efficiency.

Ancient horizontal gene transfer and the last common ancestors

Background

The genomic history of prokaryotic organismal lineages is marked by extensive horizontal gene transfer (HGT) between groups of organisms at all taxonomic levels. These HGT events have played an essential role in the origin and distribution of biological innovations. Analyses of ancient gene families show that HGT existed in the distant past, even at the time of the organismal last universal common ancestor (LUCA). Most gene transfers originated in lineages that have since gone extinct. Therefore, one cannot assume that the last common ancestors of each gene were all present in the same cell representing the cellular ancestor of all extant life.

Results

Organisms existing as part of a diverse ecosystem at the time of LUCA likely shared genetic material between lineages. If these other lineages persisted for some time, HGT with the descendants of LUCA could have continued into the bacterial and archaeal lineages. Phylogenetic analyses of aminoacyl-tRNA synthetase protein families support the hypothesis that the molecular common ancestors of the most ancient gene families did not all coincide in space and time. This is most apparent in the evolutionary histories of seryl-tRNA synthetase and threonyl-tRNA synthetase protein families, each containing highly divergent “rare” forms, as well as the sparse phylogenetic distributions of pyrrolysyl-tRNA synthetase, and the bacterial heterodimeric form of glycyl-tRNA synthetase. These topologies and phyletic distributions are consistent with horizontal transfers from ancient, likely extinct branches of the tree of life.

Conclusions

Of all the organisms that may have existed at the time of LUCA, by definition only one lineage is survived by known progeny; however, this lineage retains a genomic record of heterogeneous genetic origins. The evolutionary histories of aminoacyl-tRNA synthetases (aaRS) are especially informative in detecting this signal, as they perform primordial biological functions, have undergone several ancient HGT events, and contain many sites with low substitution rates allowing deep phylogenetic reconstruction. We conclude that some aaRS families contain groups that diverge before LUCA. We propose that these ancient gene variants be described by the term “hypnologs”, reflecting their ancient, reticulate origin from a time in life history that has been all but erased”.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0350-0) contains supplementary material, which is available to authorized users.

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.