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

The tRNA identity landscape for aminoacylation and beyond

tRNAs are key partners in ribosome-dependent protein synthesis. This process is highly dependent on the fidelity of tRNA aminoacylation by aminoacyl-tRNA synthetases and relies primarily on sets of identities within tRNA molecules composed of determinants and antideterminants preventing mischarging by non-cognate synthetases. Such identity sets were discovered in the tRNAs of a few model organisms, and their properties were generalized as universal identity rules. Since then, the panel of identity elements governing the accuracy of tRNA aminoacylation has expanded considerably, but the increasing number of reported functional idiosyncrasies has led to some confusion. In parallel, the description of other processes involving tRNAs, often well beyond aminoacylation, has progressed considerably, greatly expanding their interactome and uncovering multiple novel identities on the same tRNA molecule. This review highlights key findings on the mechanistics and evolution of tRNA and tRNA-like identities. In addition, new methods and their results for searching sets of multiple identities on a single tRNA are discussed. Taken together, this knowledge shows that a comprehensive understanding of the functional role of individual and collective nucleotide identity sets in tRNA molecules is needed for medical, biotechnological and other applications.

Neurodegenerative disorder and diffuse brain calcifications due to FARSB mutation in two siblings

Pathogenic mutations in the FARSB gene are associated with neurodevelopmental disorder involving the brain, liver, and lungs. We report genetic analysis of a family including two affected members with this disorder, which revealed a homozygous pathogenic missense variant, FARSB: NM_005687.4:c.853G > A:p.E285K in both affected patients. The parents were heterozygous for this variant.

Hierarchical groove discrimination by Class I and II aminoacyl-tRNA synthetases reveals a palimpsest of the operational RNA code in the tRNA acceptor-stem bases

Class I and II aaRS recognition of opposite grooves was likely among the earliest determinants fixed in the tRNA acceptor stem bases. A new regression model identifies those determinants in bacterial tRNAs. Integral coefficients relate digital dependent to independent variables with perfect agreement between observed and calculated grooves for all twenty isoaccepting tRNAs. Recognition is mediated by the Discriminator base 73, the first base pair, and base 2 of the acceptor stem. Subsets of these coefficients also identically compute grooves recognized by smaller numbers of aaRS. Thus, the model is hierarchical, suggesting that new rules were added to pre-existing ones as new amino acids joined the coding alphabet. A thermodynamic rationale for the simplest model implies that Class-dependent aaRS secondary structures exploited differential tendencies of the acceptor stem to form the hairpin observed in Class I aaRS•tRNA complexes, enabling the earliest groove discrimination. Curiously, groove recognition also depends explicitly on the identity of base 2 in a manner consistent with the middle bases of the codon table, confirming a hidden ancestry of codon-anticodon pairing in the acceptor stem. That, and the lack of correlation with anticodon bases support prior productive coding interaction of tRNA minihelices with proto-mRNA.

Structure-function studies of the asparaginyl-tRNA synthetase from Fasciola gigantica: understanding the role of catalytic and non-catalytic domains

The asparaginyl-tRNA synthetase (NRS) catalyzes the attachment of asparagine to its cognate tRNA during translation. NRS first catalyzes the binding of Asn and ATP to form the NRS-asparaginyl adenylate complex, followed by the esterification of Asn to its tRNA. We investigated the role of constituent domains in regulating the structure and activity of Fasciola gigantica NRS (FgNRS). We cloned the full-length FgNRS, along with its various truncated forms, expressed, and purified the corresponding proteins. Size exclusion chromatography indicated a role of the anticodon-binding domain (ABD) of FgNRS in protein dimerization. The N-terminal domain (NTD) was not essential for cognate tRNA binding, and the hinge region between the ABD and the C-terminal domain (CTD) was crucial for regulating the enzymatic activity. Molecular docking and fluorescence quenching experiments elucidated the binding affinities of the substrates to various domains. The molecular dynamics simulation of the modeled protein showed the presence of an unstructured region between the NTD and ABD that exhibited a large number of conformations over time, and further analysis indicated this region to be intrinsically disordered. The present study provides information on the structural and functional regulation, protein-substrate(s) interactions and dynamics, and the role of non-catalytic domains in regulating the activity of FgNRS.

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

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

