Going beyond base-pairs: topology-based characterization of base-multiplets in RNA.

Identification and characterization of base-multiplets, which are essentially mediated by base-pairing interactions, can provide insights into the diversity in the structure and dynamics of complex functional RNAs, and thus facilitate hypothesis driven biological research. The necessary nomenclature scheme, an extension of the geometric classification scheme for base-pairs by Leontis and Westhof, is however available only for base-triplets. In the absence of information on topology, this scheme is not applicable to quartets and higher order multiplets. Here we propose a topology-based classification scheme which, in conjunction with a graph-based algorithm, can be used for the automated identification and characterization of higher order base-multiplets in RNA structures. Here, the RNA structure is represented as a graph, where nodes represent nucleotides and edges represent base-pairing connectivity. Sets of connected components (of n nodes) within these graphs constitute subgraphs representing multiplets of "n" nucleotides. The different topological variants of the RNA multiplets thus correspond to different nonisomorphic forms of these subgraphs. To annotate RNA base-multiplets unambiguously, we propose a set of topology-based nomenclature rules for quartets, which are extendable to higher multiplets. We also demonstrate the utility of our approach toward the identification and annotation of higher order RNA multiplets, by investigating the occurrence contexts of selected examples in order to gain insights regarding their probable functional roles.

On the Nature of Nucleobase Stacking in RNA: A Comprehensive Survey of Its Structural Variability and a Systematic Classification of Associated Interactions.

The astonishing diversity in folding patterns of RNA three-dimensional (3D) structures is crafted by myriads of noncovalent contacts, of which base pairing and stacking are the most prominent. A systematic and comprehensive classification and annotation of these interactions is necessary for a molecular-level understanding of their roles. However, unlike in the case of base pairing, where a widely accepted nomenclature and classification scheme exists in the public domain, currently available classification schemes for base-base stacking need major enhancements to comprehensively capture the necessary features underlying the rich stacking diversity in RNA. Here, we extend the previous stacking classification based on nucleobase interacting faces by introducing a structurally intuitive geometry-cum topology-based scheme. Specifically, a stack is first classified in terms of the geometry described by the relative orientation of the glycosidic bonds, which generates eight basic stacking geometric families for heterodimeric stacks and six of those for homodimeric stacks. Further annotation in terms of the identity of the bases and the region of involvement of purines (five-membered, six-membered, or both rings) leads to the enumeration of 384 distinct RNA base stacks. Based on our classification scheme, we present an algorithm for automated identification of stacks in RNA crystal structures and analyze the stacking context in selected RNA structures. Overall, the work described here is expected to greatly facilitate the structure-based RNA research.

Consequences of Mg2+ binding on the geometry and stability of RNA base pairs.

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 optimization carried out in the absence of Mg2+ binding. Our density functional theory based calculations, performed using dispersion corrected M05-2X functional, suggest that, depending on its mode of binding, Mg2+ can stabilize and even fine tune a number of such base pairing geometries. These findings not only provide insights into how metal ions modulate the structure and dynamics of RNA molecules, they also provide a basis for improving the RNA structure prediction algorithms.

How Does Mg2+ Modulate the RNA Folding Mechanism: A Case Study of the G:C W:W Trans Basepair.

Reverse Watson-Crick G:C basepairs (G:C W:W Trans) occur frequently in different functional RNAs. This is one of the few basepairs whose gas-phase-optimized isolated geometry is inconsistent with the corresponding experimental geometry. Several earlier studies indicate that through post-transcriptional modification, direct protonation, or coordination with Mg2+, accumulation of positive charge near N7 of guanine can stabilize the experimental geometry. Interestingly, recent studies reveal significant variation in the position of putatively bound Mg2+. This, in conjunction with recently raised doubts regarding some of the Mg2+ assignments near the imino nitrogen of guanine, is suggestive of the existence of multiple Mg2+ binding modes for this basepair. Our detailed investigation of Mg2+-bound G:C W:W Trans pairs occurring in high-resolution RNA crystal structures shows that they are found in 14 different contexts, eight of which display Mg2+ binding at the Hoogsteen edge of guanine. Further examination of occurrences in these eight contexts led to the characterization of three different Mg2+ binding modes: 1) direct binding via N7 coordination, 2) direct binding via O6 coordination, and 3) binding via hydrogen-bonding interaction with the first-shell water molecules. In the crystal structures, the latter two modes are associated with a buckled and propeller-twisted geometry of the basepair. Interestingly, respective optimized geometries of these different Mg2+ binding modes (optimized using six different DFT functionals) are consistent with their corresponding experimental geometries. Subsequent interaction energy calculations at the MP2 level, and decomposition of its components, suggest that for G:C W:W Trans , Mg2+ binding can fine tune the basepair geometries without compromising with their stability. Our results, therefore, underline the importance of the mode of binding of Mg2+ ions in shaping RNA structure, folding and function.

