Abstract
With the discovery of regulatory roles of RNA molecules, several scientific efforts have been made to understand its structural diversity. It has been shown that, apart from the normal and modified base and noncanonical base pairs involving them, protonated bases and protonated base pairs play an important role in shaping up the complex folded structure of the RNA molecule. But the process of protonation of nucleobases is energetically unfavourable in physiological pH. It is therefore expected that, in order to compensate the cost of protonation, the protonation induced charge redistribution and geometric changes will help the protonated base to participate in different higher order interactions such as base pair, base triple, base stack, etc. Quite a few efforts has been made in recent years to understand the kinetics and thermodynamics of the process of protonation and stabilization of the protonated bases. We have observed that there is a lack of relevance with the physiological context in the methods adopted in most of those studies. Moreover, these studies mainly focus on describing the process of protonation as a function of the difference between pH of the local environment and average pKa of the nucleobase. Significance of local pKa values of different titrable sites of different RNA bases and influence of protonation on their stabilizing factors like noncovalent interactions, have not got much attention. For the sake of comparative analysis and completeness, in our study, we have considered deprotonation of RNA bases as well.
On the basis of the results obtained from advanced ab-initio quantum mechanical calculations, in this study, we therefore have proposed a physicochemically relevant method of studying the thermodynamics of the process of protonation in terms of estimation of site specific protonation propensities of RNA bases. It considers water molecule as proton donor. We have further tried to correlate the site specific protonation propensity data with the occurrence frequency of the protonated bases. We have noticed that the stability and hence the occurrence frequency of protonated bases do not depend only on their protonation propensity but the protonation induced modifications of stabilizing factors, such as, hydrogen bonding and pi-stacking potentials play a role of greater importance. We have extensively studied the case of protonated base pairs involving N7 protonated purines. Our study revealed that, although the N7 site of Guanine has a good protonation propensity, N7 protonated Guanine can not occur in the stable planar base pairs due to reasons explained by the protonation induced charge redistribution in Guanine. N7 protonated Adenine, on the other hand, can form stable and planar base pairs and we have explored two new geometries of base pairs involving N7 protonated Adenine -- A(+):G H:H Cis and Trans. Based on the charge redistribution over both the participating bases in a protonated base pair, we have also explained the stability of base triples involving protonated bases.
Our study therefore provides a systematic estimation of site specific protonation propensities of RNA bases and evaluates the influence of protonation on hydrogen bonding and pi-stacking via them. On the basis of the both it explains stability and occurrence frequency of protonated bases and base pairs involving them.