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

Occurrence and stability of anion-π interactions between phosphate and nucleobases in functional RNA molecules

We present a systematic structural and energetic characterization of phosphate(OP)-nucleobase anion…π stacking interactions in RNAs. We observed OP-nucleobase stacking contacts in a variety of structural motifs other than regular helices and spanning broadly diverse sequence distances. Apart from the stacking between a phosphate and a guanine or a uracil two-residue upstream in specific U-turns, such interactions in RNA have been scarcely characterized to date. Our QM calculations showed an energy minimum at a distance between the OP atom and the nucleobase plane centroid slightly below 3 Å for all the nucleobases. By sliding the OP atom over the nucleobase plane we localized the optimal mutual positioning of the stacked moieties, corresponding to an energy minimum below -6 kcal•mol-1, for all the nucleobases, consistently with the projections of the OP atoms over the different π-rings we observed in experimental occurrences. We also found that the strength of the interaction clearly correlates with its electrostatic component, pointing to it as the most relevant contribution. Finally, as OP-uracil and OP-guanine interactions represent together 86% of the instances we detected, we also proved their stability under dynamic conditions in model systems simulated by state-of-the art DFT-MD calculations.

RNABPDB: Molecular Modeling of RNA Structure-From Base Pair Analysis in Crystals to Structure Prediction

The stable three-dimensional structure of RNA is known to play several important biochemical roles, from post-transcriptional gene regulation to enzymatic action. These structures contain double-helical regions, which often have different types of non-canonical base pairs in addition to Watson-Crick base pairs. Hence, it is important to study their structures from experimentally obtained or even predicted ones, to understand their role, or to develop a drug against the potential targets. Molecular Modeling of RNA double helices containing non-canonical base pairs is a difficult process, particularly due to the unavailability of structural features of non-Watson-Crick base pairs. Here we show a composite web-server with an associated database that allows one to generate the structure of RNA double helix containing non-canonical base pairs using consensus parameters obtained from the database. The database classification is followed by an evaluation of the central tendency of the structural parameters as well as a quantitative estimation of interaction strengths. These parameters are used to construct three-dimensional structures of double helices composed of Watson-Crick and/or non-canonical base pairs. Our benchmark study to regenerate double-helical fragments of many experimentally derived RNA structures indicate very high accuracy. This composite server is expected to be highly useful in understanding functions of various pre-miRNA by modeling structures of the molecules and estimating binding efficiency. The database can be accessed from http://hdrnas.saha.ac.in/rnabpdb .

Alkaline pH has an unexpected effect on transcriptional pausing during synthesis of the E. coli pH-responsive riboswitch

Riboswitches are 5′-untranslated regions of mRNA that change their conformation in response to ligand binding, allowing post-transcriptional gene regulation. This ligand-based model of riboswitch function has been expanded with the discovery of a “pH-responsive element” (PRE) riboswitch in Escherichia coli. At neutral pH, the PRE folds into a translationally inactive structure with an occluded ribosome-binding sequence, whereas at alkaline pH, the PRE adopts a translationally active structure. This unique riboswitch does not rely on ligand binding in a traditional sense to modulate its alternative folding outcomes. Rather, pH controls riboswitch folding by two possible modes that are yet to be distinguished; pH either regulates the transcription rate of RNA polymerase (RNAP) or acts on the RNA itself. Previous work suggested that RNAP pausing is prolonged by alkaline pH at two sites, stimulating PRE folding into the active structure. To date, there has been no rigorous exploration into how pH influences RNAP pausing kinetics during PRE synthesis. To provide that understanding and distinguish between pH acting on RNAP versus RNA, we investigated RNAP pausing kinetics at key sites for PRE folding under different pH conditions. We find that pH influences RNAP pausing but not in the manner proposed previously. Rather, alkaline pH either decreases or has no effect on RNAP pause longevity, suggesting that the modulation of RNAP pausing is not the sole mechanism by which pH affects PRE folding. These findings invite the possibility that the RNA itself actively participates in the sensing of pH.

Theoretical investigation of the mechanism of DMAP-promoted [4 + 2]-annulation of prop-2-ynylsulfonium with isatoic anhydride

