Amino acid activation analysis of primitive aminoacyl-tRNA synthetases encoded by both strands of a single gene using the malachite green assay

Aminoacyl-tRNA synthetase urzymes optimized by deep learning behave as a quasispecies

Protein design plays a key role in our efforts to work out how genetic coding began. That effort entails urzymes. Urzymes are small, conserved excerpts from full-length aminoacyl-tRNA synthetases that remain active. Urzymes require design to connect disjoint pieces and repair naked nonpolar patches created by removing large domains. Rosetta allowed us to create the first urzymes, but those urzymes were only sparingly soluble. We could measure activity, but it was hard to concentrate those samples to levels required for structural biology. Here, we used the deep learning algorithms ProteinMPNN and AlphaFold2 to redesign a set of optimized LeuAC urzymes derived from leucyl-tRNA synthetase. We select a balanced, representative subset of eight variants for testing using principal component analysis. Most tested variants are much more soluble than the original LeuAC. They also span a range of catalytic proficiency and amino acid specificity. The data enable detailed statistical analyses of the sources of both solubility and specificity. In that way, we show how to begin to unwrap the elements of protein chemistry that were hidden within the neural networks. Deep learning networks have thus helped us surmount several vexing obstacles to further investigations into the nature of ancestral proteins. Finally, we discuss how the eight variants might resemble a sample drawn from a population similar to one subject to natural selection.

WITHDRAWN: Structural Enzymology, Phylogenetics, Differentiation, and Symbolic Reflexivity at the Dawn of Biology

This manuscript was posted without the final consent of all authors. The authors have therefore withdrawn it. The authors do not wish this work to be cited as reference for the project. If you have any questions, please contact the corresponding author, carter@med.unc.edu .

Transcriptomics, metabolomics and proteomics analyses reveal glyphosate tolerance mechanism in red swamp crayfish Procambarus clarkii

Glyphosate (Gly), the world's most widely used herbicide in agriculture, can poison the red swamp crayfish, Procambarus clarkii, via spray drift and surface runoff into surface waters. However, there is a paucity of research on the mechanisms that affect crayfish tolerance to Gly at typical environmental concentrations. To address this research gap, we investigated the effects of Gly stress (0, 6, 12, 24, and 72 h) at different concentrations (0, 1.20, 3.60, 7.20, and 10.80 mg·L-1) on antioxidant enzyme activity in crayfish hepatopancreas. Furthermore, we analyzed the species' tolerance mechanism to Gly exposure at a typical environmental concentration (3.60 mg·L-1) based on integrative transcriptomics, metabolomics, and proteomics. The Gly concentration and exposure time affected the crayfish's antioxidant system, and interacted with each other (P < 0.01). Gly concentrations higher than 7.20 mg·L-1 and exposure times longer than 48 h caused oxidative stress. When the Gly concentrations were lower than 3.60 mg·L-1, crayfish tolerated Gly exposure within 72 h by self-regulating superoxide dismutase (SOD), catalase (CAT), and malondialdehyde (MDA). A multi-omics analysis revealed that crayfish upregulated the expression of amino acid metabolites (such as glutamate, proline, and lysine) and amino acid transformation-related genes (such as GlnA and P5CS) to tolerate Gly stress by enhancing the antioxidant capacity, ammonia‑nitrogen regulation, and energy supply of the organism. Metallothionein and polyadenylate-binding proteins, which are potential markers of Gly exposure, crucially influenced crayfish tolerance to Gly by synthesizing metalloenzymes and scavenging reactive oxygen species. This study revealed the Gly tolerance mechanism in crayfish and can provide a theoretical reference for commercial eco-farming in rice-crayfish integrated aquaculture systems.

AARS Online: A collaborative database on the structure, function, and evolution of the aminoacyl-tRNA synthetases

The aminoacyl-tRNA synthetases (aaRS) are a large group of enzymes that implement the genetic code in all known biological systems. They attach amino acids to their cognate tRNAs, moonlight in various translational and non-translational activities beyond aminoacylation, and are linked to many genetic disorders. The aaRS have a subtle ontology characterized by structural and functional idiosyncrasies that vary from organism to organism, and protein to protein. Across the tree of life, the 22 coded amino acids are handled by 16 evolutionary families of Class I aaRS and 21 families of Class II aaRS. We introduce AARS Online, an interactive Wikipedia-like tool curated by an international consortium of field experts. This platform systematizes existing knowledge about the aaRS by showcasing a taxonomically diverse selection of aaRS sequences and structures. Through its graphical user interface, AARS Online facilitates a seamless exploration between protein sequence and structure, providing a friendly introduction to the material for non-experts and a useful resource for experts. Curated multiple sequence alignments can be extracted for downstream analyses. Accessible at www.aars.online, AARS Online is a free resource to delve into the world of the aaRS.