Soluble expression and purification of a full-length asparaginyl tRNA synthetase from Fasciola gigantica

We report the molecular cloning, expression, and single-step homogeneous purification of a full-length asparaginyl tRNA synthetase (NRS) from Fasciola gigantica (FgNRS). Fasciola gigantica is a parasitic liver fluke of the class Trematoda. It causes fascioliasis that infects the liver of various mammals, including humans. Aminoacyl tRNA synthetases (AARS) catalyze the first step of protein synthesis. They attach an amino acid to its cognate tRNA, forming an amino acid-tRNA complex. The gene that codes for FgNRS was generated by amplification by polymerase chain reaction. It was then inserted in the expression vector pQE30 under the transcriptional control of the bacteriophage T5 promoter and lac operator. M15 Escherichia coli strain transformed with the FgNRS expression vector pQE30-NRS accumulates large amounts of a soluble protein of about 61 kDa. The protein was purified to homogeneity using immobilized metal affinity chromatography. The recombinant protein was further confirmed by immunoblotting with anti-His antibody. Following size exclusion chromatography, the FgNRS was stable and observed to be a dimeric protein. In this study, the expression and purification procedures have provided a simple and efficient method to obtain full-length FgNRS in large quantities. This will provide an opportunity to study the structure, dynamics and function of NRS.

Identifying the ligated amino acid of archaeal tRNAs based on positions outside the anticodon

Proper recognition of tRNAs by their aminoacyl-tRNA synthetase is essential for translation accuracy. Following evidence that the enzymes can recognize the correct tRNA even when anticodon information is masked, we search for additional nucleotide positions within the tRNA molecule that potentially contain information for amino acid identification. Analyzing 3936 sequences of tRNA genes from 86 archaeal species, we show that the tRNAs' cognate amino acids can be identified by the information embedded in the tRNAs' nucleotide positions without relying on the anticodon information. We present a small set of six to 10 informative positions along the tRNA, which allow for amino acid identification accuracy of 90.6% to 97.4%, respectively. We inspected tRNAs for each of the 20 amino acid types for such informative positions and found that tRNA genes for some amino acids are distinguishable from others by as few as one or two positions. The informative nucleotide positions are in agreement with nucleotide positions that were experimentally shown to affect the loaded amino acid identity. Interestingly, the knowledge gained from the tRNA genes of one archaeal phylum does not extrapolate well to another phylum. Furthermore, each species has a unique ensemble of nucleotides in the informative tRNA positions, and the similarity between the sets of positions of two distinct species reflects their evolutionary distance. Hence, we term this set of informative positions a "tRNA cipher." It is tempting to suggest that the diverging code identified here might also serve the aminoacyl tRNA synthetase in the task of tRNA recognition.

Aminoacyl-tRNA Synthetases in the Bacterial World

Aminoacyl-tRNA synthetases (aaRSs) are modular enzymes globally conserved in the three kingdoms of life. All catalyze the same two-step reaction, i.e., the attachment of a proteinogenic amino acid on their cognate tRNAs, thereby mediating the correct expression of the genetic code. In addition, some aaRSs acquired other functions beyond this key role in translation. Genomics and X-ray crystallography have revealed great structural diversity in aaRSs (e.g., in oligomery and modularity, in ranking into two distinct groups each subdivided in 3 subgroups, by additional domains appended on the catalytic modules). AaRSs show huge structural plasticity related to function and limited idiosyncrasies that are kingdom or even species specific (e.g., the presence in many Bacteria of non discriminating aaRSs compensating for the absence of one or two specific aaRSs, notably AsnRS and/or GlnRS). Diversity, as well, occurs in the mechanisms of aaRS gene regulation that are not conserved in evolution, notably between distant groups such as Gram-positive and Gram-negative Bacteria. The review focuses on bacterial aaRSs (and their paralogs) and covers their structure, function, regulation, and evolution. Structure/function relationships are emphasized, notably the enzymology of tRNA aminoacylation and the editing mechanisms for correction of activation and charging errors. The huge amount of genomic and structural data that accumulated in last two decades is reviewed, showing how the field moved from essentially reductionist biology towards more global and integrated approaches. Likewise, the alternative functions of aaRSs and those of aaRS paralogs (e.g., during cell wall biogenesis and other metabolic processes in or outside protein synthesis) are reviewed. Since aaRS phylogenies present promiscuous bacterial, archaeal, and eukaryal features, similarities and differences in the properties of aaRSs from the three kingdoms of life are pinpointed throughout the review and distinctive characteristics of bacterium-like synthetases from organelles are outlined.

