Geometry and Interaction Energy

We have studied the geometric and energetic properties of the six water-mediated base pairs WUC, WGA, WUA, WUU, WUG, WGG (W: water) and of the related water-free complexes by quantum-chemical calculations including electron correlation at the MP2/6-31G(d,p)//HF/6-31G(d,p) and B3LYP/6-31G(d,p) levels. In the water-mediated WUC, WGA and WUU base pairs water is both donor and acceptor of H-bonds and for these complexes the calculated hydrogen bonding patterns are close to the experimental ones within the RNA helix and allow cooperative effects within the H-bond network. The geometries of these base pairs are obviously almost not affected by the nucleic acid and solvent environment. Therefore, they can be considered as structurally autonomous building blocks of RNA. In the optimized structure of the WUA base pair water still links the two bases, yet the H-bond pattern deviates somewhat from the one in the RNA crystal structure and it is not cooperative. Therefore, we classify the WUA base pair as an intermediate case. In WUG and WGG pairs water is twice acceptor or twice donor. These complexes are not structurally autonomous, but apparently have to be stabilized by additional interactions with the surrounding nucleic acid and solvent. For the UG pair a water-mediated C5-H5(U)...O6(G) contact involving an alternative water molecule (371) in the major groove is important and leads to a base pair geometry resembling the experimental structure (W 371 UG). In the GG pair short contacts to a backbone C5'-H5' donor group and to a further water molecule are likely to stabilize the experimental base pair 2 geometry. Similar to exocyclic amino groups water involved in base pairing induces nonplanar equilibrium geometries, yet except for one of the examples (WGA) the energy difference to the corresponding planar conformation is very small. Cooperative effects have been found to contribute between 9 and 13% at the MP2/HF level and between 8 and 20% at the B3LYP level of the total interaction energy to all water-mediated pairs except for WUA. WUC is the only example for which the absolute value of the total interaction energy per Hbond is significantly increasing in passing from the direct to the water-mediated pair.

Metal ions are crucial for folding and function of noncoding RNAs. The fact that RNAs have very specific metal ion binding motifs further implies that contribution of metal ions (like Mg2+) in RNA's folding is not limited to simple compensation of electrostatic repulsions. Rather, their binding to RNA is driven by very specific contextual requirements. Elucidation of such factors is necessary for a comprehensive understanding of the sequence-structure-function paradigm in RNA. In this work, we have studied the consequences of Mg2+ binding on the geometry and stability of different noncanonical base pairs that shape up the complex structural landscape of RNA. Our results show that majority of the Mg2+ bound nucleobases are also part of a base pair. Interestingly, such base pairs belong only to a specific set of base pairing geometries. Out of them, we are able to identify 14 unique cases for which the native base pairing geometries are unstable under gas phase geometry optimizati...

Noncanonical base pairs in RNA have strong structural and functional implications but are currently not considered for secondary structure predictions. We present results of comparative ab initio studies of stabilities and interaction energies for the three standard and 24 selected unusual RNA base pairs reported in the literature. Hydrogen added models of isolated base pairs, with heavy atoms frozen in their ‘away from equilibrium’ geometries, built from coordinates extracted from NDB, were geometry optimized using HF/6-31G** basis set, both before and after unfreezing the heavy atoms. Interaction energies, including BSSE and deformation energy corrections, were calculated, compared with respective single point MP2 energies, and correlated with occurrence frequencies and with types and geometries of hydrogen bonding interactions. Systems having two or more N-H...O/N hydrogen bonds had reasonable interaction energies which correlated well with respective occurrence frequencies and h...

Noncanonical base pairs in RNA have strong structural and functional implications but are currently not considered for secondary structure predictions. We present results of comparative ab initio studies of stabilities and interaction energies for the three standard and 24 selected unusual RNA base pairs reported in the literature. Hydrogen added models of isolated base pairs, with heavy atoms frozen in their 'away from equilibrium' geometries, built from coordinates extracted from NDB, were geometry optimized using HF/6-31G** basis set, both before and after unfreezing the heavy atoms. Interaction energies, including BSSE and deformation energy corrections, were calculated, compared with respective single point MP2 energies, and correlated with occurrence frequencies and with types and geometries of hydrogen bonding interactions. Systems having two or more N-H…O/N hydrogen bonds had reasonable interaction energies which correlated well with respective occurrence frequencies and highlighted the possibility of some of them playing important roles in improved secondary structure prediction methods. Several of the remaining base pairs with one N-H…O/N and/or one C-H…O/N interactions respectively, had poor interaction energies and negligible occurrences. High geometry variations on optimization of some of these were suggestive of their conformational switch like characteristics.