The role of N7 protonation of guanine in determining the structure, stability and function of RNA base pairs.

The roles of protonated nucleobases in stabilizing different structural motifs and in facilitating catalytic functions of RNA are well known. Among different polar sites of all the nucleobases, N7 of guanine has the highest protonation propensity at physiological pH. However, unlike other easily protonable sites such as N1 and N3 of adenine or N3 of cytosine, N7 protonation of guanine does not lead to the stabilization of base pairs involving its protonated Hoogsteen edge. It also does not facilitate its participation in any acid-base catalysis process. To explore the possible roles of N7 protonated guanine, we have studied its base pairing potentials involving WatsonCrick and sugar edges, which undergo major charge redistribution upon N7 protonation. We have carried out quantum chemical geometry optimization at the M05-2X/6-311G+(2d,2p) level, followed by interaction energy calculation at the MP2/aug-cc-pVDZ level, along with the analysis of the context of occurrence for selected base pairs involving the sugar edge or the WatsonCrick edge of guanine within a non-redundant set of 167 RNA crystal structures. Our results suggest that, four base pairs - G:C W:W trans, G:rC W:S cis, G:G W:H cis and G:G S:H trans may involve N7 protonated guanine. These base pairs deviate significantly from their respective experimental geometries upon QM optimization, but they retain their experimental geometries if guanine N7 protonation is considered during optimization. Our study also reveals the role of guanine N7 protonation (i) in stabilizing important RNA structural motifs, (ii) in providing a framework for designing pH driven molecular motors and (iii) in providing an alternative strategy to mimic the effect of post-transcriptional changes.

Feasibility of occurrence of different types of protonated base pairs in RNA: a quantum chemical study.

Protonated nucleobases have significant roles in facilitating catalytic functions of RNA, and in stabilizing different structural motifs. Reported pKa values of nucleobase protonation suggest that the population of neutral nucleobases is 10(3)-10(4) times higher than that of protonated nucleobases under physiological conditions (pH ∼ 7.4). Therefore, a molecular level understanding of various putative roles of protonated nucleobases cannot be achieved without addressing the question of how their occurrence propensities and stabilities are related to the free energy costs associated with the process of protonation under physiological conditions. With water as the proton donor, we use advanced QM methods to evaluate the site specific protonation propensities of nucleobases in terms of their associated free energy changes (ΔGprot). Quantitative follow up on the energetics of base pair formation and database search for evaluating their occurrence frequencies, reveal a lack of correlation between base pair stability and occurrence propensities on the one hand, and ease of protonation on the other. For example, although N7 protonated adenine (ΔGprot = 40.0 kcal mol(-1)) is found to participate in stable base pairing, base pairs involving N7 protonated guanine (ΔGprot = 36.8 kcal mol(-1)), on geometry optimization, converge to a minima where guanine transfers its extra proton to its partner base. Such observations, along with examples of weak base pairs involving N3 protonation of cytosine (ΔGprot = 37.0 kcal mol(-1)) are rationalized by analysing the protonation induced charge redistributions which are found to significantly influence, both positively and negatively, the hydrogen bonding potentials of different functional sites of individual nucleobases. Protonation induced charge redistribution is also found to strongly influence (i) the aromatic character of the rings of the participating bases and (ii) hydrogen bonding potential of the free edges of the protonated base pair. Comprehensive analysis of a non-redundant RNA crystal structure dataset further reveals that, while availability of stabilization possibilities determine the feasibility of occurrence of protonated bases, their occurrence context and specific functional roles are important factors determining their occurrence propensities.

Mg2+ Sensing by an RNA Fragment: Role of Mg2+-Coordinated Water Molecules.

RNA molecules selectively bind to specific metal ions to populate their functional active states, making it important to understand their source of ion selectivity. In large RNA systems, metal ions interact with the RNA at multiple locations, making it difficult to decipher the precise role of ions in folding. To overcome this complexity, we studied the role of different metal ions (Mg2+, Ca2+, and K+) in the folding of a small RNA hairpin motif (5'-ucCAAAga-3') using unbiased all-atom molecular dynamics simulations. The advantage of studying this system is that it requires specific binding of a single metal ion to fold to its native state. We find that even for this small RNA, the folding free energy surface (FES) is multidimensional as different metal ions present in the solution can simultaneously facilitate folding. The FES shows that specific binding of a metal ion is indispensable for its folding. We further show that in addition to the negatively charged phosphate groups, the spatial organization of electronegative nucleobase atoms drives the site-specific binding of the metal ions. Even though the binding site cannot discriminate between different metal ions, RNA folds efficiently only in a Mg2+ solution. We show that the rigid network of Mg2+-coordinated water molecules facilitates the formation of important interactions in the transition state. The other metal ions such as K+ and Ca2+ cannot facilitate the formation of such interactions. These results allow us to hypothesize possible metal-sensing mechanisms in large metalloriboswitches and also provide useful insights into the design of appropriate collective variables for studying large RNA molecules using enhanced sampling methods.