A genomic database furnishes minimal functional glycyl-tRNA synthetases homologous to other, designed class II urzymes

The hypothesis that conserved core catalytic sites could represent ancestral aminoacyl-tRNA synthetases (AARS) drove the design of functional TrpRS, LeuRS, and HisRS 'urzymes'. We describe here new urzymes detected in the genomic record of the arctic fox, Vulpes lagopus. They are homologous to the α-subunit of bacterial heterotetrameric Class II glycyl-tRNA synthetase (GlyRS-B) enzymes. AlphaFold2 predicted that the N-terminal 81 amino acids would adopt a 3D structure nearly identical to our designed HisRS urzyme (HisCA1). We expressed and purified that N-terminal segment and the spliced open reading frame GlyCA1-2. Both exhibit robust single-turnover burst sizes and ATP consumption rates higher than those previously published for HisCA urzymes and comparable to those for LeuAC and TrpAC. GlyCA is more than twice as active in glycine activation by adenosine triphosphate as the full-length GlyRS-B α2 dimer. Michaelis-Menten rate constants for all three substrates reveal significant coupling between Exon2 and both substrates. GlyCA activation favors Class II amino acids that complement those favored by HisCA and LeuAC. Structural features help explain these results. These minimalist GlyRS catalysts are thus homologous to previously described urzymes. Their properties reinforce the notion that urzymes may have the requisite catalytic activities to implement a reduced, ancestral genetic coding alphabet.

Effects of Tannin Supplementation in Diet on the Resistance to Ammonia Stress of Pacific White Shrimp Litopenaeus vannamei

Tannin (TA), as a natural phenolic compound with strong antioxidant activity, has been used as a feed additive for various animals. In this study, we fed a diet containing 800 mg/kg of tannin on Litopenaeus vannamei for 56 days and then subjected to acute ammonia stress for 48 hr to investigate the effect of dietary tannin on the ammonia stress response of L. vannamei through transcriptomic and metabolomic analysis. The transcriptome analysis indicated that ammonia stress-induced differential expression of 4,185 genes, while tannin-fed shrimp only had 964 differentially expressed genes. Compared with the TA_0 group, 59 pathways were significantly altered, and the pathways of "starch and sucrose metabolism," "retinol metabolism," "arachidonic acid metabolism," "lysosome," and "amino sugar and nucleotide sugar metabolism" were highly enriched in the TS_0 group. Compared with the TS_0 group, six pathways were significantly altered, and the pathways of "dilated cardiomyopathy," "complement and coagulation cascades," "cardiac muscle contraction," "fructose and mannose metabolism," "cGMP-PKG signaling pathway," and "beta-alanine metabolism" were significantly enriched in the TS_800 group. Metabolomic analysis showed that a total of 107 differential metabolites (DMs) were identified in the TS_0 vs. TA_0 group, while 75 DMs were identified in the TS_800 vs. TS_0 group. Based on KEGG annotation, it was found that a large amount of DM was significantly enriched in amino acid metabolism pathways in the TS_0 group, including "arginine and proline metabolism," "alanine, aspartic acid, and glutamic acid metabolism," "β-Alanine metabolism and tyrosine metabolism" indicated that tannins affect the metabolism of amino acids. The integration of DEGs and DMs indicates that dietary tannins highly alter the digestion and absorption functions of proteins, as well as the biosynthesis and metabolism of amino acids. This study provides new insights into the adaptation of Pacific white shrimp to ammonia stress and the addition of tannins to feed to enhance immune function.

Elucidation of productive alanine recognition mechanism by Escherichia coli alanyl-tRNA synthetase