Molecular level understanding of physiochemical interactions responsible for the structure and dynamics of functional RNA is a major challenge to the scientific community. RNA structural diversity is mainly due to the presence of non-canonical base-pairs, phosphate mediated interactions and interactions between two or more nucleotides involving water molecules and ions that stabilize folded RNA structures. BPFIND and INCAR are two complementary tools developed to analyse RNA 3D structure in terms of hydrogen-bonded secondary structural elements interacting through primarily non-canonical tertiary base pairing interactions. Identification of recurring combinations of such elements or motifs, important in terms of architecture enhances our understanding of the RNA tertiary structures. First tool, Base Pair Finder (BPFIND), detects and analyzes all canonical and non-canonical basepairs and base triples involving at least two direct hydrogen bonds formed between polar atoms of the bases...

Structural bioinformatics of RNA has evolved mainly in response to the rapidly accumulating evidence that non-(protein)-coding RNAs (ncRNAs) play critical roles in gene regulation and development. The structures and functions of most ncRNAs are however still unknown. Most of the available RNA structural databases rely heavily on known 3D structures, and contextually correlate base pairing geometry with actual 3D RNA structures. None of the databases provide any direct information about stabilization energies. However, the intrinsic interaction energies of constituent base pairs can provide significant insights into their roles in the overall dynamics of RNA motifs and structures. Quantum mechanical (QM) computations provide the only approach toward their accurate quantification and characterization. 'RNA Base Pair Count, Geometry and Stability' (http://bioinf.iiit.ac.in/RNABPCOGEST) brings together information, extracted from literature data, regarding occurrence frequency, experimental and quantum chemically optimized geometries, and computed interaction energies, for non-canonical base pairs observed in a non-redundant dataset of functional RNA structures. The database is designed to enable the QM community, on the one hand, to identify appropriate biologically relevant model systems and also enable the biology community to easily sift

Noncanonical base pairs in RNA have strong structural and functional implications but are currently not considered for secondary structure predictions. We present results of comparative ab initio studies of stabilities and interaction energies for the three standard and 24 selected unusual RNA base pairs reported in the literature. Hydrogen added models of isolated base pairs, with heavy atoms frozen in their 'away from equilibrium' geometries, built from coordinates extracted from NDB, were geometry optimized using HF/6-31G** basis set, both before and after unfreezing the heavy atoms. Interaction energies, including BSSE and deformation energy corrections, were calculated, compared with respective single point MP2 energies, and correlated with occurrence frequencies and with types and geometries of hydrogen bonding interactions. Systems having two or more N-H…O/N hydrogen bonds had reasonable interaction energies which correlated well with respective occurrence frequencies and highlighted the possibility of some of them playing important roles in improved secondary structure prediction methods. Several of the remaining base pairs with one N-H…O/N and/or one C-H…O/N interactions respectively, had poor interaction energies and negligible occurrences. High geometry variations on optimization of some of these were suggestive of their conformational switch like characteristics.

The trans Watson-Crick/Watson-Crick family of base pairs represent a geometric class that play important structural and possible functional roles in the ribosome, tRNA, and other functional RNA molecules. They nucleate base triplets and quartets, participate as loop closing terminal base pairs in hair pin motifs and are also responsible for several tertiary interactions that enable sequentially distant regions to interact with each other in RNA molecules. Eleven representative examples spanning nine systems belonging to this geometric family of RNA base pairs, having widely different occurrence statistics in the PDB database, were studied at the HF/6-31G (d, p) level using Morokuma decomposition, Atoms in Molecules as well as Natural Bond Orbital methods in the optimized gas phase geometries and in their crystal structure geometries, respectively. The BSSE and deformation energy corrected interaction energy values for the optimized geometries are compared with the corresponding values in the crystal geometries of the base pairs. For non protonated base pairs in their optimized geometry, these values ranged from-8.19 kcal/mol to-21.84 kcal/mol and compared favorably with those of canonical base pairs. The interaction energies of these base pairs, in their respective crystal geometries, were, however, lesser to varying extents and in one case, that of A:A W:W trans, it was actually found to be positive. The variation in RMSD between the two geometries was also large and ranged from 0.32-2.19 Å. Our analysis shows that the hydrogen bonding characteristics and interaction energies obtained, correlated with the nature and type of hydrogen bonds between base pairs; but the occurrence frequencies, interaction energies, and geometric variabilities were conspicuous by the absence of any apparent correlation. Instead, the nature of local interaction energy hyperspace of different base pairs as inferred from the degree of their respective geometric variability could be correlated with the identities of free and bound hydrogen bond donor/acceptor groups present in interacting bases in conjunction with their tertiary and neighboring group interaction potentials in the global context. It also suggests that the concept of isostericity alone may not always determine covariation potentials for base pairs, particularly for those which may be important for RNA dynamics. These considerations are more important than the absolute values of the interaction energies in their respective optimized geometries in rationalizing their occurrences in functional RNAs. They highlight the importance of revising some of the existing DNA based structure analysis approaches and may have significant implications for RNA structure and dynamics, especially in the context of structure prediction algorithms.

Using experimentally measured values of density, viscosity and refractive index at 293 K, nature, strength and variation of molecular interactions have been studied in aqueous solutions of urea. Theoretically evaluated ultrasound velocities in the mixtures have been used to predict a large number of acoustical parameters of the systems. Thermodynamical parameters such as interaction parameter, the excess free energy of mixing, and interaction energy have been calculated using experimental viscosity data. Space-filling factor, molar refraction and Ekyman parameter are evaluated using measured refractive indices of the mixtures. The variation of the excess parameters such as excess viscosity, excess refractive index, excess ultrasound velocity, excess adiabatic compressibility and excess molar volume have been studied. The results have been explained in terms of molecular interactions present in the systems under study.

Three programs have been developed to aid in the classification and visualization of RNA structure. BPViewer provides a web interface for displaying three-dimensional (3D) coordinates of individual base pairs or base pair collections. A web server, RNAview, automatically identifies and classifies the types of base pairs that are formed in nucleic acid structures by various combinations of the three edges, Watson-Crick, Hoogsteen and the Sugar edge. RNAView produces two-dimensional (2D) diagrams of secondary and tertiary structure in either Postscript, VRML or RNAML formats. The application RNAMLview can be used to rearrange various parts of the RNAView 2D diagram to generate a standard representation (like the cloverleaf structure of tRNAs) or any layout desired by the user. A 2D diagram can be rapidly reformatted using RNAMLview since all the parts of RNA (like helices and single strands) are dynamically linked while moving the selected parts. With the base pair annotation and the 2D graphic display, RNA motifs are rapidly identified and classified. A survey has been carried out for 41 unique structures selected from the NDB database. The statistics for the occurrence of each edge and of each of the 12 bp families are given for the combinations of the four bases: A, G, U and C. The program also allows for visualization of the base pair interactions by using a symbolic convention previously proposed for base pairs. The web servers for BPViewer and RNAview are available at http://ndbserver.rutgers.edu/ services/. The application RNAMLview can also be downloaded from this site. The 2D diagrams produced by RNAview are available for RNA structures in the Nucleic Acid Database (NDB) at http:// ndbserver.rutgers.edu/atlas/.

The base pairing patterns in RNA structures are more versatile and completely different as compared to DNA. We present here results of ab-initio studies of structures and interaction energies of eight selected RNA base pairs reported in literature. Interaction energies, including BSSE correction, of hydrogen added crystal geometries of base pairs have been calculated at the HF/6-31G** level. The structures and interaction energies of the base pairs in the crystal geometry are compared with those obtained after optimization of the base pairs. We find that the base pairs become more planar on full optimization. No change in the hydrogen bonding pattern is seen. It is expected that the inclusion of appropriate considerations of many of these aspects of RNA base pairing would significantly improve the accuracy of RNA secondary structure prediction.

Noncanonical base pairs in RNA have strong structural and functional implications but are currently not considered for secondary structure predictions. We present results of comparative ab initio studies of stabilities and interaction energies for the three standard and 24 selected unusual RNA base pairs reported in the literature. Hydrogen added models of isolated base pairs, with heavy atoms frozen in their 'away from equilibrium' geometries, built from coordinates extracted from NDB, were geometry optimized using HF/6-31G** basis set, both before and after unfreezing the heavy atoms. Interaction energies, including BSSE and deformation energy corrections, were calculated, compared with respective single point MP2 energies, and correlated with occurrence frequencies and with types and geometries of hydrogen bonding interactions. Systems having two or more N-H…O/N hydrogen bonds had reasonable interaction energies which correlated well with respective occurrence frequencies and highlighted the possibility of some of them playing important roles in improved secondary structure prediction methods. Several of the remaining base pairs with one N-H…O/N and/or one C-H…O/N interactions respectively, had poor interaction energies and negligible occurrences. High geometry variations on optimization of some of these were suggestive of their conformational switch like characteristics. [Bhattacharyya D, Koripella S C, Mitra A, Rajendran V B and Sinha B 2007 Theoretical analysis of noncanonical base pairing interactions in RNA molecules; J. Biosci. 32 809-825]

Base pairs involving post-transcriptionally modified nucleobases are believed to play important roles in a wide variety of functional RNAs. Here we present our attempts towards understanding the structural and functional role of naturally occurring modified base pairs by analyzing their distribution in different RNA classes, with the help of crystal structure and sequence database analyses. In addition, we quantify the variation in geometrical features of modified base pairs within RNA structures, and characterize their optimum geometries and binding energies using advanced quantum chemical methods. Further comparison of modified base pairs with their unmodified counterparts illustrates the effect of steric and electronic structure alterations due to base modifications. Analysis of specific structural contexts of modified base pairs in RNA crystal structures revealed several interesting scenarios, including those at the tRNA:rRNA interface, antibiotic-binding site and the three-way junctions within tRNA, which when analyzed in context of available experimental data, allowed us to correlate the occurrence and strength of modified base pairs with the specific functional roles they play in context of RNA macromolecules.

Structural bioinformatics of RNA has evolved mainly in response to the rapidly accumulating evidence that non-(protein)-coding RNAs (ncRNAs) play critical roles in gene regulation and development. The structures and functions of most ncRNAs are however still unknown. Most of the available RNA structural databases rely heavily on known 3D structures, and contextually correlate base pairing geometry with actual 3D RNA structures. None of the databases provide any direct information about stabilization energies. However, the intrinsic interaction energies of constituent base pairs can provide significant insights into their roles in the overall dynamics of RNA motifs and structures. Quantum mechanical (QM) computations provide the only approach toward their accurate quantification and characterization. 'RNA Base Pair Count, Geometry and Stability' (http://bioinf.iiit.ac.in/RNABPCOGEST) brings together information, extracted from literature data, regarding occurrence frequency, experimental and quantum chemically optimized geometries, and computed interaction energies, for non-canonical base pairs observed in a non-redundant dataset of functional RNA structures. The database is designed to enable the QM community, on the one hand, to identify appropriate biologically relevant model systems and also enable the biology community to easily sift V C The Author(s) Downloaded from through diverse computational results to gain theoretical insights which could promote hypothesis driven biological research.

The understanding of folding and function of RNA molecules depends on the identification and classification of interactions between ribonucleotide residues. We developed a new method named ClaRNA for computational classification of contacts in RNA 3D structures. Unique features of the program are the ability to identify imperfect contacts and to process coarse-grained models. Each doublet of spatially close ribonucleotide residues in a query structure is compared to clusters of reference doublets obtained by analysis of a large number of experimentally determined RNA structures, and assigned a score that describes its similarity to one or more known types of contacts, including pairing, stacking, base-phosphate and base-ribose interactions. The accuracy of ClaRNA is 0.997 for canonical base pairs, 0.983 for non-canonical pairs and 0.961 for stacking interactions. The generalized squared correlation coefficient (GC2) for ClaRNA is 0.969 for canonical base pairs, 0.638 for non-canonic...

Description: VARNA is a tool for the automated drawing, visualization and annotation of the secondary structure of RNA, designed as a companion software for web servers and databases. Features: VARNA implements four drawing algorithms, supports input/output using the classic formats dbn, ct, bpseq and RNAML and exports the drawing as five picture formats, either pixel-based (JPEG, PNG) or vector-based (SVG, EPS and XFIG). It also allows manual modification and structural annotation of the resulting drawing using either an interactive point and click approach, within a web server or through command-line arguments. Availability: VARNA is a free software, released under the terms of the GPLv3.0 license and available at http://varna.lri.fr Contact: ponty@lri.fr

Cis and trans sugar edge/sugar edge (SE/SE) binding patterns are essential building units of RNAs. For example, SE/SE interactions form the A-minor motifs, the most important tertiary interaction type in functional RNAs. This study provides an in-depth structure and stability analysis for these two base pair families. Gas-phaseoptimized geometries are reported for 12 cis and 7 trans SE/SE base pairs and contrasted to their X-ray counterparts. Interaction energies are computed at the RIMP2 level of theory using the density-functionaltheory-optimized geometries. There is a good overall agreement between the optimized and X-ray geometries of the cis SE/SE base pairs. In contrast, only three of the seven trans SE/SE binding patterns could be optimized without a significant distortion of the X-ray geometry. Note, however, that many SE/SE base pairs participate in broader networks of interactions; thus it is not surprising to see some of them to deviate from the X-ray geometry in a complete isolation. Computed interaction energies reveal that all 12 known cis SE/SE binding patterns are very stable. Among the trans SE/SE binding patterns, only the rG/rG, rG/rC, and rA/rG base pairs are sufficiently stable in the crystal geometry. Prediction has been made for some structures not yet detected by crystallography, namely, cis rC/rC, rG/rC, rG/rU, and rU/rU and trans rG/rA base pairs. Interestingly, the new cis SE/SE binding patterns are not necessarily isosteric with the remaining 12 members of this family. The trans rG/rA base pair represents a viable option for base pairing in RNA to be identified by future X-ray studies. In a complete lack of structural information, prediction of other unknown members of the trans SE/SE family was not attempted. Analysis of the interaction energies shows a very large electron correlation component of the interaction energy, pointing at the elevated role of dispersion energy as compared to other types of base pairs. This likely is profitable for stabilization of SE/SE binding patterns in polar environments and could be one of the reasons why the A-minor motif is the leading type of tertiary interactions in RNAs.

Complex molecular shapes of ribosomal RNA molecules are stabilized by recurrent types of tertiary interactions involving highly specific and conserved non-Watson-Crick base pairs, triplets, and quartets. We analyzed the intrinsic structure and stability of the P-motif and the four basic types of A-minor interactions (types I, II, III, and 0), which represent the most prominent RNA tertiary interaction patterns refined in the course of evolution. In the studied interactions, the electron correlation component of the stabilization usually exceeds the Hartree-Fock (HF) term, leading to a strikingly different balance of forces as compared to standard base pairing stabilized primarily by the HF term. In other words, the A-minor and P-interactions are considerably more influenced by the dispersion energy as compared to canonical base pairs, which makes them particularly suitable to zip the folded RNA structures that are substantially hydrated even in their interior. Continuum solvent COSMO calculations confirm that the stability of the canonical GC base pair is affected (reduced) by the continuum solvent screening considerably more than the stability of the A-minor interaction. Among the studied systems, the strong A-minor II and weak A-minor III interactions require water molecules to stabilize the experimental geometry. Gas-phase optimization of the canonical A-minor II A/CG triplet without water results in a geometry that is clearly inconsistent with the RNA structure. The gas-phase structure of the P-interaction and the most stable A-minor I interaction nicely agrees with the geometries occurring in the ribosome. A-minor I can also adopt an alternative water-mediated substate rather often observed in X-ray and molecular dynamics studies. The A-minor I water bridge, however, does not appear to stabilize the tertiary contact, and its role is to provide structural flexibility to this binding pattern within the context of the RNA structure. Interestingly, the insertion of a polar water molecule in the A-minor I A/CG tertiary contact occurring in the A/C tertiary pair is stabilized primarily by the HF (electrostatic) interaction energy, while the dispersioncontrolled A/G contact remains firmly bound. Thus, the intrinsic balance of forces as revealed by quantum mechanics (QM) calculations nicely correlates with many behavioral aspects of the studied interactions inside RNA. The comparison of interaction energies computed using quantum chemistry and an AMBER force field reveals that common molecular mechanics calculations perform rather well, except that the strength of the P-interaction is modestly overestimated. We also briefly discuss the non-negligible methodological differences when evaluating simple base-base nucleic acids base pairs and the complex RNA tertiary base pairing patterns using QM procedures.

The coordination of cis-[Pt(NH3)2]2+ to the N7 and O6 sites of guanine of the guanine–cytosine (GC) nucleic base pair is studied at the SCF, DFT and MP2 levels of theory, and by an ab initio BSSE-free optimization algorithm, concerning the possible mechanisms of the antitumor activity of cis-[Pt(NH3)2Cl2]. The calculations show that the cis-[Pt(NH3)2]2+ coordination results in the breakage of

Complex molecular shapes of ribosomal RNA molecules are stabilized by recurrent types of tertiary interactions involving highly specific and conserved non-Watson-Crick base pairs, triplets, and quartets. We analyzed the intrinsic structure and stability of the P-motif and the four basic types of A-minor interactions (types I, II, III, and 0), which represent the most prominent RNA tertiary interaction patterns refined in the course of evolution. In the studied interactions, the electron correlation component of the stabilization usually exceeds the Hartree-Fock (HF) term, leading to a strikingly different balance of forces as compared to standard base pairing stabilized primarily by the HF term. In other words, the A-minor and P-interactions are considerably more influenced by the dispersion energy as compared to canonical base pairs, which makes them particularly suitable to zip the folded RNA structures that are substantially hydrated even in their interior. Continuum solvent COSMO calculations confirm that the stability of the canonical GC base pair is affected (reduced) by the continuum solvent screening considerably more than the stability of the A-minor interaction. Among the studied systems, the strong A-minor II and weak A-minor III interactions require water molecules to stabilize the experimental geometry. Gas-phase optimization of the canonical A-minor II A/CG triplet without water results in a geometry that is clearly inconsistent with the RNA structure. The gas-phase structure of the P-interaction and the most stable A-minor I interaction nicely agrees with the geometries occurring in the ribosome. A-minor I can also adopt an alternative water-mediated substate rather often observed in X-ray and molecular dynamics studies. The A-minor I water bridge, however, does not appear to stabilize the tertiary contact, and its role is to provide structural flexibility to this binding pattern within the context of the RNA structure. Interestingly, the insertion of a polar water molecule in the A-minor I A/CG tertiary contact occurring in the A/C tertiary pair is stabilized primarily by the HF (electrostatic) interaction energy, while the dispersioncontrolled A/G contact remains firmly bound. Thus, the intrinsic balance of forces as revealed by quantum mechanics (QM) calculations nicely correlates with many behavioral aspects of the studied interactions inside RNA. The comparison of interaction energies computed using quantum chemistry and an AMBER force field reveals that common molecular mechanics calculations perform rather well, except that the strength of the P-interaction is modestly overestimated. We also briefly discuss the non-negligible methodological differences when evaluating simple base-base nucleic acids base pairs and the complex RNA tertiary base pairing patterns using QM procedures.

Abstract. RNA structural bioinformatics has developed, into a vibrant interdisciplinary area, in response to the emergence of noncoding RNA molecules as key regulatory players in nearly all important biological processes. Here we introduce the different conceptual tools and available resources, related to this new area. We also review its application in studying the role of RNA-protein interactions.

Base pairs involving protonated nucleobases play important roles in mediating global macromolecular conformational changes and in facilitation of catalysis in a variety of functional RNA molecules. Here we present our attempts at understanding the role of such base pairs by detecting possible protonated base pairs in the available RNA crystal structures using BPFind software, in their specific structural contexts, and by the characterization of their geometries, interaction energies, and stabilities using advanced quantum chemical computations. We report occurrences of 18 distinct protonated base pair combinations from a representative data set of RNA crystal structures and propose a theoretical model for one putative base pair combination. Optimization of base pair geometries was carried out at the B3LYP/cc-pVTZ level, and the BSSE corrected interaction energies were calculated at the MP2/aug-cc-pVDZ level of theory. The geometries for each of the base pairs were characterized in terms of H-bonding patterns observed, rmsd values observed on optimization, and base pair geometrical parameters. In addition, the intermolecular interaction in these complexes was also analyzed using Morokuma energy decomposition. The gas phase interaction energies of the base pairs range from -24 to -49 kcal/mol and reveal the dominance of Hartree-Fock component of interaction energy constituting 73% to 98% of the total interaction energy values. On the basis of our combined bioinformatics and quantum chemical analysis of different protonated base pairs, we suggest resolution of structural ambiguities and correlate their geometric and energetic features with their structural and functional roles. In addition, we also examine the suitability of specific base pairs as key elements in molecular switches and as nucleators for higher order structures such as base triplets and quartets.

The trans Watson-Crick/Watson-Crick family of base pairs represent a geometric class that play important structural and possible functional roles in the ribosome, tRNA, and other functional RNA molecules. They nucleate base triplets and quartets, participate as loop closing terminal base pairs in hair pin motifs and are also responsible for several tertiary interactions that enable sequentially distant regions to interact with each other in RNA molecules. Eleven representative examples spanning nine systems belonging to this geometric family of RNA base pairs, having widely different occurrence statistics in the PDB database, were studied at the HF/6-31G (d, p) level using Morokuma decomposition, Atoms in Molecules as well as Natural Bond Orbital methods in the optimized gas phase geometries and in their crystal structure geometries, respectively. The BSSE and deformation energy corrected interaction energy values for the optimized geometries are compared with the corresponding values in the crystal geometries of the base pairs. For non protonated base pairs in their optimized geometry, these values ranged from -8.19 kcal/mol to -21.84 kcal/mol and compared favorably with those of canonical base pairs. The interaction energies of these base pairs, in their respective crystal geometries, were, however, lesser to varying extents and in one case, that of A:A W:W trans, it was actually found to be positive. The variation in RMSD between the two geometries was also large and ranged from 0.32-2.19 Å. Our analysis shows that the hydrogen bonding characteristics and interaction energies obtained, correlated with the nature and type of hydrogen bonds between base pairs; but the occurrence frequencies, interaction energies, and geometric variabilities were conspicuous by the absence of any apparent correlation. Instead, the nature of local interaction energy hyperspace of different base pairs as inferred from the degree of their respective geometric variability could be correlated with the identities of free and bound hydrogen bond donor/acceptor groups present in interacting bases in conjunction with their tertiary and neighboring group interaction potentials in the global context. It also suggests that the concept of isostericity alone may not always determine covariation potentials for base pairs, particularly for those which may be important for RNA dynamics. These considerations are more important than the absolute values of the interaction energies in their respective optimized geometries in rationalizing their occurrences in functional RNAs. They highlight the importance of revising some of the existing DNA based structure analysis approaches and may have significant implications for RNA structure and dynamics, especially in the context of structure prediction algorithms.

The base pairing patterns in RNA structures are more versatile and completely different as compared to DNA. We present here results of ab-initio studies of structures and interaction energies of eight selected RNA base pairs reported in literature. Interaction energies, including BSSE correction, of hydrogen added crystal geometries of base pairs have been calculated at the HF/6-31G** level. The structures and interaction energies of the base pairs in the crystal geometry are compared with those obtained after optimization of the base pairs. We find that the base pairs become more planar on full optimization. No change in the hydrogen bonding pattern is seen. It is expected that the inclusion of appropriate considerations of many of these aspects of RNA base pairing would significantly improve the accuracy of RNA secondary structure prediction.

Noncanonical base pairs in RNA have strong structural and functional implications but are currently not considered for secondary structure predictions. We present results of comparative ab initio studies of stabilities and interaction energies for the three standard and 24 selected unusual RNA base pairs reported in the literature. Hydrogen added models of isolated base pairs, with heavy atoms frozen in their ‘away from equilibrium’ geometries, built from coordinates extracted from NDB, were geometry optimized using HF/6-31G** basis set, both before and after unfreezing the heavy atoms. Interaction energies, including BSSE and deformation energy corrections, were calculated, compared with respective single point MP2 energies, and correlated with occurrence frequencies and with types and geometries of hydrogen bonding interactions. Systems having two or more N-H...O/N hydrogen bonds had reasonable interaction energies which correlated well with respective occurrence frequencies and highlighted the possibility of some of them playing important roles in improved secondary structure prediction methods. Several of the remaining base pairs with one N-H...O/N and/or one C-H...O/N interactions respectively, had poor interaction energies and negligible occurrences. High geometry variations on optimization of some of these were suggestive of their conformational switch like characteristics.

The base pairing patterns in RNA structures are more versatile and completely different as compared to DNA. We present here results of ab-initio studies of structures and interaction energies of eight selected RNA base pairs reported in literature. Interaction energies, including BSSE correction, of hydrogen added crystal geometries of base pairs have been calculated at the HF/6-31G** level. The structures and interaction energies of the base pairs in the crystal geometry are compared with those obtained after optimization of the base pairs. We find that the base pairs become more planar on full optimization. No change in the hydrogen bonding pattern is seen. It is expected that the inclusion of appropriate considerations of many of these aspects of RNA base pairing would significantly improve the accuracy of RNA secondary structure prediction.