Fitness Landscapes of Functional RNAs

The notion of fitness landscapes, a map between genotype and fitness, was proposed more than 80 years ago. For most of this time data was only available for a few alleles, and thus we had only a restricted view of the whole fitness landscape. Recently, advances in genetics and molecular biology allow a more detailed view of them. Here we review experimental and theoretical studies of fitness landscapes of functional RNAs, especially aptamers and ribozymes. We find that RNA structures can be divided into critical structures, connecting structures, neutral structures and forbidden structures. Such characterisation, coupled with theoretical sequence-to-structure predictions, allows us to construct the whole fitness landscape. Fitness landscapes then can be used to study evolution, and in our case the development of the RNA world.

Does the genetic code have a eukaryotic origin?

In the RNA world, RNA is assumed to be the dominant macromolecule performing most, if not all, core “house-keeping” functions. The ribo-cell hypothesis suggests that the genetic code and the translation machinery may both be born of the RNA world, and the introduction of DNA to ribo-cells may take over the informational role of RNA gradually, such as a mature set of genetic code and mechanism enabling stable inheritance of sequence and its variation. In this context, we modeled the genetic code in two content variables—GC and purine contents—of protein-coding sequences and measured the purine content sensitivities for each codon when the sensitivity (% usage) is plotted as a function of GC content variation. The analysis leads to a new pattern—the symmetric pattern—where the sensitivity of purine content variation shows diagonally symmetry in the codon table more significantly in the two GC content invariable quarters in addition to the two existing patterns where the table is divided into either four GC content sensitivity quarters or two amino acid diversity halves. The most insensitive codon sets are GUN (valine) and CAN (CAR for asparagine and CAY for aspartic acid) and the most biased amino acid is valine (always over-estimated) followed by alanine (always under-estimated). The unique position of valine and its codons suggests its key roles in the final recruitment of the complete codon set of the canonical table. The distinct choice may only be attributable to sequence signatures or signals of splice sites for spliceosomal introns shared by all extant eukaryotes.

Aminoacyl-tRNA Synthetases in the Bacterial World

Aminoacyl-tRNAsynthetases (aaRSs) are modular enzymesglobally conserved in the three kingdoms of life. All catalyze the same two-step reaction, i.e., the attachment of a proteinogenic amino acid on their cognate tRNAs, thereby mediating the correct expression of the genetic code. In addition, some aaRSs acquired other functions beyond this key role in translation.Genomics and X-ray crystallography have revealed great structural diversity in aaRSs (e.g.,in oligomery and modularity, in ranking into two distinct groups each subdivided in 3 subgroups, by additional domains appended on the catalytic modules). AaRSs show hugestructural plasticity related to function andlimited idiosyncrasies that are kingdom or even speciesspecific (e.g.,the presence in many Bacteria of non discriminating aaRSs compensating for the absence of one or two specific aaRSs, notably AsnRS and/or GlnRS).Diversity, as well, occurs in the mechanisms of aaRS gene regulation that are not conserved in evolution, notably betweendistant groups such as Gram-positive and Gram-negative Bacteria.Thereview focuses on bacterial aaRSs (and their paralogs) and covers their structure, function, regulation,and evolution. Structure/function relationships are emphasized, notably the enzymology of tRNA aminoacylation and the editing mechanisms for correction of activation and charging errors. The huge amount of genomic and structural data that accumulatedin last two decades is reviewed,showing how thefield moved from essentially reductionist biologytowards more global and integrated approaches. Likewise, the alternative functions of aaRSs and those of aaRSparalogs (e.g., during cellwall biogenesis and other metabolic processes in or outside protein synthesis) are reviewed. Since aaRS phylogenies present promiscuous bacterial, archaeal, and eukaryal features, similarities and differences in the properties of aaRSs from the three kingdoms of life are pinpointedthroughout the reviewand distinctive characteristics of bacterium-like synthetases from organelles are outlined.

Mapping hidden potential identity elements by computing the average discriminating power of individual tRNA positions