Alanyl-tRNA synthetase (AlaRS) incorrectly recognizes both a slightly smaller glycine and a slightly larger serine in addition to alanine, and the probability of incorrect identification is extremely low at 1/300 and 1/170, respectively. Alanine is the second smallest amino acid after glycine; however, the mechanism by which AlaRS specifically identifies small differences in side chains with high accuracy remains unknown. In this study, using a malachite green assay, we aimed to elucidate the alanine recognition mechanism of a fragment (AlaRS368N) containing only the amino acid activation domain of Escherichia coli AlaRS. This method quantifies monophosphate by decomposing pyrophosphate generated during aminoacyl-AMP production. AlaRS368N produced far more pyrophosphate when glycine or serine was used as a substrate than when alanine was used. Among several mutants tested, an AlaRS mutant in which the widely conserved aspartic acid at the 235th position (D235) near the active center was replaced with glutamic acid (D235E) increased pyrophosphate release for the alanine substrate, compared to that from glycine and serine. These results suggested that D235 is optimal for AlaRS to specifically recognize alanine. Alanylation activities of an RNA minihelix by the mutants of valine at the 214th position (V214) of another fragment (AlaRS442N), which is the smallest AlaRS with alanine charging activity, suggest the existence of the van der Waals-like interaction between the side chain of V214 and the methyl group of the alanine substrate.

Base Pairing Promoted the Self-Organization of Genetic Coding, Catalysis, and Free-Energy Transduction

How Nature discovered genetic coding is a largely ignored question, yet the answer is key to explaining the transition from biochemical building blocks to life. Other, related puzzles also fall inside the aegis enclosing the codes themselves. The peptide bond is unstable with respect to hydrolysis. So, it requires some form of chemical free energy to drive it. Amino acid activation and acyl transfer are also slow and must be catalyzed. All living things must thus also convert free energy and synchronize cellular chemistry. Most importantly, functional proteins occupy only small, isolated regions of sequence space. Nature evolved heritable symbolic data processing to seek out and use those sequences. That system has three parts: a memory of how amino acids behave in solution and inside proteins, a set of code keys to access that memory, and a scoring function. The code keys themselves are the genes for cognate pairs of tRNA and aminoacyl-tRNA synthetases, AARSs. The scoring function is the enzymatic specificity constant, kcat/kM, which measures both catalysis and specificity. The work described here deepens the evidence for and understanding of an unexpected consequence of ancestral bidirectional coding. Secondary structures occur in approximately the same places within antiparallel alignments of their gene products. However, the polar amino acids that define the molecular surface of one are reflected into core-defining non-polar side chains on the other. Proteins translated from base-paired coding strands fold up inside out. Bidirectional genes thus project an inverted structural duality into the proteome. I review how experimental data root the scoring functions responsible for the origins of coding and catalyzed activation of unfavorable chemical reactions in that duality.

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

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

Genomic database furnishes a spontaneous example of a functional Class II glycyl-tRNA synthetase urzyme

The chief barrier to studies of how genetic coding emerged is the lack of experimental models for ancestral aminoacyl-tRNA synthetases (AARS). We hypothesized that conserved core catalytic sites could represent such ancestors. That hypothesis enabled engineering functional "urzymes" from TrpRS, LeuRS, and HisRS. We describe here a fourth urzyme, GlyCA, detected in an open reading frame from the genomic record of the arctic fox, Vulpes lagopus. GlyCA is homologous to a bacterial heterotetrameric Class II GlyRS-B. Alphafold2 predicted that the N-terminal 81 amino acids would adopt a 3D structure nearly identical to the HisRS urzyme (HisCA1). We expressed and purified that N-terminal segment. Enzymatic characterization revealed a robust single-turnover burst size and a catalytic rate for ATP consumption well in excess of that previously published for HisCA1. Time-dependent aminoacylation of tRNAGly proceeds at a rate consistent with that observed for amino acid activation. In fact, GlyCA is actually 35 times more active in glycine activation by ATP than the full-length GlyRS-B α-subunit dimer. ATP-dependent activation of the 20 canonical amino acids favors Class II amino acids that complement those favored by HisCA and LeuAC. These properties reinforce the notion that urzymes represent the requisite ancestral catalytic activities to implement a reduced genetic coding alphabet.

Domain acquisition by class I aminoacyl-tRNA synthetase urzymes coordinated the catalytic functions of HVGH and KMSKS motifs

