ePrints@IIScePrints@IISc Home | About | Browse | Latest Additions | Advanced Search | Contact | Help

Estimating Strengths of Individual Hydrogen Bonds in RNA Base Pairs: Toward a Consensus between Different Computational Approaches

Haider, Antarip and Data, Dhruv and Seelam, Preethi P and Bhattacharyya, Dhananjay and Mitra, Abhijit (2019) Estimating Strengths of Individual Hydrogen Bonds in RNA Base Pairs: Toward a Consensus between Different Computational Approaches. In: ACS OMEGA, 4 (4). pp. 7354-7368.

Acs_Omg_4_4_7354-7368_2019.pdf - Published Version

Download (2MB) | Preview
ao8b03689_si_001.pdf - Published Supplemental Material

Download (1MB) | Preview
[img] Archive (ZIP)
Dhruv-paper-optimized-geometries.rar - Published Supplemental Material

Download (455kB)
Official URL: https://doi.org/10.1021/acsomega.8b03689


Noncoding RNA molecules are composed of a large variety of noncanonical base pairs that shape up their functionally competent folded structures. Each base pair is composed of at least two interbase hydrogen bonds (H-bonds). It is expected that the characteristic geometry and stability of different noncanonical base pairs are determined collectively by the properties of these interbase H-bonds. We have studied the ground-state electronic properties using density functional theory (DFT) and DFT-D3-based methods] of all the 118 normal base pairs and 36 modified base pairs, belonging to 12 different geometric families (cis and trans of WW, WH, HH, WS, HS, and SS) that occur in a nonredundant set of high-resolution RNA crystal structures. Having addressed some of the limitations of the earlier approaches, we provide here a comprehensive compilation of the average energies of different types of interbase H-bonds (E-HB). We have also characterized each interbase H-bond using 13 different parameters that describe its geometry, charge distribution at its bond critical point (BCP), and n -> sigma*-type charge transfer from filled pi orbitals of the H-bond acceptor to the empty antibonding orbital of the H-bond donor. On the basis of the extent of their linear correlation with the H-bonding energy, we have shortlisted five parameters to model linear equations for predicting E-HB values. They are (i) electron density at the BCP: rho, (ii) its Laplacian: del(2)rho, (iii) stabilization energy due to n -> sigma*-type charge transfer: E(2), (iv) donor-hydrogen distance, and (v) hydrogen-acceptor distance. We have performed single variable and multivariable linear regression analysis over the normal base pairs and have modeled sets of linear relationships between these five parameters and E-HB. Performance testing of our model over the set of modified base pairs shows promising results, at least for the moderately strong H-bonds.

Item Type: Journal Article
Additional Information: copyright for this article belongs to American Chemical Society
Department/Centre: Division of Chemical Sciences > Solid State & Structural Chemistry Unit
Depositing User: Id for Latest eprints
Date Deposited: 06 Aug 2019 10:31
Last Modified: 06 Aug 2019 10:31
URI: http://eprints.iisc.ac.in/id/eprint/62889

Actions (login required)

View Item View Item