The recently published discrete mathematical method, extended consensus partition (ECP), identifies nucleotide types at each position that are strictly absent from a given sequence set, while occur in other sets. These are defined as discriminating elements (DEs). In this study using the ECP approach, we mapped potential hidden identity elements that discriminate the 20 different tRNA identities. We filtered the tDNA data set for the obligatory presence of well-established tRNA features, and then separately for each identity set, the presence of already experimentally identified strictly present identity elements. The analysis was performed on the three kingdoms of life. We determined the number of DE, e.g. the number of sets discriminated by the given position, for each tRNA position of each tRNA identity set. Then, from the positional DE numbers obtained from the 380 pairwise comparisons of the 20 identity sets, we calculated the average excluding value (AEV) for each tRNA position. The AEV provides a measure on the overall discriminating power of each position. Using a statistical analysis, we show that positional AEVs correlate with the number of already identified identity elements. Positions having high AEV but lacking published identity elements predict hitherto undiscovered tRNA identity elements.

Undetected antisense tRNAs in mitochondrial genomes?

BACKGROUND

The hypothesis that both mitochondrial (mt) complementary DNA strands of tRNA genes code for tRNAs (sense-antisense coding) is explored. This could explain why mt tRNA mutations are 6.5 times more frequently pathogenic than in other mt sequences. Antisense tRNA expression is plausible because tRNA punctuation signals mt sense RNA maturation: both sense and antisense tRNAs form secondary structures potentially signalling processing. Sense RNA maturation processes by default 11 antisense tRNAs neighbouring sense genes. If antisense tRNAs are expressed, processed antisense tRNAs should have adapted more for translational activity than unprocessed ones. Four tRNA properties are examined: antisense tRNA 5' and 3' end processing by sense RNA maturation and its accuracy, cloverleaf stability and misacylation potential.

RESULTS

Processed antisense tRNAs align better with standard tRNA sequences with the same cognate than unprocessed antisense tRNAs, suggesting less misacylations. Misacylation increases with cloverleaf fragility and processing inaccuracy. Cloverleaf fragility, misacylation and processing accuracy of antisense tRNAs decrease with genome-wide usage of their predicted cognate amino acid.

CONCLUSIONS

These properties correlate as if they adaptively coevolved for translational activity by some antisense tRNAs, and to avoid such activity by other antisense tRNAs. Analyses also suggest previously unsuspected particularities of aminoacylation specificity in mt tRNAs: combinations of competition between tRNAs on tRNA synthetases with competition between tRNA synthetases on tRNAs determine specificities of tRNA amino acylations. The latter analyses show that alignment methods used to detect tRNA cognates yield relatively robust results, even when they apparently fail to detect the tRNA's cognate amino acid and indicate high misacylation potential.

Computational analysis of tRNA identity

I review recent developments in computational analysis of tRNA identity. I suggest that the tRNA-protein interaction network is hierarchically organized, and coevolutionarily flexible. Its functional specificity of recognition and discrimination persists despite generic structural constraints and perturbative evolutionary forces. This flexibility comes from its arbitrary nature as a self-recognizing shape code. A revisualization of predicted Proteobacterial tRNA identity highlights open research problems. tRNA identity elements and their coevolution with proteins must be mapped structurally over the Tree of Life. These traits can also resolve deep roots in the Tree. I show that histidylation identity elements phylogenetically reposition Pelagibacter ubique within alpha-Proteobacteria.

Toward a more complete view of tRNA biology

Transfer RNAs are ancient molecules present in all domains of life. In addition to translating the genetic code into protein and defining the second genetic code together with aminoacyl-tRNA synthetases, tRNAs act in many other cellular functions. Robust phenomenological observations on the role of tRNAs in translation, together with massive sequence and crystallographic data, have led to a deeper physicochemical understanding of tRNA architecture, dynamics and identity. In vitro studies complemented by cell biology data already indicate how tRNA behaves in cellular environments, in particular in higher Eukarya. From an opposite approach, reverse evolution considerations suggest how tRNAs emerged as simplified structures from the RNA world. This perspective discusses what basic questions remain unanswered, how these answers can be obtained and how a more rational understanding of the function and dysfunction of tRNA can have applications in medicine and biotechnology.