We propose an improved transfer entropy approach called the dynamic version of the force constant fitted Gaussian network model based on molecular dynamics ensemble (dfcfGNMMD) to explore the allosteric mechanism of human mitochondrial phenylalanyl-tRNA synthetase (hmPheRS), one of the aminoacyl-tRNA synthetases that play a crucial role in translation of the genetic code. The dfcfGNMMD method can provide reliable estimates of the transfer entropy and give new insights into the role of the anticodon binding domain in driving the catalytic domain in aminoacylation activity and into the effects of tRNA binding and residue mutation on the enzyme activity, revealing the causal mechanism of the allosteric communication in hmPheRS. In addition, we incorporate the residue dynamic and co-evolutionary information to further investigate the key residues in hmPheRS allostery. This study sheds light on the mechanisms of hmPheRS allostery and can provide important information for related drug design.

Aminoacyl-tRNA synthetases (aaRSs), a family of ubiquitous and essential enzymes, can bind target tRNAs and catalyze the aminoacylation reaction in genetic code translation. In this work, we explore the dynamic properties and allosteric communication of human mitochondrial phenylalanyl-tRNA synthetase (hmPheRS) in free and bound states to understand the mechanisms of its tRNAPhe recognition and allostery using molecular dynamics simulations combined with the torsional mutual information-based network model. Our results reveal that hmPheRS\'s residue mobility and inter-residue motional coupling are significantly enhanced by tRNAPhe binding, and there occurs a strong allosteric communication which is critical for the aminoacylation reaction, suggesting the vital role of tRNAPhe binding in the enzyme\'s function. The identified signaling pathways mainly make the connections between the anticodon binding domain (ABD) and catalytic domain (CAD), as well as within the CAD composed of many functional fragments and active sites, revealing the co-regulation role of them to act coordinately and achieve hmPheRS\'s aminoacylation function. Besides, several key residues along the communication pathways are identified to be involved in mediating the coordinated coupling between anticodon recognition at the ABD and activation process at the CAD, showing their pivotal role in the allosteric network, which are well consistent with the experimental observation. This study sheds light on the allosteric communication mechanism in hmPheRS and can provide important information for the structure-based drug design targeting aaRSs.

The GA codon box incorporates the two-fold degeneracy of aspartic acid and of glutamic acid. Using the molecular mechanics approach of the AMBER suite, the four codons of the GA box are paired via H-bonding with two aspartic acid anticodons and two glutamic acid anticodons to yield 8 cognate and 11 non-cognate codon-anticodon duplexes. In addition four select non-cognate duplexes between the GA box codons and three alanine anticodons are also studied. These 23 duplexes display a variety of base-pairing possibilities at the wobble position. Cognate duplexes are differentiated from non-cognate duplexes on the grounds of structure and stability (chiefly the former). The results are in line with Crick\'s wobble hypothesis, and corroborate the observed reading properties of the aspartic acid anticodons GUC and QUC and of the glutamic acid anticodons CUC and SmnUC.

Tautomerism is important in many biomolecular interactions, not least in RNA biology. Crystallographic studies show the possible presence of minor tautomer forms of transfer-RNA (tRNA) anticodon bases in the ribosome. The hydrogen positions are not resolved in the X-ray studies, and we have used ab initio calculations and molecular dynamics simulations to understand if and how the minor enol form of uracil (U), or the modified uracil 5-oxyacetic acid (cmo5U), can be accommodated in the tRNA-messenger-RNA interactions in the ribosome decoding center. Ab initio calculations on isolated bases show that the modification affects the keto-enol equilibrium of the uracil base only slightly; the keto form is dominant (>99.99%) in both U and cmo5U. Other factors such as interactions with the surrounding nucleotides or ions would be required to shift the equilibrium toward the enol tautomer. Classical molecular simulations show a better agreement with the X-ray structures for the enol form, but free energy calculations indicate that the most stable form is the keto. In the ribosome, the enol tautomers of U and cmo5U pair with a guanine forming two hydrogen bonds, which do not involve the enol group. The oxyacetic acid modification has a minor effect on the keto-enol equilibrium.

Eukaryotic and archaeal translation initiation complexes have a common structural core comprising e/aIF1, e/aIF1A, the ternary complex (TC, e/aIF2-GTP-Met-tRNAiMet) and mRNA bound to the small ribosomal subunit. e/aIF2 plays a crucial role in this process but how this factor controls start codon selection remains unclear. Here, we present cryo-EM structures of the full archaeal 30S initiation complex showing two conformational states of the TC. In the first state, the TC is bound to the ribosome in a relaxed conformation with the tRNA oriented out of the P site. In the second state, the tRNA is accommodated within the peptidyl (P) site and the TC becomes constrained. This constraint is compensated by codon/anticodon base pairing, whereas in the absence of a start codon, aIF2 contributes to swing out the tRNA. This spring force concept highlights a mechanism of codon/anticodon probing by the initiator tRNA directly assisted by aIF2.

