Determining the Identity of Nucleotides and the Energy of Binding of tRNAs to Their Aminoacyl-tRNA Synthetases Using a Simple Logistic Model

Abstract

This study showed that the predictor in logistic regression can be applied to estimating the Gibbs free energy of tRNAs’ recognition of and binding to their aminoacyl-tRNA synthetases. Then, 24 linear logistic regression models predicting different classes of tRNAs loaded with a corresponding amino acid were trained in a machine learning classification method, reducing the misclassification error to zero. The models were based on minimal subsets of Boolean explanatory variables describing the favorite presence of nucleotides or nucleosides localized in the different parts of the tRNA. In 90% of cases, they agree with the components of the consensus strand in a class of tRNAs loaded by a given amino acid. According to the proposed theoretical model, the values of the free energy for the entry of the recognition state in the process of tRNA charging were obtained, and the inputs from identity nucleotides and the tRNA strand backbone were distinguished. Almost all the resulting models indicated leading anticodon tandems defining the first and second positions of the anticodon (positions 35 and 36 of the tRNA strand) and the small sets (up to six positions) of the other nucleotides as the natural identity nucleotides most influential in the free energy balance. The magnitude of their input to this energy depends on the position in the strand, favoring positions −1, 35, and 36. The role of position 34 is relatively smaller. These identity attributes may not always be fully arranged in a real single adaptor molecule but were comprehensively present in a given tRNA class. A detailed analysis of the resulting models showed that the absolute value of the energy of binding the tandem 35–36 decreases with the number of identity positions, as well as with the decreasing number of possible hydrogen bonds. On the other hand, in these conditions, the absolute value of the energy of binding of other identity nucleotides increases. All the models indicate that the nucleotide-independent energy of the repulsion tRNA backbone decreases with the number of identity nucleotides. It was also shown that the total free energy change in entering the recognition state increases with the amino acid mass, making this process less spontaneous, which may have an evolutionary reference.