Radial scan of the molecular electrostatic potential of RNA double helices: An application to the enzyme-tRNA recognition

We introduced a method to characterize quantitatively the molecular electrostatic potential (MEP) of the minor and major grooves of base pairs located at nucleic acid double helices. By means of a radial MEP scan, we obtained a n-tuple of potential values corresponding to each groove, which can be analyzed by plotting the MEP values as a function of the angle in the radial scan. We studied base pairs of two different tRNAs, relevant in the recognition process with their cognate aminoacyl tRNA synthetases (aaRSs), in order to correlate their electrostatic behavior with the corresponding aminoacylation activity. We analyzed the first three base pairs of the Escherichia coli tRNAAla acceptor stem, finding several cases where the MEP profiles obtained from the plots are in agreement with the reported aminoacylation activities. Additionally, a non-hierarchical clustering performed over the MEP n-tuples resulted in meaningful classifications that correlate with the activity and with the predicted stereochemistry of the reaction. We also studied the first two base pairs of the E. coli tRNAThr acceptor stem but constraining the analysis to the angle intervals that seem relevant for the binding sites of the enzyme. These intervals were deduced from the ThrRS-tRNAThr complex crystal structure. In this case, we also found a good agreement between the MEP profiles and the activity, supporting the idea that the tRNA identity elements function is to allow an optimal electrostatic complementarity between the aminoacyl-tRNA synthetase and the tRNA.