Leucyl-tRNA synthetase (LeuRS) is a Class I aminoacyl-tRNA synthetase (aaRS) that synthesizes leucyl-tRNAleu for codon-directed protein synthesis. Two signature sequences, HxGH and KMSKS help stabilize transition-states for amino acid activation and tRNA aminoacylation by all Class I aaRS. Separate alanine mutants of each signature, together with the double mutant, behave in opposite ways in Pyrococcus horikoshii LeuRS and the 129-residue urzyme ancestral model generated from it (LeuAC). Free energy coupling terms, Δ(ΔG‡), for both reactions are large and favourable for LeuRS, but unfavourable for LeuAC. Single turnover assays with 32Pα-ATP show correspondingly different internal products. These results implicate domain motion in catalysis by full-length LeuRS. The distributed thermodynamic cycle of mutational changes authenticates LeuAC urzyme catalysis far more convincingly than do single point mutations. Most importantly, the evolutionary gain of function induced by acquiring the anticodon-binding (ABD) and multiple insertion modules in the catalytic domain appears to be to coordinate the catalytic function of the HxGH and KMSKS signature sequences. The implication that backbone elements of secondary structures achieve a major portion of the overall transition-state stabilization by LeuAC is also consistent with coevolution of the genetic code and metabolic pathways necessary to produce histidine and lysine sidechains.

Integration of transcriptomic and metabolomic reveals carbonate alkalinity stress responses in the hepatopancreas of Litopenaeus vannamei

Carbonate alkalinity (CA) is one of the environmental factors affecting the survival and growth of aquatic animals. However, the toxic effects of CA stress on Pacific white shrimp Litopenaeus vannamei at the molecular level are completely unclear. In this study, we investigated the changes of the survival and growth, and hepatopancreas histology of L. vannamei under different levels of CA stress, and integrated transcriptomics and metabolomics to explore major functional changes in the hepatopancreas and identify biomarkers. After CA exposure for 14 days, the survival and growth of the shrimp were reduced, and the hepatopancreas showed obvious histological damage. A total of 253 genes were differentially expressed in the three CA stress groups, and immune-related genes such as pattern recognition receptors, phenoloxidase system and detoxification metabolism were affected; substance transport-related regulators and transporters were mostly downregulated. Furthermore, the metabolic pattern of the shrimp was also altered by CA stress, especially amino acids, arachidonic acid and B-vitamin metabolites. The integration analysis of differential metabolites and genes further showed that the functions of ABC transporters, protein digestion and absorption, and amino acid biosynthesis and metabolism were highly altered by CA stress. The results of this study revealed that CA stress caused immune, substance transport, and amino acid metabolic variations in L. vannamei, and identified several potential biomarkers related to stress response.

Toxicological response of Pacific white shrimp Litopenaeus vannamei to a hazardous cyanotoxin nodularin exposure

Nodularin (NOD) is a harmful cyanotoxin that affects shrimp farming. The hepatopancreas and intestine of shrimp are the main target organs of cyanotoxins. In this study, we exposed Litopenaeus vannamei to NOD at 0.1 and 1 μg/L for 72 h, respectively, and changes in histology, oxidative stress, gene transcription, metabolism, and intestinal microbiota were investigated. After NOD exposure, the hepatopancreas and intestine showed obvious histopathological damage and elevated oxidative stress response. Transcription patterns of immune genes related to detoxification, prophenoloxidase and coagulation system were altered in the hepatopancreas. Furthermore, metabolic patterns, especially amino acid metabolism and arachidonic acid related metabolites, were also disturbed. The integration of differential genes and metabolites revealed that the functions of "alanine, aspartic acid and glutamate metabolism" and "aminoacyl-tRNA biosynthesis" were highly affected. Alternatively, NOD exposure induced the variation of the diversity and composition of intestinal microbiota, especially the abundance of potentially beneficial bacteria (Demequina, Phyllobacterium and Pseudoalteromonas) and pathogenic bacteria (Photobacterium and Vibrio). Several intestinal bacteria were correlated with the changes of host the metabolic function and immune factors. These results revealed the toxic effects of NOD on shrimp, and identified some biomarkers.

Fragment screening and structural analyses highlight the ATP-assisted ligand binding for inhibitor discovery against type 1 methionyl-tRNA synthetase