Role of Guanidinium-Carboxylate Ion Interaction in Enzyme Inhibition with Implications for Drug Design.

Guanidinium cation (Gdm+) interacts strongly with amino acids of different polarities modulating protein structure and function. Using density functional theory calculations and molecular dynamics simulations, we studied the interaction of Gdm+ with carboxylate ions mimicking its interaction with acidic amino acids and explored its effect in enzymatic folding and activity. We show that, in low concentrations, Gdm+ stabilizes carboxylate ion dimers by acting as a bridge between them, thereby reducing the electrostatic repulsion. We further show that this carboxylate-Gdm+-carboxylate interaction can have an effect on the structure-activity relationship in enzymes with active sites containing two acidic residues. Using five enzymes (hen egg white lysozyme, T4 lysozyme, HIV-1 protease, pepsin, and creatine kinase), which have two acidic amino acids in their active sites, we show that, in low concentrations (<0.5 M), Gdm+ strongly binds to the enzyme active site, thereby potentially inhibiting its activity without unfolding it. This can lead to misleading conclusions in experiments, which infer the extent of enzyme unfolding from activity measurements. However, the carboxylate-Gdm+-carboxylate specific interaction can be exploited in drug discovery as drugs based on guanidinium derivatives are already being used to treat various maladies related to muscle weakness, cancer, diabetes etc. Guanidinium derivatives can be designed as potential drug molecules to inhibit activity or functioning of enzymes, which have binding pockets with two acidic residues in close vicinity.

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

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 → σ*-type charge transfer from filled π 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: ρ, (ii) its Laplacian: ∇2ρ, (iii) stabilization energy due to n → σ*-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.

Evidence for Hidden Involvement of N3-Protonated Guanine in RNA Structure and Function.

Charged nucleobases have been found to occur in several known RNA molecules and are considered essential for their structure and function. The mechanism of their involvement is however not yet fully understood. Revelation of the role of N7-protonated guanine, in modulating the geometry and stability of noncanonical base pairs formed through its unprotonated edges [Watson-Crick (WC) and sugar], has triggered the need to evaluate the feasibility of similar roles of other protonated nucleobases [Halder et al., Phys Chem Chem Phys, 2015, 17, 26249]. In this context, N3 protonation of guanine makes an interesting case as its influence on the charge distribution of the WC edge is similar to that of N7 protonation, though its thermodynamic cost of protonation is significantly higher. In this work, we have carried out structural bioinformatics analyses and quantum mechanics-based calculations to show that N3 protonation of guanine may take place in a cellular environment, at least in the G:C W:W Trans and G:G W:H Cis base pairs. Our results provide a reasonable starting point for future investigations in order to address the larger mechanistic question.

RNABP COGEST: a resource for investigating functional RNAs.

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 through diverse computational results to gain theoretical insights which could promote hypothesis driven biological research.

Why does substitution of thymine by 6-ethynylpyridone increase the thermostability of DNA double helices?

Efficiency of 6-ethynylpyridone (E), a potential thymine (T) analogue, which forms high-fidelity base pairs with adenine (A) and gives rise to stabler DNA duplexes, with stability comparable to those containing canonical cytosine(C):guanine(G) base pairs, has been reported recently. Estimates of the interaction energies, involving geometry optimization at the DFT level (including middle range dispersion interactions) followed by single point energy calculation at MP2 level, in excellent correlation with the experimentally observed trends, show that E binds more strongly and more discriminately with A than T does. Detailed analysis reveals that the increase in base-base interaction arises out of conjugation of acetylenic π electrons with the ring π system of E, which results in not only an extra stabilizing C-H···π interaction in the EA pair, but also a strengthening of the conventional hydrogen bonds. However, the computed base-base interaction energy for the EA pair was found to be much less than that of the canonical CG pair, implying that the difference in the TA versus EA base pairing interaction alone cannot explain the large experimentally observed increase in the thermostability of DNA duplexes, where a TA pair is replaced with an EA pair. Our computations show that the conjugation of acetylenic π electrons with the ring π system also possibly plays a role in increasing the stacking potential of the EA pair, which in turn can explain its marked influence in the enhancement of duplex stability.