We report the first systematic evolution and study of tRNA variants that are able to read a set of UAGN (N = A, G, U, C) codons in a genomically recoded E. coli strain that lacks any endogenous in-frame UAGN sequences and release factor 1. Through randomizing bases in anticodon stem-loop followed by a functional selection, we identified tRNA mutants with significantly improved UAGN decoding efficiency, which will augment the current efforts on genetic code expansion through quadruplet decoding. We found that an extended anticodon loop with an extra nucleotide was required for a detectable efficiency in UAGN decoding. We also observed that this crucial extra nucleotide was converged to a U (position 33.5) in all of the top tRNA hits no matter which UAGN codon they suppress. The insertion of U33.5 in the anticodon loop likely causes tRNA distortion and affects anticodon-codon interaction, which induces +1 frameshift in the P site of ribosome. A new model was proposed to explain the observed features of UAGN decoding. Overall, our findings elevate our understanding of the +1 frameshift mechanism and provide a useful guidance for further efforts on the genetic code expansion using a non-canonical quadruplet reading frame.

Hypermodified nucleosides lysidine (L) and N(6)-threonylcarbamoyladenosine (t(6)A) influence codon-anticodon interactions during the protein biosynthesis process. Lysidine prevents the misrecognition of the AUG codon as isoleucine and that of AUA as methionine. The structural significance of these modified bases has not been studied in detail at the atomic level. Hence, in the present study we performed multiple molecular dynamics (MD) simulations of anticodon stem loop (ASL) of tRNA(Ile) in the presence and absence of modified bases \'L\' and \'t(6)A\' at the 34th and 37th positions respectively along with trinucleotide \'AUA\' and \'AUG\' codons. Hydrogen bonding interactions formed by the tautomeric form of lysidine may assist in reading the third base adenine of the \'AUA\' codon, unlike the guanine of the \'AUG\' codon. Such interactions might be useful to restrict codon specificity to recognize isoleucine tRNA instead of methionine tRNA. The t(6)A side chain interacts with the purine ring of the first codon nucleotide adenine, which might provide base stacking interactions and could be responsible for restricting extended codon-anticodon recognition. We found that ASL tRNA(Ile) in the absence of modifications at the 34th and 37th positions cannot establish proper hydrogen bonding interactions to recognize the isoleucine codon \'AUA\' and subsequently disturbs the anticodon loop structure. The binding free energy calculations revealed that tRNA(Ile) ASL with modified nucleosides prefers the codon AUA over AUG. Thus, these findings might be useful to understand the role of modified bases L and t(6)A to recognize the AUA codon instead of AUG.

Over the past several years, structural studies have led to the unexpected discovery of iron-sulfur clusters in enzymes that are involved in DNA replication/repair and protein biosynthesis. Although these clusters are generally well-studied cofactors, their significance in the new contexts often remains elusive. One fascinating example is a tryptophanyl-tRNA synthetase from the thermophilic bacterium Thermotoga maritima, TmTrpRS, that has recently been structurally characterized. It represents an unprecedented connection among a primordial iron-sulfur cofactor, RNA and protein biosynthesis. Here, a possible role of the [Fe4S4] cluster in tRNA anticodon-loop recognition is investigated by means of density functional theory and comparison with the structure of a human tryptophanyl-tRNA synthetase/tRNA complex. It turns out that a cluster-coordinating cysteine residue, R224, and polar main chain atoms form a characteristic structural motif for recognizing a putative 5\' cytosine or 5\' 2-thiocytosine moiety in the anticodon loop of the tRNA molecule. This motif provides not only affinity but also specificity by creating a structural and energetical penalty for the binding of other bases, such as uracil.

Two alternative hypotheses aim to predict the wobble nucleotide of tRNA anticodons in mitochondrion. The codon-anticodon adaptation hypothesis predicts that the wobble nucleotide of tRNA anticodon should evolve toward maximizing the Watson-Crick base pairing with the most frequently used codon within each synonymous codon family. In contrast, the wobble versatility hypothesis argues that the nucleotide at the wobble site should be occupied by a nucleotide most versatile in wobble pairing, i.e., the wobble site of the tRNA anticodon should be G for NNY codon families and U for NNR and NNN codon families (where Y stands for C or U, R for A or G, and N for any nucleotide). We examined codon usage and anticodon wobble sites in 36 fungal genomes to evaluate these two alternative hypotheses and identify exceptional cases that deserve new explanations. While the wobble versatility hypothesis is generally supported, there are interesting exceptions involving tRNA(Arg) translating the CGN codon family, tRNA(Trp) translating the UGR codon family, and tRNA(Met) translating the AUR codon family. Our results suggest that the potential to suppress stop codons, the historical inertia, and the conflict between translation initiation and elongation can all contribute to determining the wobble nucleotide of tRNA anticodons.

The Crick wobble hypothesis attributes the phenomenon of codon degeneracy to a certain impreciseness of pairing between the third base of the codon and the first base of the anticodon. This theoretical study investigates the pairing properties of some wobble bases, including both, observed and unobserved pairs. Some wobble base-pairs are predicted to follow the Watson-Crick pairs in configuration and pairing facility, while others deviate from this norm. The observed U:V pair is unique in that a pairing configuration may be suggested for it wherein the hydrogen-bonding involves the exocyclic 5-carboxymethoxy group of V. By comparing the theoretical data on the configurations of these pairs with the evidence for their existence/non-existence in nature, some guidelines emerge for differentiating between observed and unobserved base pairs on the basis of the pairing configuration.