Methionyl-tRNA synthetase (MetRS) charges tRNAMet with l-methionine (L-Met) to decode the ATG codon for protein translation, making it indispensable for all cellular lives. Many gram-positive bacteria use a type 1 MetRS (MetRS1), which is considered a promising antimicrobial drug target due to its low sequence identity with human cytosolic MetRS (HcMetRS, which belongs to MetRS2). Here, we report crystal structures of a representative MetRS1 from Staphylococcus aureus (SaMetRS) in its apo and substrate-binding forms. The connecting peptide (CP) domain of SaMetRS differs from HcMetRS in structural organization and dynamic movement. We screened 1049 chemical fragments against SaMetRS preincubated with or without substrate ATP, and ten hits were identified. Four cocrystal structures revealed that the fragments bound to either the L-Met binding site or an auxiliary pocket near the tRNA CCA end binding site of SaMetRS. Interestingly, fragment binding was enhanced by ATP in most cases, suggesting a potential ATP-assisted ligand binding mechanism in MetRS1. Moreover, co-binding with ATP was also observed in our cocrystal structure of SaMetRS with a class of newly reported inhibitors that simultaneously occupied the auxiliary pocket, tRNA site and L-Met site. Our findings will inspire the development of new MetRS1 inhibitors for fighting microbial infections.

A Leucyl-tRNA Synthetase Urzyme: Authenticity of tRNA Synthetase Catalytic Activities and Promiscuous Phosphorylation of Leucyl-5'AMP

Aminoacyl-tRNA synthetase (aaRS)/tRNA cognate pairs translate the genetic code by synthesizing specific aminoacyl-tRNAs that are assembled on messenger RNA by the ribosome. Deconstruction of the two distinct aaRS superfamilies (Classes) has provided conceptual and experimental models for their early evolution. Urzymes, containing ~120–130 amino acids excerpted from regions where genetic coding sequence complementarities have been identified, are key experimental models motivated by the proposal of a single bidirectional ancestral gene. Previous reports that Class I and Class II urzymes accelerate both amino acid activation and tRNA aminoacylation have not been extended to other synthetases. We describe a third urzyme (LeuAC) prepared from the Class IA Pyrococcus horikoshii leucyl-tRNA synthetase. We adduce multiple lines of evidence for the authenticity of its catalysis of both canonical reactions, amino acid activation and tRNA^Leu aminoacylation. Mutation of the three active-site lysine residues to alanine causes significant, but modest reduction in both amino acid activation and aminoacylation. LeuAC also catalyzes production of ADP, a non-canonical enzymatic function that has been overlooked since it first was described for several full-length aaRS in the 1970s. Structural data suggest that the LeuAC active site accommodates two ATP conformations that are prominent in water but rarely seen bound to proteins, accounting for successive, in situ phosphorylation of the bound leucyl-5′AMP phosphate, accounting for ADP production. This unusual ATP consumption regenerates the transition state for amino acid activation and suggests, in turn, that in the absence of the editing and anticodon-binding domains, LeuAC releases leu-5′AMP unusually slowly, relative to the two phosphorylation reactions.

Multidimensional Phylogenetic Metrics Identify Class I Aminoacyl-tRNA Synthetase Evolutionary Mosaicity and Inter-Modular Coupling

The role of aminoacyl-tRNA synthetases (aaRS) in the emergence and evolution of genetic coding poses challenging questions concerning their provenance. We seek evidence about their ancestry from curated structure-based multiple sequence alignments of a structurally invariant “scaffold” shared by all 10 canonical Class I aaRS. Three uncorrelated phylogenetic metrics—mutation frequency, its uniformity, and row-by-row cladistic congruence—imply that the Class I scaffold is a mosaic assembled from successive genetic sources. Metrics for different modules vary in accordance with their presumed functionality. Sequences derived from the ATP– and amino acid– binding sites exhibit specific two-way coupling to those derived from Connecting Peptide 1, a third module whose metrics suggest later acquisition. The data help validate: (i) experimental fragmentations of the canonical Class I structure into three partitions that retain catalytic activities in proportion to their length; and (ii) evidence that the ancestral Class I aaRS gene also encoded a Class II ancestor in frame on the opposite strand. A 46-residue Class I “protozyme” roots the Class I tree prior to the adaptive radiation of the Rossmann dinucleotide binding fold that refined substrate discrimination. Such rooting implies near simultaneous emergence of genetic coding and the origin of the proteome, resolving a conundrum posed by previous inferences that Class I aaRS evolved after the genetic code had been implemented in an RNA world. Further, pinpointing discontinuous enhancements of aaRS fidelity establishes a timeline for the growth of coding from a binary amino acid alphabet.