The mechanism for DMAP-promoted [4 + 2]-annulation of prop-2-ynylsulfonium with isatoic anhydride is investigated using the M06-2X functional. The reaction comprises isomerization of prop-2-ynylsulfonium in stage 1. Stage 2 includes DMAP-promoted deprotonation, nucleophilic addition, ring opening, and decarboxylation. Three steps of intramolecular cycloaddition, DMAP-promoted protonation, and dealkylation occur in stage 3, generating methylated DMAP and neutral thioether, which undergo double-bond isomerization to yield 3-methylthio-4-quinolone. The ability of DMAP to promote the reaction lies in the barrier decrease for alkyne isomerization, deprotonation/protonation of allenes, and dealkylation as effective bases for transferring protons and methyl groups. The roles of prop-2-ynylsulfonium and isatoic anhydride were demonstrated to be C2 and C4 synthons via Multiwfn analysis on the frontier molecular orbital. An alternative path was also confirmed by the Mayer bond order of the vital transition states.

AC-motif: a DNA motif containing adenine and cytosine repeat plays a role in gene regulation

I-motif or C4 is a four-stranded DNA structure with a protonated cytosine:cytosine base pair (C+:C) found in cytosine-rich sequences. We have found that oligodeoxynucleotides containing adenine and cytosine repeats form a stable secondary structure at a physiological pH with magnesium ion, which is similar to i-motif structure, and have named this structure ‘adenine:cytosine-motif (AC-motif)’. AC-motif contains C⁺:C base pairs intercalated with putative A⁺:C base pairs between protonated adenine and cytosine. By investigation of the AC-motif present in the CDKL3 promoter (AC-motif_CDKL3), one of AC-motifs found in the genome, we confirmed that AC-motif_CDKL3 has a key role in regulating CDKL3 gene expression in response to magnesium. This is further supported by confirming that genome-edited mutant cell lines, lacking the AC-motif formation, lost this regulation effect. Our results verify that adenine-cytosine repeats commonly present in the genome can form a stable non-canonical secondary structure with a non-Watson–Crick base pair and have regulatory roles in cells, which expand non-canonical DNA repertoires.

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 (Mg²⁺, Ca²⁺, 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 Mg²⁺ solution. We show that the rigid network of Mg²⁺-coordinated water molecules facilitates the formation of important interactions in the transition state. The other metal ions such as K⁺ and Ca²⁺ 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.

Frequency and hydrogen bonding of nucleobase homopairs in small molecule crystals

We used the high resolution and accuracy of the Cambridge Structural Database (CSD) to provide detailed information regarding base pairing interactions of selected nucleobases. We searched for base pairs in which nucleobases interact with each other through two or more hydrogen bonds and form more or less planar structures. The investigated compounds were either free forms or derivatives of adenine, guanine, hypoxanthine, thymine, uracil and cytosine. We divided our findings into categories including types of pairs, protonation patterns and whether they are formed by free bases or substituted ones. We found base pair types that are exclusive to small molecule crystal structures, some that can be found only in RNA containing crystal structures and many that are native to both environments. With a few exceptions, nucleobase protonation generally followed a standard pattern governed by pKa values. The lengths of hydrogen bonds did not depend on whether the nucleobases forming a base pair were charged or not. The reasons why particular nucleobases formed base pairs in a certain way varied significantly.

An IRMPD spectroscopic and computational study of protonated guanine-containing mismatched base pairs in the gas phase

Infrared multiple photon dissociation spectroscopy has been used to probe the structures of the three protonated base-pair mismatches containing 9-ethylguanine (9eG) in the gas phase. Some of these protonated base-pairs have been identified in RNA.

The effect of adenine protonation on RNA phosphodiester backbone bond cleavage elucidated by deaza-nucleobase modifications and mass spectrometry

The catalytic strategies of small self-cleaving ribozymes often involve interactions between nucleobases and the ribonucleic acid (RNA) backbone. Here we show that multiply protonated, gaseous RNA has an intrinsic preference for the formation of ionic hydrogen bonds between adenine protonated at N3 and the phosphodiester backbone moiety on its 5'-side that facilitates preferential phosphodiester backbone bond cleavage upon vibrational excitation by low-energy collisionally activated dissociation. Removal of the basic N3 site by deaza-modification of adenine was found to abrogate preferential phosphodiester backbone bond cleavage. No such effects were observed for N1 or N7 of adenine. Importantly, we found that the pH of the solution used for generation of the multiply protonated, gaseous RNA ions by electrospray ionization affects phosphodiester backbone bond cleavage next to adenine, which implies that the protonation patterns in solution are at least in part preserved during and after transfer into the gas phase. Our study suggests that interactions between protonated adenine and phosphodiester moieties of RNA may play a more important mechanistic role in biological processes than considered until now.

Single-stranded DNA structural diversity: TAGGGT from monomers to dimers to tetramer formation

RATIONALE

DNA quadruplex structures have emerged as novel drug targets due to their role in preventing abnormal gene transcription and maintaining telomere stability. Trapped Ion Mobility Spectrometry-Mass Spectrometry (TIMS-MS), combined with theoretical modeling, is a powerful tool for studying the kinetic intermediates of DNA complexes formed in solution and interrogated in the gas phase after desolvation.

METHODS

A TAGGGT ssDNA sequence was purchased and studied in 10 mM ammonium acetate using nanospray electrospray ionization (nESI)-TIMS-MS in positive and negative ion mode. Collisional cross section (CCS) profiles were measured using internal calibration (Tune Mix). Theoretical structures were proposed based on molecular dynamics, charge location and geometry optimization for the most intense IMS bands based on the number of TAGGGT units, adduct form and charge states.

RESULTS

A distribution of monomeric, dimeric and tetrameric TAGGGT structures were formed in solution and separated in the gas phase based on their mobility and m/z value (e.g., [M + 2H]⁺² , [2M + 3H]⁺³ , [M - 2H]^-2 , [2M - 3H]^-3 , [4M + 4H]⁺⁴ , [4M + 3H + NH₄ ]⁺⁴ , [4M + 2H + 2NH₄ ]⁺⁴ and [4M + H + 3NH₄ ]⁺⁴ ). The high mobility resolution of the TIMS-MS analyzer permitted the observation of multiple CCS bands per molecular ion form. Comparison with theoretical candidate structures suggests that monomeric TAGGGT species are stabilized by A-T and G⁺ -G interactions, with the size of the conformer influenced by the proton location. In the case of the TAGGGT quadruplex, the protonated species displayed a broad CCS distribution, while six discrete conformers were stabilized by the presence of ammonium ions (n = 1-3).

CONCLUSIONS

This is the first observation of multiple conformations of TAGGGT complexes (n = 1, 2 and 4) in 10 mM ammonium acetate. Candidate structures with intramolecular interactions of the form of G⁺ -G and traditional A-T base pairing agreed with the experimental trends. Our results demonstrate the structural diversity of TAGGGT monomers, dimers and tetramers in the gas phase beyond the previously reported solution structure, using 10 mM ammonium acetate to replicate biological conditions.

Alternated Branching Ratios by Anomaly in Collision-Induced Dissociation of Proton-Bound Hoogsteen Base Pairs of 1-Methylcytosine with 1-Methylguanine and 9-Methylguanine

A comparative study on the proton-bound complexes of 1-methylcytosine (1-mC) with 1-methylguanine (1-mG) and 9-methylguanine (9-mG), [1-mC:1-mG:H]⁺ and [1-mC:9-mG:H]⁺, respectively, was carried out using energy-resolved collision-induced dissociation (ER-CID) experiments in combination with quantum chemical calculations. In ER-CID experiments, the measured survival yields indicated an essentially identical stability for the two proton-bound complexes. In comparison with the lowest-energy structures and base-pairing energetics predicted at the B3LYP/6-311+G(2d,2p) theory level, both complexes produced in this study were suggested to be proton-bound Hoogsteen base pairs. Curiously, despite the similarity in structures, binding energetics, and potential energy surfaces predicted by the B3LYP theory, the fragment branching ratios exhibited an intriguing alternation between the two proton-bound Hoogsteen base pairs. The CID of [1-mC:1-mG:H]⁺ produced protonated cytosines, [1-mC:H]⁺, more abundantly than [1-mG:H]⁺, whereas that of [1-mC:9-mG:H]⁺ gave rise to a more pronounced production of protonated guanines, [9-mG:H]⁺. However, using the proton affinities of moieties predicted by the high-accuracy methods, including CBS-QB3 and the Guassian-4 theory, the anomaly known for [Cytosine:Guanine:H]⁺ (J. Am. Soc. Mass Spectrom. 29, 2368-2379 (2018)) successfully accounted for the alternated branching ratios. Thereby, the anomaly, more specifically, the production of proton-transferred fragments of O-protonated cytosines in the CID of proton-bound Hoogsteen base pairs, is indeed real, which is disclosed as the alternated branching ratios in the CID spectra of [1-mC:1-mG:H]⁺ and [1-mC:9-mG:H]⁺ in this study. Graphical Abstract .

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.

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: ∇²ρ, (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.

Binding preference of nitroimidazolic radiosensitizers to nucleobases and nucleosides probed by electrospray ionization mass spectrometry and density functional theory

Nitroimidazolic radiosensitizers are used in radiation therapy to selectively sensitize cancer cells deprived of oxygen, and the actual mechanism of radiosensitization is still not understood. Selecting five radiosensitizers (1-methyl-5-nitroimidazole, ronidazole, ornidazole, metronidazole, and nimorazole) with a common 5-nitroimidazolic ring with different substitutions at N1 and C2 positions of the imidazole moiety, we investigate here their binding to nucleobases (A, T, G, and C) and nucleosides (As, Td, Gs, and Cd) via the positive electrospray ionization mass spectrometry experiments. In addition, quantum chemical calculations at the M062x/6-311+G(d,p) level of theory and basis set were used to determine binding energies of the proton bound dimers of a radiosensitizer and a nucleobase. The positive electrospray ionization leads to the formation of proton bound dimers of all radiosensitizers except 1-methyl-5-nitroimidazole in high abundance with C and smaller abundance with G. Ronidazole and metronidazole formed less abundant dimers also with A, while no dimers were observed to be formed at all with T. In contrast to the case of the nucleoside Td, the dimer intensity is as high as that with Cd, while the abundance of the dimer with Gs is smaller than that of the former. The experimental results are consistent with the calculations of binding energies suggesting proton bound dimers with C and G to be the strongest bound ones. Finally, a barrier-free proton transfer is observed when protonated G or C approaches the nitroimidazole ring.

Proton Transfer Accounting for Anomalous Collision-Induced Dissociation of Proton-Bound Hoogsteen Base Pair of Cytosine and Guanine

To understand the anomalous collision-induced dissociation (CID) behavior of the proton-bound Hoogsteen base pair of cytosine (C) and guanine (G), C:H⁺∙∙∙G, we investigated CID of a homologue series of proton-bound heterodimers of C, 1-methylcytosine, and 5-methylcytosine with G as a common base partner. The CID experiments were performed in an energy-resolved way (ER-CID) under both multiple and near-single collision conditions. The relative stabilities of the protonated complexes examined by ER-CID suggested that the proton-bound complexes produced by electrospray ionization in this study are proton-bound Hoogsteen base pairs. On the other hand, in contrast to the other base pairs, CID of C:H⁺∙∙∙G exhibited more abundant productions of C:H⁺, the fragment protonated on the moiety with a smaller proton affinity, than that of G:H⁺. This appeared to contradict general prediction based on the kinetic method. However, further theoretical exploration of potential energy surfaces found that there can be facile proton transfers in the proton-bound Hoogsteen base pairs during the CID process, which makes the process accessible to an additional product state of O-protonated C for C:H⁺ fragments. The presence of an additional dissociation channel, which in other words corresponds to twofold degeneracy in the transition state leading to C:H⁺ fragments, effectively doubles the apparent reaction rate for production of C:H⁺. In this way, the process gives rise to the anomaly, the observed pronounced formation of C:H⁺ in the CID of the proton-bound Hoogsteen base pair, C:H⁺∙∙∙G. Graphical Abstract ᅟ.

Protonated nucleobases are not fully ionized in their chloride salt crystals and form metastable base pairs further stabilized by the surrounding anions

This paper presents experimental charge-density studies of cytosinium chloride, adeninium chloride hemihydrate and guaninium dichloride crystals based on ultra-high-resolution X-ray diffraction data and extensive theoretical calculations. The results confirm that the cohesive energies of the studied systems are dominated by contributions from intermolecular electrostatic interactions, as expected for ionic crystals. Electrostatic interaction energies (Ees) usually constitute 95% of the total interaction energy. The Ees energies in this study were several times larger in absolute value when compared, for example, with dimers of neutral nucleobases. However, they were not as large as some theoretical calculations have predicted. This was because the molecules appeared not to be fully ionized in the studied crystals. Apart from charge transfer from chlorine to the protonated nucleobases, small but visible charge redistribution within the nucleobase cations was observed. Some dimers of singly protonated bases in the studied crystals, namely a cytosinium-cytosinium trans sugar/sugar edge pair and an adeninium-adeninium trans Hoogsteen/Hoogsteen edge pair, exhibited attractive interactions (negative values of Ees) or unusually low repulsion despite identical molecular charges. The pairs are metastable as a result of strong hydrogen bonding between bases which overcompensates the overall cation-cation repulsion, the latter being weakened due to charge transfer and molecular charge-density polarization.

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

Quantum chemical calculations reveal the role of magnesium in stabilizing the geometries of intrinsically unstable RNA base pairs.

Occurrence and stability of lone pair-π stacking interactions between ribose and nucleobases in functional RNAs

The specific folding pattern and function of RNA molecules lies in various weak interactions, in addition to the strong base-base pairing and stacking. One of these relatively weak interactions, characterized by the stacking of the O4' atom of a ribose on top of the heterocycle ring of a nucleobase, has been known to occur but has largely been ignored in the description of RNA structures. We identified 2015 ribose-base stacking interactions in a high-resolution set of non-redundant RNA crystal structures. They are widespread in structured RNA molecules and are located in structural motifs other than regular stems. Over 50% of them involve an adenine, as we found ribose-adenine contacts to be recurring elements in A-minor motifs. Fewer than 50% of the interactions involve a ribose and a base of neighboring residues, while approximately 30% of them involve a ribose and a nucleobase at least four residues apart. Some of them establish inter-domain or inter-molecular contacts and often implicate functionally relevant nucleotides. In vacuo ribose-nucleobase stacking interaction energies were calculated by quantum mechanics methods. Finally, we found that lone pair-π stacking interactions also occur between ribose and aromatic amino acids in RNA-protein complexes.

Fast coarse-grained model for RNA titration

A new numerical scheme for RNA (ribonucleic acid) titration based on the Debye-Hückel framework for the salt description is proposed in an effort to reduce the computational costs for further applications to study protein-RNA systems. By means of different sets of Monte Carlo simulations, we demonstrated that this new scheme is able to correctly reproduce the experimental titration behavior and salt pK_a shifts. In comparison with other theoretical approaches, similar or even better outcomes are achieved at much lower computational costs. The model was tested on the lead-dependent ribozyme, the branch-point helix, and the domain 5 from Azotobacter vinelandii Intron 5.

Structural landscape of base pairs containing post-transcriptional modifications in RNA

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 toward understanding the structural and functional role of naturally occurring modified base pairs using a combination of X-ray crystal structure database analysis, sequence analysis, and advanced quantum chemical methods. Our bioinformatics analysis reveals that despite their presence in all major secondary structural elements, modified base pairs are most prevalent in tRNA crystal structures and most commonly involve guanine or uridine modifications. Further, analysis of tRNA sequences reveals additional examples of modified base pairs at structurally conserved tRNA regions and highlights the conservation patterns of these base pairs in three domains of life. Comparison of structures and binding energies of modified base pairs with their unmodified counterparts, using quantum chemical methods, allowed us to classify the base modifications in terms of the nature of their electronic structure effects on base-pairing. 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 sites on the ribosome, and the three-way junctions within tRNA. These scenarios, when analyzed in the context of available experimental data, allowed us to correlate the occurrence and strength of modified base pairs with their specific functional roles. Overall, our study highlights the structural importance of modified base pairs in RNA and points toward the need for greater appreciation of the role of modified bases and their interactions, in the context of many biological processes involving RNA.

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 Mg²⁺, 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 Mg²⁺. This, in conjunction with recently raised doubts regarding some of the Mg²⁺ assignments near the imino nitrogen of guanine, is suggestive of the existence of multiple Mg²⁺ binding modes for this basepair. Our detailed investigation of Mg²⁺-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 Mg²⁺ binding at the Hoogsteen edge of guanine. Further examination of occurrences in these eight contexts led to the characterization of three different Mg²⁺ 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 Mg²⁺ 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 , Mg²⁺ binding can fine tune the basepair geometries without compromising with their stability. Our results, therefore, underline the importance of the mode of binding of Mg²⁺ ions in shaping RNA structure, folding and function.

Crystal structure of a poly(rA) staggered zipper at acidic pH: evidence that adenine N1 protonation mediates parallel double helix formation

We have solved at 1.07 Å resolution the X-ray crystal structure of a polyriboadenylic acid (poly(rA)) parallel and continuous double helix. Fifty-nine years ago, double helices of poly(rA) were first proposed to form at acidic pH. Here, we show that 7-mer oligo(rA), i.e. rA₇, hybridizes and overlaps in all registers at pH 3.5 to form stacked double helices that span the crystal. Under these conditions, rA₇ forms well-ordered crystals, whereas rA₆ forms fragile crystalline-like structures, and rA₅, rA₈ and rA₁₁ fail to crystallize. Our findings support studies from ∼50 years ago: one showed using spectroscopic methods that duplex formation at pH 4.5 largely starts with rA₇ and begins to plateau with rA₈; another proposed a so-called ‘staggered zipper’ model in which oligo(rA) strands overlap in multiple registers to extend the helical duplex. While never shown, protonation of adenines at position N1 has been hypothesized to be critical for helix formation. Bond angles in our structure suggest that N1 is protonated on the adenines of every other rAMP−rAMP helix base pair. Our data offer new insights into poly(rA) duplex formation that may be useful in developing a pH sensor.

A comparison of four different conformations adopted by human telomeric G-quadruplex using computer simulations

The telomeric G-quadruplexes for their unique structural features are considered as potential anticancer drug targets. These, however, exhibit structural polymorphism as different topology types for the intra-molecular G-quadruplexes from human telomeric G-rich sequences have been reported based on NMR spectroscopy and X-ray crystallography. These techniques provide detailed atomic-level information about the molecule but relative conformational stability of the different topologies remains unsolved. Therefore, to understand the conformational preference, we have carried out quantum chemical calculations on G-quartets; used all-atom molecular dynamics (MD) simulations and steered molecular dynamics (SMD) simulations to characterize the four human telomeric G-quadruplex topologies based on its G-tetrad core-types, viz., parallel, anti-parallel, mixed-(3 + 1)-form1 and mixed-(3 + 1)-form2. We have also studied a non-telomeric sequence along with these telomeric forms giving a comparison between the two G-rich forms. The structural properties such as base pairing, stacking geometry and backbone conformations have been analyzed. The quantum calculations indicate that presence of a sodium ion inside the G-tetrad plane or two potassium ions on both sides of the plane give it an overall planarity which is much needed for good stacking to form a helix. MD simulations indicate that capping of the G-tetrad core by the TTA loops keep the terminal guanine bases away from water. The SMD simulations along with equilibrium MD studies indicate that the parallel and non-telomeric forms are comparatively less stable. We could come to the conclusion that the anti-parallel form and also the mixed-(3 + 1)-form1 topology are most likely to represent the major conformation., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 83-99, 2016.

The effects of tautomerization and protonation on the adenine-cytosine mismatches: a density functional theory study

In the present work, we demonstrate the results of a theoretical study concerned with the question how tautomerization and protonation of adenine affect the various properties of adenine-cytosine mismatches. The calculations, in gas phase and in water, are performed at B3LYP/6-311++G(d,p) level. In gas phase, it is observed that any tautomeric form of investigated mismatches is more stabilized when adenine is protonated. As for the neutral mismatches, the mismatches containing amino form of cytosine and imino form of protonated adenine are more stable. The role of aromaticity on the stability of tautomeric forms of mismatches is investigated by NICS(1)ZZ index. The stability of mispairs decreases by going from gas phase to water. It can be explained using dipole moment parameter. The influence of hydrogen bonds on the stability of mismatches is examined by atoms in molecules and natural bond orbital analyses. In addition to geometrical parameters and binding energies, the study of the topological properties of electron charge density aids in better understanding of these mispairs.

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.

Database URL:http://bioinf.iiit.ac.in/RNABPCOGEST

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

Ab initio computations and bioinformatics studies reveal that stabilization of some important RNA structural motifs might involve N7 protonation of guanine.